savePrediction_simple.C

Go to the documentation of this file.

#ifdef __CINT__
void savePrediction_simple(){
  std::cout << "Sorry, you must run in compiled mode" << std::endl;
}
#else

#include "../WorkshopIncludes.h"

void savePrediction_simple()
{
  std::cout << "Drawing plots..." << std::endl;

  // Production Concats
  std::string directoryD   = "/pnfs/nova/persistent/production/concat/R16-03-03-prod2reco.d/";
  std::string directoryF   = "/pnfs/nova/persistent/production/concat/R16-03-03-prod2reco.f/";
  // Data
  std::string nd_data      = directoryD+"prod_decaf_R16-03-03-prod2reco.d_nd_numi_fhc_full_numu_contain_v1_goodruns_prod2-snapshot.root";

  // Baseline MC
  // Only one of the files. In true life, use all MC!
  std::string nd_nonswap_few  = directoryD+"prod_decaf_R16-03-03-prod2reco.d_nd_genie_nonswap_genierw_fhc_nova_v08_epoch1-3c_numu_contain_v1_prod2-snapshot/*_1_of_16.root";
  // Epoch 3b only for FD to run quickly. In practice, use a PredictionCombinePeriods 
  std::string fd_nonswap_3b  = directoryF+"/epoch3b/prod_decaf_R16-03-03-prod2reco.f_fd_genie_nonswap_fhc_nova_v08_epoch3b_numu_contain_v1_prod2-snapshot.root";
  std::string fd_fluxswap_3b = directoryF+"/epoch3b/prod_decaf_R16-03-03-prod2reco.f_fd_genie_fluxswap_fhc_nova_v08_epoch3b_numu_contain_v1_prod2-snapshot.root";
  std::string fd_tau_3b      = directoryF+"/epoch3b/prod_decaf_R16-03-03-prod2reco.f_fd_genie_tau_fhc_nova_v08_epoch3b_numu_contain_v1_prod2-snapshot.root";

  Loaders loaders;
  loaders.SetLoaderPath( nd_data,        caf::kNEARDET, Loaders::kData, ana::kBeam, Loaders::kNonSwap );
  loaders.SetLoaderPath( nd_nonswap_few, caf::kNEARDET, Loaders::kMC,   ana::kBeam, Loaders::kNonSwap );
  loaders.SetLoaderPath( fd_nonswap_3b,  caf::kFARDET,  Loaders::kMC,   ana::kBeam, Loaders::kNonSwap );
  loaders.SetLoaderPath( fd_fluxswap_3b, caf::kFARDET,  Loaders::kMC,   ana::kBeam, Loaders::kFluxSwap);
  loaders.SetLoaderPath( fd_tau_3b,      caf::kFARDET,  Loaders::kMC,   ana::kBeam, Loaders::kTauSwap );

  loaders.SetSpillCut(kStandardSpillCuts);

  osc::OscCalculatorPMNSOpt calc;
  ana::ResetOscCalcToDefault(&calc); // This seeds to default (in particular max-mixing)

  calc.SetTh23(TMath::ASin(sqrt(0.403))); // Second analysis best fit                             
  calc.SetDmsq32(2.67e-3); // Same

  TFile* f = new TFile("../Predictions/Prediction_epoch3b_simple_kCCE.root","RECREATE"); 

  HistAxis axis( ";Reconstructed neutrino energy (GeV);Events / (0.25 GeV)", kNumuEnergyBinning, kCCE );

  PredictionNoExtrap* prediction = new PredictionNoExtrap(loaders, axis, kNumuFD, kNoShift, kTuftsWeightCC);

  // GO!
  loaders.Go();

  prediction->SaveTo(f->mkdir( "pred_3b" ));

  f->Close();

  std::cout << "Prediction saved!" << std::endl;

} // End of function

#endif