data_fit.fcl
Go to the documentation of this file.
1 #include "FitterPlotStyle.fcl"
2 #include "ProtoDUNECalibration.fcl"
3 #include "FitSystematics.fcl"
4 
5 #MCFileName: "/dune/data/users/calcuttj/Prod4a_6_8_21_with_data/eventSelection_mc_all.root"
6 #DataFileName: "/dune/data/users/calcuttj/Prod4a_6_8_21_with_data/eventSelection_data_BeamQuality.root"
7 #MCFileName: "/dune/data/users/calcuttj/Prod4a_6_29_21_with_data_5809/eventSelection_mc_all.root"
8 #DataFileName: "/dune/data/users/calcuttj/Prod4a_6_29_21_with_data_5809/eventSelection_data_BeamQuality.root"
9 MCFileName: "/dune/data/users/calcuttj/Prod4a_7_12_21_with_data_58XX/eventSelection_mc_all.root"
10 DataFileName: "/dune/data/users/calcuttj/Prod4a_7_12_21_with_data_58XX/eventSelection_data_BeamQuality.root"
11 TreeName: "pduneana/beamana"
12 
13 DriverName: "protoana::AbsCexDriver"
14 AnalysisOptions: {
15  #FakeDataRoutine: "SampleScales"
16  #FakeDataRoutine: "BinnedScales"
17  #FakeDataRoutine: "G4RW"
18  #FakeDataRoutine: "dEdX"
19  #FakeDataRoutine: "EffVar"
20  #FakeDataRoutine: "G4RWGrid"
21  #FakeDataRoutine: "PionAngle"
22  FakeDataRoutine: "AngleVar"
23  #FakeDataRoutine: "Asimov"
24  #FakeDataRoutine: "Toy"
25 
26  FakeDataScales: [
27  #[1, 2.]
28  [2, 2.]
29  ]
30 
31  FakeDataBinnedScales: [
32  [2, [1., 1.488, 1.394, 1.469, 1.]],
33  [1, [1., .748, .782, .832, .912, 1.]]
34  ]
35 
36  FakeDataG4RW: {
37  Position: 2
38  Shift: 1
39  Full: true
40  }
41 
42  FakeDataG4RWGrid: {
43  Position: [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 ]
44  Shift: [12, 12, 12, 12, 12, 12, 12, 08, 08, 08, 08, 08, 08, 08 ]
45  Branch: "g4rw_full_grid_piplus_weights_fake_data"
46 
47  ###(reac, 0), (abs, 1), (cex, 2), (dcex, 3), (prod, 4), (inel, 5)
48  ### x7 each (200 MeV bins) -- total range: 0, 41 (0, 1400)
49  }
50 
51  FakeDatadEdX: {
52  Cal_set: @local::CalorimetryParameters_SPProd4_MC_SCE
53  VariedCCal: 1.051e-3
54  }
55 
56  FakeDataEffVar: {
57  F: .333
58  }
59 
60  FakeDataPionAngle: {
61  #RatioFile: "angle_weight.root"
62  #RatioFile: "flat_pion_weights.root"
63  RatioFile: "half_shift.root"
64  Limits: [0., .4, .6, .8, 1.]
65  RatioNames: ["r1", "r2", "r3", "r4", "r5"]
66  }
67 
68  FakeDataAngleVar: {
69  #RatioFile: "angle_weight.root"
70  RatioFile: "gaus_pi0_weights_080.root"
71  #Limits: [0., .4, .6, .8, 1.]
72  Limits: [
73  [1, [0., 400., 500., 600., 800., 1000.]],
74  [2, [0., 450., 600., 800., 1200.]],
75  [3, [0., 400., 500. , 600., 700., 800., 900., 1000.]]
76  ]
77 
78  RatioNames: [
79  [1, ["r1", "r2", "r3", "r4", "r5", "r6"]],
80  [2, ["r1", "r2", "r3", "r4", "r5"]],
81  [3, ["r1", "r2", "r3", "r4", "r5", "r6", "r7", "r8"]]
82  ]
83  PDG: 111
84  }
85 
86  EndZCut: 222.10561
87  SliceCut: 464
88  WirePitch: .47974
89  Z0: 0.56035
90 
91  SliceMethod: "Traj" #E, Alt, Default
92 
93  TrajZStart: -.49375
94 
95  #Fixing reco
96  DoEnergyFix: true
97  EnergyFix: 80.
98 
99  EventSelection: {
100 
101  }
102 }
103 
104 Systematics: [
105  @local::SystematicsList.dEdX_Cal_Spline,
106  @local::G4RWSysts.g4rw_QE1,
107  @local::G4RWSysts.g4rw_QE2,
108  @local::G4RWSysts.g4rw_Prod,
109  @local::G4RWSysts.g4rw_DCex,
110  @local::G4RWSysts.g4rw_proton,
111  @local::SystematicsList.beam_shift,
112  @local::SystematicsList.ediv_weight,
113  @local::SystematicsList.no_track_weight,
114  @local::SystematicsList.beam_cut_weight
115 ]
116 
117 AddSystTerm: true
118 CovarianceBins: [
119  ["dEdX_Cal_Spline", 0],
120  ["g4rw_QE1", 1],
121  ["g4rw_QE2", 2],
122  ["g4rw_Prod", 3],
123  ["g4rw_DCex", 4],
124  ["g4rw_proton", 5],
125  ["ediv_weight", 6],
126  ["beam_shift", 7],
127  ["no_track_weight", 8],
128  ["beam_cut_weight", 9]
129 ]
130 CovarianceFile: "full_with_dcex_beamcut.root"
131 CovarianceMatrix: "m"
132 
133 FixVariables: true
134 SystsToFix: [
135  ["beam_shift", 0.]
136 ]
137 
138 Selections: [
139  {
140  Name: "Abs"
141  ID: 1
142  #RecoBins: [[-5000., 0., 200., 400., 600., 800., 1200.]]
143  #BinLabels: [["< 0.", "0 - 200", "200 - 400", "400 - 600", "600 - 800", "800 - 1200"]]
144  #RecoBins: [[0., 150., 300., 350., 400., 450., 500., 550., 600., 650., 700., 750., 800., 850., 900., 950., 1000., 1200.]]
145  #BinLabels: [["0 - 150", "150 - 300", "300 - 350", "350 - 400", "400 - 450", "450 - 500",
146  # "500 - 550", "550 - 600", "600 - 650", "650 - 700", "700 - 750",
147  # "750 - 800", "800 - 850", "850 - 900", "900 - 950", "950 - 1000", "1000 - 1200"]]
148  RecoBins: [[0., 400., 500., 600., 800., 1200.]]
149  BinLabels: [["0 - 400", "400 - 500", "500 - 600", "600 - 800", "800 - 1200"]]
150  AxisTitles: ["Reconstructed KE (MeV)"]
151  },
152  {
153  Name: "Cex"
154  ID: 2
155  #RecoBins: [[-5000., 0., 200., 400., 600., 800., 1200.]]
156  #BinLabels: [["< 0.", "0 - 200", "200 - 400", "400 - 600", "600 - 800", "800 - 1200"]]
157  #RecoBins: [[0., 300., 350., 400., 450., 500., 550., 600., 650., 700., 750., 800., 850., 900., 950., 1000., 1200.]]
158  #BinLabels: [["0 - 300", "300 - 350", "350 - 400", "400 - 450", "450 - 500",
159  # "500 - 550", "550 - 600", "600 - 650", "650 - 700", "700 - 750",
160  # "750 - 800", "800 - 850", "850 - 900", "900 - 950", "950 - 1000", "1000 - 1200"]]
161  RecoBins: [[0., 500., 600., 800., 1200.]]
162  BinLabels: [["0 - 500", "500 - 600", "600 - 800", "800 - 1200"]]
163  AxisTitles: ["Reconstructed KE (MeV)"]
164  } ,
165  {
166  Name: "RejectedInt"
167  ID: 3
168  #RecoBins: [[0., 300., 350., 400., 450., 500., 550., 600., 650., 700., 750., 800., 850., 900., 950., 1000., 1200.]]
169  #BinLabels: [["0 - 300", "300 - 350", "350 - 400", "400 - 450", "450 - 500",
170  # "500 - 550", "550 - 600", "600 - 650", "650 - 700", "700 - 750",
171  # "750 - 800", "800 - 850", "850 - 900", "900 - 950", "950 - 1000", "1000 - 1200"]]
172  RecoBins: [[0., 400., 500., 600., 800., 1200.]]
173  BinLabels: [["0 - 400", "400 - 500", "500 - 600", "600 - 800", "800 - 1200"]]
174  AxisTitles: ["Reconstructed KE (MeV)"]
175  },
176  {
177  Name: "APA2"
178  ID: 4
179  #RecoBins: [[222, 275, 325, 375, 425, 580]]
180  #RecoBins: [[222, 275, 580]]
181  RecoBins: [[222, 234, 580]]
182  AxisTitles: ["Reconstructed End Z (cm)"]
183  },
184  {
185  Name: "FailedBeamCuts"
186  ID: 5
187  RecoBins: [[0, 1]]
188  AxisTitles: [""]
189  },
190  {
191  Name: "NoBeamTrack"
192  ID: 6
193  RecoBins: [[0, 1]]
194  AxisTitles: [""]
195  }
196 ]
197 
198 
199 IncidentRecoBins: [-5000., 0., 200., 400., 600., 800., 1000., 1200.]
200 SelectedRecoBins: [-5000., 0., 200., 400., 600., 800., 1200.]
201 
202 FluxTypes: [
203  [2, "Muons"],
204  [1, "Pions"]
205  ]
206 
207 
208 #DefaultSignalBins: [450., 600., 800, 1200.]
209 DefaultSignalBins: [0., 400., 800., 1000, 1200.]
210 TrueIncidentBins: @local::DefaultSignalBins
211 
212 BeamEnergyBins: [0., 800., 900., 1000., 1100., 1200., 4000.]
213 #BeamEnergyBins: [0., 800., 1200., 4000.]
214 #BeamEnergyBins: [0., 4000.]
215 
216 #For use identifying the true incident samples
217 #used to make the final total pion incident histogram
218 IncidentSamples: [1, 2, 3, 6, 7]
219 MeasurementSamples: [1, 2]
220 
221 Samples: [
222  {
223  Name: "Abs"
224  ID: 1
225  IsSignal: true
226  SignalBins: [400., 500., 600., 800., 1000.]
227  #SignalBins: [200., 500., 600., 800., 1200.]
228  #SignalBins: [000., 50.,
229  # 100., 150.,
230  # 200., 250.,
231  # 300., 350.,
232  # 400., 450.,
233  # 500., 550.,
234  # 600., 650.,
235  # 700., 750.,
236  # 800., 850.,
237  # 900., 950.,
238  # 1000.]
239  FluxType: 1
240  },
241  {
242  Name: "Cex"
243  ID: 2
244  IsSignal: true
245  SignalBins: [450., 600., 800., 1200.]
246  #SignalBins: [000., 50.,
247  # 100., 150.,
248  # 200., 250.,
249  # 300., 350.,
250  # 400., 450.,
251  # 500., 550.,
252  # 600., 650.,
253  # 700., 750.,
254  # 800., 850.,
255  # 900., 950.,
256  # 1000.]
257  FluxType: 1
258  },
259  {
260  Name: "OtherInel"
261  ID: 3
262  IsSignal: false
263  SignalBins: @local::DefaultSignalBins
264  FluxType: 1
265  },
266  {
267  Name: "UpstreamInt"
268  ID: 4
269  IsSignal: false
270  SignalBins: []
271  FluxType: 1
272  },
273  {
274  Name: "Muons"
275  ID: 5
276  IsSignal: false
277  SignalBins: []
278  FluxType: 2
279  },
280  {
281  Name: "PionPastFV"
282  ID: 6
283  IsSignal: false
284  SignalBins: []
285  FluxType: 1
286  },
287  {
288  Name: "Other"
289  ID: 7
290  IsSignal: false
291  SignalBins: []
292  FluxType: 1
293  }
294 ]
295 
296 ## Minimizer setup
297 #MaxCalls: 2000
298 NScanSteps: 100
299 Tolerance: 0.001
300 UpperLimit: 10.0
301 LowerLimit: 0.05
302 ReducedIncidentChi2: false
303 FitFlux: true
304 ScaleDataToNorm: false
305 DataNorm: 20500.
306 #RNGSeed: 0
307 
308 FitFunctionType: 2
309 #####################
310 
311 ## Plotting
312 #PlotStyle: @local::ReducedColorsStyle
313 PlotStyle: @local::DefaultColorsStyle
314 PlotRebinned: true
315 DrawXSecUnderflow: false
316 #####################
317 
318 ## Validation
319 RandomStart: false
320 SplitMC: false
321 DoFakeData: false
322 FluctuateStats: false
323 DoThrows: true
324 NThrows: 10 #00
325 DoScans: false
326 RunHesse: true
327 Do1DShifts: true
328 DoSysts: true
329 MaxRethrows: 1000
330 #####################
331 
332 FitType: "Normal"
333 #FitType: "None"
334 #FitType: "Pulls"
335 #FitType: "Toy"
336 
337 NPulls: 10
Generated by   doxygen 1.8.11