Ana2017_box_opening_macro.C

Go to the documentation of this file.

#include "CAFAna/Cuts/Cuts.h"
#include "3FlavorAna/Cuts/NueCuts2017.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Cuts/TimingCuts.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 "CAFAna/Vars/GenieWeights.h"

#include "OscLib/OscCalcDumb.h"
#include "CAFAna/Prediction/PredictionNoExtrap.h"

#include "TH1D.h"

using namespace ana;


void Ana2017_box_opening_macro()
{
  //Make time axis for 2017 analysis period
  const Var kSliceTime = SIMPLEVAR(hdr.unixtime);
  const Var kSliceMeanTime = SIMPLEVAR(hdr.subevtmeantime);
  const HistAxis kTimeAxis("unix time (s)",Binning::Simple(10000, 1391000000, 1491000000), kSliceTime);
  const HistAxis kSliceTimeAxis("slice time (ns)",Binning::Simple(600, -50000, 550000), kSliceMeanTime);

  const HistAxis kNue2017EnergyAxis( "Nue Energy Estimator (GeV)",
                                     Binning::Simple( 10, 0, 5 ), kNueEnergy2017 );

  const std::string fnamecos  =
    // !!! Nue 2017 retricted dataset !!!
    "prod_sumrestricteddecaf_R17-03-01-prod3reco.k_fd_numi_fhc_full_v1_goodruns_nue2017";

  SpectrumLoader loadercos(fnamecos);
  
  loadercos.SetSpillCut(kStandardSpillCuts);

  // out of time spectra
  Spectrum spec_cvn_outtime(loadercos, kNue2017Axis,
                            kNue2017FDAllSamples && kInTimingSideband,
                            kNoShift);
  
  Spectrum spec_cvn_nueE_core_outtime(loadercos, kNue2017EnergyAxis,
                            kNue2017FD && kInTimingSideband,
                            kNoShift);

  Spectrum spec_cvn_nueE_peripheral_outtime(loadercos, kNue2017EnergyAxis,
                            kNue2017FDPeripheral && kInTimingSideband,
                            kNoShift);
  
  Spectrum spec_cvn_time_outtime(loadercos, kTimeAxis,
                            kNue2017FDAllSamples && kInTimingSideband,
                            kNoShift);
  
  Spectrum spec_cvn_slctime_outtime(loadercos, kSliceTimeAxis,
                            kNue2017FDAllSamples && kInTimingSideband,
                            kNoShift);

  // in time spectra
  Spectrum spec_cvn_intime(loadercos, kNue2017Axis,
                            kNue2017FDAllSamples && kInBeamSpill,
                            kNoShift);
  
  Spectrum spec_cvn_nueE_core_intime(loadercos, kNue2017EnergyAxis,
                            kNue2017FD && kInBeamSpill,
                            kNoShift);

  Spectrum spec_cvn_nueE_peripheral_intime(loadercos, kNue2017EnergyAxis,
                            kNue2017FDPeripheral && kInBeamSpill,
                            kNoShift);
  
  Spectrum spec_cvn_time_intime(loadercos, kTimeAxis,
                            kNue2017FDAllSamples && kInBeamSpill,
                            kNoShift);
  
  Spectrum spec_cvn_slctime_intime(loadercos, kSliceTimeAxis,
                            kNue2017FDAllSamples && kInBeamSpill,
                            kNoShift);

  // let's go!
  loadercos.Go();
  
  std::cout<<"Making output file ";
  TFile fout("Ana2017_numi_data.root", "RECREATE");
  std::cout<< fout.GetName() << std::endl;

  spec_cvn_intime.SaveTo(&fout, "spec_cvn_cvn2d_numi");    
  spec_cvn_outtime.SaveTo(&fout, "spec_cvn_cvn2d_cosmic");    

  spec_cvn_nueE_core_intime.SaveTo(&fout, "spec_cvn_nueE_core_numi");    
  spec_cvn_nueE_core_outtime.SaveTo(&fout, "spec_cvn_nueE_core_cosmic");    

  spec_cvn_nueE_peripheral_intime.SaveTo(&fout, "spec_cvn_nueE_peripheral_numi");    
  spec_cvn_nueE_peripheral_outtime.SaveTo(&fout, "spec_cvn_nueE_peripheral_cosmic");    

  spec_cvn_time_intime.SaveTo(&fout, "spec_cvn_time_numi");    
  spec_cvn_time_outtime.SaveTo(&fout, "spec_cvn_time_cosmic");    

  spec_cvn_slctime_intime.SaveTo(&fout, "spec_cvn_slctime_numi");    
  spec_cvn_slctime_outtime.SaveTo(&fout, "spec_cvn_slctime_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,
                                  &kNue2017AnaBin,
                                  &kCVNe, 
                                  &kNueEnergy2017                               
  };

  MakeTextListFile(fnamecos,
                   {kInBeamSpill && kNue2017FDAllSamples,
                       kInTimingSideband && kNue2017FDAllSamples},
                   {"nue2017_selected_numi.txt",
                       "nue2017_selected_cosmic.txt"},
                   vars,
                   &kStandardSpillCuts );

}