NumuCCIncBins.h
Go to the documentation of this file.
1 #pragma once
2 
3 #include "CAFAna/Core/Binning.h"
4 
5 namespace ana
6 {
7 
8  // binning for scattering LL in remid
9  const Binning scatbins = Binning::Simple(100,-0.45, 0.45);
10 
11  // binning for dedx LL in remid
12  const Binning dEdXbins = Binning::Simple(100,-3.0, 1.0);
13 
14  // binning for energy per length
15  const Binning ebylbins = Binning::Simple(100,0,0.006);
16 
17  // binning for x or y coordinate
18  const Binning xybins = Binning::Simple(20,-200,200);
19 
20  // binning for z coordinate
21  const Binning zbins = Binning::Simple(80,0 ,1600);
22 
23  // binning for energy
24  const Binning ebins = Binning::Simple(50,0,5);
25 
26  // binning for muon transverse momentum
27  const Binning ptbins = Binning::Simple(30,0,1.5);
28 
29 
30  // binning for resolution
31  const Binning resbins = Binning::Simple(101,-1.005,1.005);
32 
33  // binning for zoom costheta
34  const Binning angzoombins = Binning::Simple(101,0.895,1.005);
35 
36  // binning for costheta
37  const Binning angbins = Binning::Simple(101,-1.005,1.005);
38 
39  // binning for efficiency
40  const Binning effbins = Binning::Simple(101,-0.005,1.005);
41 
42  // binning for energy resolution vs energy, 2d
43  const Binning resebins = Binning::Simple(ebins.NBins()*resbins.NBins(),
44  0,
45  ebins.NBins()*resbins.NBins());
46 
47  // binning for costheta resolution vs costheta, 2d
48  const Binning rescosthetabins = Binning::Simple(angzoombins.NBins()*resbins.NBins(),
49  0,
50  angzoombins.NBins()*resbins.NBins());
51 
52  // binning for remid
53  const Binning remidbins = Binning::Simple(81, -1/160., 1+1/160.);
54 
55  // binning for momentum vs costheta, 2d
56  const Binning pvsthetabins = Binning::Simple(ebins.NBins()*angbins.NBins(),
57  0, ebins.NBins()*angbins.NBins());
58 
59  // binning for pt vs pz, 2d
60  const Binning ptvspzbins = Binning::Simple(ebins.NBins()*ebins.NBins(),
61  0, ebins.NBins()*ebins.NBins());
62 
63 
64  //-----------------------------------------------------------------------------------
65  // Standard analysis binning
66 
67  // bin edges for custom binning in costheta, allowing unphysical and outside our phasespace
68  const std::vector<double> angedges {-1.0, 0.5, 0.56, 0.62, 0.68, 0.74, 0.80, 0.85, 0.88, 0.91, 0.94, 0.96, 0.98, 0.99, 1.0};
69  const std::vector<double> angedges_display {0.5, 0.56, 0.62, 0.68, 0.74, 0.80, 0.85, 0.88, 0.91, 0.94, 0.96, 0.98, 0.99, 1.0};
70 
71  // custom binning for costheta
74 
75  // bin edges for custom binning in muon KE, allowing unphysical and outside our phasespace
76  const std::vector<double> mukeedges {-10.0,
77  0.5, 0.6, 0.7, 0.8, 0.9,
78  1.0, 1.1, 1.2, 1.3, 1.4,
79  1.5, 1.6, 1.7, 1.8, 1.9,
80  2.0, 2.1, 2.2, 2.3, 2.4,
81  2.5, 120};
82  const std::vector<double> mukeedges_display {
83  0.5, 0.6, 0.7, 0.8, 0.9,
84  1.0, 1.1, 1.2, 1.3, 1.4,
85  1.5, 1.6, 1.7, 1.8, 1.9,
86  2.0, 2.1, 2.2, 2.3, 2.4,
87  2.5};
88 
89  // Custom binning for muon Kinetic Energy
92 
93  // bin edges for custom binning in neutrino energy
94  const std::vector<double> enuedges { -10,
95  0.0, 0.25, 0.50, 0.75,
96  1.0, 1.25, 1.50, 1.75,
97  2.0, 2.25, 2.50, 2.75,
98  3.0, 3.25, 3.50, 3.75,
99  4.0, 4.25, 4.50, 4.75,
100  5.0, 120};
101  const std::vector<double> enuedges_display {
102  0.0, 0.25, 0.50, 0.75,
103  1.0, 1.25, 1.50, 1.75,
104  2.0, 2.25, 2.50, 2.75,
105  3.0, 3.25, 3.50, 3.75,
106  4.0, 4.25, 4.50, 4.75,
107  5.0};
108  // Custom binning for Neutrino Energy
111 
112 
113  // bin edges for custom binning in neutrino energy
114  const std::vector<double> eavailedges {
115  0.0, 0.25, 0.50, 0.75,
116  1.0, 1.25, 1.50, 1.75,
117  2.0, 2.25, 2.50, 120};
118  const std::vector<double> eavailedges_display {
119  0.0, 0.25, 0.50, 0.75,
120  1.0, 1.25, 1.50, 1.75,
121  2.0, 2.25, 2.50};
122  // Custom binning for Neutrino Energy
125 
126  // custom binning for costheta vs momentum, 2d
128  angbinsCustom.NBins(),
129  0,
130  mukebins.NBins()*
131  angbinsCustom.NBins());
132 
134  angbinsCustom.NBins()*
135  enubins.NBins(),
136  0,
137  mukebins.NBins()*
138  angbinsCustom.NBins()*
139  enubins.NBins());
140 
142  angbinsCustom.NBins()*
143  eavailbins.NBins(),
144  0,
145  mukebins.NBins()*
146  angbinsCustom.NBins()*
147  eavailbins.NBins());
148 
149  const std::vector<double> q2edges {
150  0.0, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.40, 0.50, 0.60, 0.75, 0.90, 1.10, 1.40, 1.80, 2.8, 120.
151  };
152 
154 
155  // Binning for Energy Estimation
162 }
const Binning kTrkLenCatBinning
const Binning angvsmukevseavailbins
const Binning rescosthetabins
Definition: NumuCCIncBins.h:48
const Binning effbins
Definition: NumuCCIncBins.h:40
Represent the binning of a Spectrum&#39;s x-axis.
Definition: Binning.h:16
Cuts and Vars for the 2020 FD DiF Study.
Definition: vars.h:6
const Binning eavailbinsDisplay
const Binning zbins
Definition: NumuCCIncBins.h:21
const Binning resbins
Definition: NumuCCIncBins.h:31
const Binning xybins
Definition: NumuCCIncBins.h:18
const std::vector< double > mukeedges_display
Definition: NumuCCIncBins.h:82
const Binning scatbins
Definition: NumuCCIncBins.h:9
const Binning ebins
Definition: NumuCCIncBins.h:24
const Binning eavailbins
const Binning remidbins
Definition: NumuCCIncBins.h:53
const Binning ebylbins
Definition: NumuCCIncBins.h:15
const Binning kHadEnergyBinning
const Binning q2bins
const Binning ptvspzbins
Definition: NumuCCIncBins.h:60
const Binning kMuonECatBinning
const Binning angbins
Definition: NumuCCIncBins.h:37
const Binning enubins
const std::vector< double > q2edges
const std::vector< double > angedges
Definition: NumuCCIncBins.h:68
const std::vector< double > eavailedges_display
const Binning angvsmukebins
const Binning angzoombins
Definition: NumuCCIncBins.h:34
const Binning resebins
Definition: NumuCCIncBins.h:43
const Binning kMuonEActBinning
const Binning mukebinsDisplay
Definition: NumuCCIncBins.h:91
static Binning Custom(const std::vector< double > &edges)
Definition: Binning.cxx:145
const std::vector< double > eavailedges
const Binning mukebins
Definition: NumuCCIncBins.h:90
const Binning kVisHadEBinning
const Binning dEdXbins
Definition: NumuCCIncBins.h:12
int NBins() const
Definition: Binning.h:29
const Binning angvsmukevsebins
const Binning angbinsDisplay
Definition: NumuCCIncBins.h:73
const std::vector< double > enuedges
Definition: NumuCCIncBins.h:94
const std::vector< double > mukeedges
Definition: NumuCCIncBins.h:76
const Binning enubinsDisplay
const Binning pvsthetabins
Definition: NumuCCIncBins.h:56
const Binning angbinsCustom
Definition: NumuCCIncBins.h:72
const std::vector< double > angedges_display
Definition: NumuCCIncBins.h:69
static Binning Simple(int n, double lo, double hi, const std::vector< std::string > &labels={})
Definition: Binning.cxx:107
const Binning ptbins
Definition: NumuCCIncBins.h:27
const Binning kTrkLenActBinning
const std::vector< double > enuedges_display