futureSig_getNueNumuSyst.C
Go to the documentation of this file.
1 // make three predictions, one for each period
2 
3 #include "CAFAna/Cuts/Cuts.h"
4 #include "3FlavorAna/Cuts/NumuCuts.h"
5 #include "CAFAna/Cuts/SpillCuts.h"
6 #include "CAFAna/Core/Loaders.h"
7 #include "CAFAna/Core/Binning.h"
8 #include "3FlavorAna/Prediction/PredictionGenerator3Flavor.h"
9 #include "CAFAna/Prediction/PredictionNoExtrap.h"
10 #include "CAFAna/Prediction/PredictionInterp.h"
11 #include "CAFAna/Systs/Systs.h"
12 #include "3FlavorAna/Vars/HistAxes.h"
13 #include "CAFAna/Vars/GenieWeights.h"
14 #include "OscLib/OscCalcDumb.h"
15 #include "TFile.h"
16 #include "TString.h"
17 
18 #include "futureSig_systs.h"
19 
20 using namespace ana;
21 
22 void setPathsFNALP123b( Loaders& loaders);
23 void setPathsFNALP3c( Loaders& loaders);
24 void setPathsRHC( Loaders& loadersRHC);
25 void SavePrediction (const TString part, const TString period, PredictionInterp * pred,
26  const TString predName, const TString version);
27 
28 void futureSig_getNueNumuSyst(const TString version="v01", const TString options = "fhc1 fhc2 rhc")
29 {
30  bool getFHC1 = options.Contains("fhc1");
31  bool getFHC2 = options.Contains("fhc2");
32  bool getRHC = options.Contains("rhc");
33 
34  ////////////////////////////////////////////////////////////
35  /// Setup loaders
36  ////////////////////////////////////////////////////////////
37  Loaders loadersp123b, loadersp3c;
38  loadersp123b.SetSpillCut(kStandardSpillCuts);
39  loadersp3c.SetSpillCut(kStandardSpillCuts);
40  setPathsFNALP123b(loadersp123b);
41  setPathsFNALP3c(loadersp3c);
42 
43  Loaders loadersRHC;
44  setPathsRHC(loadersRHC);
45 
46  osc::OscCalcDumb dumb;
47 
48  ////////////////////////////////////////////////////////////
49  /// Predictions: nue, numu; FHC P123b, FHC P3c,RHC
50  ////////////////////////////////////////////////////////////
51  // Nue FHC
52  const NuePropExtrapGenerator kNueFHCExtrap (kNueAxisSecondAnaCVN, kNumuCCAxis, kNueSAFDCVNSsb,kNueNDCVNSsb, kNumuND, kNoShift, kTuftsWeightCC);
53  PredictionInterp predNue_123b (systsNueFHC,dumb.Copy(),kNueFHCExtrap,loadersp123b);
54  PredictionInterp predNue_3c (systsNueFHC,dumb.Copy(),kNueFHCExtrap,loadersp3c);
55 
56  // Nue RHC
57  const NoExtrapGenerator kNueRHCExtrap(kNueAxisSecondAnaCVN, kNueSAFDCVNSsb);
58  PredictionInterp predNue_RHC (systsNueRHC,dumb.Copy(),kNueRHCExtrap,loadersRHC);
59 
60  // Numu FHC
61  const NumuExtrapGenerator kNumuFHCExtrap_full(kNumuCCAxis, kNumuFD, kNumuND, kNoShift, kTuftsWeightCC);
62  PredictionInterp predNumu_123b_full (systsNumuFHC,dumb.Copy(),kNumuFHCExtrap_full,loadersp123b);
63  PredictionInterp predNumu_3c_full (systsNumuFHC,dumb.Copy(),kNumuFHCExtrap_full,loadersp3c);
64 
65  //Numu RHC
66  const NoExtrapGenerator kNumuRHCExtrap_full(kNumuCCAxis, kNumuFD);
67  PredictionInterp predNumu_RHC_full (systsNumuRHC,dumb.Copy(),kNumuRHCExtrap_full,loadersRHC);
68 
69  // Test: Numu FHC RHC with half the background
70  const NumuExtrapGenerator kNumuFHCExtrap_half(kNumuCCAxis, kNumuFD, kNumuND, kNoShift, kTuftsWeightCC * kkNumuBkgHalf);
71  PredictionInterp predNumu_123b_half (systsNumuFHC,dumb.Copy(),kNumuFHCExtrap_half,loadersp123b);
72  PredictionInterp predNumu_3c_half (systsNumuFHC,dumb.Copy(),kNumuFHCExtrap_half,loadersp3c);
73 
74  const NoExtrapGenerator kNumuRHCExtrap_half(kNumuCCAxis, kNumuFD, kkNumuBkgHalf);
75  PredictionInterp predNumu_RHC_half (systsNumuRHC,dumb.Copy(),kNumuRHCExtrap_half,loadersRHC);
76 
77  ////////////////////////////////////////////////////////////
78  /// Go! Save everything
79  ////////////////////////////////////////////////////////////
80 
81  if(getFHC1) {
82  loadersp123b.Go();
83  SavePrediction("nue", "FHC1", &predNue_123b, "nue_fhcPrediction_p123b", version);
84  SavePrediction("numu", "FHC1_full", &predNumu_123b_full, "numu_fhcPrediction_p123b_full", version);
85  SavePrediction("numu", "FHC1_half", &predNumu_123b_half, "numu_fhcPrediction_p123b_half", version);
86  }
87 
88  if(getFHC2) {
89  loadersp3c.Go();
90  SavePrediction("nue", "FHC2", &predNue_3c, "nue_fhcPrediction_p3c", version);
91  SavePrediction("numu", "FHC2_full", &predNumu_3c_full, "numu_fhcPrediction_p3c_full", version);
92  SavePrediction("numu", "FHC2_half", &predNumu_3c_half, "numu_fhcPrediction_p3c_half", version);
93  }
94  if(getRHC) {
95  loadersRHC.Go();
96  SavePrediction("nue", "RHC", &predNue_RHC, "nue_rhcPrediction", version);
97  SavePrediction("numu", "RHC_full", &predNumu_RHC_full, "numu_rhcPrediction_full", version);
98  SavePrediction("numu", "RHC_half", &predNumu_RHC_half, "numu_rhcPrediction_half", version);
99  }
100 
101 }
102 
103 
104 void setPathsFNALP123b( Loaders& loaders)
105 {
106 
107  std::string path = "/pnfs/nova/persistent/production/concat/R16-03-03-prod2reco.f/";
108  std::string pathd= "/pnfs/nova/persistent/production/concat/R16-03-03-prod2reco.d/";
109 
110  std::string fname1(path+"period1/prod_decaf_R16-03-03-prod2reco.f_fd_genie_nonswap_fhc_nova_v08_period1_nue_or_numu_contain_v1_prod2-snapshot.root");
111  std::string fnameSwap1(path+"period1/prod_decaf_R16-03-03-prod2reco.f_fd_genie_fluxswap_fhc_nova_v08_period1_nue_or_numu_contain_v1_prod2-snapshot.root");
112  std::string fnameTau1(path+"period1/prod_decaf_R16-03-03-prod2reco.f_fd_genie_tau_fhc_nova_v08_period1_nue_or_numu_contain_v1_prod2-snapshot.root");
113  std::string fnameNDData(pathd+"prod_decaf_R16-03-03-prod2reco.d_nd_numi_fhc_full_nue_or_numu_contain_v1_goodruns_prod2-snapshot.root");
114  std::string fnameNDMC(pathd+"prod_decaf_R16-03-03-prod2reco.d_nd_genie_nonswap_genierw_fhc_nova_v08_epoch1-3c_nue_or_numu_contain_v1_prod2-snapshot/*.root");
115 
116  std::string fname2(path+"period2/prod_decaf_R16-03-03-prod2reco.f_fd_genie_nonswap_fhc_nova_v08_period2_nue_or_numu_contain_v1_prod2-snapshot.root");
117  std::string fnameSwap2(path+"period2/prod_decaf_R16-03-03-prod2reco.f_fd_genie_fluxswap_fhc_nova_v08_period2_nue_or_numu_contain_v1_prod2-snapshot.root");
118  std::string fnameTau2(path+"period2/prod_decaf_R16-03-03-prod2reco.f_fd_genie_tau_fhc_nova_v08_period2_nue_or_numu_contain_v1_prod2-snapshot.root");
119 
120  std::string fname3b(path+"epoch3b/prod_decaf_R16-03-03-prod2reco.f_fd_genie_nonswap_fhc_nova_v08_epoch3b_nue_or_numu_contain_v1_prod2-snapshot.root");
121  std::string fnameSwap3b(path+"epoch3b/prod_decaf_R16-03-03-prod2reco.f_fd_genie_fluxswap_fhc_nova_v08_epoch3b_nue_or_numu_contain_v1_prod2-snapshot.root");
122  std::string fnameTau3b(path+"epoch3b/prod_decaf_R16-03-03-prod2reco.f_fd_genie_tau_fhc_nova_v08_epoch3b_nue_or_numu_contain_v1_prod2-snapshot.root");
123 
124  loaders.SetLoaderFiles({fname1, fname2, fname3b}, caf::kFARDET, Loaders::kMC,
125  ana::kBeam, Loaders::kNonSwap);
126  loaders.SetLoaderFiles({fnameSwap1,fnameSwap2, fnameSwap3b}, caf::kFARDET, Loaders::kMC,
127  ana::kBeam, Loaders::kFluxSwap);
128  loaders.SetLoaderFiles({fnameTau1, fnameTau2, fnameTau3b}, caf::kFARDET, Loaders::kMC,
129  ana::kBeam, Loaders::kTauSwap);
130  loaders.SetLoaderPath(fnameNDData, caf::kNEARDET, Loaders::kData);
131  loaders.SetLoaderPath(fnameNDMC, caf::kNEARDET, Loaders::kMC);
132 }
133 
134 void setPathsFNALP3c( Loaders& loaders)
135 {
136  std::string path = "/pnfs/nova/persistent/production/concat/R16-03-03-prod2reco.f/";
137  std::string pathd= "/pnfs/nova/persistent/production/concat/R16-03-03-prod2reco.d/";
138 
139  std::string fname3c(path+"epoch3c/prod_decaf_R16-03-03-prod2reco.f_fd_genie_nonswap_fhc_nova_v08_epoch3c_nue_or_numu_contain_v1_prod2-snapshot.root");
140  std::string fnameSwap3c(path+"epoch3c/prod_decaf_R16-03-03-prod2reco.f_fd_genie_fluxswap_fhc_nova_v08_epoch3c_nue_or_numu_contain_v1_prod2-snapshot.root");
141  std::string fnameTau3c(path+"epoch3c/prod_decaf_R16-03-03-prod2reco.f_fd_genie_tau_fhc_nova_v08_epoch3c_nue_or_numu_contain_v1_prod2-snapshot.root");
142  std::string fnameNDData(pathd+"prod_decaf_R16-03-03-prod2reco.d_nd_numi_fhc_full_nue_or_numu_contain_v1_goodruns_prod2-snapshot.root");
143  std::string fnameNDMC(pathd+"prod_decaf_R16-03-03-prod2reco.d_nd_genie_nonswap_genierw_fhc_nova_v08_epoch1-3c_nue_or_numu_contain_v1_prod2-snapshot/*.root");
144 
145  loaders.SetLoaderPath(fname3c, caf::kFARDET, Loaders::kMC,
146  ana::kBeam, Loaders::kNonSwap);
147  loaders.SetLoaderPath(fnameSwap3c, caf::kFARDET, Loaders::kMC,
148  ana::kBeam, Loaders::kFluxSwap);
149  loaders.SetLoaderPath(fnameTau3c, caf::kFARDET, Loaders::kMC,
150  ana::kBeam, Loaders::kTauSwap);
151  loaders.SetLoaderPath(fnameNDData, caf::kNEARDET, Loaders::kData);
152  loaders.SetLoaderPath(fnameNDMC, caf::kNEARDET, Loaders::kMC);
153 }
154 
155 void setPathsRHC( Loaders& loadersRHC){
156  std::string fname = "prod_decaf_R16-03-03-prod2reco.e_fd_genie_nonswap_nogenierw_rhc_nova_v08_full_nue_or_numu_contain_v1";
157  std::string fswap ="prod_decaf_R16-03-03-prod2reco.e_fd_genie_fluxswap_nogenierw_rhc_nova_v08_full_nue_or_numu_contain_v1";
158 
159  //Unofficial concats
160  fname = "/pnfs/nova/users/ecatanom/concat/prod_decaf_R16-03-03-prod2reco.e_fd_genie_nonswap_nogenierw_rhc_nova_v08_full_nue_or_numu_contain_v1.root";
161  fswap = "/pnfs/nova/users/ecatanom/concat/prod_decaf_R16-03-03-prod2reco.e_fd_genie_fluxswap_nogenierw_rhc_nova_v08_full_nue_or_numu_contain_v1.root";
162 
163  loadersRHC.SetLoaderPath(fname,caf::kFARDET, Loaders::kMC, ana::kBeam, Loaders::kNonSwap);
164  loadersRHC.SetLoaderPath(fswap,caf::kFARDET, Loaders::kMC, ana::kBeam, Loaders::kFluxSwap);
165  loadersRHC.DisableLoader(caf::kFARDET, Loaders::kMC, ana::kBeam, Loaders::kTauSwap);
166 }
167 
168 void SavePrediction (const TString part, const TString period, PredictionInterp * pred, const TString predName, const TString version){
169  TString dir = "root_predictions/";
170  gSystem->mkdir(dir);
171  TFile * output = new TFile (dir+"predInterp_"+part+"_" + period + "_"+ version +".root", "RECREATE");
172  pred->SaveTo(output, predName);
173  std::cout << "Wrote prediction to " << output->GetName() << std::endl;
174  delete output;
175 }
caf::kNEARDET
Near Detector underground.
Definition: SREnums.h:10
plot_validation_datamc.fname
fname
Definition: plot_validation_datamc.py:180
output
ofstream output
Definition: plot_prediction_systs.C:32
plot_validation_datamc.pred
pred
Definition: plot_validation_datamc.py:175
ana::PredictionInterp
Implements systematic errors by interpolation between shifted templates.
Definition: PredictionInterp.h:23
caf::kFARDET
Far Detector at Ash River.
Definition: SREnums.h:11
ana
Cuts and Vars for the 2020 FD DiF Study.
Definition: vars.h:6
PandAna.cut.analysis_cuts.kNueNDCVNSsb
kNueNDCVNSsb
Definition: analysis_cuts.py:143
ana::Loaders
Collection of SpectrumLoaders for many configurations.
Definition: Loaders.h:23
systsNumuFHC
std::vector< const ISyst * > systsNumuFHC
Definition: futureSig_systs.h:88
ana::Loaders::DisableLoader
void DisableLoader(caf::Det_t det, DataMC datamc, DataSource src=kBeam, SwappingConfig swap=kNonSwap)
Definition: Loaders.cxx:65
ana::kNumuFD
const Cut kNumuFD
Definition: NumuCuts.h:53
ana::NoExtrapGenerator
Generates FD-only predictions (no extrapolation)
Definition: PredictionGenerator.h:32
kkNumuBkgHalf
const Var kkNumuBkgHalf([](const caf::SRProxy *sr){if((kkIsFD &&!kIsNumuCC)(sr)) return 0.5;return 1.;})
ana::NumuExtrapGenerator
Generates extrapolated Numu predictions.
Definition: PredictionGenerator3Flavor.h:20
ana::Loaders::Go
void Go()
Call Go() on all the loaders.
Definition: Loaders.cxx:162
ana::Loaders::SetSpillCut
void SetSpillCut(const SpillCut &cut)
Definition: Loaders.cxx:121
ana::PredictionInterp::SaveTo
virtual void SaveTo(TDirectory *dir, const std::string &name) const override
Definition: PredictionInterp.cxx:589
ana::NuePropExtrapGenerator
Generates extrapolated Nue predictions using ProportionalDecomp.
Definition: PredictionGenerator3Flavor.h:48
part
TString part[npart]
Definition: Style.C:32
ana::kNumuCCAxis
const HistAxis kNumuCCAxis("Reconstructed Neutrino Energy (GeV)", kNumuEnergyBinning, kCCE)
Definition: HistAxes.h:8
systsNueFHC
std::vector< const ISyst * > systsNueFHC
Definition: futureSig_systs.h:33
ana::kNumuND
const Cut kNumuND
Definition: NumuCuts.h:55
osc::OscCalcDumb::Copy
virtual IOscCalc * Copy() const
Definition: OscCalcDumb.h:19
futureSig_systs.h
osc::OscCalcDumb
Simple oscillation probability calculator that has no solar term or mass hierarchy or delta so it&#39;s s...
Definition: OscCalcDumb.h:16
setPathsRHC
void setPathsRHC(Loaders &loadersRHC)
Definition: futureSig_getNueNumuSyst.C:155
setPathsFNALP3c
void setPathsFNALP3c(Loaders &loaders)
Definition: futureSig_getNueNumuSyst.C:134
systsNumuRHC
std::vector< const ISyst * > systsNumuRHC
Definition: futureSig_systs.h:94
ana::kNoShift
const SystShifts kNoShift
Definition: SystShifts.h:115
om::cout
OStream cout
Definition: OStream.cxx:6
systsNueRHC
std::vector< const ISyst * > systsNueRHC
Definition: futureSig_systs.h:36
version
const XML_Char * version
Definition: expat.h:187
path
const std::string path
Definition: plot_BEN.C:43
ana::Loaders::SetLoaderFiles
void SetLoaderFiles(const std::vector< std::string > &files, caf::Det_t det, DataMC datamc, DataSource src=kBeam, SwappingConfig swap=kNonSwap)
Configure loader via explicit file list.
Definition: Loaders.cxx:45
SavePrediction
void SavePrediction(const TString part, const TString period, PredictionInterp *pred, const TString predName, const TString version)
Definition: futureSig_getNueNumuSyst.C:168
xml_schema::string
::xsd::cxx::tree::string< char, simple_type > string
Definition: Database.h:154
ana::kTuftsWeightCC
const Var kTuftsWeightCC
Definition: XsecTunes.h:31
dir
TDirectory * dir
Definition: macro.C:5
loaders
std::vector< Loaders * > loaders
Definition: syst_header.h:386
DCS_db_parser.period
int period
Definition: DCS_db_parser.py:19
futureSig_getNueNumuSyst
void futureSig_getNueNumuSyst(const TString version="v01", const TString options="fhc1 fhc2 rhc")
Definition: futureSig_getNueNumuSyst.C:28
ana::kStandardSpillCuts
const SpillCut kStandardSpillCuts
Apply this unless you&#39;re doing something special.
Definition: SpillCuts.h:46
setPathsFNALP123b
void setPathsFNALP123b(Loaders &loaders)
Definition: futureSig_getNueNumuSyst.C:104
ana::Loaders::SetLoaderPath
void SetLoaderPath(const std::string &path, caf::Det_t det, DataMC datamc, DataSource src=kBeam, SwappingConfig swap=kNonSwap)
Configure loader via wildcard path.
Definition: Loaders.cxx:25
Generated on Wed Sep 23 2020 05:13:29 for NOvA by   doxygen 1.8.11