test_numu_summed_genie_systs.C

Go to the documentation of this file.

// Fills spectra to test the summed small genie systs 

#include "CAFAna/Cuts/Cuts.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "CAFAna/Prediction/PredictionGenerator.h"
#include "CAFAna/Core/Loaders.h"
#include "CAFAna/Core/Binning.h"
#include "CAFAna/Analysis/Calcs.h"
#include "CAFAna/Prediction/PredictionInterp.h"
#include "CAFAna/Systs/Systs.h"
#include "3FlavorAna/Systs/NumuSysts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Analysis/Plots.h"
#include "3FlavorAna/Vars/Binnings.h"


#include "OscLib/IOscCalc.h"
#include "OscLib/OscCalcGeneral.h"
#include "OscLib/OscCalcPMNS.h"
#include "OscLib/OscCalc.h"


#include "TCanvas.h"
#include "TFile.h"
#include "TGraph.h"
#include "TH1.h"
#include "TH2.h"
#include "TMath.h"
#include "TGaxis.h"
#include "TMultiGraph.h"
#include "TLegend.h"
#include "TLatex.h"
#include "TStyle.h"
#include "THStack.h"
#include "TPaveText.h"

#include <cmath>
#include <iostream>
#include <vector>
#include <sstream>
#include <string>
#include <fstream>
#include <iomanip>



using namespace ana;







void test_numu_summed_genie_systs()
{


  Loaders loaders;
  std::string fnameNDMC("defname: ");
  fnameNDMC+="rocco_decaf_S15-05-22_nd_genie_fhc_nonswap_genierw_numu_contain";
  fnameNDMC+=" with limit 200";
  

  loaders.SetLoaderPath(fnameNDMC, caf::kFARDET, Loaders::kMC, 
                              ana::kBeam, Loaders::kNonSwap);

  
  loaders.SetSpillCut(kStandardSpillCuts);

 
  const HistAxis axisCC("Reconstructed Neutrino Energy (GeV)",
                                     kNumuEnergyBinning, kCCE);


  NoExtrapGenerator gen(axisCC, kNumuND);


  std::vector<const ISyst*> systs;
  systs.push_back(&kNumuSummedSmallGENIESyst);
  //systs.push_back(new GenieRwgtTableSyst( rwgt::fReweightMaCCQE));
  osc::NoOscillations calc;

  //Predictions with extrap
  PredictionInterp pred(systs,
        &calc ,
        gen,
        loaders);



  loaders.Go();

  PlotWithSystErrorBand(&pred,systs,&calc, 1e18);


  
}