nus_ana01_box_opening.C

Go to the documentation of this file.

#include "CAFAna/Cuts/Cuts.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 "NuXAna/macros/NuSLoadMacroDefs.h"
#include "TH1D.h"

using namespace ana;

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

void nus_ana01_box_opening()
{

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

  //Make time axis for second analysis period
  const HistAxis kTimeAxis("unix time (s)",Binning::Simple(7200, 1391000000, 1463000000), 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,
                                kNusFD && kInTimingSideband,
                                kNoShift);
  
  Spectrum spec_nus_time_outtime(loader, kTimeAxis,
                                 kNusFD && kInTimingSideband,
                                 kNoShift);
  
  Spectrum spec_nus_slctime_outtime(loader, kSliceTimeAxis,
                                    kNusFD && kInTimingSideband,
                                    kNoShift);
  
  Spectrum spec_nus_cale_outtime(loader, kNCAxis,
                                 kNusFD && kInTimingSideband,
                                 kNoShift);
    
  // in time spectra
  Spectrum spec_nus_cvn_intime(loader, kNusAxisCVN,
                               kNusFD && kInBeamSpill,
                               kNoShift);
  
  Spectrum spec_nus_time_intime(loader, kTimeAxis,
                                kNusFD && kInBeamSpill,
                                kNoShift);
  
  Spectrum spec_nus_slctime_intime(loader, kSliceTimeAxis,
                                kNusFD && kInBeamSpill,
                                kNoShift);
  
  Spectrum spec_nus_cale_intime(loader, kNCAxis,
                                kNusFD && kInBeamSpill,
                                kNoShift);
  
  
  // let's go!
  loader.Go();
  
  std::cout<<"Making output file, the money number is in here! \n";
  
  TFile fout("fout_nus_ana01_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_cale_intime.SaveTo(&fout, "spec_nus_cale_numi");    
  spec_nus_cale_outtime.SaveTo(&fout, "spec_nus_cale_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, &kCaloE};

  
  MakeTextListFile(fnameunblind,
                   {kInBeamSpill && kNusFD,
                       kInTimingSideband && kNusFD},
                   {"nus_selected_numi.txt","nus_selected_cosmic.txt"},
                   vars,
                   &kStandardSpillCuts );
  
}

std::vector<std::string> MakeUnblindList()
{
  std::vector<std::string> fname1 = Wildcard(
    diskdirFD + "by_period/prod_restricteddecaf_S16-05-20_fd_numi_fhc_period1_nus_contain_v1_goodruns_prod2-snapshot/"
    + "prod_restricteddecaf_S16-05-20_fd_numi_fhc_period1_nus_contain_v1_goodruns_prod2-snapshot*.root"
  );
  std::vector<std::string> fname2 = Wildcard(
    diskdirFD + "by_period/prod_restricteddecaf_S16-05-20_fd_numi_fhc_period2_nus_contain_v1_goodruns_prod2-snapshot/"
    + "prod_restricteddecaf_S16-05-20_fd_numi_fhc_period2_nus_contain_v1_goodruns_prod2-snapshot*.root"
  );
  std::vector<std::string> fname3b = Wildcard(
    diskdirFD + "by_period/prod_restricteddecaf_S16-05-20_fd_numi_fhc_epoch3*_nus_contain_v1_goodruns_prod2-snapshot.root"
  );
  std::vector<std::string> fname3c = Wildcard(
    diskdirFD + "by_period/prod_restricteddecaf_S16-05-20_fd_numi_fhc_epoch3c_nus_contain_v1_goodruns_prod2-snapshot/"
    + "prod_restricteddecaf_S16-05-20_fd_numi_fhc_epoch3c_nus_contain_v1_goodruns_prod2-snapshot*.root"
  );
  std::vector<std::string> fname3d = Wildcard(
    diskdirFD + "by_period/prod_restricteddecaf_S16-05-20_fd_numi_fhc_epoch3d_nus_contain_v1_goodruns_prod2-snapshot/"
    + "prod_restricteddecaf_S16-05-20_fd_numi_fhc_epoch3d_nus_contain_v1_goodruns_prod2-snapshot*.root"
  );
  std::vector<std::string> ret;
  for(const auto& file : fname1)  { ret.push_back(file); }
  for(const auto& file : fname2)  { ret.push_back(file); }
  for(const auto& file : fname3b) { ret.push_back(file); }
  for(const auto& file : fname3c) { ret.push_back(file); }
  for(const auto& file : fname3d) { ret.push_back(file); }
  return ret;
}