PredictionFilePath.h
Go to the documentation of this file.
1 #pragma once
2 
3 namespace ana {
4 
5  std::string getNue2018PredFile(std::string beam, bool IsRealNDData = true){
7  bool isFHC = false;
8  if(beam == "FHC") isFHC = true;
9 
10  if(isFHC && IsRealNDData)
11  ret = "/nova/ana/nu_e_ana/Ana2018/Predictions/provisional_singles/NueProd4SystematicsOnRealNDData_2018-05-22/pred_allxp_Nue2018Axis_full_FHCAllSyst_nueconcat_realNDData_2018-05-22.root";
12  else if(isFHC && !IsRealNDData)
13  ret = "/nova/ana/nu_e_ana/Ana2018/Predictions/provisional_singles/NueProd4SystematicsOnFakeNDData_2018-05-23/pred_allxp_Nue2018Axis_full_FHCAllSysts_nueconcat_fakeNDData_05_22_2018_pass1_merged.root";
14 
15  else if(!isFHC && IsRealNDData)
16  ret = "/nova/ana/nu_e_ana/Ana2018/Predictions/provisional_singles/NueProd4SystematicsOnRealNDData_2018-05-22/pred_allxp_Nue2018Axis_full_RHCAllSyst_nueconcat_realNDData_2018-05-22.root";
17  else if(!isFHC && !IsRealNDData)
18  ret = "/nova/ana/nu_e_ana/Ana2018/Predictions/provisional_singles/NueProd4SystematicsOnFakeNDData_2018-05-23/pred_allxp_Nue2018Axis_full_RHCAllSysts_nueconcat_fakeNDData_05_22_2018_pass1_merged.root";
19 
20  else
21  std::cout<<"please, leave me along :)"<<std::endl;
22 
23  return ret;
24  }
25 
26  std::string getNumu2018PredFile(std::string beam, bool IsRealNDData = true){
28  bool isFHC = false;
29  if(beam == "FHC") isFHC = true;
30 
31  if(!IsRealNDData){
32  std::cout<<" Sorry there is no fake ND Data prediction for numu analysis, why do you need them? Ask systematics guys for fakd ND Data prediction files."<<std::endl;
33  std::cout<<" Returning REAL REAL REAL ND Data prediction path of numu prediction files..."<<std::endl;
34  }
35 
36  if(isFHC)
37  ret = "/nova/ana/nu_mu_ana/Ana2018/Predictions/pred_numuconcat_fhc__numu2018.root";
38  else
39  ret = "/nova/ana/nu_mu_ana/Ana2018/Predictions/pred_numuconcat_rhc__numu2018.root";
40  return ret;
41  }
42 
43 /* extrap = {nxp, xp_prop} , sigma = { -2, -1, 0, 1, 2}, not done yet*/
45 
47 
48  ret = "pred_xp_numu_Nominal/noShift/";
49 
50  if(extrap == "nxp")
51  ret = "pred_nxp_Nominal/noShift/";
52 
53  if( quant == -1)
54  ret += "numu_pred_NumuEnergy";
55  else ret += "numu_pred_Energy_Quant" + std::to_string(quant) ;
56  std::cout<<ret<<std::endl;
57  return ret;
58  }
59 
60 /* extrap = {nxp, xp_prop, xp_combo} , sigma = { -2, -1, 0, 1, 2}, not done yet*/
63 
64  if(sigma == 0)
65  ret = "pred_" +extrap+ "_Nominal/noShift/nue_pred_Nue2018Axis";
66  else{
67  std::cout<<"please implement for different systematics yourself :)"<<std::endl;
68  abort();
69  }
70  return ret;
71  }
72 
73 }
std::string getNumu2018PredVar(std::string extrap, int quant)
std::string getNumu2018PredFile(std::string beam, bool IsRealNDData=true)
Cuts and Vars for the 2020 FD DiF Study.
Definition: vars.h:6
static bool isFHC
double sigma(TH1F *hist, double percentile)
OStream cout
Definition: OStream.cxx:6
std::string getNue2018PredFile(std::string beam, bool IsRealNDData=true)
std::string getNue2018PredVar(std::string extrap, int sigma=0)
std::string to_string(ModuleType mt)
Definition: ModuleType.h:32
enum BeamMode string