Skip to content

Commit 8bf571b

Browse files
sjgardinersjgardiner
committed
Add SBN25_22i_00_000 CMC that is identical to SBN25_20i_00_000 except that
CRPA is used for CCQE instead of the Valencia model.
1 parent 2cb68bd commit 8bf571b

File tree

6 files changed

+639
-0
lines changed

6 files changed

+639
-0
lines changed

‎config/SBN25_22i/CommonParam.xml‎

Lines changed: 337 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,337 @@
1+
<?xml version="1.0" encoding="ISO-8859-1"?>
2+
3+
<!--
4+
***************************************************************************************************
5+
This file defines parameters sets for parameters which are common between algorithms and therefore
6+
has to be the same for those parameters to be meaningful.
7+
This is file is read from the AlgConfigPool and the parameters set are expected to be used by the
8+
algorithms to configure these common parameters.
9+
10+
Author:
11+
Costas Andreopoulos <c.andreopoulos \at cern.ch>
12+
University of Liverpool
13+
14+
and
15+
Marco Roda <marco.roda \at liverpool.ac.uk>
16+
University of Liverpool
17+
***************************************************************************************************
18+
-->
19+
20+
21+
22+
<common_Param_list>
23+
24+
25+
<param_set name="Tunable">
26+
<!-- Reserved register for tuning -->
27+
28+
<param type="double" name="RES-Ma"> 1.088962 </param>
29+
<param type="double" name="DIS-XSecScale"> 1.062213 </param>
30+
<param type="double" name="RES-CC-XSecScale"> 0.838730 </param>
31+
32+
</param_set>
33+
34+
35+
<param_set name="WeakInt">
36+
<param type="double" name="WeinbergAngle"> 0.501716712132 </param> <!-- 0.501568 -->
37+
<param type="double" name="FermiConstant"> 1.16639E-5 </param> <!-- GeV^-2 --> <!-- historical GENIE value -->
38+
<param type="double" name="FineStructureConstant"> 7.2973525332858855E-3 </param> <!-- historical GENIE value -->
39+
40+
</param_set>
41+
42+
<param_set name="StrongInt">
43+
<param type="double" name="SU3-D"> 0.804 </param>
44+
<param type="double" name="SU3-F"> 0.463 </param>
45+
46+
</param_set>
47+
48+
49+
<param_set name="CKM">
50+
<!-- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
51+
CKM quark mixing parameters
52+
Chinese Physics C Vol.40, No. 10 (2016) 100001
53+
Review of Particle Physics
54+
-->
55+
<param type="double" name="CKM-Vud"> 0.97417 </param>
56+
<param type="double" name="CKM-Vus"> 0.2248 </param>
57+
<param type="double" name="CKM-Vcd"> 0.220 </param> <!-- 0.2239 -->
58+
<param type="double" name="CKM-Vcs"> 0.995 </param>
59+
60+
<param type="double" name="CabibboAngle"> 0.227780466682 </param> <!-- 0.226893 -->
61+
62+
</param_set>
63+
64+
<param_set name="NUCL">
65+
66+
<!-- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
67+
NUCL-Ro (in fm) is a scale parameter driving the effective nuclear sizes (Ro in R=Ro*A^1/3)
68+
-->
69+
<param type="double" name="NUCL-R0"> 1.4 </param>
70+
<param type="double" name="NUCL-NR"> 3.0 </param>
71+
72+
</param_set>
73+
74+
<param_set name="Masses">
75+
<param type="double" name="Charm-Mass"> 1.430 </param> <!-- GeV -->
76+
</param_set>
77+
78+
79+
<param_set name="NonResBackground">
80+
<!--
81+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
82+
RES/DIS joining scheme
83+
At this point there is a single scheme, the one originally developed in neugen3.
84+
For W > Wcut : RES -> 0, + DIS -> full
85+
For W <= Wcut : RES -> full + `DIS' (non-RES bkg) -> modified by DIS-HMultWgt-* params
86+
-
87+
- ...
88+
89+
-->
90+
<param type="bool" name="UseDRJoinScheme"> true </param>
91+
<param type="double" name="Wcut"> 1.809000 </param>
92+
93+
<!--
94+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
95+
NEUGEN parameters applied to DIS hadronic multiplicity distributions
96+
to avoid 2-ble counting with RES
97+
-->
98+
99+
<param type="double" name="DIS-HMultWgt-vp-CC-m2"> 0.008 </param>
100+
<param type="double" name="DIS-HMultWgt-vp-CC-m3"> 0.943153 </param>
101+
<param type="double" name="DIS-HMultWgt-vp-NC-m2"> 0.100 </param>
102+
<param type="double" name="DIS-HMultWgt-vp-NC-m3"> 1.000 </param>
103+
<param type="double" name="DIS-HMultWgt-vp-EM-m2"> 0.100 </param>
104+
<param type="double" name="DIS-HMultWgt-vp-EM-m3"> 1.000 </param>
105+
<param type="double" name="DIS-HMultWgt-vn-CC-m2"> 0.030802 </param>
106+
<param type="double" name="DIS-HMultWgt-vn-CC-m3"> 2.338338 </param>
107+
<param type="double" name="DIS-HMultWgt-vn-NC-m2"> 0.300 </param>
108+
<param type="double" name="DIS-HMultWgt-vn-NC-m3"> 1.000 </param>
109+
<param type="double" name="DIS-HMultWgt-vn-EM-m2"> 0.300 </param>
110+
<param type="double" name="DIS-HMultWgt-vn-EM-m3"> 1.000 </param>
111+
<param type="double" name="DIS-HMultWgt-vbp-CC-m2"> 0.030802 </param>
112+
<param type="double" name="DIS-HMultWgt-vbp-CC-m3"> 2.338338 </param>
113+
<param type="double" name="DIS-HMultWgt-vbp-NC-m2"> 0.300 </param>
114+
<param type="double" name="DIS-HMultWgt-vbp-NC-m3"> 1.000 </param>
115+
<param type="double" name="DIS-HMultWgt-vbp-EM-m2"> 0.300 </param>
116+
<param type="double" name="DIS-HMultWgt-vbp-EM-m3"> 1.000 </param>
117+
<param type="double" name="DIS-HMultWgt-vbn-CC-m2"> 0.008 </param>
118+
<param type="double" name="DIS-HMultWgt-vbn-CC-m3"> 0.943153 </param>
119+
<param type="double" name="DIS-HMultWgt-vbn-NC-m2"> 0.100 </param>
120+
<param type="double" name="DIS-HMultWgt-vbn-NC-m3"> 1.000 </param>
121+
<param type="double" name="DIS-HMultWgt-vbn-EM-m2"> 0.100 </param>
122+
<param type="double" name="DIS-HMultWgt-vbn-EM-m3"> 1.000 </param>
123+
</param_set>
124+
125+
126+
<param_set name="FermiGas">
127+
<!--
128+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
129+
Parameters related with GENIE's Fermi Gas (RFG) model implementation
130+
- RFG-Momentum-CutOff is a momentum cut-off for the NN correlation tail
131+
- RFG-NucRemovalE@Pdg=10LZZZAAAI is the removal energy for the nucleus with the specified pdg
132+
code. If none is used then the average binding energy will be computed from Wapstra's semi-
133+
empirical formula.
134+
Currently, if you explicitly specify a binding energy for a nucleus then the same value will
135+
be used for all isotopes.
136+
-->
137+
<param type="double" name="RFG-MomentumCutOff"> 0.0001 </param> <!-- No BR tail -->
138+
<param type="double" name="RFG-NucRemovalE@Pdg=1000020040"> 0.0150 </param> <!-- He4 -->
139+
<param type="double" name="RFG-NucRemovalE@Pdg=1000030060"> 0.0150 </param> <!-- Li6 -->
140+
<param type="double" name="RFG-NucRemovalE@Pdg=1000060120"> 0.0200 </param> <!-- C12 -->
141+
<param type="double" name="RFG-NucRemovalE@Pdg=1000080160"> 0.0160 </param> <!-- O16 -->
142+
<param type="double" name="RFG-NucRemovalE@Pdg=1000120240"> 0.0250 </param> <!-- Mg24 -->
143+
<param type="double" name="RFG-NucRemovalE@Pdg=1000180400"> 0.0280 </param> <!-- Ar40 -->
144+
<param type="double" name="RFG-NucRemovalE@Pdg=1000200400"> 0.0280 </param> <!-- Ca40 -->
145+
<param type="double" name="RFG-NucRemovalE@Pdg=1000260560"> 0.0230 </param> <!-- Fe56 -->
146+
<param type="double" name="RFG-NucRemovalE@Pdg=1000280580"> 0.0300 </param> <!-- Ni58 -->
147+
<param type="double" name="RFG-NucRemovalE@Pdg=1000501190"> 0.0280 </param> <!-- Sn119 -->
148+
<param type="double" name="RFG-NucRemovalE@Pdg=1000791970"> 0.0310 </param> <!-- Au197 -->
149+
<param type="double" name="RFG-NucRemovalE@Pdg=1000822080"> 0.0310 </param> <!-- Pb208 -->
150+
<param type="bool" name="RFG-UseParametrization"> false </param>
151+
152+
<!--
153+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
154+
Table of Fermi momentum (kF) constants for various nuclei
155+
The tables can be found in $GENIE/config/FermiMomentumTables.xml
156+
-->
157+
<param type="string" name="FermiMomentumTable"> SuSA </param>
158+
159+
</param_set>
160+
161+
162+
<param_set name="Coherent">
163+
164+
<param type="double" name="COH-Ro"> 1.000 </param>
165+
166+
<!--
167+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
168+
Minimum and maximum considered Q^2 for Berger-Sehgal coherent reactions when estimating
169+
the max cross section.
170+
Units in GeV^2.
171+
-->
172+
<param type="double" name="COH-Q2-min"> 0.000 </param>
173+
<param type="double" name="COH-Q2-max"> 1.000 </param>
174+
175+
<!--
176+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
177+
Maximum considered t for Berger-Sehgal coherent reactions when estimating the
178+
max cross section.
179+
Units in GeV^2.
180+
-->
181+
<param type="double" name="COH-t-max"> 0.250 </param>
182+
183+
</param_set>
184+
185+
186+
<param_set name="Diffractive">
187+
188+
<!--
189+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
190+
Slope parameter beta for diffractive scattering (GeV^-2) (b ~ 0.333 * nucleon_size^2)
191+
-->
192+
<param type="double" name="DFR-Beta"> 7.0 </param>
193+
194+
<!--
195+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
196+
Maximum considered t for diffractive scattering reactions when estimating the max cross section.
197+
Units in GeV^2.
198+
199+
This value is read in KpHase space!!! It's mental, but we have to be careful before removing this
200+
Or changing the name of this parameter set
201+
-->
202+
<param type="double" name="DFR-t-max"> 0.35 </param>
203+
204+
</param_set>
205+
206+
<param_set name="QuasiElastic">
207+
<!--
208+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
209+
Value of axial form factor at Q2=0
210+
-->
211+
<param type="double" name="QEL-FA0"> -1.2670 </param>
212+
213+
<!--
214+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
215+
Axial and vector masses for quasi-elastic scattering
216+
-->
217+
<param type="double" name="QEL-Ma"> 0.994989 </param>
218+
<param type="double" name="QEL-Mv"> 0.840 </param>
219+
220+
</param_set>
221+
222+
<param_set name="MagnMoments">
223+
<!--
224+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
225+
Proton and neutron anomalous magnetic moments
226+
-->
227+
<param type="double" name="AnomMagnMoment-P"> 2.7930 </param> <!-- 2.7928473 -->
228+
<param type="double" name="AnomMagnMoment-N"> -1.913042 </param>
229+
230+
</param_set>
231+
232+
<param_set name="MultiNucleons">
233+
<!-- Q3 max for 2p2h model -->
234+
<param type="double" name="NSV-Q3Max"> 2.0 </param>
235+
</param_set>
236+
237+
<param_set name="ElasticFF">
238+
<!--
239+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
240+
Elastic form factors used for QEL CC cross section calculation.
241+
Options are:
242+
- genie::DipoleELFormFactorsModel
243+
- genie::BBA03ELFormFactorsModel, H.Budd, NuINT-02 proceedings
244+
- genie::BBA05ELFormFactorsModel, R.Bradford, A.Bodek, H.Budd and J.Arrington, hep-ex/0602017
245+
- genie::BBA07ELFormFactorsModel, R.Bradford, A.Bodek, H.Budd and S.Avvakumov, Euro.Phys.J.C53 (2008);[arXiv:0708.1946 [hep-ex]]
246+
-->
247+
<param type="alg" name="ElasticFormFactorsModel"> genie::BBA07ELFormFactorsModel/Default </param>
248+
249+
<!-- Option for turning on Transverse Enhancement by Elastic Form Factor adjustment.
250+
See http://arxiv.org/abs/1405.0583 and http://arxiv.org/abs/1106.0340
251+
-->
252+
<param type="bool" name="UseElFFTransverseEnhancement"> false </param>
253+
<param type="alg" name="TransverseEnhancement"> genie::TransverseEnhancementFFModel/Default </param>
254+
255+
</param_set>
256+
257+
<param_set name="Resonances">
258+
<!--
259+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
260+
Default list of baryon resonances included in cross section models and generation threads
261+
-->
262+
<param type="string" name="ResonanceNameList">
263+
P33(1232),S11(1535),D13(1520),S11(1650),D13(1700),D15(1675),
264+
S31(1620),D33(1700),P11(1440),P33(1600),P13(1720),F15(1680),
265+
P31(1910),P33(1920),F35(1905),F37(1950),P11(1710)
266+
</param>
267+
268+
</param_set>
269+
270+
<param_set name="Lepton">
271+
<param type="bool" name="ApplyCoulombCorrection"> false </param>
272+
273+
</param_set>
274+
275+
276+
<param_set name="BoostedDarkMatter">
277+
<!--
278+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
279+
Set the coupling of dark matter to nucleons
280+
-->
281+
<param type="double" name="ZpCoupling"> 0.1 </param>
282+
<param type="double" name="DarkLeftCharge"> 1.0 </param>
283+
<param type="double" name="DarkRightCharge"> 1.0 </param>
284+
<param type="double" name="DarkScalarCharge"> 1.0 </param>
285+
<param type="double" name="UpLeftCharge"> 1.0 </param>
286+
<param type="double" name="UpRightCharge"> 1.0 </param>
287+
<param type="double" name="DownLeftCharge"> 1.0 </param>
288+
<param type="double" name="DownRightCharge"> 1.0 </param>
289+
<param type="double" name="StrangeLeftCharge"> 1.0 </param>
290+
<param type="double" name="StrangeRightCharge"> 1.0 </param>
291+
<param type="double" name="CharmLeftCharge"> 1.0 </param>
292+
<param type="double" name="CharmRightCharge"> 1.0 </param>
293+
<param type="double" name="ElectronLeftCharge"> 1.0 </param>
294+
<param type="double" name="ElectronRightCharge"> 1.0 </param>
295+
<param type="double" name="DMEL-Mp"> 1.441 </param>
296+
<param type="double" name="DMEL-Mpi"> 0.1349766 </param>
297+
<param type="double" name="DMEL-Meta"> 0.547862 </param>
298+
<param type="double" name="AxialVectorSpin-u"> 0.827 </param>
299+
<param type="double" name="AxialVectorSpin-d"> -0.38 </param>
300+
<param type="double" name="AxialVectorSpin-s"> -0.0427 </param>
301+
</param_set>
302+
303+
304+
<param_set name="KNO2Pythia">
305+
<!--
306+
307+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
308+
Invariant mass window for the transition from KNO model to PYTHIA
309+
-->
310+
311+
<param type="double" name="KNO2PYTHIA-Wmin"> 2.30 </param>
312+
<param type="double" name="KNO2PYTHIA-Wmax"> 3.00 </param>
313+
</param_set>
314+
315+
<param_set name="PDG">
316+
<!-- when using the PDG table provideed by genie,
317+
that are contained in $GENIE/data/evgen/catalogues/pdg
318+
this string specifies the file to be used.
319+
Please be aware that changing the PDG table might invalidated the selected tune -->
320+
<param type="string" name="PDG-TableName"> genie_pdg_table.txt </param>
321+
</param_set>
322+
323+
324+
325+
<param_set name="Validation">
326+
<!--
327+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
328+
GENIE validity range
329+
-->
330+
<param type="double" name="GVLD-Emin"> 0.010 </param>
331+
<param type="double" name="GVLD-Emax"> 1000.000 </param>
332+
333+
</param_set>
334+
335+
336+
337+
</common_Param_list>

‎config/SBN25_22i/MECInteractionListGenerator.xml‎

Lines changed: 47 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,47 @@
1+
<?xml version="1.0" encoding="ISO-8859-1"?>
2+
3+
<alg_conf>
4+
5+
<!--
6+
Configuration for the MECInteractionListGenerator InteractionListGeneratorI
7+
8+
Configurable Parameters:
9+
.......................................................................................................
10+
Name Type Optional Comment Default
11+
.......................................................................................................
12+
-->
13+
14+
<param_set name="Default">
15+
16+
<!--
17+
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
18+
SetDiNucleonCode: Include the dinucleon code in the interaction list generator.
19+
Want `true` for EmpiricalMEC, `false` for NSV
20+
-->
21+
<param type="bool" name="SetDiNucleonCode"> false </param>
22+
<param type="bool" name="SetDiNucleonCodeEM"> true </param>
23+
24+
</param_set>
25+
26+
27+
28+
<param_set name="CC-Default">
29+
<param type="bool" name="is-CC"> true </param>
30+
<param type="bool" name="is-NC"> false </param>
31+
<param type="bool" name="is-EM"> false </param>
32+
</param_set>
33+
34+
<param_set name="NC-Default">
35+
<param type="bool" name="is-CC"> false </param>
36+
<param type="bool" name="is-NC"> true </param>
37+
<param type="bool" name="is-EM"> false </param>
38+
</param_set>
39+
40+
<param_set name="EM-Default">
41+
<param type="bool" name="is-CC"> false </param>
42+
<param type="bool" name="is-NC"> false </param>
43+
<param type="bool" name="is-EM"> true </param>
44+
</param_set>
45+
46+
</alg_conf>
47+

0 commit comments

Comments
 (0)