make_quantiles_histogram_2020.C

Go to the documentation of this file.

#include "3FlavorAna/Cuts/NumuCuts2020.h"
#include "3FlavorAna/Cuts/QuantileCuts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Cuts/TruthCuts.h"

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

#include "3FlavorAna/Vars/HistAxes.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "CAFAna/Vars/XsecTunes.h"

#include "TCanvas.h"                                                                                                                                                      
#include "TH2.h"
#include "TFile.h"
#include <iostream>

using namespace ana;


void make_quantiles_histogram_2020( bool isFHC, std::string period   = "full" ) {
  
  std::string polarity = "";
  if ( isFHC) polarity = "fhc";
  if (!isFHC) polarity = "rhc";

  std::cout << "\n================================= \n"
            << " make distribution for quantiles"
            << "\n================================= \n"
            << " polarity: " << polarity << ", period: " << period
            << "\n================================= \n" 
            << std::endl;
  
  // Which cut am I using?
  const Cut kThisCut  = kNumu2020FD  && kIsNumuCC;

  // What XSec Weight am I using?
  const Var kXSecWeight = kXSecCVWgt2020;

  // Figure out my definition, and then declare my spectrum loader.
  std::string MyDef = "prod_caf_R19-11-18-prod5reco.f_fd_genie_N1810j0211a_nonswap_"+polarity+"_nova_v08_"+period+"_v1";
  SpectrumLoader loader( MyDef );
  loader.SetSpillCut(kStandardSpillCuts);
  
  // Declare my spectra
  Spectrum sQuantiles (loader, kNumuCCOptimisedAxis2020, kHadEFracAxis2020, kThisCut , kNoShift, kPPFXFluxCVWgt*kXSecWeight );
  
  // Set my loader going.
  loader.Go();
  
  // Call the MakeQuantileHistogram function.
  TH2* hQuantiles  = sQuantiles .ToTH2( sQuantiles .POT() );

  // Make my OldPresel Output file
  std::string Out_OldPresel = "quantiles_" + polarity + "_" + period + "_numu2020.root";
  TFile *OutFile = TFile::Open(Out_OldPresel.c_str(), "RECREATE");

  hQuantiles -> Write("FDSpec2D" );

  OutFile -> Close();

  new TCanvas;
  hQuantiles->Draw("colz");

}