7 #include "../WorkshopIncludes.h" 14 std::string directoryD =
"/pnfs/nova/persistent/production/concat/R16-03-03-prod2reco.d/";
15 std::string directoryF =
"/pnfs/nova/persistent/production/concat/R16-03-03-prod2reco.f/";
17 std::string nd_data = directoryD+
"prod_decaf_R16-03-03-prod2reco.d_nd_numi_fhc_full_numu_contain_v1_goodruns_prod2-snapshot.root";
21 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";
23 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";
24 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";
25 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";
36 osc::OscCalculatorPMNSOpt
calc;
39 calc.SetTh23(TMath::ASin(
sqrt(0.403)));
40 calc.SetDmsq32(2.67
e-3);
42 TFile*
f =
new TFile(
"../Predictions/Prediction_epoch3b_simple_kCCE.root",
"RECREATE");
51 prediction->
SaveTo(f->mkdir(
"pred_3b" ));
Near Detector underground.
Far Detector at Ash River.
void ResetOscCalcToDefault(osc::IOscCalcAdjustable *calc)
Reset calculator to default assumptions for all parameters.
Collection of SpectrumLoaders for many configurations.
void savePrediction_simple()
void Go()
Call Go() on all the loaders.
void SetSpillCut(const SpillCut &cut)
const Binning kNumuEnergyBinning
const SystShifts kNoShift
const SpillCut kStandardSpillCuts
Apply this unless you're doing something special.
void SetLoaderPath(const std::string &path, caf::Det_t det, DataMC datamc, DataSource src=kBeam, SwappingConfig swap=kNonSwap)
Configure loader via wildcard path.