rock.C

Go to the documentation of this file.

#include "OscLib/OscCalcDumb.h"

#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"

#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "3FlavorAna/Cuts/NueCuts2017.h"
#include "CAFAna/Cuts/TimingCuts.h"

#include "CAFAna/Vars/Vars.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "CAFAna/Vars/GenieWeights.h"

#include "CAFAna/Prediction/PredictionNoExtrap.h"

#include "3FlavorAna/Ana2017/nue/sidebands/vars.h"
//#include "./vars.h"

#include "TFile.h"

#include <iostream>
#include <cmath>

using namespace ana;


void rock(std::string suffix = "full")
{

//  SpectrumLoader lNonS("prod_caf_R17-03-01-prod3reco.m_fd_rock+cosmic-overlay_genie_nonswap_fhc_nova_v08_full_v2");
//  SpectrumLoader lSwap("prod_caf_R17-03-01-prod3reco.m_fd_rock+cosmic-overlay_genie_fluxswap_fhc_nova_v08_full_v2");

  SpectrumLoader lNonS("/pnfs/nova/scratch/users/bnayak/nue2017_fd_mc/prod_sumdecaf_R17-03-01-prod3reco.m_fd_rock+cosmic-overlay_genie_nonswap_fhc_nova_v08_full_v2_nue2017/*.root");
  SpectrumLoader lSwap("/pnfs/nova/scratch/users/bnayak/nue2017_fd_mc/prod_sumdecaf_R17-03-01-prod3reco.m_fd_rock+cosmic-overlay_genie_fluxswap_fhc_nova_v08_full_v2_nue2017/*.root");

  lNonS.SetSpillCut(kStandardSpillCuts);
  lSwap.SetSpillCut(kStandardSpillCuts);

  PredictionNoExtrap* pred[kNumSels][kNumVars];

  cout << "ATTENTION! THERE IS A DECAF CUT BELOW, BE CAREFUL!" << endl;

  for(int selIdx = 0; selIdx < kNumSels; ++selIdx){
     for (int varIdx=0; varIdx<kNumVars;varIdx++){
         const HistAxis& axis = defs[varIdx].axis;
         pred[selIdx][varIdx] = new PredictionNoExtrap(lNonS, lSwap, 
                                                       axis.GetLabels()[0], axis.GetBinnings()[0], axis.GetVars()[0], 
                                                       sels[selIdx] && kInBeamSpill && kNueFD2017DecafCut, //don't forget about decaf cut here
                                                       kNoShift, 
                                                       kXSecCVWgt2017 * kPPFXFluxCVWgt); 

     }
  }

  lNonS.Go();
  lSwap.Go();

  TFile* file = new TFile(("rock_bkg_"+suffix+".root").c_str(),"RECREATE");
  for(int selIdx = 0; selIdx < kNumSels; ++selIdx){
     TDirectory* d = file->mkdir(selNames[selIdx].c_str());
     for (int varIdx=0; varIdx<kNumVars;varIdx++){
         const char* name = defs[varIdx].name.c_str();
         pred[selIdx][varIdx]->SaveTo(d, TString::Format("rock_%s", name));
     }
  }

}