nus17_box_opening.C

Go to the documentation of this file.

#include "CAFAna/Cuts/Cuts.h"
#include "NuXAna/Cuts/NusCuts17.h"
#include "CAFAna/Cuts/SpillCuts.h"


#include "CAFAna/Core/EventList.h"
#include "CAFAna/Core/Utilities.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Vars/Vars.h"
#include "CAFAna/Vars/HistAxes.h"
#include "NuXAna/Vars/HistAxes.h"
#include "CAFAna/Vars/GenieWeights.h"
#include "CAFAna/Cuts/TimingCuts.h"
#include "OscLib/OscCalcDumb.h"
#include "CAFAna/Prediction/PredictionNoExtrap.h"

#include "TH1D.h"

using namespace ana;

std::vector<std::string> MakeUnblindList();

void nus17_box_opening()
{

  const Var kSliceTime     = SIMPLEVAR(hdr.unixtime);
  const Var kSliceMeanTime = SIMPLEVAR(hdr.subevtmeantime);

  //Make time axis for summer 2017 (nus17) period
  const HistAxis kTimeAxis("Unix Time (s)",Binning::Simple(9900, 1391000000, 1490000000), kSliceTime);
  const HistAxis kSliceTimeAxis("Slice Time (ns)",Binning::Simple(600, -50000, 555000), kSliceMeanTime);

  const HistAxis kNusAxisCVN("CVN NCID", Binning::Simple(21,0,1.05), kCVNnc);

  const std::vector<std::string> fnameunblind(MakeUnblindList());
  SpectrumLoader loader(fnameunblind);  
  loader.SetSpillCut(kStandardSpillCuts);

  // out of time spectra
  Spectrum spec_nus_cvn_outtime(loader, kNusAxisCVN,
                                kNus17FDAlt && kInTimingSideband,
                                kNoShift);
  
  Spectrum spec_nus_time_outtime(loader, kTimeAxis,
                                 kNus17FDAlt && kInTimingSideband,
                                 kNoShift);
  
  Spectrum spec_nus_slctime_outtime(loader, kSliceTimeAxis,
                                    kNus17FDAlt && kInTimingSideband,
                                    kNoShift);
  
  Spectrum spec_nus_vise_outtime(loader, kNus17Axis,
                                 kNus17FDAlt && kInTimingSideband,
                                 kNoShift);
    
  // in time spectra
  Spectrum spec_nus_cvn_intime(loader, kNusAxisCVN,
                               kNus17FDAlt && kInBeamSpill,
                               kNoShift);
  
  Spectrum spec_nus_time_intime(loader, kTimeAxis,
                                kNus17FDAlt && kInBeamSpill,
                                kNoShift);
  
  Spectrum spec_nus_slctime_intime(loader, kSliceTimeAxis,
                                kNus17FDAlt && kInBeamSpill,
                                kNoShift);
  
  Spectrum spec_nus_vise_intime(loader, kNus17Axis,
                                kNus17FDAlt && kInBeamSpill,
                                kNoShift);
  
  
  // let's go!
  loader.Go();
  
  std::cout << "Making output file, the money number is in here!" << std::endl;
  
  TFile fout("fout_nus17_box_opening_restricted.root", "RECREATE");

  spec_nus_cvn_intime.SaveTo(&fout, "spec_nus_cvn_numi");    
  spec_nus_cvn_outtime.SaveTo(&fout, "spec_nus_cvn_cosmic");    

  spec_nus_time_intime.SaveTo(&fout, "spec_nus_time_numi");    
  spec_nus_time_outtime.SaveTo(&fout, "spec_nus_time_cosmic");    

  spec_nus_slctime_intime.SaveTo(&fout, "spec_nus_slctime_numi");    
  spec_nus_slctime_outtime.SaveTo(&fout, "spec_nus_slctime_cosmic");    

  spec_nus_vise_intime.SaveTo(&fout, "spec_nus_vise_numi");    
  spec_nus_vise_outtime.SaveTo(&fout, "spec_nus_vise_cosmic");    


  Var kRun    = SIMPLEVAR(hdr.run);
  Var kSubrun = SIMPLEVAR(hdr.subrun);
  Var kEvt    = SIMPLEVAR(hdr.evt);
  Var kSlice  = SIMPLEVAR(hdr.subevt);
  std::vector<const Var*> vars = {&kRun, &kSubrun, &kEvt, &kSlice, 
                                  &kSliceMeanTime, &kCVNnc, &kNus17Energy};

  
  MakeTextListFile(fnameunblind,
                   {kInBeamSpill && kNus17FDAlt,
                       kInTimingSideband && kNus17FDAlt},
                   {"nus17_selected_numi.txt","nus17_selected_cosmic.txt"},
                   vars,
                   &kStandardSpillCuts );
  
}

std::vector<std::string> MakeUnblindList()
{
  std::vector<std::string> ret;
  const std::string pnfsdir = "/pnfs/nova/persistent/users/gsdavies/nus/prod3/concat/";

  std::vector<std::string> fnumi_period1 = Wildcard(
    pnfsdir + "prod_sumrestricteddecaf_R17-03-01-prod3reco.k_fd_numi_fhc_period1_v1_goodruns_nus2017/*.root"
  );
  std::vector<std::string> fnumi_period2 = Wildcard(
    pnfsdir + "prod_sumrestricteddecaf_R17-03-01-prod3reco.k_fd_numi_fhc_period2_v1_goodruns_nus2017/*.root"
  );
  std::vector<std::string> fnumi_period3 = Wildcard(
    pnfsdir + "prod_sumrestricteddecaf_R17-03-01-prod3reco.k_fd_numi_fhc_period3_v1_goodruns_nus2017/*.root"
  );
  std::vector<std::string> fnumi_period5 = Wildcard(
    pnfsdir + "prod_sumrestricteddecaf_R17-03-01-prod3reco.k_fd_numi_fhc_period5_v1_goodruns_nus2017/*.root"
  );
  
  for(const auto& file : fnumi_period1) {ret.push_back(file);}
  for(const auto& file : fnumi_period2) {ret.push_back(file);}
  for(const auto& file : fnumi_period3) {ret.push_back(file);}
  for(const auto& file : fnumi_period5) {ret.push_back(file);}
  return ret;
  
}