Nus18Systs.h
Go to the documentation of this file.
1 #pragma once
2 
3 #include "CAFAna/Core/ISyst.h"
5 
6 #include "NuXAna/Systs/NusSysts.h" // for NusChannel
7 
8 #include "TString.h"
9 #include "TObjArray.h"
10 #include "TObjString.h"
11 
12 #include <vector>
13 
14 class TH1D;
15 
17 
18 namespace ana
19 {
20  /// Function to load systematics
21  std::string LoadSystematics(std::string polarity, bool extrap = false) {
22 
23  DontAddDirectory guard;
24 
25  char* path = std::getenv("NUSDATA_NUS18_SYSTS");
26  if (!path)
27  throw std::runtime_error("NUSDATA_NUS18_SYSTS environment variable not set! Did you forget to set up nusdata?");
28  else if (polarity != "fhc" && polarity != "rhc")
29  throw std::runtime_error(Form("Polarity %s invalid: systematics can only be loaded for \"fhc\" or \"rhc\".", polarity.c_str()));
30  // else
31  // std::cout << "Loading " << polarity << " systematics from " << path << std::endl;
32 
33  std::string systFile = Form("%s/Nus18%sSystRatios%s.root", path, polarity == "fhc"? "FHC" : "RHC", extrap? "Extrap" : "");
34  return systFile;
35  }
36 
37  //----------------------------------------------------------------------
38  /// Nus group systematic applied as a shape shift, with shape information
39  /// loaded from input histograms
40  class Nus18SystFromHist: public ISyst
41  {
42  public:
43  Nus18SystFromHist(const std::string &fname, const std::string &det, const std::string &polarity,
44  const std::string &systname, const std::string &latexname);
46 
47  void Shift(double sigma,
48  caf::SRProxy* sr, double &weight) const override;
49 
50  double WeightFor(NusChannel chan, double sigma, double nueenergy) const;
51 
52  private:
56 
57  void LoadHists() const;
58 
59  // get oscillation channel for current event
61  // Store hist for sig and bkg for each syst, store a pair of a
62  // sigma shift value and a histogram evaluated at that shift
63  mutable std::vector< std::vector< std::pair<int, TH1D*> > > fHists;
64  };
65 
66  //----------------------------------------------------------------------
67  /// Nus group systematic applied as a flat shift only dependent on NC or CC
68  class Nus18FlatSyst: public ISyst
69  {
70  public:
71  Nus18FlatSyst(const std::string &systname, const std::string &latexname,const std::string &det,const std::string &polarity,
72  const double &ncWei, const double &bgWei);
73  ~Nus18FlatSyst();
74 
75  void Shift(double sigma,
76  caf::SRProxy* sr, double &weight) const override;
77 
78  private:
79  double fNCWei;
80  double fBGWei;
83  };
84 
85  //----------------------------------------------------------------------
86  class Nus18BaseSyst: public ISyst
87  {
88  public:
89  Nus18BaseSyst(const std::string& shortname,
90  const std::string& latexname);
91  ~Nus18BaseSyst();
92 
93  void Shift(double sigma,
94  caf::SRProxy* sr, double &weight) const override;
95  };
96 
100 
101  //----------------------------------------------------------------------
102  // Systematics based on docdb-15330
103 
104  // Systematics handled using shape information
105 
106  // Beam
113 
114  // Kaon
121 
122  // Neutron
129 
130  // Tau
137 
138  // Cherenkov
145 
146  // CalibFlat
153 
154  // CalibRel
161 
162  // CalibShape
169 
170  // GENIEGen
177 
178  // GENIENOvA
185 
186  // LightL
193 
194  // NueCC
195  extern const Nus18BaseSyst kNus18BaseNueCCSyst;
201 
202  // NumuCC
203  extern const Nus18BaseSyst kNus18BaseNumuCCSyst;
209 
210  // NC
211  extern const Nus18BaseSyst kNus18BaseNCSyst;
217 
218  // Nue ND Univ
222 // extern const Nus18SystFromHist kNus18RHCNDNueUnivSyst;
223 // extern const Nus18FlatSyst kNus18RHCFDNueUnivSyst;
224 // extern const Nus18SystFromHist kNus18RHCEXNueUnivSyst;
225 
226  // Tau ND Univ
230 // extern const Nus18SystFromHist kNus18RHCNDTauUnivSyst;
231 // extern const Nus18FlatSyst kNus18RHCFDTauUnivSyst;
232 // extern const Nus18SystFromHist kNus18RHCEXTauUnivSyst;
233 
234  // FD OscP Univ
241 
242  // Acc
243  extern const Nus18BaseSyst kNus18BaseAccSyst;
249 
250  // Systematics applied as flat shifts
251  extern const Nus18FlatSyst kNus18FHCNormSyst;
253 
254  extern const Nus18FlatSyst kNus18RHCNormSyst;
256 
257  // PPFX
264 
265  // ND normalisation
269  // extern const Nus18FlatSyst kNus18RHCNDNormNDSyst;
270  // extern const Nus18FlatSyst kNus18RHCFDNormNDSyst;
271  // extern const Nus18FlatSyst kNus18RHCEXNormSyst;
272 
273  // FD normalisation
277  // extern const Nus18FlatSyst kNus18RHCNDNormSyst;
278  // extern const Nus18FlatSyst kNus18RHCFDNormSyst;
279  // extern const Nus18FlatSyst kNus18RHCEXNormSyst;
280 
281  // XSec
288 
289  // GENIE individual
369  //GENIE bg dominated syst
370 
401 
402 
403  //----------------------------------------------------------------------
404 
405  // Get a vector of all the nus group systs
406  void GetSystType(std::string syst_type, bool &use_xsec,
407  bool &use_nonxsec);
408  // Base
409  std::vector<const ISyst*> GetNus18Systs(bool rhc, std::string det_type,
410  std::string syst_type);
411 
412  std::vector<const ISyst*> GetNus18BaseSysts(bool rhc,
413  std::string syst_type);
414 
415  // FHC
416  std::vector<const ISyst*> GetNus18FHCNDSysts(std::string syst_type);
417  std::vector<const ISyst*> GetNus18FHCFDSysts(std::string syst_type);
418 
419  // RHC
420  std::vector<const ISyst*> GetNus18RHCNDSysts(std::string syst_type);
421  std::vector<const ISyst*> GetNus18RHCFDSysts(std::string syst_type);
422  std::vector<const ISyst*> GetNus18RHCExtrapSysts();
423 
424 } // namespace
const Nus18BaseSyst kNus18BaseGENIEDISvnCC0piSyst("DISvnCC0pi18","GENIE DISvnCC0pi")
Definition: Nus18Systs.h:297
std::vector< const ISyst * > GetNus18BaseSysts(bool rhc, std::string syst_type)
Definition: Nus18Systs.cxx:460
const Nus18BaseSyst kNus18BaseGENIEGenSyst("GENIEg18","GENIE general syst")
Definition: Nus18Systs.h:171
const Nus18SystFromHist kNus18RHCEXCalibShapeSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","CalShape18","Calibration shape")
Definition: Nus18Systs.h:168
const Nus18SystFromHist kNus18RHCNDCalibRelSyst
const Nus18SystFromHist kNus18FHCNDNCSyst
const Nus18SystFromHist kNus18FHCFDGENIEMECEShapeNuSyst
double WeightFor(NusChannel chan, double sigma, double nueenergy) const
Definition: Nus18Systs.cxx:304
const Nus18BaseSyst kNus18BaseGENIEMFPpiSyst("MFPpi18","GENIE MFPpi")
Definition: Nus18Systs.h:375
const Nus18SystFromHist kNus18RHCFDNeutronSyst(kNusAna18RHCSystFile,"FD","RHC","Neutron18","Neutron")
Definition: Nus18Systs.h:127
std::vector< const ISyst * > GetNus18FHCFDSysts(std::string syst_type)
Definition: Nus18Systs.cxx:630
const Nus18SystFromHist kNus18RHCEXGENIEGenSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","GENIEg18","GENIE General Systs")
Definition: Nus18Systs.h:176
const Nus18SystFromHist kNus18FHCNDGENIENOvASyst(kNusAna18FHCSystFile,"ND","FHC","GENIEn18","GENIE NOvA Specific Systs (ND)")
Definition: Nus18Systs.h:180
const Nus18SystFromHist kNus18RHCEXGENIEDISvbarpNC2piSyst
const Nus18SystFromHist kNus18RHCEXNumuCCSyst
const Nus18FlatSyst kNus18FHCFDNormFDSyst("NormFD18","Normalisation FD","FD","FHC", 1.5, 1.5)
Definition: Nus18Systs.h:276
const Nus18SystFromHist kNus18FHCNDNueCCSyst
Cuts and Vars for the 2020 FD DiF Study.
Definition: vars.h:6
const Nus18BaseSyst kNus18BaseGENIEMECEShapeNuSyst("MECEShapeNu18","GENIE MECEShapeNu")
Definition: Nus18Systs.h:374
const Nus18SystFromHist kNus18RHCNDGENIEDISvbarpNC2piSyst(kNusAna18RHCSystFile,"ND","RHC","DISvbarpNC2pi18","GENIE DISvbarpNC2pi")
Definition: Nus18Systs.h:351
const Nus18BaseSyst kNus18BaseGENIEMaCOHpiSyst("MaCOHpi18","GENIE MaCOHpi")
Definition: Nus18Systs.h:292
const Nus18SystFromHist kNus18FHCNDNueUnivSyst(kNusAna18FHCSystFile,"ND","FHC","NueNDUniv18","Nue ND")
Definition: Nus18Systs.h:220
Nus18SystFromHist(const std::string &fname, const std::string &det, const std::string &polarity, const std::string &systname, const std::string &latexname)
Definition: Nus18Systs.cxx:227
std::string LoadSystematics(std::string polarity, bool extrap=false)
Function to load systematics.
Definition: Nus18Systs.h:21
void LoadHists() const
Definition: Nus18Systs.cxx:243
const Nus18BaseSyst kNus18BaseNueCCSyst
const Nus18SystFromHist kNus18RHCFDGENIERPARESSyst(kNusAna18RHCSystFile,"FD","RHC","RPARES18","GENIE RPARES")
Definition: Nus18Systs.h:388
const Nus18SystFromHist kNus18RHCFDGENIEMECEShapeNuSyst(kNusAna18RHCSystFile,"FD","RHC","MECEShapeNu18","GENIE MECEShapeNu")
Definition: Nus18Systs.h:392
const Nus18SystFromHist kNus18FHCFDGENIEMaCCRESSyst(kNusAna18FHCSystFile,"FD","FHC","MaCCRES18","GENIE MaCCRES")
Definition: Nus18Systs.h:380
const Nus18BaseSyst kNus18BaseGENIEDISvbarpNC2piSyst("DISvbarpNC2pi18","GENIE DISvbarpNC2pi")
Definition: Nus18Systs.h:296
const Nus18SystFromHist kNus18FHCFDGENIERPARESSyst
const Nus18BaseSyst kNus18BaseGENIEDISvpNC1piSyst("DISvpNC1pi18","GENIE DISvpNC1pi")
Definition: Nus18Systs.h:307
const Nus18SystFromHist kNus18RHCFDGENIER0COHpiSyst(kNusAna18RHCSystFile,"FD","RHC","R0COHpi18","GENIE R0COHpi")
Definition: Nus18Systs.h:336
const Nus18SystFromHist kNus18FHCNDCherenkovSyst(kNusAna18FHCSystFile,"ND","FHC","Cherenkov18","Cherenkov (ND)")
Definition: Nus18Systs.h:140
const Var weight
const Nus18BaseSyst kNus18BaseGENIEDISvbarpNC3piSyst("DISvbarpNC3pi18","GENIE DISvbarpNC3pi")
Definition: Nus18Systs.h:299
const Nus18SystFromHist kNus18FHCNDCalibRelSyst
Definition: Nus18Systs.h:157
const Nus18SystFromHist kNus18RHCNDGENIEDISvbarpNC3piSyst(kNusAna18RHCSystFile,"ND","RHC","DISvbarpNC3pi18","GENIE DISvbarpNC3pi")
Definition: Nus18Systs.h:354
const Nus18SystFromHist kNus18FHCNDGENIEDISvnNC1piSyst(kNusAna18FHCSystFile,"ND","FHC","DISvnNC1pi18","GENIE DISvnNC1pi")
Definition: Nus18Systs.h:326
const Nus18FlatSyst kNus18FHCNDNormFDSyst("NormFD18","Normalisation FD","ND","FHC", 0, 0)
Definition: Nus18Systs.h:275
const Nus18SystFromHist kNus18RHCNDNueCCSyst
const Nus18SystFromHist kNus18RHCEXTauSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","TauScale18","Tau")
Definition: Nus18Systs.h:136
const Nus18SystFromHist kNus18FHCFDGENIEDISvpNC3piSyst(kNusAna18FHCSystFile,"FD","FHC","DISvpNC3pi18","GENIE DISvpNC3pi")
Definition: Nus18Systs.h:320
const Nus18SystFromHist kNus18RHCEXGENIEDISvbarpNC3piSyst
const Nus18BaseSyst kNus18BaseGENIEDISvnNC1piSyst("DISvnNC1pi18","GENIE DISvnNC1pi")
Definition: Nus18Systs.h:302
Proxy for caf::StandardRecord.
Definition: SRProxy.h:2126
const Nus18SystFromHist kNus18RHCEXGENIEMaCOHpiSyst
const Nus18BaseSyst kNus18BaseNormFDSyst("NormFD18","FD normalisation syst")
Definition: Nus18Systs.h:274
const Nus18BaseSyst kNus18BaseXSecOnOffSyst("XSecOnOff18","Cross-section On/Off")
Definition: Nus18Systs.h:282
const Nus18SystFromHist kNus18FHCFDAccSyst
const Nus18SystFromHist kNus18RHCNDGENIER0COHpiSyst(kNusAna18RHCSystFile,"ND","RHC","R0COHpi18","GENIE R0COHpi")
Definition: Nus18Systs.h:348
const Nus18FlatSyst kNus18RHCFDPPFXSyst("PPFX18","Fluxweight FD","FD","RHC", 11.3, 11.2)
Definition: Nus18Systs.h:262
const Nus18SystFromHist kNus18FHCNDGENIEDISvnNC2piSyst(kNusAna18FHCSystFile,"ND","FHC","DISvnNC2pi18","GENIE DISvnNC2pi")
Definition: Nus18Systs.h:327
const Nus18SystFromHist kNus18FHCFDGENIEDISvpNC1piSyst(kNusAna18FHCSystFile,"FD","FHC","DISvpNC1pi18","GENIE DISvpNC1pi")
Definition: Nus18Systs.h:319
const Nus18SystFromHist kNus18RHCFDGENIEMaCCRESSyst(kNusAna18RHCSystFile,"FD","RHC","MaCCRES18","GENIE MaCCRES")
Definition: Nus18Systs.h:378
const Nus18SystFromHist kNus18RHCEXNCSyst
const Nus18FlatSyst kNus18FHCFDTauUnivSyst("TauNDUniv18","Tau ND","FD","FHC", 0., 0.)
Definition: Nus18Systs.h:229
const Nus18SystFromHist kNus18RHCNDKaonSyst(kNusAna18RHCSystFile,"ND","RHC","Kaon18","Kaon (ND)")
Definition: Nus18Systs.h:118
const Nus18SystFromHist kNus18RHCEXNueCCSyst
const Nus18SystFromHist kNus18RHCEXCherenkovSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","Cherenkov18","Cherenkov")
Definition: Nus18Systs.h:144
const Nus18SystFromHist kNus18RHCNDGENIEMaCOHpiSyst(kNusAna18RHCSystFile,"ND","RHC","MaCOHpi18","GENIE MaCOHpi")
Definition: Nus18Systs.h:347
const Nus18BaseSyst kNus18BaseGENIER0COHpiSyst("R0COHpi18","GENIE R0COHpi")
Definition: Nus18Systs.h:293
const Nus18SystFromHist kNus18FHCNDAccSyst
const Nus18SystFromHist kNus18FHCNDTauUnivSyst(kNusAna18FHCSystFile,"ND","FHC","TauNDUniv18","Tau ND")
Definition: Nus18Systs.h:228
const Nus18SystFromHist kNus18FHCNDGENIEMaNCRESSyst(kNusAna18FHCSystFile,"ND","FHC","MaNCRES18","GENIE MaNCRES")
Definition: Nus18Systs.h:321
const Nus18SystFromHist kNus18RHCEXGENIEMANCELSyst
const Nus18SystFromHist kNus18FHCNDCalibFlatSyst(kNusAna18FHCSystFile,"ND","FHC","CalFlat18","Flat Miscalibration (ND)")
Definition: Nus18Systs.h:148
const Nus18SystFromHist kNus18FHCFDGENIEMaCOHpiSyst(kNusAna18FHCSystFile,"FD","FHC","MaCOHpi18","GENIE MaCOHpi")
Definition: Nus18Systs.h:311
const Nus18BaseSyst kNus18BaseCalibShapeSyst("CalShape18","Calibration shape")
Definition: Nus18Systs.h:163
const std::string kNusAna18FHCSystFile
Definition: Nus18Systs.h:97
std::vector< const ISyst * > GetNus18RHCFDSysts(std::string syst_type)
Definition: Nus18Systs.cxx:680
const Nus18SystFromHist kNus18RHCNDGENIERPARESSyst(kNusAna18RHCSystFile,"ND","RHC","RPARES18","GENIE RPARES")
Definition: Nus18Systs.h:387
std::vector< std::vector< std::pair< int, TH1D * > > > fHists
Definition: Nus18Systs.h:63
Encapsulate code to systematically shift a caf::SRProxy.
Definition: ISyst.h:14
const Nus18BaseSyst kNus18BaseCalibRelSyst("CalRelFD18","Relative detector calibration")
Definition: Nus18Systs.h:155
const Nus18SystFromHist kNus18RHCNDLightLSyst(kNusAna18RHCSystFile,"ND","RHC","CalLight18","Light Level syst (ND)")
Definition: Nus18Systs.h:190
const Nus18SystFromHist kNus18RHCFDCalibShapeSyst(kNusAna18RHCSystFile,"FD","RHC","CalShape18","Calibration shape")
Definition: Nus18Systs.h:167
const Nus18SystFromHist kNus18FHCNDKaonSyst(kNusAna18FHCSystFile,"ND","FHC","Kaon18","Kaon (ND)")
Definition: Nus18Systs.h:116
const Nus18SystFromHist kNus18RHCFDBeamSyst(kNusAna18RHCSystFile,"FD","RHC","Beam18","All Beam")
Definition: Nus18Systs.h:111
const Nus18SystFromHist kNus18FHCFDGENIEMaNCRESSyst(kNusAna18FHCSystFile,"FD","FHC","MaNCRES18","GENIE MaNCRES")
Definition: Nus18Systs.h:309
const Nus18BaseSyst kNus18BaseAccSyst
const Nus18BaseSyst kNus18BaseGENIEMvNCRESSyst("MvNCRES18","GENIE MvNCRES")
Definition: Nus18Systs.h:305
const Nus18SystFromHist kNus18RHCFDGENIEGenSyst(kNusAna18RHCSystFile,"FD","RHC","GENIEg18","GENIE General Systs")
Definition: Nus18Systs.h:175
const Nus18SystFromHist kNus18RHCFDGENIERadCorrNueSyst(kNusAna18RHCSystFile,"FD","RHC","RadCorrNue18","GENIE RadCorrNue")
Definition: Nus18Systs.h:341
const Nus18SystFromHist kNus18RHCFDGENIEMvCCRESSyst(kNusAna18RHCSystFile,"FD","RHC","MvCCRES18","GENIE MvCCRES")
Definition: Nus18Systs.h:382
std::vector< const ISyst * > GetNus18RHCNDSysts(std::string syst_type)
Definition: Nus18Systs.cxx:584
const Nus18SystFromHist kNus18RHCNDGENIEDISvbarpCC0piSyst(kNusAna18RHCSystFile,"ND","RHC","DISvbarpCC0pi18","GENIE DISvbarpCC0pi")
Definition: Nus18Systs.h:356
const Nus18BaseSyst kNus18BaseNCSyst
const Nus18SystFromHist kNus18FHCFDGENIEMFPpiSyst(kNusAna18FHCSystFile,"FD","FHC","MFPpi18","GENIE MFPpi")
Definition: Nus18Systs.h:394
const Nus18FlatSyst kNus18RHCOscParamSyst
const Nus18SystFromHist kNus18RHCFDGENIEDISvbarnNC3piSyst(kNusAna18RHCSystFile,"FD","RHC","DISvbarnNC3pi18","GENIE DISvbarnNC3pi")
Definition: Nus18Systs.h:338
const Nus18FlatSyst kNus18FHCFDPPFXSyst("PPFX18","Fluxweight FD","FD","FHC", 9.04, 11.8)
Definition: Nus18Systs.h:260
const Nus18BaseSyst kNus18BaseGENIEMANCELSyst("MANCEL18","GENIE MANCEL")
Definition: Nus18Systs.h:294
const Nus18SystFromHist kNus18RHCNDXSecOnOffSyst(kNusAna18RHCSystFile,"ND","RHC","XSecOnOff18","Cross-section On/Off")
Definition: Nus18Systs.h:285
const Nus18FlatSyst kNus18FHCNormSyst
const Nus18SystFromHist kNus18FHCFDTauSyst(kNusAna18FHCSystFile,"FD","FHC","TauScale18","Tau scale")
Definition: Nus18Systs.h:133
const Nus18SystFromHist kNus18FHCNDGENIEDISvpNC0piSyst(kNusAna18FHCSystFile,"ND","FHC","DISvpNC0pi18","GENIE DISvpNC0pi")
Definition: Nus18Systs.h:330
const Nus18SystFromHist kNus18RHCFDOscPUnivSyst(kNusAna18RHCSystFile,"FD","RHC","FDOscPUniv18","Oscillation parameters")
Definition: Nus18Systs.h:239
const Nus18BaseSyst kNus18BaseTauSyst("TauScale18","Tau")
Definition: Nus18Systs.h:131
const Nus18SystFromHist kNus18RHCFDGENIEDISvbarpNC3piSyst(kNusAna18RHCSystFile,"FD","RHC","DISvbarpNC3pi18","GENIE DISvbarpNC3pi")
Definition: Nus18Systs.h:342
const Nus18SystFromHist kNus18RHCFDCalibRelSyst
const Nus18BaseSyst kNus18BaseGENIEDISvbarnNC3piSyst("DISvbarnNC3pi18","GENIE DISvbarnNC3pi")
Definition: Nus18Systs.h:295
const Nus18SystFromHist kNus18FHCNDGENIEMFPpiSyst(kNusAna18FHCSystFile,"ND","FHC","MFPpi18","GENIE MFPpi")
Definition: Nus18Systs.h:393
const Nus18SystFromHist kNus18RHCFDGENIEMANCELSyst(kNusAna18RHCSystFile,"FD","RHC","MANCEL18","GENIE MANCEL")
Definition: Nus18Systs.h:337
const Nus18SystFromHist kNus18RHCFDGENIEDISvbarpNC2piSyst(kNusAna18RHCSystFile,"FD","RHC","DISvbarpNC2pi18","GENIE DISvbarpNC2pi")
Definition: Nus18Systs.h:339
const Nus18SystFromHist kNus18RHCNDGENIENOvASyst(kNusAna18RHCSystFile,"ND","RHC","GENIEn18","GENIE NOvA Specific Systs (ND)")
Definition: Nus18Systs.h:182
const Nus18SystFromHist kNus18RHCNDGENIEMFPpiSyst
const Nus18SystFromHist kNus18FHCNDNumuCCSyst
const Nus18FlatSyst kNus18RHCNDOscPUnivSyst("FDOscPUniv18","FD oscillation parameters","ND","RHC", 0., 0.)
Definition: Nus18Systs.h:238
const Nus18BaseSyst kNus18BaseNeutronSyst("Neutron18","Neutron")
Definition: Nus18Systs.h:123
const Nus18SystFromHist kNus18FHCNDGENIEDISvnNC3piSyst(kNusAna18FHCSystFile,"ND","FHC","DISvnNC3pi18","GENIE DISvnNC3pi")
Definition: Nus18Systs.h:328
const Nus18SystFromHist kNus18RHCNDGENIEMECEShapeNuSyst(kNusAna18RHCSystFile,"ND","RHC","MECEShapeNu18","GENIE MECEShapeNu")
Definition: Nus18Systs.h:391
const Nus18SystFromHist kNus18FHCFDCherenkovSyst(kNusAna18FHCSystFile,"FD","FHC","Cherenkov18","Cherenkov")
Definition: Nus18Systs.h:141
const Nus18BaseSyst kNus18BaseCherenkovSyst("Cherenkov18","Cherenkov")
Definition: Nus18Systs.h:139
const Nus18SystFromHist kNus18FHCFDLightLSyst(kNusAna18FHCSystFile,"FD","FHC","CalLight18","Light Level syst")
Definition: Nus18Systs.h:189
const Nus18SystFromHist kNus18FHCNDGENIEDISvpNC1piSyst(kNusAna18FHCSystFile,"ND","FHC","DISvpNC1pi18","GENIE DISvpNC1pi")
Definition: Nus18Systs.h:331
const Nus18SystFromHist kNus18RHCEXGENIEMaNCRESSyst
const Nus18SystFromHist kNus18RHCNDNeutronSyst(kNusAna18RHCSystFile,"ND","RHC","Neutron18","Neutron (ND)")
Definition: Nus18Systs.h:126
const Nus18SystFromHist kNus18FHCFDNeutronSyst(kNusAna18FHCSystFile,"FD","FHC","Neutron18","Neutron")
Definition: Nus18Systs.h:125
std::string getenv(std::string const &name)
const Nus18SystFromHist kNus18RHCNDGENIEFormZoneSyst(kNusAna18RHCSystFile,"ND","RHC","FormZone18","GENIE FormZone")
Definition: Nus18Systs.h:346
const Nus18SystFromHist kNus18FHCFDOscPUnivSyst(kNusAna18FHCSystFile,"FD","FHC","FDOscPUniv18","FD oscillation parameters")
Definition: Nus18Systs.h:237
const Nus18SystFromHist kNus18RHCNDGENIEMANCELSyst(kNusAna18RHCSystFile,"ND","RHC","MANCEL18","GENIE MANCEL")
Definition: Nus18Systs.h:349
const Nus18SystFromHist kNus18RHCEXGENIERadCorrNueSyst
const Nus18SystFromHist kNus18FHCFDGENIEFormZoneSyst(kNusAna18FHCSystFile,"FD","FHC","FormZone18","GENIE FormZone")
Definition: Nus18Systs.h:310
const Nus18SystFromHist kNus18RHCFDGENIEDISvbarpNC1piSyst(kNusAna18RHCSystFile,"FD","RHC","DISvbarpNC1pi18","GENIE DISvbarpNC1pi")
Definition: Nus18Systs.h:343
const Nus18BaseSyst kNus18BaseGENIEBhBYSyst("BhtBY18","GENIE BhtBY")
Definition: Nus18Systs.h:376
const Nus18BaseSyst kNus18BaseGENIEDISvnNC3piSyst("DISvnNC3pi18","GENIE DISvnNC3pi")
Definition: Nus18Systs.h:304
const Nus18BaseSyst kNus18BaseGENIEDISvbarpCC0piSyst("DISvbarpCC0pi18","GENIE DISvbarpCC0pi")
Definition: Nus18Systs.h:301
const Nus18BaseSyst kNus18BaseGENIERPARESSyst("RPARES18","GENIE RPARES")
Definition: Nus18Systs.h:373
void GetSystType(std::string syst_type, bool &use_xsec, bool &use_nonxsec)
Definition: Nus18Systs.cxx:415
const Nus18SystFromHist kNus18FHCFDGENIENOvASyst(kNusAna18FHCSystFile,"FD","FHC","GENIEn18","GENIE NOvA Specific Systs")
Definition: Nus18Systs.h:181
const Nus18SystFromHist kNus18FHCNDGENIEMvNCRESSyst(kNusAna18FHCSystFile,"ND","FHC","MvNCRES18","GENIE MvNCRES")
Definition: Nus18Systs.h:329
const Nus18FlatSyst kNus18FHCNDOscPUnivSyst("FDOscPUniv18","FD oscillation parameters","ND","FHC", 0., 0.)
Definition: Nus18Systs.h:236
const Nus18SystFromHist kNus18FHCFDGENIEDISvnNC2piSyst(kNusAna18FHCSystFile,"FD","FHC","DISvnNC2pi18","GENIE DISvnNC2pi")
Definition: Nus18Systs.h:315
const Nus18SystFromHist kNus18RHCEXGENIEFormZoneSyst
std::string fDet
Definition: Nus18Systs.h:81
caf::StandardRecord * sr
const Nus18SystFromHist kNus18RHCFDGENIEMaNCRESSyst(kNusAna18RHCSystFile,"FD","RHC","MaNCRES18","GENIE MaNCRES")
Definition: Nus18Systs.h:333
const Nus18SystFromHist kNus18RHCEXBeamSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","Beam18","All Beam")
Definition: Nus18Systs.h:112
const Nus18SystFromHist kNus18FHCFDGENIEMvNCRESSyst(kNusAna18FHCSystFile,"FD","FHC","MvNCRES18","GENIE MvNCRES")
Definition: Nus18Systs.h:317
const Nus18FlatSyst kNus18FHCFDNormNDSyst("NormND18","Normalisation ND","FD","FHC", 0, 0)
Definition: Nus18Systs.h:268
const Nus18SystFromHist kNus18RHCFDGENIENOvASyst(kNusAna18RHCSystFile,"FD","RHC","GENIEn18","GENIE NOvA Specific Systs")
Definition: Nus18Systs.h:183
const Nus18SystFromHist kNus18RHCNDNumuCCSyst
const Nus18SystFromHist kNus18RHCFDCalibFlatSyst(kNusAna18RHCSystFile,"FD","RHC","CalFlat18","Flat Miscalibration")
Definition: Nus18Systs.h:151
const Nus18BaseSyst kNus18BaseGENIERadCorrNueSyst("RadCorrNue18","GENIE RadCorrNue")
Definition: Nus18Systs.h:298
const Nus18SystFromHist kNus18FHCFDCalibFlatSyst(kNusAna18FHCSystFile,"FD","FHC","CalFlat18","Flat Miscalibration")
Definition: Nus18Systs.h:149
const Nus18SystFromHist kNus18FHCFDNCSyst
const Nus18SystFromHist kNus18FHCNDGENIEMaCCRESSyst(kNusAna18FHCSystFile,"ND","FHC","MaCCRES18","GENIE MaCCRES")
Definition: Nus18Systs.h:379
const std::string kNusAna18RHCExtrapSystFile
Definition: Nus18Systs.h:99
const Nus18BaseSyst kNus18BaseNumuCCSyst
const Nus18SystFromHist kNus18FHCFDXSecOnOffSyst(kNusAna18FHCSystFile,"FD","FHC","XSecOnOff18","Cross-section On/Off")
Definition: Nus18Systs.h:284
const Nus18SystFromHist kNus18RHCEXCalibRelSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","CalRelFD18","Relative detector calibration")
Definition: Nus18Systs.h:160
const Nus18SystFromHist kNus18RHCEXGENIEDISvnCC0piSyst
const Nus18BaseSyst kNus18BaseGENIEDISvpNC0piSyst("DISvpNC0pi18","GENIE DISvpNC0pi")
Definition: Nus18Systs.h:306
const Nus18SystFromHist kNus18RHCNDCalibShapeSyst(kNusAna18RHCSystFile,"ND","RHC","CalShape18","Calibration shape (ND)")
Definition: Nus18Systs.h:166
const Nus18SystFromHist kNus18FHCNDGENIERPARESSyst
const Nus18SystFromHist kNus18FHCFDGENIEMvCCRESSyst(kNusAna18FHCSystFile,"FD","FHC","MvCCRES18","GENIE MvCCRES")
Definition: Nus18Systs.h:384
const Nus18BaseSyst kNus18BaseNormNDSyst("NormND18","ND normalisation syst")
Definition: Nus18Systs.h:266
double sigma(TH1F *hist, double percentile)
const Nus18BaseSyst kNus18BaseGENIEDISvpNC3piSyst("DISvpNC3pi18","GENIE DISvpNC3pi")
Definition: Nus18Systs.h:308
const Nus18SystFromHist kNus18FHCNDBeamSyst(kNusAna18FHCSystFile,"ND","FHC","Beam18","All Beam (ND)")
Definition: Nus18Systs.h:108
const Nus18BaseSyst kNus18BaseKaonSyst("Kaon18","Kaon")
Definition: Nus18Systs.h:115
const Nus18SystFromHist kNus18RHCFDGENIEMFPpiSyst
det_type
Definition: Constants.h:21
const Nus18SystFromHist kNus18FHCFDGENIEMANCELSyst(kNusAna18FHCSystFile,"FD","FHC","MANCEL18","GENIE MANCEL")
Definition: Nus18Systs.h:313
const Nus18FlatSyst kNus18FHCOscParamSyst
const Nus18FlatSyst kNus18RHCEXPPFXSyst("PPFX18","Fluxweight EX","EX","RHC", 1.33, 4.63)
Definition: Nus18Systs.h:263
const Nus18SystFromHist kNus18RHCFDGENIEMaCOHpiSyst(kNusAna18RHCSystFile,"FD","RHC","MaCOHpi18","GENIE MaCOHpi")
Definition: Nus18Systs.h:335
const Nus18FlatSyst kNus18FHCFDNueUnivSyst("NueNDUniv18","Nue ND","FD","FHC", 0., 0.)
Definition: Nus18Systs.h:221
const Nus18SystFromHist kNus18FHCNDNeutronSyst(kNusAna18FHCSystFile,"ND","FHC","Neutron18","Neutron (ND)")
Definition: Nus18Systs.h:124
const std::string path
Definition: plot_BEN.C:43
const Nus18SystFromHist kNus18RHCEXGENIEDISvbarpNC1piSyst
const Nus18SystFromHist kNus18RHCNDGENIEBhBYSyst
const Nus18BaseSyst kNus18BaseGENIEMaCCRESSyst("MaCCRES18","GENIE MaCCRES")
Definition: Nus18Systs.h:371
const Nus18SystFromHist kNus18FHCFDGENIEDISvnNC3piSyst(kNusAna18FHCSystFile,"FD","FHC","DISvnNC3pi18","GENIE DISvnNC3pi")
Definition: Nus18Systs.h:316
const Nus18SystFromHist kNus18RHCEXNeutronSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","Neutron18","Neutron")
Definition: Nus18Systs.h:128
const Nus18FlatSyst kNus18FHCNDTauSyst("TauScale18","Tau Scale (ND)","ND","FHC", 0., 0.)
Definition: Nus18Systs.h:132
Nus group systematic applied as a flat shift only dependent on NC or CC.
Definition: Nus18Systs.h:68
const Nus18SystFromHist kNus18RHCFDAccSyst
std::string fFileName
Definition: Nus18Systs.h:53
const Nus18BaseSyst kNus18BaseBeamSyst("Beam18","All Beam")
Definition: Nus18Systs.h:107
const Nus18SystFromHist kNus18FHCNDGENIEMECEShapeNuSyst
const Nus18SystFromHist kNus18FHCFDGENIEGenSyst(kNusAna18FHCSystFile,"FD","FHC","GENIEg18","GENIE General Systs")
Definition: Nus18Systs.h:173
std::string fPolarity
Definition: Nus18Systs.h:82
const Nus18SystFromHist kNus18RHCNDGENIERadCorrNueSyst(kNusAna18RHCSystFile,"ND","RHC","RadCorrNue18","GENIE RadCorrNue")
Definition: Nus18Systs.h:353
const Nus18SystFromHist kNus18FHCNDXSecOnOffSyst(kNusAna18FHCSystFile,"ND","FHC","XSecOnOff18","Cross-section On/Off")
Definition: Nus18Systs.h:283
const Nus18SystFromHist kNus18RHCEXGENIEDISvbarnNC3piSyst
const Nus18SystFromHist kNus18RHCEXXSecOnOffSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","XSecOnOff18","Cross-section On/Off")
Definition: Nus18Systs.h:287
std::vector< const ISyst * > GetNus18RHCExtrapSysts()
Definition: Nus18Systs.cxx:727
const Nus18SystFromHist kNus18FHCFDGENIER0COHpiSyst(kNusAna18FHCSystFile,"FD","FHC","R0COHpi18","GENIE R0COHpi")
Definition: Nus18Systs.h:312
const Nus18SystFromHist kNus18RHCFDNCSyst
const Nus18BaseSyst kNus18BaseNueUnivSyst("NueNDUniv18","Nue ND")
Definition: Nus18Systs.h:219
const Nus18SystFromHist kNus18FHCFDCalibShapeSyst(kNusAna18FHCSystFile,"FD","FHC","CalShape18","Calibration shape")
Definition: Nus18Systs.h:165
const Nus18BaseSyst kNus18BaseGENIEDISvbarpNC1piSyst("DISvbarpNC1pi18","GENIE DISvbarpNC1pi")
Definition: Nus18Systs.h:300
NusChannel
Definition: NusSysts.h:13
const Nus18SystFromHist kNus18FHCFDNumuCCSyst
const Nus18SystFromHist kNus18RHCFDGENIEBhBYSyst
const Nus18SystFromHist kNus18FHCNDGENIEMvCCRESSyst(kNusAna18FHCSystFile,"ND","FHC","MvCCRES18","GENIE MvCCRES")
Definition: Nus18Systs.h:383
const Nus18SystFromHist kNus18RHCNDGENIEDISvnCC0piSyst(kNusAna18RHCSystFile,"ND","RHC","DISvnCC0pi18","GENIE DISvnCC0pi")
Definition: Nus18Systs.h:352
const Nus18SystFromHist kNus18FHCNDGENIEBhBYSyst(kNusAna18FHCSystFile,"ND","FHC","BhtBY18","GENIE BhtBY")
Definition: Nus18Systs.h:397
const Nus18BaseSyst kNus18BaseGENIEMaNCRESSyst("MaNCRES18","GENIE MaNCRES")
Definition: Nus18Systs.h:290
const Nus18BaseSyst kNus18BasePPFXSyst("PPFX18","PPFX syst")
Definition: Nus18Systs.h:258
const Nus18BaseSyst kNus18BaseGENIEFormZoneSyst("FormZone18","GENIE FormZone")
Definition: Nus18Systs.h:291
const Nus18SystFromHist kNus18FHCFDKaonSyst(kNusAna18FHCSystFile,"FD","FHC","Kaon18","Kaon")
Definition: Nus18Systs.h:117
const Nus18SystFromHist kNus18RHCNDCherenkovSyst(kNusAna18RHCSystFile,"ND","RHC","Cherenkov18","Cherenkov (ND)")
Definition: Nus18Systs.h:142
const Nus18SystFromHist kNus18RHCEXLightLSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","CalLight18","Light Level syst")
Definition: Nus18Systs.h:192
const Nus18SystFromHist kNus18RHCNDAccSyst
const Nus18BaseSyst kNus18BaseCalibFlatSyst("CalFlat18","Flat miscalibration")
Definition: Nus18Systs.h:147
const Nus18SystFromHist kNus18RHCFDGENIEDISvbarpCC0piSyst(kNusAna18RHCSystFile,"FD","RHC","DISvbarpCC0pi18","GENIE DISvbarpCC0pi")
Definition: Nus18Systs.h:344
const Nus18FlatSyst kNus18FHCNDPPFXSyst("PPFX18","Fluxweight ND","ND","FHC", 8.56, 8.78)
Definition: Nus18Systs.h:259
const Nus18BaseSyst kNus18BaseLightLSyst("CalLight18","Light level syst")
Definition: Nus18Systs.h:187
const Nus18SystFromHist kNus18RHCEXAccSyst
const Nus18SystFromHist kNus18FHCNDGENIEFormZoneSyst(kNusAna18FHCSystFile,"ND","FHC","FormZone18","GENIE FormZone")
Definition: Nus18Systs.h:322
const Nus18SystFromHist kNus18RHCFDTauSyst(kNusAna18RHCSystFile,"FD","RHC","TauScale18","Tau scale")
Definition: Nus18Systs.h:135
const Nus18SystFromHist kNus18FHCFDNueCCSyst
void Shift(double sigma, caf::SRProxy *sr, double &weight) const override
Perform the systematic shift.
Definition: Nus18Systs.cxx:337
const Nus18FlatSyst kNus18RHCNormSyst
const Nus18SystFromHist kNus18FHCFDBeamSyst(kNusAna18FHCSystFile,"FD","FHC","Beam18","All Beam")
Definition: Nus18Systs.h:109
const Nus18FlatSyst kNus18FHCNDNormNDSyst("NormND18","Normalisation ND","ND","FHC", 1.5, 1.5)
Definition: Nus18Systs.h:267
NusChannel GetNusChannel(caf::SRProxy *sr) const
Definition: Nus18Systs.cxx:293
const Nus18SystFromHist kNus18RHCFDNumuCCSyst
const Nus18SystFromHist kNus18FHCNDGENIEDISvpNC3piSyst(kNusAna18FHCSystFile,"ND","FHC","DISvpNC3pi18","GENIE DISvpNC3pi")
Definition: Nus18Systs.h:332
const Nus18SystFromHist kNus18RHCFDNueCCSyst
const Nus18BaseSyst kNus18BaseTauUnivSyst("TauNDUniv18","Tau ND")
Definition: Nus18Systs.h:227
const Nus18SystFromHist kNus18RHCEXCalibFlatSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","CalFlat18","Flat Miscalibration")
Definition: Nus18Systs.h:152
const Nus18SystFromHist kNus18RHCNDGENIEMvCCRESSyst(kNusAna18RHCSystFile,"ND","RHC","MvCCRES18","GENIE MvCCRES")
Definition: Nus18Systs.h:381
const Nus18SystFromHist kNus18RHCFDGENIEDISvnCC0piSyst(kNusAna18RHCSystFile,"FD","RHC","DISvnCC0pi18","GENIE DISvnCC0pi")
Definition: Nus18Systs.h:340
const Nus18BaseSyst kNus18BaseOscPUnivSyst("FDOscPUniv18","FD oscillation parameters")
Definition: Nus18Systs.h:235
const Nus18SystFromHist kNus18RHCEXOscPUnivSyst
const Nus18SystFromHist kNus18RHCNDGENIEDISvbarpNC1piSyst(kNusAna18RHCSystFile,"ND","RHC","DISvbarpNC1pi18","GENIE DISvbarpNC1pi")
Definition: Nus18Systs.h:355
const Nus18SystFromHist kNus18FHCNDGENIEMaCOHpiSyst(kNusAna18FHCSystFile,"ND","FHC","MaCOHpi18","GENIE MaCOHpi")
Definition: Nus18Systs.h:323
const Nus18BaseSyst kNus18BaseGENIEMvCCRESSyst("MvCCRES18","GENIE MvCCRES")
Definition: Nus18Systs.h:372
const Nus18SystFromHist kNus18FHCFDGENIEDISvpNC0piSyst(kNusAna18FHCSystFile,"FD","FHC","DISvpNC0pi18","GENIE DISvpNC0pi")
Definition: Nus18Systs.h:318
const Nus18SystFromHist kNus18FHCFDGENIEDISvnNC1piSyst(kNusAna18FHCSystFile,"FD","FHC","DISvnNC1pi18","GENIE DISvnNC1pi")
Definition: Nus18Systs.h:314
const Nus18SystFromHist kNus18RHCNDNCSyst
const Nus18SystFromHist kNus18RHCFDKaonSyst(kNusAna18RHCSystFile,"FD","RHC","Kaon18","Kaon")
Definition: Nus18Systs.h:119
Prevent histograms being added to the current directory.
Definition: UtilsExt.h:46
const Nus18SystFromHist kNus18RHCNDCalibFlatSyst(kNusAna18RHCSystFile,"ND","RHC","CalFlat18","Flat Miscalibration (ND)")
Definition: Nus18Systs.h:150
const Nus18FlatSyst kNus18RHCNDPPFXSyst("PPFX18","Fluxweight ND","ND","RHC", 11.7, 11.4)
Definition: Nus18Systs.h:261
const Nus18SystFromHist kNus18FHCNDGENIER0COHpiSyst(kNusAna18FHCSystFile,"ND","FHC","R0COHpi18","GENIE R0COHpi")
Definition: Nus18Systs.h:324
const Nus18SystFromHist kNus18RHCFDLightLSyst(kNusAna18RHCSystFile,"FD","RHC","CalLight18","Light Level syst")
Definition: Nus18Systs.h:191
const Nus18SystFromHist kNus18RHCNDBeamSyst(kNusAna18RHCSystFile,"ND","RHC","Beam18","All Beam (ND)")
Definition: Nus18Systs.h:110
std::vector< const ISyst * > GetNus18Systs(bool rhc, std::string det_type, std::string syst_type)
Definition: Nus18Systs.cxx:439
const Nus18SystFromHist kNus18RHCFDXSecOnOffSyst(kNusAna18RHCSystFile,"FD","RHC","XSecOnOff18","Cross-section On/Off")
Definition: Nus18Systs.h:286
const Nus18SystFromHist kNus18FHCNDGENIEMANCELSyst(kNusAna18FHCSystFile,"ND","FHC","MANCEL18","GENIE MANCEL")
Definition: Nus18Systs.h:325
const Nus18SystFromHist kNus18RHCNDGENIEMaCCRESSyst(kNusAna18RHCSystFile,"ND","RHC","MaCCRES18","GENIE MaCCRES")
Definition: Nus18Systs.h:377
const Nus18SystFromHist kNus18FHCNDCalibShapeSyst(kNusAna18FHCSystFile,"ND","FHC","CalShape18","Calibration shape (ND)")
Definition: Nus18Systs.h:164
const Nus18SystFromHist kNus18FHCFDGENIEBhBYSyst(kNusAna18FHCSystFile,"FD","FHC","BhtBY18","GENIE BhtBY")
Definition: Nus18Systs.h:398
const Nus18SystFromHist kNus18RHCNDGENIEGenSyst(kNusAna18RHCSystFile,"ND","RHC","GENIEg18","GENIE General Systs (ND)")
Definition: Nus18Systs.h:174
std::string fPolarity
Definition: Nus18Systs.h:55
const Nus18SystFromHist kNus18RHCFDCherenkovSyst(kNusAna18RHCSystFile,"FD","RHC","Cherenkov18","Cherenkov")
Definition: Nus18Systs.h:143
const Nus18SystFromHist kNus18RHCEXKaonSyst(kNusAna18RHCExtrapSystFile,"EX","RHC","Kaon18","Kaon")
Definition: Nus18Systs.h:120
const Nus18SystFromHist kNus18RHCEXGENIER0COHpiSyst
const Nus18SystFromHist kNus18FHCNDLightLSyst(kNusAna18FHCSystFile,"ND","FHC","CalLight18","Light Level syst (ND)")
Definition: Nus18Systs.h:188
const Nus18SystFromHist kNus18RHCFDGENIEFormZoneSyst(kNusAna18RHCSystFile,"FD","RHC","FormZone18","GENIE FormZone")
Definition: Nus18Systs.h:334
std::vector< const ISyst * > GetNus18FHCNDSysts(std::string syst_type)
Definition: Nus18Systs.cxx:534
const Nus18SystFromHist kNus18RHCNDGENIEMaNCRESSyst(kNusAna18RHCSystFile,"ND","RHC","MaNCRES18","GENIE MaNCRES")
Definition: Nus18Systs.h:345
const Nus18BaseSyst kNus18BaseGENIENOvASyst("GENIEn18","GENIE NOvA specific syst")
Definition: Nus18Systs.h:179
const Nus18SystFromHist kNus18RHCEXGENIEDISvbarpCC0piSyst
const Nus18SystFromHist kNus18RHCEXGENIENOvASyst(kNusAna18RHCExtrapSystFile,"EX","RHC","GENIEn18","GENIE NOvA Specific Systs")
Definition: Nus18Systs.h:184
const Nus18FlatSyst kNus18RHCNDTauSyst("TauScale18","Tau Scale (ND)","ND","RHC", 0., 0.)
Definition: Nus18Systs.h:134
const std::string kNusAna18RHCSystFile
Definition: Nus18Systs.h:98
const Nus18SystFromHist kNus18RHCNDGENIEDISvbarnNC3piSyst(kNusAna18RHCSystFile,"ND","RHC","DISvbarnNC3pi18","GENIE DISvbarnNC3pi")
Definition: Nus18Systs.h:350
const Nus18SystFromHist kNus18FHCNDGENIEGenSyst(kNusAna18FHCSystFile,"ND","FHC","GENIEg18","GENIE General Systs (ND)")
Definition: Nus18Systs.h:172
const Nus18BaseSyst kNus18BaseGENIEDISvnNC2piSyst("DISvnNC2pi18","GENIE DISvnNC2pi")
Definition: Nus18Systs.h:303
enum BeamMode string