-
Notifications
You must be signed in to change notification settings - Fork 0
/
config.txt
167 lines (128 loc) · 5.48 KB
/
config.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
//------------------------------------------------------------------------
// Configuration file for PhysicsBookPlots.C
//------------------------------------------------------------------------
// You should only have to modify this file.
// Input parameters go to the right of ":".
// Any number of spaces is fine, but please don't use tabs.
// Some inputs accept multiple parameters separated by ","'s.
// Send bugs/suggestions to [email protected]
// Authors: Tom Brooks, Rhiannon Jones, Ed Tyley, Natalie Wright
//------------------------------------------------------------------------
//------------------------------------------------------------------------
// File configuration
//------------------------------------------------------------------------
// Name of input file
InputFile: Trees/xsectree_v2rwt.root
// Name of model/tune in input file (1 per input file)
TuneName: GENIEv2_
// Name of output file (.root for modifiable canvas)
OutputFile: temp/2Dxsecv2.png
//------------------------------------------------------------------------
// Select neutrino interaction (commas for multiple)
//------------------------------------------------------------------------
// Muon ([-]14) or electron ([-]12) neutrino
NuPdg: 14
// Charged current (1), neutral current (0), not selected (-1, reco only)
IsCC: 1
// Require lepton to be contained (put "true, false" for contained and exiting)
ContainedLepton: true, false
// Require secondary particles to be contained
ContainedParticles: true, false
// Apply a fiducial volume cut to the neutrino vertex (-1 for AV)
// Cuts +ve values from face: min x, min y, min z, max x, max y, max z(, CPA cut, APA cut)
Fiducial: 10, 10, 15, 10, 20, 50, 5, 2.5
// FSI (true) or true interaction (false)
PlotByFsi: true
//- Determine FSI (ignored if plotting by interaction)
//- Number of protons (-1 for all)
NumProtons: -1
//- Number of charged pions
NumPiPM: -1
//- Number of neutral pions
NumPi0: -1
//- Choose true interaction (ignored if plotting by fsi)
//- QE (0), RES (1), DIS (2), COH (3), MEC (10)
InteractionType: -1
//------------------------------------------------------------------------
// Choose what to plot
//------------------------------------------------------------------------
// Choose what stage to plot
//- true (selection and kinematic variables all in truth)
//- smeareff (Applies efficiency and smears plotting variables only)
//- reco (Efficiency + smearing + PID + selection) TODO: only for numuCC
Stage: reco
// Variable to plot in
// Bools
//- lep_contained (containment of primary lepton)
//- particles_contained (containment of secondary particles)
//- cc (charged or neutral current)
// Integers
//- nu_pdg (neutrino PDG code)
//- int_type (interaction type, see above)
//- n_pr (number of protons)
//- n_pipm (number of charged pions)
//- n_pi0 (number of neutral pions)
// Doubles - can add cos_ to the front of angular variables
//- nu_energy (neutrino energy) [GeV]
//- lep_mom (lepton momentum) [GeV]
//- lep_theta (lepton angle wrt beam) [rad]
//- pr1_mom (leading proton momentum) [GeV]
//- pr1_theta (leading proton angle wrt beam) [rad]
//- lep_pr1_angle (lepton angle wrt leading proton) [rad]
//- pipm1_mom (leading charged pion momentum) [GeV]
//- pipm1_theta (leading charged pion angle wrt beam) [rad]
//- lep_pipm1_angle (angle between lepton and charged pion) [rad]
//- delta_pt (transverse kinematic variable) [GeV]
//- delta_alphat (transverse kinematic variable) [rad]
//- delta_phit (transverse kinematic variable) [rad]
PlotVariable: cos_lep_theta, lep_mom
// Save the plots for each variable
PlotSlices: false
// Scale to a certain POT (-1 for no scaling) (6.6e20 is nominal)
PotScale: 6.6e20
//------------------------------------------------------------------------
// Choose how to plot
//------------------------------------------------------------------------
// Produce a stacked histogram (true)
PlotStacked: false
// Stack in true FSI (fsi), true interaction (int), or neutrino type (nu)
StackBy: int
// Minimum bin (-1 for min tree value)
MinValue: -1, 0
// Maximum bin (-1 for max tree value)
MaxValue: 1, 2
// Number of bins (-1 for automatic)
NumBins: 25, 20
// Bin edges (separate edges with , and plot variables in [] brackets)
BinEdges: [-1], [-1]
// Maximum percentage error per bin (0->1, -1 for none) (merges bins)
MaxError: 0.01
// Plot cross-section (true) or rate (false)
PlotXSec: true
// Plot filled histograms (true)
PlotFilled: false
//------------------------------------------------------------------------
// Choose plot options
//------------------------------------------------------------------------
// Show selection information (true)
ShowInfo: false
// Show error band at bottom of plot (true)
ShowErrorBand: true
// Show percentage statistical error (true)
ShowStatError: true
// Show systematic error (true)
ShowSystError: true
// List of systematics to include (genie, flux, detector, background, constant)
Systematics: genie, flux, detector, background, constant
// Value of constant percentage systematic error (0->1)
ConstantSyst: 0.0224
// Show error bars on histograms (true)
ShowErrorBars: true
// Produce covariance and correlation matrices for systematic errors (true)
PlotCorrelation: true
// Plot all of the different universes for systematic uncertainties
PlotUniverses: true
// Produce efficiency and purity plots (true)
PlotEffPur: false
// Produce response matrix (true)
PlotResponse: true