NuebarCCIncVars.h
Go to the documentation of this file.
1 #pragma once
2 
3 #include "CAFAna/Core/Var.h"
4 #include "CAFAna/Core/Cut.h"
5 #include "CAFAna/Core/HistAxis.h"
6 #include "CAFAna/Core/Utilities.h"
7 
8 #include "CAFAna/Cuts/SpillCuts.h"
9 #include "3FlavorAna/Vars/NueVars.h"
10 #include "3FlavorAna/Vars/NueVarsExtra.h"
11 
12 #include "NDAna/numucc_inc/NDXSecMuonPID.h"
13 #include "NDAna/nuebarcc_inc/NuebarCCIncBins.h"
14 
15 namespace nuebarccinc
16 {
17 
18  // some basic kinematic vars
19  // prongs
20  extern const ana::Var kNProngs;
21 
22  // slice hits
23  extern const ana::Var kNHits;
24 
25  // vertex
26  const ana::Var kRecoVtxX = ana::kVtxX;
27  const ana::Var kRecoVtxY = ana::kVtxY;
28  const ana::Var kRecoVtxZ = ana::kVtxZ;
29 
30  extern const ana::NuTruthVar kTrueVtxXST;
31  extern const ana::NuTruthVar kTrueVtxYST;
32  extern const ana::NuTruthVar kTrueVtxZST;
33 
34  const ana::Var kTrueVtxX = ana::VarFromNuTruthVar(kTrueVtxXST);
35  const ana::Var kTrueVtxY = ana::VarFromNuTruthVar(kTrueVtxYST);
36  const ana::Var kTrueVtxZ = ana::VarFromNuTruthVar(kTrueVtxZST);
37 
38  // prong distance from sides
39  const ana::Var kAllPngDistTop = ana::kDistAllTop;
40  const ana::Var kAllPngDistBottom = ana::kDistAllBottom;
41  const ana::Var kAllPngDistWest = ana::kDistAllWest;
42  const ana::Var kAllPngDistEast = ana::kDistAllEast;
43  const ana::Var kAllPngDistBack = ana::kDistAllBack;
44  const ana::Var kAllPngDistFront = ana::kDistAllFront;
45 
46 
47 
48  int true_electron_png(const caf::SRProxy * sr);
49  ////////////////////////////////////////////////////////////////////
50  //Analysis Variables
51  ////////////////////////////////////////////////////////////////////
52 
53  extern const ana::Var kMuonID;
54  extern const ana::Var kCVNNueID;
55 
56  const ana::Var kRecoNeutrinoE = ana::kNueEnergy2018; //2018 energy estimator
57 
58  extern const ana::Var kRecoElectronCosTheta;
59  extern const ana::Var kRecoElectronTheta;
60 
61  extern const ana::Var kTrueElectronRecoCosTheta;
62  extern const ana::Var kTrueElectronRecoTheta;
63 
64  //Electron Energy w/ Bias Correction
65  extern const ana::Var kElectronEUntuned;
66  extern const ana::Var kElectronEFHCTuned;
67  extern const ana::Var kElectronERHCTuned;
68  const ana::Var kElectronE = kElectronERHCTuned;
69  extern const ana::Var kTrueElectronRecoEUntuned;
70  extern const ana::Var kTrueElectronRecoEFHCTuned;
71  extern const ana::Var kTrueElectronRecoERHCTuned;
72  const ana::Var kTrueElectronRecoE = kTrueElectronRecoERHCTuned;
73 
74  const ana::Var kRecoElectronEVsCos = ana::Var2D(kRecoElectronCosTheta,
75  costhetabins,
76  kElectronE,
77  eelecbins);
78 
79  const ana::Var kRecoElectronEVsCosVsEnu = ana::Var3D(kRecoElectronCosTheta,
80  costhetabins,
81  kElectronE,
82  eelecbins,
83  kRecoNeutrinoE,
84  enubins);
85  extern const ana::Var kRecoQ2;
86  ////////////////////////////////////////////////////////////////////
87  //Prong CVN stuff
88  ////////////////////////////////////////////////////////////////////
89  extern const ana::Var kProngCVNElectronScore;
90  extern const ana::Var kProngCVNElectronE;
91  extern const ana::Var kProngCVNElectronCosTheta;
92  extern const ana::Var kProngCVNElectronTheta;
93 
94 
95  ////////////////////////////////////////////////////////////////////
96  //MC Truth Variables
97  ////////////////////////////////////////////////////////////////////
98  extern const ana::NuTruthVar kTrueQ2ST;
99  extern const ana::NuTruthVar kTrueElectronEST;
100  extern const ana::NuTruthVar kTrueElectronCosThetaST;
101  extern const ana::NuTruthVar kTrueElectronThetaST;
102 
103  extern const ana::NuTruthVar kTrueNeutrinoEST;
104  extern const ana::NuTruthVar kTrueHadronicYST;
105  const ana::NuTruthVar kTrueHadEST = kTrueNeutrinoEST * kTrueHadronicYST;
106  const ana::NuTruthVar kTrueEMEST = kTrueNeutrinoEST - kTrueHadEST;
107 
108  const ana::Var kTrueQ2 = ana::VarFromNuTruthVar(kTrueQ2ST);
109  const ana::Var kTrueElectronE = ana::VarFromNuTruthVar(kTrueElectronEST);
110  const ana::Var kTrueNeutrinoE = ana::VarFromNuTruthVar(kTrueNeutrinoEST);
111  const ana::Var kTrueElectronCosTheta = ana::VarFromNuTruthVar(kTrueElectronCosThetaST);
112  const ana::Var kTrueElectronTheta = ana::VarFromNuTruthVar(kTrueElectronThetaST);
113 
114  const ana::NuTruthVar kTrueElectronEVsCosST = ana::Var2D(kTrueElectronCosThetaST,
115  costhetabins,
116  kTrueElectronEST,
117  eelecbins);
118 
119  const ana::Var kTrueElectronEVsCos = ana::VarFromNuTruthVar(kTrueElectronEVsCosST);
120 
121  const ana::NuTruthVar kTrueElectronEVsCosVsEnuST = ana::Var3D(kTrueElectronCosThetaST,
122  costhetabins,
123  kTrueElectronEST,
124  eelecbins,
125  kTrueNeutrinoEST,
126  enubins);
127 
128 
129  ////////////////////////////////////////////////////////////////////
130  //Resolution Variables
131  ////////////////////////////////////////////////////////////////////
132  const ana::Var kFracResNeutrinoE = ( kRecoNeutrinoE - kTrueNeutrinoE ) / kTrueNeutrinoE;
133  const ana::Var kFracResElectronE = ( kElectronE - kTrueElectronE ) / kTrueElectronE;
134  const ana::Var kFracResElectronCosTheta = ( kRecoElectronCosTheta - kTrueElectronCosTheta ) / kTrueElectronCosTheta;
135  const ana::Var kFracResElectronTheta = ( kRecoElectronTheta - kTrueElectronTheta ) / kTrueElectronTheta;
136 
137  const ana::Var kFracResTrueElectronE = ( kTrueElectronRecoE - kTrueElectronE ) / kTrueElectronE;
138  const ana::Var kFracResTrueElectronCosTheta = ( kTrueElectronRecoCosTheta - kTrueElectronCosTheta ) / kTrueElectronCosTheta;
139  const ana::Var kFracResTrueElectronTheta = ( kTrueElectronRecoTheta - kTrueElectronTheta ) / kTrueElectronTheta;
140 
141  const ana::Var kFracResProngCVNElectronE = ( kProngCVNElectronE - kTrueElectronE ) / kTrueElectronE;
142  const ana::Var kFracResProngCVNElectronCosTheta = ( kProngCVNElectronCosTheta - kTrueElectronCosTheta ) / kTrueElectronCosTheta;
143  const ana::Var kFracResProngCVNElectronTheta = ( kProngCVNElectronTheta - kTrueElectronTheta ) / kTrueElectronTheta;
144 
145  const ana::Var kAbsResNeutrinoE = kRecoNeutrinoE - kTrueNeutrinoE ;
146  const ana::Var kAbsResElectronE = kElectronE - kTrueElectronE ;
147  const ana::Var kAbsResElectronCosTheta = kRecoElectronCosTheta - kTrueElectronCosTheta ;
148  const ana::Var kAbsResElectronTheta = kRecoElectronTheta - kTrueElectronTheta ;
149 
150  const ana::Var kAbsResTrueElectronE = kTrueElectronRecoE - kTrueElectronE ;
151  const ana::Var kAbsResTrueElectronCosTheta = kTrueElectronRecoCosTheta - kTrueElectronCosTheta ;
152  const ana::Var kAbsResTrueElectronTheta = kTrueElectronRecoTheta - kTrueElectronTheta ;
153 
154  const ana::Var kAbsResProngCVNElectronE = kProngCVNElectronE - kTrueElectronE ;
155  const ana::Var kAbsResProngCVNElectronCosTheta = kProngCVNElectronCosTheta - kTrueElectronCosTheta;
156  const ana::Var kAbsResProngCVNElectronTheta = kProngCVNElectronTheta - kTrueElectronTheta ;
157 
158  ////////////////////////////////////////////////////////////////////
159  //Defining Hist Axis
160  ////////////////////////////////////////////////////////////////////
161  struct mHistAxisDef
162  {
163  std::string name;
164  ana::HistAxis axis;
165  };
166 
167  struct mHistAxisSTDef
168  {
169  std::string name;
170  ana::NuTruthHistAxis axis;
171  };
172  //Multidimension Hist Axes
173  // electron kinematics 2d
174  const ana::NuTruthHistAxis kTrueElectronEVsCosStandardAxisST("True E_{e} vs cos #theta; Elec_{e} (GeV); cos #{theta}", double_diff_bins, kTrueElectronEVsCosST);
175  const ana::HistAxis kTrueElectronEVsCosStandardAxis = ana::HistAxisFromNuTruthHistAxis(kTrueElectronEVsCosStandardAxisST);
176  const ana::HistAxis kRecoElectronEVsCosStandardAxis("Reconstructed E_{e} vs cos #theta; Elec_{e} (GeV); cos #{theta}", double_diff_bins, kRecoElectronEVsCos);
177 
178  // 3d analysis binning
179  const ana::NuTruthHistAxis kTrueElectronEVsCosVsEnuStandardAxisST("True E_{e} vs cos #{theta} vs Neutrino Energy (GeV)", three_d_bins,kTrueElectronEVsCosVsEnuST);
180  const ana::HistAxis kTrueElectronEVsCosVsEnuStandardAxis =
181  ana::HistAxisFromNuTruthHistAxis(kTrueElectronEVsCosVsEnuStandardAxisST);
182  const ana::HistAxis kRecoElectronEVsCosVsEnuStandardAxis("Reco. E_{e} vs cos #{theta} vs Neutrino Energy (GeV)", three_d_bins,kRecoElectronEVsCosVsEnu);
183 
184  // electron energy 1d
185  const ana::NuTruthHistAxis kTrueElectronEStandardAxisST("True E_{e} (GeV)", eelecbins, kTrueElectronEST);
186  const ana::HistAxis kTrueElectronEStandardAxis = ana::HistAxisFromNuTruthHistAxis(kTrueElectronEStandardAxisST);
187  const ana::HistAxis kRecoElectronEStandardAxis("Reco E_{e} (GeV)", eelecbins, kElectronE);
188 
189  // electron angle 1d
190  const ana::NuTruthHistAxis kTrueElectronCosThetaStandardAxisST("True cos #theta_{e}", costhetabins, kTrueElectronCosThetaST);
191  const ana::HistAxis kTrueElectronCosThetaStandardAxis = ana::HistAxisFromNuTruthHistAxis(kTrueElectronCosThetaStandardAxisST);
192  const ana::HistAxis kRecoElectronCosThetaStandardAxis("Reco cos #theta_{e}", costhetabins, kRecoElectronCosTheta);
193 
194  // neutrino energy
195  const ana::NuTruthHistAxis kTrueNeutrinoEnergyStandardAxisST("True Neutrino Energy (GeV)", enubins, kTrueNeutrinoEST);
196  const ana::HistAxis kTrueNeutrinoEnergyStandardAxis = ana::HistAxisFromNuTruthHistAxis(kTrueNeutrinoEnergyStandardAxisST);
197  const ana::HistAxis kRecoNeutrinoEnergyStandardAxis("Reconstructed Neutrino Energy (GeV)", enubins, kRecoNeutrinoE);
198 
199  // q2
200  const ana::NuTruthHistAxis kTrueQ2StandardAxisST("True Q^{2} (GeV)",
201  q2bins,
202  kTrueQ2ST);
203  const ana::HistAxis kTrueQ2StandardAxis = ana::HistAxisFromNuTruthHistAxis(kTrueQ2StandardAxisST);
204  const ana::HistAxis kRecoQ2StandardAxis("Reconstructed Q^{2} (GeV)",
205  q2bins,
206  kRecoQ2);
207 } //end of namespace
nuebarccinc::kAllPngDistBack
const ana::Var kAllPngDistBack
Definition: NuebarCCIncVars.h:43
nuebarccinc::kCVNNueID
const ana::Var kCVNNueID
Definition: NuebarCCIncVars.cxx:58
nuebarccinc::kAbsResElectronTheta
const ana::Var kAbsResElectronTheta
Definition: NuebarCCIncVars.h:148
nuebarccinc::kAllPngDistBottom
const ana::Var kAllPngDistBottom
Definition: NuebarCCIncVars.h:40
nuebarccinc::kNHits
const ana::Var kNHits([](const caf::SRProxy *sr){return sr->slc.nhit;})
Definition: NuebarCCIncVars.h:23
nuebarccinc::kAbsResTrueElectronE
const ana::Var kAbsResTrueElectronE
Definition: NuebarCCIncVars.h:150
nuebarccinc::three_d_bins
const ana::Binning three_d_bins
Definition: NuebarCCIncBins.h:46
ana::kDistAllBottom
const Var kDistAllBottom([](const caf::SRProxy *sr){if(std::isnan(1.*sr->sel.nuecosrej.distallpngbottom)) return-1000.0f;return float(sr->sel.nuecosrej.distallpngbottom);})
Distance of all showers in slice from the bottom edge of detector.
Definition: NueVars.h:33
nuebarccinc::kProngCVNElectronCosTheta
const ana::Var kProngCVNElectronCosTheta([](const caf::SRProxy *sr){if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;int nProngs=sr->vtx.elastic.fuzzyk.npng;float maxProngCVNe=-5.;TVector3 elecdir={0, 0, 0};for(int iProng=0;iProng< nProngs;iProng++){if(sr->vtx.elastic.fuzzyk.png[iProng].cvnpart.electronid > maxProngCVNe){maxProngCVNe=sr->vtx.elastic.fuzzyk.png[iProng].cvnpart.electronid;elecdir=(TVector3) sr->vtx.elastic.fuzzyk.png[iProng].dir;}}if(maxProngCVNe==-5) return-5.f;TVector3 beamdir=ana::NuMIBeamDirection(caf::kNEARDET);return(float) elecdir.Dot(beamdir);})
Definition: NuebarCCIncVars.h:91
nuebarccinc::eelecbins
const ana::Binning eelecbins
Definition: NuebarCCIncBins.h:19
nuebarccinc::kTrueElectronRecoERHCTuned
const ana::Var kTrueElectronRecoERHCTuned([](const caf::SRProxy *sr) ->float{ double p0=-0.0861793;double p1=0.235812;double p2=-0.0348964;float energy=-1000.f;if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;int true_electron_idx=true_electron_png(sr);if(true_electron_idx< 0) return-5.f;float x=(sr->vtx.elastic.fuzzyk.png[true_electron_idx].shwlid.shwE);if(x< 0) return 0.0f;if(x<=6) energy=x-(p0+x *p1+p2 *pow(x, 2));else energy=x;if(energy< 0) energy=x;return energy;})
Definition: NuebarCCIncVars.h:71
ana::kDistAllWest
const Var kDistAllWest([](const caf::SRProxy *sr){if(std::isnan(1.*sr->sel.nuecosrej.distallpngwest)) return-1000.0f;return float(sr->sel.nuecosrej.distallpngwest);})
Distance of all showers in slice from the west edge of detector.
Definition: NueVars.h:36
nuebarccinc::kTrueHadEST
const ana::NuTruthVar kTrueHadEST
Definition: NuebarCCIncVars.h:105
ana::kDistAllTop
const Var kDistAllTop([](const caf::SRProxy *sr){if(std::isnan(1.*sr->sel.nuecosrej.distallpngtop)) return-1000.0f;return float(sr->sel.nuecosrej.distallpngtop);})
Distance of all showers in slice from the top edge of detector.
Definition: NueVars.h:30
nuebarccinc::kTrueElectronEStandardAxis
const ana::HistAxis kTrueElectronEStandardAxis
Definition: NuebarCCIncVars.h:186
nuebarccinc::kFracResNeutrinoE
const ana::Var kFracResNeutrinoE
Definition: NuebarCCIncVars.h:132
nuebarccinc::kTrueEMEST
const ana::NuTruthVar kTrueEMEST
Definition: NuebarCCIncVars.h:106
nuebarccinc::kElectronERHCTuned
const ana::Var kElectronERHCTuned([](const caf::SRProxy *sr) ->float{ double p0=-0.0861793;double p1=0.235812;double p2=-0.0348964;float energy=-1000.f;if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;float x=(sr->vtx.elastic.fuzzyk.png[0].shwlid.shwE);if(x< 0) return 0.0f;if(x<=6) energy=x-(p0+x *p1+p2 *pow(x, 2));else energy=x;if(energy< 0) energy=x;return energy;})
Definition: NuebarCCIncVars.h:67
nuebarccinc::kTrueVtxYST
const ana::NuTruthVar kTrueVtxYST
Definition: NuebarCCIncVars.cxx:35
nuebarccinc::kAllPngDistFront
const ana::Var kAllPngDistFront
Definition: NuebarCCIncVars.h:44
nuebarccinc::costhetabins
const ana::Binning costhetabins
Definition: NuebarCCIncBins.h:23
ana::kDistAllBack
const Var kDistAllBack([](const caf::SRProxy *sr){if(std::isnan(1.*sr->sel.nuecosrej.distallpngback)) return-1000.0f;return float(sr->sel.nuecosrej.distallpngback);})
Distance of all showers in slice from the back edge of detector.
Definition: NueVars.h:42
nuebarccinc::kRecoElectronEVsCosVsEnu
const ana::Var kRecoElectronEVsCosVsEnu
Definition: NuebarCCIncVars.h:79
caf::Proxy< caf::StandardRecord >
Proxy for caf::StandardRecord.
Definition: SRProxy.h:2126
ana::Var2D
_Var< T > Var2D(const _Var< T > &a, const Binning &binsa, const _Var< T > &b, const Binning &binsb)
Variable formed from two input variables.
Definition: Var.cxx:247
nuebarccinc::kTrueElectronEVsCosVsEnuStandardAxisST
const ana::NuTruthHistAxis kTrueElectronEVsCosVsEnuStandardAxisST("True E_{e} vs cos #{theta} vs Neutrino Energy (GeV)", three_d_bins, kTrueElectronEVsCosVsEnuST)
nuebarccinc::kAbsResElectronE
const ana::Var kAbsResElectronE
Definition: NuebarCCIncVars.h:146
nuebarccinc::kAbsResProngCVNElectronE
const ana::Var kAbsResProngCVNElectronE
Definition: NuebarCCIncVars.h:154
nuebarccinc::kTrueElectronRecoCosTheta
const ana::Var kTrueElectronRecoCosTheta([](const caf::SRProxy *sr) ->float{if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;int true_electron_idx=true_electron_png(sr);if(true_electron_idx< 0) return-5.f;TVector3 elecdir=(TVector3) sr->vtx.elastic.fuzzyk.png[true_electron_idx].dir;TVector3 beamdir=ana::NuMIBeamDirection(caf::kNEARDET);return elecdir.Dot(beamdir);})
Definition: NuebarCCIncVars.h:61
nuebarccinc::kRecoQ2StandardAxis
const ana::HistAxis kRecoQ2StandardAxis("Reconstructed Q^{2} (GeV)", q2bins, kRecoQ2)
nuebarccinc::kTrueVtxXST
const ana::NuTruthVar kTrueVtxXST
Definition: NuebarCCIncVars.cxx:34
nuebarccinc::kAbsResTrueElectronTheta
const ana::Var kAbsResTrueElectronTheta
Definition: NuebarCCIncVars.h:152
nuebarccinc::kAbsResElectronCosTheta
const ana::Var kAbsResElectronCosTheta
Definition: NuebarCCIncVars.h:147
nuebarccinc::kTrueElectronEStandardAxisST
const ana::NuTruthHistAxis kTrueElectronEStandardAxisST("True E_{e} (GeV)", eelecbins, kTrueElectronEST)
ana::kVtxX
const Var kVtxX
Definition: NueVarsExtra.cxx:31
nuebarccinc::kTrueElectronEVsCosStandardAxis
const ana::HistAxis kTrueElectronEVsCosStandardAxis
Definition: NuebarCCIncVars.h:175
nuebarccinc::kRecoElectronEVsCos
const ana::Var kRecoElectronEVsCos
Definition: NuebarCCIncVars.h:74
nuebarccinc::kFracResElectronCosTheta
const ana::Var kFracResElectronCosTheta
Definition: NuebarCCIncVars.h:134
nuebarccinc::kRecoElectronCosThetaStandardAxis
const ana::HistAxis kRecoElectronCosThetaStandardAxis("Reco cos #theta_{e}", costhetabins, kRecoElectronCosTheta)
ana::Var3D
_Var< T > Var3D(const _Var< T > &a, const Binning &binsa, const _Var< T > &b, const Binning &binsb, const _Var< T > &c, const Binning &binsc)
This is just like a Var2D, but useful for 3D Spectra.
Definition: Var.cxx:273
nuebarccinc::kTrueElectronRecoE
const ana::Var kTrueElectronRecoE
Definition: NuebarCCIncVars.h:72
nuebarccinc::kAbsResNeutrinoE
const ana::Var kAbsResNeutrinoE
Definition: NuebarCCIncVars.h:145
nuebarccinc::kTrueNeutrinoEST
const ana::NuTruthVar kTrueNeutrinoEST
Definition: NuebarCCIncVars.cxx:350
nuebarccinc::kFracResTrueElectronCosTheta
const ana::Var kFracResTrueElectronCosTheta
Definition: NuebarCCIncVars.h:138
nuebarccinc::kTrueElectronEVsCosVsEnuStandardAxis
const ana::HistAxis kTrueElectronEVsCosVsEnuStandardAxis
Definition: NuebarCCIncVars.h:180
ana::kDistAllEast
const Var kDistAllEast([](const caf::SRProxy *sr){if(std::isnan(1.*sr->sel.nuecosrej.distallpngeast)) return-1000.0f;return float(sr->sel.nuecosrej.distallpngeast);})
Distance of all showers in slice from the east edge of detector.
Definition: NueVars.h:39
nuebarccinc::kTrueHadronicYST
const ana::NuTruthVar kTrueHadronicYST
Definition: NuebarCCIncVars.cxx:351
nuebarccinc::kFracResElectronTheta
const ana::Var kFracResElectronTheta
Definition: NuebarCCIncVars.h:135
nuebarccinc::kElectronE
const ana::Var kElectronE
Definition: NuebarCCIncVars.h:68
NDXSecMuonPID.h
nuebarccinc::double_diff_bins
const ana::Binning double_diff_bins
Definition: NuebarCCIncBins.h:39
nuebarccinc::kFracResProngCVNElectronTheta
const ana::Var kFracResProngCVNElectronTheta
Definition: NuebarCCIncVars.h:143
nuebarccinc::kTrueElectronCosThetaStandardAxis
const ana::HistAxis kTrueElectronCosThetaStandardAxis
Definition: NuebarCCIncVars.h:191
nuebarccinc::kTrueElectronCosTheta
const ana::Var kTrueElectronCosTheta
Definition: NuebarCCIncVars.h:111
nuebarccinc::kAllPngDistEast
const ana::Var kAllPngDistEast
Definition: NuebarCCIncVars.h:42
ana::kNueEnergy2018
const Var kNueEnergy2018([](const caf::SRProxy *sr){if(!sr->vtx.elastic.IsValid) return-1.0;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1.0;if(kIsRHC(sr)) return kNueEnergyRHC(sr);else return kNueEnergyFHC(sr);})
Definition: NueEnergy2018.h:25
nuebarccinc::kAbsResProngCVNElectronCosTheta
const ana::Var kAbsResProngCVNElectronCosTheta
Definition: NuebarCCIncVars.h:155
nuebarccinc::kTrueElectronEVsCos
const ana::Var kTrueElectronEVsCos
Definition: NuebarCCIncVars.h:119
nuebarccinc::kRecoElectronEVsCosVsEnuStandardAxis
const ana::HistAxis kRecoElectronEVsCosVsEnuStandardAxis("Reco. E_{e} vs cos #{theta} vs Neutrino Energy (GeV)", three_d_bins, kRecoElectronEVsCosVsEnu)
nuebarccinc::kRecoVtxY
const ana::Var kRecoVtxY
Definition: NuebarCCIncVars.h:27
nuebarccinc::kRecoElectronTheta
const ana::Var kRecoElectronTheta([](const caf::SRProxy *sr) ->float{if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;TVector3 elecdir=(TVector3) sr->vtx.elastic.fuzzyk.png[0].dir;TVector3 beamdir=ana::NuMIBeamDirection(caf::kNEARDET);float theta=elecdir.Dot(beamdir);return(float) TMath::ACos(theta)*180/TMath::Pi();})
Definition: NuebarCCIncVars.h:59
nuebarccinc::kFracResProngCVNElectronE
const ana::Var kFracResProngCVNElectronE
Definition: NuebarCCIncVars.h:141
nuebarccinc::kTrueNeutrinoEnergyStandardAxisST
const ana::NuTruthHistAxis kTrueNeutrinoEnergyStandardAxisST("True Neutrino Energy (GeV)", enubins, kTrueNeutrinoEST)
sr
caf::StandardRecord * sr
Definition: FillTrainingTrees.h:18
ana::HistAxisFromNuTruthHistAxis
HistAxis HistAxisFromNuTruthHistAxis(NuTruthHistAxis ntha, double _default)
Definition: HistAxis.cxx:9
nuebarccinc::kTrueElectronEVsCosVsEnuST
const ana::NuTruthVar kTrueElectronEVsCosVsEnuST
Definition: NuebarCCIncVars.h:121
nuebarccinc::kFracResTrueElectronTheta
const ana::Var kFracResTrueElectronTheta
Definition: NuebarCCIncVars.h:139
nuebarccinc::kAllPngDistTop
const ana::Var kAllPngDistTop
Definition: NuebarCCIncVars.h:39
nuebarccinc::kTrueElectronTheta
const ana::Var kTrueElectronTheta
Definition: NuebarCCIncVars.h:112
nuebarccinc::kAbsResProngCVNElectronTheta
const ana::Var kAbsResProngCVNElectronTheta
Definition: NuebarCCIncVars.h:156
nuebarccinc::kTrueElectronCosThetaStandardAxisST
const ana::NuTruthHistAxis kTrueElectronCosThetaStandardAxisST("True cos #theta_{e}", costhetabins, kTrueElectronCosThetaST)
nuebarccinc::kRecoElectronEStandardAxis
const ana::HistAxis kRecoElectronEStandardAxis("Reco E_{e} (GeV)", eelecbins, kElectronE)
nuebarccinc::mHistAxisSTDef::axis
ana::NuTruthHistAxis axis
Definition: NuebarCCIncVars.h:170
nuebarccinc::kTrueQ2StandardAxis
const ana::HistAxis kTrueQ2StandardAxis
Definition: NuebarCCIncVars.h:203
nuebarccinc::kFracResElectronE
const ana::Var kFracResElectronE
Definition: NuebarCCIncVars.h:133
nuebarccinc::enubins
const ana::Binning enubins
Definition: NuebarCCIncBins.h:14
nuebarccinc::kTrueElectronRecoEFHCTuned
const ana::Var kTrueElectronRecoEFHCTuned([](const caf::SRProxy *sr) ->float{double p0=-0.10588;double p1=0.57885;double p2=-0.408767;double p3=0.153795;double p4=-0.0270565;double p5=0.00171687;float energy=-1000.f;if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;int true_electron_idx=true_electron_png(sr);if(true_electron_idx< 0) return-5.f;float x=(sr->vtx.elastic.fuzzyk.png[true_electron_idx].shwlid.shwE);if(x< 0) return 0.0f;if(x<=6) energy=x-(p0+x *p1+p2 *pow(x, 2)+p3 *pow(x, 3)+p4 *pow(x, 4)+p5 *pow(x, 5));else energy=x;if(energy< 0) energy=x;return energy;})
Definition: NuebarCCIncVars.h:70
ana::VarFromNuTruthVar
Var VarFromNuTruthVar(const NuTruthVar &stv, double _default)
Definition: Var.cxx:7
ana::kDistAllFront
const Var kDistAllFront([](const caf::SRProxy *sr){if(std::isnan(1.*sr->sel.nuecosrej.distallpngfront)) return-1000.0f;return float(sr->sel.nuecosrej.distallpngfront);})
Distance of all showers in slice from the front edge of detector.
Definition: NueVars.h:45
nuebarccinc::kFracResProngCVNElectronCosTheta
const ana::Var kFracResProngCVNElectronCosTheta
Definition: NuebarCCIncVars.h:142
ana::kVtxY
const Var kVtxY
Definition: NueVarsExtra.cxx:32
ana::_HistAxis< Var >
nuebarccinc::kRecoNeutrinoEnergyStandardAxis
const ana::HistAxis kRecoNeutrinoEnergyStandardAxis("Reconstructed Neutrino Energy (GeV)", enubins, kRecoNeutrinoE)
nuebarccinc::kTrueQ2StandardAxisST
const ana::NuTruthHistAxis kTrueQ2StandardAxisST("True Q^{2} (GeV)", q2bins, kTrueQ2ST)
nuebarccinc::kTrueVtxY
const ana::Var kTrueVtxY
Definition: NuebarCCIncVars.h:35
nuebarccinc::kTrueVtxX
const ana::Var kTrueVtxX
Definition: NuebarCCIncVars.h:34
nuebarccinc::kProngCVNElectronTheta
const ana::Var kProngCVNElectronTheta([](const caf::SRProxy *sr){if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;int nProngs=sr->vtx.elastic.fuzzyk.npng;float maxProngCVNe=-5.;TVector3 elecdir={0, 0, 0};for(int iProng=0;iProng< nProngs;iProng++){if(sr->vtx.elastic.fuzzyk.png[iProng].cvnpart.electronid > maxProngCVNe){maxProngCVNe=sr->vtx.elastic.fuzzyk.png[iProng].cvnpart.electronid;elecdir=(TVector3) sr->vtx.elastic.fuzzyk.png[iProng].dir;}}if(maxProngCVNe==-5) return-5.f;TVector3 beamdir=ana::NuMIBeamDirection(caf::kNEARDET);float theta=elecdir.Dot(beamdir);return(float)(TMath::ACos(theta)*180/TMath::Pi());})
Definition: NuebarCCIncVars.h:92
nuebarccinc::kTrueNeutrinoE
const ana::Var kTrueNeutrinoE
Definition: NuebarCCIncVars.h:110
ana::kVtxZ
const Var kVtxZ
Definition: NueVarsExtra.cxx:33
nuebarccinc::kElectronEUntuned
const ana::Var kElectronEUntuned([](const caf::SRProxy *sr) ->float{if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;float x=(sr->vtx.elastic.fuzzyk.png[0].shwlid.shwE);return x;})
Definition: NuebarCCIncVars.h:65
nuebarccinc::kTrueElectronE
const ana::Var kTrueElectronE
Definition: NuebarCCIncVars.h:109
nuebarccinc::kAllPngDistWest
const ana::Var kAllPngDistWest
Definition: NuebarCCIncVars.h:41
nuebarccinc::kTrueVtxZST
const ana::NuTruthVar kTrueVtxZST
Definition: NuebarCCIncVars.cxx:36
nuebarccinc::kRecoVtxZ
const ana::Var kRecoVtxZ
Definition: NuebarCCIncVars.h:28
nuebarccinc::kNProngs
const ana::Var kNProngs([](const caf::SRProxy *sr) ->int{if(!sr->vtx.elastic.IsValid) return-5;return sr->vtx.elastic.fuzzyk.npng;})
Definition: NuebarCCIncVars.h:20
nuebarccinc::kProngCVNElectronScore
const ana::Var kProngCVNElectronScore([](const caf::SRProxy *sr){if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;int nProngs=sr->vtx.elastic.fuzzyk.npng;float maxProngCVNe=-5.;for(int iProng=0;iProng< nProngs;iProng++){if(sr->vtx.elastic.fuzzyk.png[iProng].cvnpart.electronid > maxProngCVNe){maxProngCVNe=sr->vtx.elastic.fuzzyk.png[iProng].cvnpart.electronid;}}return maxProngCVNe;})
Definition: NuebarCCIncVars.h:89
nuebarccinc::kTrueElectronRecoEUntuned
const ana::Var kTrueElectronRecoEUntuned([](const caf::SRProxy *sr) ->float{if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;int true_electron_idx=true_electron_png(sr);if(true_electron_idx< 0) return-5.f;float x=(sr->vtx.elastic.fuzzyk.png[true_electron_idx].shwlid.shwE);return x;})
Definition: NuebarCCIncVars.h:69
nuebarccinc::kProngCVNElectronE
const ana::Var kProngCVNElectronE([](const caf::SRProxy *sr){if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;int nProngs=sr->vtx.elastic.fuzzyk.npng;float maxProngCVNe=-5.;float maxProngCVNe_ElecE=-5.;for(int iProng=0;iProng< nProngs;iProng++){if(sr->vtx.elastic.fuzzyk.png[iProng].cvnpart.electronid > maxProngCVNe){maxProngCVNe=sr->vtx.elastic.fuzzyk.png[iProng].cvnpart.electronid;maxProngCVNe_ElecE=sr->vtx.elastic.fuzzyk.png[iProng].calE;}}return maxProngCVNe_ElecE;})
Definition: NuebarCCIncVars.h:90
nuebarccinc::true_electron_png
int true_electron_png(const caf::SRProxy *sr)
Definition: NuebarCCIncVars.cxx:38
nuebarccinc::kTrueVtxZ
const ana::Var kTrueVtxZ
Definition: NuebarCCIncVars.h:36
nuebarccinc::kElectronEFHCTuned
const ana::Var kElectronEFHCTuned([](const caf::SRProxy *sr) ->float{double p0=-0.10588;double p1=0.57885;double p2=-0.408767;double p3=0.153795;double p4=-0.0270565;double p5=0.00171687;float energy=-1000.f;if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;float x=(sr->vtx.elastic.fuzzyk.png[0].shwlid.shwE);if(x< 0) return 0.0f;if(x<=6) energy=x-(p0+x *p1+p2 *pow(x, 2)+p3 *pow(x, 3)+p4 *pow(x, 4)+p5 *pow(x, 5));else energy=x;if(energy< 0) energy=x;return energy;})
Definition: NuebarCCIncVars.h:66
ana::_Var< caf::SRProxy >
nuebarccinc::kAbsResTrueElectronCosTheta
const ana::Var kAbsResTrueElectronCosTheta
Definition: NuebarCCIncVars.h:151
nuebarccinc::q2bins
const ana::Binning q2bins
Definition: NuebarCCIncBins.h:28
nuebarccinc::kTrueNeutrinoEnergyStandardAxis
const ana::HistAxis kTrueNeutrinoEnergyStandardAxis
Definition: NuebarCCIncVars.h:196
nuebarccinc::kTrueElectronThetaST
const ana::NuTruthVar kTrueElectronThetaST([](const caf::SRNeutrinoProxy *nu){if(abs(nu->pdg)!=12|| !nu->iscc) return-5.0;int nprims=nu->prim.size();for(int iprim=0;iprim< nprims;iprim++){if(abs(nu->prim[iprim].pdg)==11){TVector3 edir=nu->prim[iprim].p.Vect();TVector3 beamdir=ana::NuMIBeamDirection(caf::kNEARDET);float theta=edir.Unit().Dot(beamdir.Unit());return(float) TMath::ACos(theta)*180/TMath::Pi();}}return-5.0;})
Definition: NuebarCCIncVars.h:101
nuebarccinc::kMuonID
const ana::Var kMuonID
Definition: NuebarCCIncVars.cxx:50
nuebarccinc::kRecoVtxX
const ana::Var kRecoVtxX
Definition: NuebarCCIncVars.h:26
nuebarccinc::kTrueQ2ST
const ana::NuTruthVar kTrueQ2ST
Definition: NuebarCCIncVars.cxx:353
nuebarccinc::kRecoElectronEVsCosStandardAxis
const ana::HistAxis kRecoElectronEVsCosStandardAxis("Reconstructed E_{e} vs cos #theta; Elec_{e} (GeV); cos #{theta}", double_diff_bins, kRecoElectronEVsCos)
nuebarccinc::kRecoNeutrinoE
const ana::Var kRecoNeutrinoE
Definition: NuebarCCIncVars.h:56
nuebarccinc::kTrueElectronEVsCosStandardAxisST
const ana::NuTruthHistAxis kTrueElectronEVsCosStandardAxisST("True E_{e} vs cos #theta; Elec_{e} (GeV); cos #{theta}", double_diff_bins, kTrueElectronEVsCosST)
nuebarccinc::kTrueElectronRecoTheta
const ana::Var kTrueElectronRecoTheta([](const caf::SRProxy *sr) ->float{if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;int true_electron_idx=true_electron_png(sr);if(true_electron_idx< 0) return-5.f;TVector3 elecdir=(TVector3) sr->vtx.elastic.fuzzyk.png[true_electron_idx].dir;TVector3 beamdir=ana::NuMIBeamDirection(caf::kNEARDET);float theta=elecdir.Dot(beamdir);return(float) TMath::ACos(theta)*180/TMath::Pi();})
Definition: NuebarCCIncVars.h:62
nuebarccinc::kRecoElectronCosTheta
const ana::Var kRecoElectronCosTheta([](const caf::SRProxy *sr) ->float{if(!sr->vtx.elastic.IsValid) return-1000.f;if(sr->vtx.elastic.fuzzyk.npng< 1) return-1000.f;TVector3 elecdir=(TVector3) sr->vtx.elastic.fuzzyk.png[0].dir;TVector3 beamdir=ana::NuMIBeamDirection(caf::kNEARDET);return elecdir.Dot(beamdir);})
Definition: NuebarCCIncVars.h:58
nuebarccinc::kTrueElectronEVsCosST
const ana::NuTruthVar kTrueElectronEVsCosST
Definition: NuebarCCIncVars.h:114
nuebarccinc::kTrueQ2
const ana::Var kTrueQ2
Definition: NuebarCCIncVars.h:108
nuebarccinc::kRecoQ2
const ana::Var kRecoQ2([](const caf::SRProxy *sr){const double M_e_sqrd=std::pow(0.000511, 2);double E_e=kElectronE(sr);double E_nu=kRecoNeutrinoE(sr);double cos_e=kRecoElectronCosTheta(sr);double p_e=std::sqrt(std::pow(E_e, 2)-M_e_sqrd);return 2 *E_nu *(E_e-p_e *cos_e)-M_e_sqrd;})
Definition: NuebarCCIncVars.h:85
nuebarccinc::kTrueElectronCosThetaST
const ana::NuTruthVar kTrueElectronCosThetaST([](const caf::SRNeutrinoProxy *nu){int nprims=nu->prim.size();for(int iprim=0;iprim< nprims;iprim++){if(abs(nu->prim[iprim].pdg)==11){TVector3 edir=nu->prim[iprim].p.Vect();TVector3 beamdir=ana::NuMIBeamDirection(caf::kNEARDET);return edir.Unit().Dot(beamdir.Unit());}}return-5.0;})
Definition: NuebarCCIncVars.h:100
nuebarccinc::kFracResTrueElectronE
const ana::Var kFracResTrueElectronE
Definition: NuebarCCIncVars.h:137
string
enum BeamMode string
nuebarccinc::kTrueElectronEST
const ana::NuTruthVar kTrueElectronEST([](const caf::SRNeutrinoProxy *nu){float ee=-5;int nprims=nu->prim.size();for(int iprim=0;iprim< nprims;iprim++){if(abs(nu->prim[iprim].pdg)==11){double E=nu->prim[iprim].p.E;ee=E;}}return ee;})
Definition: NuebarCCIncVars.h:99
Generated on Sun Jan 24 2021 14:12:47 for NOvA by   doxygen 1.8.11