LArSoft  v07_13_02
Liquid Argon Software toolkit - http://larsoft.org/
genie_eventweight.fcl
Go to the documentation of this file.
1 BEGIN_PROLOG
2 
3 genie_eventweight:
4 {
5  module_type: "EventWeight"
6 
7  #Use this to reweight one parameter at a time (to study the effect of varying one single parameter)
8  #weight_functions: [ genie_qema, genie_ncelAxial, genie_ncelEta, genie_qevec, genie_ccresAxial, genie_ccresVector, genie_ncresAxial, genie_ncresVector, genie_cohMA, genie_cohR0, genie_NonResRvp1pi, genie_NonResRvbarp1pi, genie_NonResRvp2pi, genie_NonResRvbarp2pi, genie_ResDecayGamma, genie_ResDecayEta, genie_ResDecayTheta, genie_NC, genie_DISAth, genie_DISBth, genie_DISCv1u, genie_DISCv2u, genie_AGKYxF, genie_AGKYpT, genie_FormZone, genie_FermiGasModelKf, genie_FermiGasModelSf, genie_IntraNukeNmfp, genie_IntraNukeNcex, genie_IntraNukeNel, genie_IntraNukeNinel, genie_IntraNukeNabs, genie_IntraNukeNpi, genie_IntraNukePImfp, genie_IntraNukePIcex, genie_IntraNukePIel, genie_IntraNukePIinel, genie_IntraNukePIabs, genie_IntraNukePIpi ]
9 
10  #This was used to use the same MaCCQE value as MiniBooNE best fit
11  weight_functions: [ genie_qema_miniboone ]
12 
13  #Use to reweight all GENIE parameters at the same time (see below for definition)
14  #weight_functions: [ genie_all ]
15 
16  #throw exception if weight outside limits
17  min_weight: 0
18  max_weight: 100
19 
20  genie_module_label: generator
21 
22 #########################################
23 #
24 # -----------------------------------------
25 # Function definition for GENIE reweighting
26 # -----------------------------------------
27 #
28 # One function is defined for each GENIE parameter
29 # *pm1sigma* is specifed: two weights will be evaluated for each parameter,
30 # one with +1 sigma the other with -1 sigma,
31 # it doesn't matter what you write in parameter_sigma and number _of_multisims
32 # if pm1sigma is specifed
33 # Defining one function for each parameter means that the reweighting procedure is done
34 # separately for each parameter
35 # We can define one function with more parameters in parameter_list and in this case the reweight will
36 # be done simultaneously for all the parameters specified
37 #
38 # A list of available parameters is in GenieWeightCalc.cxx
39 #
40 #########################################
41 
42 # MaCCQE to get same result as MiniBooNE
43  genie_qema_miniboone : {
44  type: Genie
45  random_seed: 1
46  parameter_list: ["QEMA"]
47  parameter_sigma: [1.4545] # corresponds to Ma = 1.35 GeV
48  mode: reweight
49  number_of_multisims: 1
50  }
51 
52 # MaCCQE
53  genie_qema: {
54  type: Genie
55  random_seed: 1
56  parameter_list: ["QEMA"]
57  parameter_sigma: [1 ]
58  mode: pm1sigma
59  number_of_multisims: 15
60  }
61 # MaNCEL
62  genie_ncelAxial: {
63  type: Genie
64  random_seed: 1
65  parameter_list: ["NCELaxial"]
66  parameter_sigma: [1 ]
67  mode: pm1sigma
68  number_of_multisims: 15
69  }
70 # EtaNCEL
71  genie_ncelEta: {
72  type: Genie
73  random_seed: 1
74  parameter_list: ["NCELeta"]
75  parameter_sigma: [1 ]
76  mode: pm1sigma
77  number_of_multisims: 15
78  }
79 # CCQE_VecFF
80  genie_qevec: {
81  type: Genie
82  random_seed: 1
83  parameter_list: ["QEVec"]
84  parameter_sigma: [1 ]
85  mode: pm1sigma
86  number_of_multisims: 15
87  }
88 # MaCCRES
89  genie_ccresAxial: {
90  type: Genie
91  random_seed: 1
92  parameter_list: ["CCResAxial"]
93  parameter_sigma: [1 ]
94  mode: pm1sigma
95  number_of_multisims: 15
96  }
97 # MvCCRES
98  genie_ccresVector: {
99  type: Genie
100  random_seed: 1
101  parameter_list: ["CCResVector"]
102  parameter_sigma: [1 ]
103  mode: pm1sigma
104  number_of_multisims: 15
105  }
106 # MaNCRES
107  genie_ncresAxial: {
108  type: Genie
109  random_seed: 1
110  parameter_list: ["NCResAxial"]
111  parameter_sigma: [1 ]
112  mode: pm1sigma
113  number_of_multisims: 15
114  }
115 # MvNCRES
116  genie_ncresVector: {
117  type: Genie
118  random_seed: 1
119  parameter_list: ["NCResVector"]
120  parameter_sigma: [1 ]
121  mode: pm1sigma
122  number_of_multisims: 15
123  }
124 
125 # MaCCRES && MvCCRES && MaNCRES && MvNCRES
126 # genie_resganged: {
127 # type: Genie
128 # random_seed: 1
129 # parameter_list: ["ResGanged"]
130 # parameter_sigma: [1 ]
131 # mode: pm1sigma
132 # number_of_multisims: 15
133 # }
134 
135 # MaCOHPi
136  genie_cohMA: {
137  type: Genie
138  random_seed: 1
139  parameter_list: ["CohMA"]
140  parameter_sigma: [1 ]
141  mode: pm1sigma
142  number_of_multisims: 15
143  }
144 # R0COHPi
145  genie_cohR0: {
146  type: Genie
147  random_seed: 1
148  parameter_list: ["CohR0"]
149  parameter_sigma: [1 ]
150  mode: pm1sigma
151  number_of_multisims: 15
152  }
153 # ReweightNonResRvp1pi
154  genie_NonResRvp1pi: {
155  type: Genie
156  random_seed: 1
157  parameter_list: ["NonResRvp1pi"]
158  parameter_sigma: [1 ]
159  mode: pm1sigma
160  number_of_multisims: 15
161  }
162 # NonResRvbarp1pi
163  genie_NonResRvbarp1pi: {
164  type: Genie
165  random_seed: 1
166  parameter_list: ["NonResRvbarp1pi"]
167  parameter_sigma: [1 ]
168  mode: pm1sigma
169  number_of_multisims: 15
170  }
171 # NonResRvp2pi
172  genie_NonResRvp2pi: {
173  type: Genie
174  random_seed: 1
175  parameter_list: ["NonResRvp2pi"]
176  parameter_sigma: [1 ]
177  mode: pm1sigma
178  number_of_multisims: 15
179  }
180 # NonResRvbarp2pi
181  genie_NonResRvbarp2pi: {
182  type: Genie
183  random_seed: 1
184  parameter_list: ["NonResRvbarp2pi"]
185  parameter_sigma: [1 ]
186  mode: pm1sigma
187  number_of_multisims: 15
188  }
189 # ResDecayGamma
190  genie_ResDecayGamma: {
191  type: Genie
192  random_seed: 1
193  parameter_list: ["ResDecayGamma"]
194  parameter_sigma: [1 ]
195  mode: pm1sigma
196  number_of_multisims: 15
197  }
198 # ResDecayEta
199  genie_ResDecayEta: {
200  type: Genie
201  random_seed: 1
202  parameter_list: ["ResDecayEta"]
203  parameter_sigma: [1 ]
204  mode: pm1sigma
205  number_of_multisims: 15
206  }
207 # ResDecayTheta
208  genie_ResDecayTheta: {
209  type: Genie
210  random_seed: 1
211  parameter_list: ["ResDecayTheta"]
212  parameter_sigma: [1 ]
213  mode: pm1sigma
214  number_of_multisims: 15
215  }
216 # NC
217  genie_NC: {
218  type: Genie
219  random_seed: 1
220  parameter_list: ["NC"]
221  parameter_sigma: [1 ]
222  mode: pm1sigma
223  number_of_multisims: 15
224  }
225 # DISAth
226  genie_DISAth: {
227  type: Genie
228  random_seed: 1
229  parameter_list: ["DISAth"]
230  parameter_sigma: [1 ]
231  mode: pm1sigma
232  number_of_multisims: 15
233  }
234 # DISBth
235  genie_DISBth: {
236  type: Genie
237  random_seed: 1
238  parameter_list: ["DISBth"]
239  parameter_sigma: [1 ]
240  mode: pm1sigma
241  number_of_multisims: 15
242  }
243 # DISCv1u
244  genie_DISCv1u: {
245  type: Genie
246  random_seed: 1
247  parameter_list: ["DISCv1u"]
248  parameter_sigma: [1 ]
249  mode: pm1sigma
250  number_of_multisims: 15
251  }
252 # DISCv2u
253  genie_DISCv2u: {
254  type: Genie
255  random_seed: 1
256  parameter_list: ["DISCv2u"]
257  parameter_sigma: [1 ]
258  mode: pm1sigma
259  number_of_multisims: 15
260  }
261 
262 # DISnucl
263 # genie_DISnucl: {
264 # type: Genie
265 # random_seed: 1
266 # parameter_list: ["DISnucl"]
267 # parameter_sigma: [1 ]
268 # mode: pm1sigma
269 # number_of_multisims: 15
270 # }
271 
272 # AGKYxF
273  genie_AGKYxF: {
274  type: Genie
275  random_seed: 1
276  parameter_list: ["AGKYxF"]
277  parameter_sigma: [1 ]
278  mode: pm1sigma
279  number_of_multisims: 15
280  }
281 # AGKYpT
282  genie_AGKYpT: {
283  type: Genie
284  random_seed: 1
285  parameter_list: ["AGKYpT"]
286  parameter_sigma: [1 ]
287  mode: pm1sigma
288  number_of_multisims: 15
289  }
290 # FormZone
291  genie_FormZone: {
292  type: Genie
293  random_seed: 1
294  parameter_list: ["FormZone"]
295  parameter_sigma: [1 ]
296  mode: pm1sigma
297  number_of_multisims: 15
298  }
299 # FermiGasModelKf
300  genie_FermiGasModelKf: {
301  type: Genie
302  random_seed: 1
303  parameter_list: ["FermiGasModelKf"]
304  parameter_sigma: [1 ]
305  mode: pm1sigma
306  number_of_multisims: 15
307  }
308 # FermiGasModelSf
309  genie_FermiGasModelSf: {
310  type: Genie
311  random_seed: 1
312  parameter_list: ["FermiGasModelSf"]
313  parameter_sigma: [1 ]
314  mode: pm1sigma
315  number_of_multisims: 15
316  }
317 # IntraNukeNmfp
318  genie_IntraNukeNmfp: {
319  type: Genie
320  random_seed: 1
321  parameter_list: ["IntraNukeNmfp"]
322  parameter_sigma: [1 ]
323  mode: pm1sigma
324  number_of_multisims: 15
325  }
326 # IntraNukeNcex
327  genie_IntraNukeNcex: {
328  type: Genie
329  random_seed: 1
330  parameter_list: ["IntraNukeNcex"]
331  parameter_sigma: [1 ]
332  mode: pm1sigma
333  number_of_multisims: 15
334  }
335 # IntraNukeNel
336  genie_IntraNukeNel: {
337  type: Genie
338  random_seed: 1
339  parameter_list: ["IntraNukeNel"]
340  parameter_sigma: [1 ]
341  mode: pm1sigma
342  number_of_multisims: 15
343  }
344 # IntraNukeNinel
345  genie_IntraNukeNinel: {
346  type: Genie
347  random_seed: 1
348  parameter_list: ["IntraNukeNinel"]
349  parameter_sigma: [1 ]
350  mode: pm1sigma
351  number_of_multisims: 15
352  }
353 # IntraNukeNabs
354  genie_IntraNukeNabs: {
355  type: Genie
356  random_seed: 1
357  parameter_list: ["IntraNukeNabs"]
358  parameter_sigma: [1 ]
359  mode: pm1sigma
360  number_of_multisims: 15
361  }
362 # IntraNukeNpi
363  genie_IntraNukeNpi: {
364  type: Genie
365  random_seed: 1
366  parameter_list: ["IntraNukeNpi"]
367  parameter_sigma: [1 ]
368  mode: pm1sigma
369  number_of_multisims: 15
370  }
371 # IntraNukePImfp
372  genie_IntraNukePImfp: {
373  type: Genie
374  random_seed: 1
375  parameter_list: ["IntraNukePImfp"]
376  parameter_sigma: [1 ]
377  mode: pm1sigma
378  number_of_multisims: 15
379  }
380 # IntraNukePIcex
381  genie_IntraNukePIcex: {
382  type: Genie
383  random_seed: 1
384  parameter_list: ["IntraNukePIcex"]
385  parameter_sigma: [1 ]
386  mode: pm1sigma
387  number_of_multisims: 15
388  }
389 # IntraNukePIel
390  genie_IntraNukePIel: {
391  type: Genie
392  random_seed: 1
393  parameter_list: ["IntraNukePIel"]
394  parameter_sigma: [1 ]
395  mode: pm1sigma
396  number_of_multisims: 15
397  }
398 # IntraNukePIinel
399  genie_IntraNukePIinel: {
400  type: Genie
401  random_seed: 1
402  parameter_list: ["IntraNukePIinel"]
403  parameter_sigma: [1 ]
404  mode: pm1sigma
405  number_of_multisims: 15
406  }
407 # IntraNukePIabs
408  genie_IntraNukePIabs: {
409  type: Genie
410  random_seed: 1
411  parameter_list: ["IntraNukePIabs"]
412  parameter_sigma: [1 ]
413  mode: pm1sigma
414  number_of_multisims: 15
415  }
416 # IntraNukePIpi
417  genie_IntraNukePIpi: {
418  type: Genie
419  random_seed: 1
420  parameter_list: ["IntraNukePIpi"]
421  parameter_sigma: [1 ]
422  mode: pm1sigma
423  number_of_multisims: 15
424  }
425 
426 # Here the GENIE parameters are reweighted all at once (40 parameters)
427 genie_all: {
428  type: Genie
429  random_seed: 1
430  parameter_list: ["QEMA", "NCELaxial", "NCELeta", "CCResAxial", "CCResVector", "NCResAxial", "NCResVector", "CohMA", "CohR0", "NonResRvp1pi", "NonResRvbarp1pi", "NonResRvp2pi", "NonResRvbarp2pi", "ResDecayGamma", "ResDecayTheta", "NC", "DISAth", "DISBth", "DISCv1u", "DISCv2u", "AGKYxF", "AGKYpT", "FormZone", "FermiGasModelKf", "IntraNukeNmfp", "IntraNukeNcex", "IntraNukeNel", "IntraNukeNinel", "IntraNukeNabs", "IntraNukeNpi", "IntraNukePImfp", "IntraNukePIcex", "IntraNukePIel", "IntraNukePIinel", "IntraNukePIabs", "IntraNukePIpi"]
431  parameter_sigma: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
432  mode: multisim
433  number_of_multisims: 1000
434  }
435 
436 
437 ####################################################
438 
439 
440 
441 # Old function: here the GENIE parameters are reweighted all at once
442 # genie_all: {
443 # type: Genie
444 # random_seed: 1
445 # parameter_list: ["QEMA", "ResGanged", "CCRes",
446 # "NCRes", "NonResRvp1pi", "NonResRvbarp1pi",
447 # "NonResRvp2pi", "NonResRvbarp2pi", "NC" ]
448 # parameter_sigma: [1 , 1 , 1 ,
449 # 1 , 1 , 1 ,
450 # 1 , 1 , 1 ]
451 # mode: multisim
452 # number_of_multisims: 15
453 # }
454 #
455 }
456 
457 END_PROLOG