make_pid.C File Reference

#include "TCanvas.h"
#include "TH1.h"
#include "TF1.h"
#include "TFile.h"
#include "TLegend.h"
#include "TLegendEntry.h"
#include "TArrayD.h"
#include "TAxis.h"
#include "TLatex.h"
#include "CAFAna/Core/Cut.h"
#include "CAFAna/Cuts/NueCutsFirstAna.h"
#include "CAFAna/Analysis/Plots.h"
#include "CAFAna/Analysis/Style.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Vars/Vars.h"
#include "CAFAna/Vars/GenieWeights.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "StandardRecord/Proxy/SRProxy.h"
#include "CAFAna/reco/prongID/make_pid.h"

Go to the source code of this file.

Functions
void	make_pid (bool isFHC=1)

Function Documentation

void make_pid

(

bool

isFHC = 1

)

Definition at line 35 of file make_pid.C.

References ana::bins, om::cout, cut, allTimeWatchdog::endl, submit_syst::fout, ana::SpectrumLoaderBase::Go(), MECModelEnuComparisons::i, isFHC, calib::j, ana::kNoShift, kNPlots, kNSels, ana::kPPFXFluxCVWgt, ana::kStandardDQCuts, ana::kTightBeamQualityCuts, ana::kXSecCVWgt2017, label, demo0::loader, plots, ana::Spectrum::SaveTo(), ana::sels, ana::SpectrumLoaderBase::SetSpillCut(), string, and PandAna.Demos.tute_pid_validation::var.

                             {
  const std::string fData = isFHC ? 
    "prod_decaf_R17-03-01-prod3reco.d_nd_numi_fhc_full_nue_or_numu_or_nus_contain_v1_goodruns" :
    "prod_decaf_R17-03-01-prod3reco.l_nd_numi_rhc_full_nue_or_numu_or_nus_contain_v1_goodruns" ;

  const std::string fMC = isFHC ?
    "prod_decaf_R17-03-01-prod3reco.d_nd_genie_nonswap_fhc_nova_v08_full_nue_or_numu_or_nus_contain_v1" :
    "prod_decaf_R17-03-01-prod3reco.l_nd_genie_nonswap_rhc_nova_v08_full_nue_or_numu_or_nus_contain_v1" ;

  SpectrumLoaderBase* loader = new SpectrumLoader(fData);
  SpectrumLoaderBase* loaderMC = new SpectrumLoader(fMC);

  std::cout << "MC:   " << fMC << std::endl;
  std::cout << "Data: " << fData << std::endl;

  loader->SetSpillCut(kTightBeamQualityCuts && kStandardDQCuts);
  loaderMC->SetSpillCut(kStandardDQCuts);

  // One spectrum per selection, per Var (see .h)

  Spectrum* spect[kNSels][kNPlots];
  
  for(unsigned int i = 0;i < kNSels;++i){
    for(unsigned int j = 0;j < kNPlots;++j){
      if(j < 5) 
        spect[i][j] = new Spectrum(plots[j].label,plots[j].bins,
                                   *loader,plots[j].var,sels[i].cut);
      else      
        spect[i][j] = new Spectrum(plots[j].label,plots[j].bins,
                                   *loaderMC,plots[j].var,sels[i].cut,
                                   kNoShift,kXSecCVWgt2017*kPPFXFluxCVWgt);
    }
  }

  // Go, loaders, go!

  loader->Go();
  loaderMC->Go();
  

  // Save output

  std::string cStr = isFHC ? "FHC" : "RHC";
  TFile fout(("make_pid_output_"+cStr+".root").c_str(),"RECREATE");

  for(unsigned int i = 0;i < kNSels;++i){
    for(unsigned int j = 0;j < kNPlots;++j){
      spect[i][j]->SaveTo(&fout, sels[i].suffix+"_"+plots[j].suffix);
    }
  }


}