make_ehade_histogram.C

Go to the documentation of this file.

#include "3FlavorAna/Ana2018/numu/includes.h"
//#include "/nova/app/users/dmendez/Releases/dev/CAFAna/numu/Analysis2018/varsandcuts.h"
#include "3FlavorAna/Ana2018/numu/varsandcuts.h"

using namespace ana;


void make_ehade_histogram(std::string polarity = "rhc",
                          std::string period   = "full")
{
  
  std::cout << "\n=================================" << std::endl;
  std::cout << " make distribution for quantiles   " << std::endl;
  std::cout << "=================================\n" << std::endl;

  std::string production = "sumdecaf_R17-11-14-prod4reco";
  if(polarity == "fhc") production += ".d";
  if(polarity == "rhc") production += ".e";


  // define selection
  const Cut kNumuNDCut = kNumuCutND2018;
  const Cut kNumuFDCut = kNumuCutFD2018;
  const Var kXSecWeight = kXSecCVWgt2018Fix;

  std::cout << "polarity: " << polarity << ", period: " << period << ", files: " << production << std::endl;

  std::string OutDir = "";
  // Get a vector of file names to pass to the spectrum loader
  std::string fd_genie_nonswap  =  "prod_" + production + "_fd_genie_nonswap_" + polarity + "_nova_v08_" + period + "_v1_numu2018";
  SpectrumLoader loader(fd_genie_nonswap);
  TH2* hQuantiles = MakeQuantileHistogram(loader, kNumuCCOptimisedAxis, kHadEFracAxis, kNumuFDCut && kIsNumuCC, kStandardSpillCuts, kNoShift, kPPFXFluxCVWgt*kXSecWeight);
  std::string OutFileName = OutDir + "quantiles__" + polarity + "_" + period + "__numu2018.root";
  std::cerr << "Saving distribution for quantiles to " << OutFileName << std::endl;
  TFile *f = TFile::Open(OutFileName.c_str(), "RECREATE");
  hQuantiles->Write("FDSpec2D");
  f->Close();
    
}