savePrediction_addExpt.C File Reference

#include "../WorkshopIncludes.h"

Go to the source code of this file.

Functions
void	savePrediction_addExpt ()

Function Documentation

void savePrediction_addExpt

(

)

Definition at line 9 of file savePrediction_addExpt.C.

References allInOneTrainingPlots::axis, calc, om::cout, e, allTimeWatchdog::endl, MakeMiniprodValidationCuts::f, ana::Loaders::Go(), ana::kBeam, ana::kCCE, cmf::kData, caf::kFARDET, kMC, caf::kNEARDET, ana::kNoShift, ana::kNumuContainFD, ana::kNumuCosmicRej, ana::kNumuEnergyBinning, ana::kNumuQuality, ana::kStandardSpillCuts, cmf::kTauSwap, ana::kTuftsWeightCC, loaders, caf::SRRemid::pid, caf::SRIDBranch::remid, ana::ResetOscCalcToDefault(), ana::PredictionNoExtrap::SaveTo(), caf::StandardRecord::sel, ana::Loaders::SetLoaderPath(), ana::Loaders::SetSpillCut(), std::sqrt(), sr, and string.

{
  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_addExpt_kCCE.root","RECREATE"); // different name

  const Cut kUglyRemID(
    {"sel.remid.pid"},
    [](const caf::StandardRecord* sr)
      { return (sr->sel.remid.pid > 0.25 && sr->sel.remid.pid < 0.75 );
      }
  );

  const Cut kBadRemID(
    {"sel.remid.pid"},
    [](const caf::StandardRecord* sr)
      { return (sr->sel.remid.pid >= 0.75 && sr->sel.remid.pid < 1.0 );
      }
  );

  const Cut kGoodRemID(
    {"sel.remid.pid"},
    [](const caf::StandardRecord* sr)
      { return (sr->sel.remid.pid == 1.0 );
      }
  );

  const Cut kUglyFD = kNumuQuality && kNumuContainFD && kUglyRemID && kNumuCosmicRej;
  const Cut kBadFD  = kNumuQuality && kNumuContainFD && kBadRemID  && kNumuCosmicRej;
  const Cut kGoodFD = kNumuQuality && kNumuContainFD && kGoodRemID && kNumuCosmicRej;

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

  PredictionNoExtrap* prediction_ugly = new PredictionNoExtrap(loaders, axis, kUglyFD, kNoShift, kTuftsWeightCC);
  PredictionNoExtrap* prediction_bad  = new PredictionNoExtrap(loaders, axis, kBadFD,  kNoShift, kTuftsWeightCC);
  PredictionNoExtrap* prediction_good = new PredictionNoExtrap(loaders, axis, kGoodFD, kNoShift, kTuftsWeightCC);

  // GO!
  loaders.Go();

  prediction_ugly->SaveTo(f->mkdir( "pred_ugly" ));
  prediction_bad->SaveTo(f->mkdir( "pred_bad" ));
  prediction_good->SaveTo(f->mkdir( "pred_good" ));

  f->Close();

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

} // End of function