makePred.C File Reference

#include "CAFAna/Core/Loaders.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Cuts/TimingCuts.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "3FlavorAna/Cuts/NumuCuts2017.h"
#include "3FlavorAna/Cuts/QuantileCuts.h"
#include "3FlavorAna/Vars/HistAxes.h"
#include "CAFAna/Vars/Vars.h"
#include "3FlavorAna/Vars/NumuVars.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "3FlavorAna/Systs/NumuSysts.h"
#include "CAFAna/Analysis/Calcs.h"
#include "OscLib/OscCalcPMNSOpt.h"
#include "CAFAna/Prediction/PredictionExtrap.h"
#include "3FlavorAna/Prediction/PredictionGenerator3Flavor.h"
#include "CAFAna/Prediction/PredictionInterp.h"
#include "3FlavorAna/Prediction/PredictionSystNumu2017.h"
#include "TFile.h"

Go to the source code of this file.

Functions
void	setLoadersPathNumu2017 (Loaders &loaders, bool fakenddata)
void	makePred (bool usesysts=true, bool useFakeNDData=false)

Function Documentation

void makePred	(	bool	usesysts = true,
		bool	useFakeNDData = false
	)

Definition at line 64 of file makePred.C.

References std::asin(), om::cerr, om::cout, e, allTimeWatchdog::endl, MakeMiniprodValidationCuts::f, ana::getAllDirectNumuSysts2017(), ana::Loaders::Go(), ana::kHadEFracAxis, ana::kNoShift, ana::kNumuCCOptimisedAxis, ana::kNumuCutFD2017, ana::kNumuCutND2017, ana::kPPFXFluxCVWgt, ana::kStandardSpillCuts, ana::kXSecCVWgt2018, loaders, ana::pnfs2xrootd(), ana::QuantileCutsFromTH2(), ana::ResetOscCalcToDefault(), ana::PredictionInterp::SaveTo(), ana::SaveTo(), osc::_OscCalcPMNSOpt< T >::SetDmsq32(), setLoadersPathNumu2017(), ana::Loaders::SetSpillCut(), osc::_OscCalcPMNSOpt< T >::SetTh23(), std::sqrt(), string, and systs.

Referenced by getPredictions(), and getTimePeakPlots().

                                                               {                        
  
  std::cout << "\n use full numu systs?: " << usesysts
           << "\n " <<std::endl;

  // custom osc. parameters (non-max mixing) Get the numu result parameters!
  osc::OscCalcPMNSOpt calc_nonMax;
  ana::ResetOscCalcToDefault(&calc_nonMax);
  calc_nonMax.SetDmsq32(2.67e-3);
  calc_nonMax.SetTh23(asin(sqrt(0.402)));  // PRL-published non-max mixing

  std::vector<const ISyst*> systs = {};
  if(usesysts) systs = getAllDirectNumuSysts2017();

  // Get FD files for HadEFrac vs. neutrino energy distribution
  std::string quantpath = "/pnfs/nova/persistent/analysis/numu/Ana2017/all_FD_histo_for_quantile_cuts.root";
  TFile* quantfile = TFile::Open(pnfs2xrootd(quantpath).c_str());
  TH2* quanthist = (TH2*)quantfile->Get("dir_FDSpec2D/FDSpec2D");
  std::vector<Cut> QuantCuts = QuantileCutsFromTH2(quanthist, kNumuCCOptimisedAxis, kHadEFracAxis, 4);

  Loaders loaders;
  loaders.SetSpillCut(kStandardSpillCuts);
  setLoadersPathNumu2017(loaders, useFakeNDData);

  //NumuExtrapGenerator NumuCCExtrapAllQuants(kNumuCCOptimisedAxis, kNumuCutFD2017, kNumuCutND2017, kNoShift, kPPFXFluxCVWgt * kXSecCVWgt2017);
  NumuExtrapGenerator NumuCCExtrapAllQuants(kNumuCCOptimisedAxis, kNumuCutFD2017, kNumuCutND2017, kNoShift, kPPFXFluxCVWgt * kXSecCVWgt2018);
  PredictionInterp predictionAllQuants(systs, &calc_nonMax, NumuCCExtrapAllQuants, loaders);
  
  // make a vector with an entry for each quantiles
  std::vector<NumuExtrapGenerator*> NumuCCExtrapVec;
  std::vector<PredictionInterp*> predictionVec;

  for(auto thisCut : QuantCuts){
    //NumuCCExtrapVec.push_back(new NumuExtrapGenerator(kNumuCCOptimisedAxis, thisCut && kNumuCutFD2017, thisCut && kNumuCutND2017, kNoShift, kPPFXFluxCVWgt * kXSecCVWgt2017));
    NumuCCExtrapVec.push_back(new NumuExtrapGenerator(kNumuCCOptimisedAxis, thisCut && kNumuCutFD2017, thisCut && kNumuCutND2017, kNoShift, kPPFXFluxCVWgt * kXSecCVWgt2018));
    predictionVec.push_back(new PredictionInterp(systs, &calc_nonMax, *(NumuCCExtrapVec.back()), loaders));

  } // loop over quantiles

  ////////////////////
  // GOOOOOOOO!!!!!
  loaders.Go();
  ////////////////////

  // Now we should save the predictions
  std::cerr << "Saving Predictions" << std::endl;
  std::string ndDataType = "";
  if(useFakeNDData) ndDataType = "fakeNDData_";
  std::string OutFileName = "pred_3A_prod3_SAor2017decafs_2018Wgt_AllQuants_" + ndDataType + ".root";
  TFile *fAQ = TFile::Open(OutFileName.c_str(), "RECREATE");
  predictionAllQuants.SaveTo(fAQ, "prediction");
  fAQ->Close();
  delete fAQ;

  for(size_t quant=1; quant <= predictionVec.size();quant++){
    auto prediction = predictionVec[quant-1];
    OutFileName = "pred_3A_prod3_SAor2017decafs_2018Wgt_" + (std::string)Form("Quant%d", quant) + "_" + ndDataType + ".root";
    std::cerr << "Saving Quantile " << quant << " to Prediction file: " << OutFileName << std::endl;
    TFile *f = TFile::Open(OutFileName.c_str(), "RECREATE");
    prediction -> SaveTo(f, "prediction");
    f->Close();
    delete f;
  } 

} 

