CutFlow_MC.C File Reference

#include "CAFAna/Analysis/Prod4Loaders.h"
#include "CAFAna/Cuts/AnalysisMasks.h"
#include "CAFAna/Cuts/Cuts.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "3FlavorAna/Cuts/NumuCuts2018.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Cuts/TimingCuts.h"
#include "CAFAna/Core/Loaders.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Vars/Vars.h"
#include "3FlavorAna/Vars/Binnings.h"
#include "3FlavorAna/Vars/NumuVars.h"
#include "CAFAna/Vars/GenieWeights.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "CAFAna/Extrap/ModularExtrap.h"
#include "CAFAna/Prediction/PredictionNoExtrap.h"
#include <fstream>

Go to the source code of this file.

Functions
void	CutFlow_MC ()

Function Documentation

void CutFlow_MC

(

)

Definition at line 28 of file CutFlow_MC.C.

References cd(), caf::Proxy< caf::SRCosRej >::cdirscore, caf::Proxy< caf::SRIDBranch >::cosrej, cut, CutNames, ana::Loaders::Go(), ana::kBoxCut, ana::kCCE, ana::Loaders::kFHC, ana::kFullCAF, ana::kNoCut, ana::kNoShift, ana::kNumuCCEOptimisedBinning, ana::kNumuContainFD2017, ana::kNumuCosmicRej2018, ana::kNumuPID2018, ana::kNumuQuality, ana::kPPFXFluxCVWgt, ana::Loaders::kRHC, ana::kStandardSpillCuts, ana::kXSecCVWgt2018, NCuts, OutFile, ana::RatioBins, ana::SaveTo(), caf::Proxy< caf::StandardRecord >::sel, ana::Loaders::SetSpillCut(), ana::Binning::Simple(), sr, and string.

{
  // --- Set my loader
  Loaders *FHC_Loader = new Prod4NomLoaders( ECAFType::kFullCAF, Loaders::kFHC, "full");
  Loaders *RHC_Loader = new Prod4NomLoaders( ECAFType::kFullCAF, Loaders::kRHC, "full");
  FHC_Loader->SetSpillCut(kStandardSpillCuts);
  RHC_Loader->SetSpillCut(kStandardSpillCuts);
  
  // --- Define my cuts and their names.
  std::vector<Cut> TieredCuts;
  std::vector<std::string> CutNames;
  // --- A set without true containment.
  CutNames.emplace_back("NoCut");       TieredCuts.emplace_back(kNoCut);
  CutNames.emplace_back("Quality");     TieredCuts.emplace_back(kNumuQuality);
  CutNames.emplace_back("Containment"); TieredCuts.emplace_back(kNumuQuality && kNumuContainFD2017);
  CutNames.emplace_back("Cosmic_Rej");  TieredCuts.emplace_back(kNumuQuality && kNumuContainFD2017 && kNumuCosmicRej2018);
  CutNames.emplace_back("Particle_ID"); TieredCuts.emplace_back(kNumuQuality && kNumuContainFD2017 && kNumuCosmicRej2018 && kNumuPID2018 );
  // --- A set with true containment.
  CutNames.emplace_back("True-NoCut");       TieredCuts.emplace_back(kBoxCut && kNoCut);
  CutNames.emplace_back("True-Quality");     TieredCuts.emplace_back(kBoxCut && kNumuQuality);
  CutNames.emplace_back("True-Containment"); TieredCuts.emplace_back(kBoxCut && kNumuQuality && kNumuContainFD2017);
  CutNames.emplace_back("True-Cosmic_Rej");  TieredCuts.emplace_back(kBoxCut && kNumuQuality && kNumuContainFD2017 && kNumuCosmicRej2018);
  CutNames.emplace_back("True-Particle_ID"); TieredCuts.emplace_back(kBoxCut && kNumuQuality && kNumuContainFD2017 && kNumuCosmicRej2018 && kNumuPID2018 );

  size_t NCuts = TieredCuts.size();

  const Binning RatioBins  = Binning::Simple(100, -0.05, 0.05 );
  const Var kCDirScore([](const caf::SRProxy* sr){
      return ((float)sr->sel.cosrej.cdirscore);
    });
  
  // --- Define my objects.
  PredictionNoExtrap* ReconEnergy_FHC[NCuts];
  PredictionNoExtrap* ReconEnergy_RHC[NCuts];

  PredictionNoExtrap* CDir_FHC[NCuts];
  PredictionNoExtrap* CDir_RHC[NCuts];
  
  // --- Set my objects.
  for(size_t cut = 0; cut < NCuts; ++cut) {
    const Cut ThisCut = TieredCuts[cut];    
    ReconEnergy_FHC[cut] = new PredictionNoExtrap( *FHC_Loader, "Reco Energy (GeV)", kNumuCCEOptimisedBinning, kCCE, ThisCut, kNoShift, kPPFXFluxCVWgt*kXSecCVWgt2018 );
    ReconEnergy_RHC[cut] = new PredictionNoExtrap( *RHC_Loader, "Reco Energy (GeV)", kNumuCCEOptimisedBinning, kCCE, ThisCut, kNoShift, kPPFXFluxCVWgt*kXSecCVWgt2018 );

    CDir_FHC[cut] = new PredictionNoExtrap( *FHC_Loader, "CosRej CDir", RatioBins, kCDirScore, ThisCut, kNoShift, kPPFXFluxCVWgt*kXSecCVWgt2018 );
    CDir_RHC[cut] = new PredictionNoExtrap( *RHC_Loader, "CosRej CDir", RatioBins, kCDirScore, ThisCut, kNoShift, kPPFXFluxCVWgt*kXSecCVWgt2018 );
  }

  // --- Set the loader off!
  FHC_Loader->Go();
  RHC_Loader->Go();

  // --- Define an outfile.
  std::string OutName = "MonteCarloCutFlow.root";
  TFile *OutFile = new TFile(OutName.c_str(), "RECREATE");
  OutFile  -> cd();
  for(size_t cut = 0; cut < NCuts; ++cut){
    ReconEnergy_FHC[cut] -> SaveTo( OutFile,  TString(CutNames[cut])+TString("_FullMC_FHC"  ) ) ;
    ReconEnergy_RHC[cut] -> SaveTo( OutFile,  TString(CutNames[cut])+TString("_FullMC_RHC"  ) ) ;

    CDir_FHC[cut] -> SaveTo( OutFile,  TString(CutNames[cut])+TString("_FullMC_CDir_FHC"  ) ) ;
    CDir_RHC[cut] -> SaveTo( OutFile,  TString(CutNames[cut])+TString("_FullMC_CDir_RHC"  ) ) ;
  }
} // End of function