void setLoadersPathNumu2017	(	Loaders &	loaders,
		bool	fakenddata
	)

Definition at line 33 of file makePred.C.

References ana::kBeam, ana::Loaders::kData, caf::kFARDET, ana::Loaders::kFluxSwap, ana::Loaders::kMC, caf::kNEARDET, ana::Loaders::kNonSwap, ana::Loaders::kTauSwap, ana::Loaders::SetLoaderPath(), and string.

Referenced by makePred().

                                                              {

  // std::string nd_data = "prod_sumdecaf_R17-03-01-prod3reco.d_nd_numi_fhc_full_v1_goodruns_numu2017";
  // std::string nd_nonswap = "prod_sumdecaf_R17-03-01-prod3reco.d_nd_genie_nonswap_fhc_nova_v08_full_v1_numu2017";
  // std::string fd_nonswap  = "prod_sumdecaf_R17-03-01-prod3reco.l_fd_genie_nonswap_fhc_nova_v08_full_v1_numu2017";
  std::string nd_data    = "prod_sumdecaf_R17-03-01-prod3reco.d_nd_numi_fhc_full_v1_goodruns_numu_sa_2017";
  std::string nd_nonswap = "prod_sumdecaf_R17-03-01-prod3reco.d_nd_genie_nonswap_fhc_nova_v08_full_v1_numu_sa_2017";
  std::string fd_nonswap = "prod_sumdecaf_R17-03-01-prod3reco.l_fd_genie_nonswap_fhc_nova_v08_full_v1_numu_sa_2017";

  std::string fd_fluxswap = "prod_sumdecaf_R17-03-01-prod3reco.l_fd_genie_fluxswap_fhc_nova_v08_full_v1_numu2017";
  std::string fd_tau      = "prod_sumdecaf_R17-03-01-prod3reco.l_fd_genie_tau_fhc_nova_v08_full_v1_numu2017";

  //////////////////// set loaders ////////////////////
  if(fakenddata) loaders.SetLoaderPath( nd_nonswap, caf::kNEARDET,   
                                        Loaders::kData, ana::kBeam, 
                                        Loaders::kNonSwap );
  if(!fakenddata) loaders.SetLoaderPath( nd_data, caf::kNEARDET,   
                                        Loaders::kData, ana::kBeam, 
                                        Loaders::kNonSwap );

  loaders.SetLoaderPath( nd_nonswap,  caf::kNEARDET,    Loaders::kMC,   ana::kBeam,
                         Loaders::kNonSwap );
  loaders.SetLoaderPath( fd_nonswap,  caf::kFARDET,     Loaders::kMC,   ana::kBeam,
                         Loaders::kNonSwap );
  loaders.SetLoaderPath( fd_fluxswap, caf::kFARDET,     Loaders::kMC,   ana::kBeam,
                         Loaders::kFluxSwap);
  loaders.SetLoaderPath( fd_tau,      caf::kFARDET,     Loaders::kMC,   ana::kBeam,
                         Loaders::kTauSwap );

}