lars-p GENIE(Reweight) development #478
There are no files selected for viewing
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
|
|
@@ -14,3 +14,4 @@ src/Framework/Conventions/GVersion.h | |
| *.rootmap | ||
| *_ROOT_DICT_* | ||
| *.swp | ||
| .DS_Store | ||
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -1 +1 @@ | ||
| 3.06.02 |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,312 @@ | ||
| <?xml version="1.0" encoding="ISO-8859-1"?> | ||
|
|
||
| <!-- | ||
| *************************************************************************************************** | ||
| This file defines parameters sets for parameters which are common between algorithms and therefore | ||
| has to be the same for those parameters to be meaningful. | ||
| This is file is read from the AlgConfigPool and the parameters set are expected to be used by the | ||
| algorithms to configure these common parameters. | ||
|
|
||
| Author: | ||
| Costas Andreopoulos <c.andreopoulos \at cern.ch> | ||
| University of Liverpool | ||
|
|
||
| and | ||
| Marco Roda <marco.roda \at liverpool.ac.uk> | ||
| University of Liverpool | ||
| *************************************************************************************************** | ||
| --> | ||
|
|
||
|
|
||
|
|
||
| <common_Param_list> | ||
|
|
||
|
|
||
| <param_set name="Tunable"> | ||
| <!-- Reserved register for tuning --> | ||
| <param type="double" name="RES-Ma"> 1.088962 </param> | ||
| <param type="double" name="DIS-CC-XSecScale"> 1.062213 </param> | ||
| <param type="double" name="RES-CC-XSecScale"> 0.838730 </param> | ||
| </param_set> | ||
|
|
||
| <param_set name="KNO2Pythia"> | ||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Invariant mass window for the transition from KNO model to PYTHIA | ||
| --> | ||
|
|
||
| <param type="double" name="KNO2PYTHIA-Wmin"> 2.30 </param> | ||
| <param type="double" name="KNO2PYTHIA-Wmax"> 3.00 </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="WeakInt"> | ||
| <param type="double" name="WeinbergAngle"> 0.501716712132 </param> <!-- 0.501568 --> | ||
| <param type="double" name="FermiConstant"> 1.16639E-5 </param> <!-- GeV^-2 --> <!-- historical GENIE value --> | ||
| <param type="double" name="FineStructureConstant"> 7.2973525332858855E-3 </param> <!-- historical GENIE value --> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="StrongInt"> | ||
| <param type="double" name="SU3-D"> 0.804 </param> | ||
| <param type="double" name="SU3-F"> 0.463 </param> | ||
|
|
||
| </param_set> | ||
|
|
||
|
|
||
| <param_set name="CKM"> | ||
| <!-- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| CKM quark mixing parameters | ||
| Chinese Physics C Vol.40, No. 10 (2016) 100001 | ||
| Review of Particle Physics | ||
| --> | ||
| <param type="double" name="CKM-Vud"> 0.97417 </param> | ||
| <param type="double" name="CKM-Vus"> 0.2248 </param> | ||
| <param type="double" name="CKM-Vcd"> 0.220 </param> <!-- 0.2239 --> | ||
| <param type="double" name="CKM-Vcs"> 0.995 </param> | ||
|
|
||
| <param type="double" name="CabibboAngle"> 0.227780466682 </param> <!-- 0.226893 --> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="NUCL"> | ||
|
|
||
| <!-- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| NUCL-Ro (in fm) is a scale parameter driving the effective nuclear sizes (Ro in R=Ro*A^1/3) | ||
| --> | ||
| <param type="double" name="NUCL-R0"> 1.4 </param> | ||
| <param type="double" name="NUCL-NR"> 3.0 </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="Masses"> | ||
| <param type="double" name="Charm-Mass"> 1.430 </param> <!-- GeV --> | ||
| </param_set> | ||
|
|
||
|
|
||
| <param_set name="NonResBackground"> | ||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| RES/DIS joining scheme | ||
| At this point there is a single scheme, the one originally developed in neugen3. | ||
| For W > Wcut : RES -> 0, + DIS -> full | ||
| For W <= Wcut : RES -> full + `DIS' (non-RES bkg) -> modified by DIS-HMultWgt-* params | ||
| - | ||
| - ... | ||
|
|
||
| --> | ||
| <param type="bool" name="UseDRJoinScheme"> true </param> | ||
| <param type="double" name="Wcut"> 1.809000 </param> | ||
|
|
||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| NEUGEN parameters applied to DIS hadronic multiplicity distributions | ||
| to avoid 2-ble counting with RES | ||
| --> | ||
| <param type="double" name="DIS-HMultWgt-vp-CC-m2"> 0.008 </param> | ||
| <param type="double" name="DIS-HMultWgt-vp-CC-m3"> 0.943153 </param> | ||
| <param type="double" name="DIS-HMultWgt-vp-NC-m2"> 0.100 </param> | ||
| <param type="double" name="DIS-HMultWgt-vp-NC-m3"> 1.000 </param> | ||
| <param type="double" name="DIS-HMultWgt-vp-EM-m2"> 0.100 </param> | ||
| <param type="double" name="DIS-HMultWgt-vp-EM-m3"> 1.000 </param> | ||
| <param type="double" name="DIS-HMultWgt-vn-CC-m2"> 0.030802 </param> | ||
| <param type="double" name="DIS-HMultWgt-vn-CC-m3"> 2.338338 </param> | ||
| <param type="double" name="DIS-HMultWgt-vn-NC-m2"> 0.300 </param> | ||
| <param type="double" name="DIS-HMultWgt-vn-NC-m3"> 1.000 </param> | ||
| <param type="double" name="DIS-HMultWgt-vn-EM-m2"> 0.300 </param> | ||
| <param type="double" name="DIS-HMultWgt-vn-EM-m3"> 1.000 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbp-CC-m2"> 0.030802 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbp-CC-m3"> 2.338338 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbp-NC-m2"> 0.300 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbp-NC-m3"> 1.000 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbp-EM-m2"> 0.300 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbp-EM-m3"> 1.000 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbn-CC-m2"> 0.008 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbn-CC-m3"> 0.943153 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbn-NC-m2"> 0.100 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbn-NC-m3"> 1.000 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbn-EM-m2"> 0.100 </param> | ||
| <param type="double" name="DIS-HMultWgt-vbn-EM-m3"> 1.000 </param> | ||
|
|
||
| </param_set> | ||
|
|
||
|
|
||
| <param_set name="FermiGas"> | ||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Parameters related with GENIE's Fermi Gas (RFG) model implementation | ||
| - RFG-Momentum-CutOff is a momentum cut-off for the NN correlation tail | ||
| - RFG-NucRemovalE@Pdg=10LZZZAAAI is the removal energy for the nucleus with the specified pdg | ||
| code. If none is used then the average binding energy will be computed from Wapstra's semi- | ||
| empirical formula. | ||
| Currently, if you explicitly specify a binding energy for a nucleus then the same value will | ||
| be used for all isotopes. | ||
| --> | ||
| <param type="double" name="RFG-MomentumCutOff"> 0.5000 </param> <!-- BR tail --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000020040"> 0.0150 </param> <!-- He4 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000030060"> 0.0170 </param> <!-- Li6 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000060120"> 0.0250 </param> <!-- C12 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000080160"> 0.0270 </param> <!-- O16 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000120240"> 0.0320 </param> <!-- Mg24 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000180400"> 0.0295 </param> <!-- Ar40 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000200400"> 0.0280 </param> <!-- Ca40 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000220480"> 0.0300 </param> <!-- Ti48 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000260560"> 0.0360 </param> <!-- Fe56 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000280580"> 0.0360 </param> <!-- Ni58 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000501190"> 0.0280 </param> <!-- Sn119 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000791970"> 0.0310 </param> <!-- Au197 --> | ||
| <param type="double" name="RFG-NucRemovalE@Pdg=1000822080"> 0.0440 </param> <!-- Pb208 --> | ||
| <param type="bool" name="RFG-UseParametrization"> false </param> | ||
|
|
||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Table of Fermi momentum (kF) constants for various nuclei | ||
| The tables can be found in $GENIE/config/FermiMomentumTables.xml | ||
| --> | ||
| <param type="string" name="FermiMomentumTable"> Default </param> | ||
|
|
||
| </param_set> | ||
|
|
||
|
|
||
| <param_set name="Coherent"> | ||
|
|
||
| <param type="double" name="COH-Ro"> 1.000 </param> | ||
|
|
||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Minimum and maximum considered Q^2 for Berger-Sehgal coherent reactions when estimating | ||
| the max cross section. | ||
| Units in GeV^2. | ||
| --> | ||
| <param type="double" name="COH-Q2-min"> 0.000 </param> | ||
| <param type="double" name="COH-Q2-max"> 1.000 </param> | ||
|
|
||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Maximum considered t for Berger-Sehgal coherent reactions when estimating the | ||
| max cross section. | ||
| Units in GeV^2. | ||
| --> | ||
| <param type="double" name="COH-t-max"> 0.250 </param> | ||
|
|
||
| </param_set> | ||
|
|
||
|
|
||
| <param_set name="Diffractive"> | ||
|
|
||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Slope parameter beta for diffractive scattering (GeV^-2) (b ~ 0.333 * nucleon_size^2) | ||
| --> | ||
| <param type="double" name="DFR-Beta"> 7.0 </param> | ||
|
|
||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Maximum considered t for diffractive scattering reactions when estimating the max cross section. | ||
| Units in GeV^2. | ||
|
|
||
| This value is read in KpHase space!!! It's mental, but we have to be careful before removing this | ||
| Or changing the name of this parameter set | ||
| --> | ||
| <param type="double" name="DFR-t-max"> 0.35 </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="QuasiElastic"> | ||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Value of axial form factor at Q2=0 | ||
| --> | ||
| <param type="double" name="QEL-FA0"> -1.2670 </param> | ||
|
|
||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Axial and vector masses for quasi-elastic scattering | ||
| --> | ||
| <param type="double" name="QEL-Ma"> 0.994989 </param> | ||
| <param type="double" name="QEL-Mv"> 0.840 </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="MagnMoments"> | ||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Proton and neutron anomalous magnetic moments | ||
| --> | ||
| <param type="double" name="AnomMagnMoment-P"> 2.7930 </param> <!-- 2.7928473 --> | ||
| <param type="double" name="AnomMagnMoment-N"> -1.913042 </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="MultiNucleons"> | ||
| <!-- Q3 max for 2p2h model --> | ||
| <!-- <param type="double" name="NSV-Q3Max"> 1.2 </param> --> \\ Nieves | ||
| <param type="double" name="NSV-Q3Max"> 2.0 </param> \\ SuSAv2, Martini | ||
| </param_set> | ||
|
|
||
| <param_set name="ElasticFF"> | ||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Elastic form factors used for QEL CC cross section calculation. | ||
| Options are: | ||
| - genie::DipoleELFormFactorsModel | ||
| - genie::BBA03ELFormFactorsModel, H.Budd, NuINT-02 proceedings | ||
| - genie::BBA05ELFormFactorsModel, R.Bradford, A.Bodek, H.Budd and J.Arrington, hep-ex/0602017 | ||
| - genie::BBA07ELFormFactorsModel, R.Bradford, A.Bodek, H.Budd and S.Avvakumov, Euro.Phys.J.C53 (2008);[arXiv:0708.1946 [hep-ex]] | ||
| --> | ||
| <param type="alg" name="ElasticFormFactorsModel"> genie::BBA07ELFormFactorsModel/Default </param> | ||
|
|
||
| <!-- Option for turning on Transverse Enhancement by Elastic Form Factor adjustment. | ||
| See http://arxiv.org/abs/1405.0583 and http://arxiv.org/abs/1106.0340 | ||
| --> | ||
| <param type="bool" name="UseElFFTransverseEnhancement"> false </param> | ||
| <param type="alg" name="TransverseEnhancement"> genie::TransverseEnhancementFFModel/Default </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="Resonances"> | ||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| Default list of baryon resonances included in cross section models and generation threads | ||
| --> | ||
| <param type="string" name="ResonanceNameList"> | ||
| P33(1232),S11(1535),D13(1520),S11(1650),D13(1700),D15(1675), | ||
| S31(1620),D33(1700),P11(1440),P33(1600),P13(1720),F15(1680), | ||
| P31(1910),P33(1920),F35(1905),F37(1950),P11(1710) | ||
| </param> | ||
|
|
||
|
|
||
| <!-- | ||
| In this configuration set we include only the first two resonances | ||
| <param type="string" name="ResonanceNameList"> P33(1232),P11(1440) </param> | ||
| --> | ||
|
|
||
|
|
||
| </param_set> | ||
|
|
||
| <param_set name="Lepton"> | ||
| <param type="bool" name="ApplyCoulombCorrection"> false </param> | ||
|
|
||
| </param_set> | ||
|
|
||
|
|
||
| <param_set name="PDG"> | ||
| <!-- when using the PDG table provideed by genie, | ||
| that are contained in $GENIE/data/evgen/catalogues/pdg | ||
| this string specifies the file to be used. | ||
| Please be aware that changing the PDG table might invalidated the selected tune --> | ||
| <param type="string" name="PDG-TableName"> genie_pdg_table.txt </param> | ||
| </param_set> | ||
|
|
||
|
|
||
| <param_set name="Validation"> | ||
| <!-- | ||
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| GENIE validity range | ||
| --> | ||
| <param type="double" name="GVLD-Emin"> 0.010 </param> | ||
| <param type="double" name="GVLD-Emax"> 1000.000 </param> | ||
|
|
||
| </param_set> | ||
|
|
||
| </common_Param_list> |
Oops, something went wrong.
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
An erroneous removal of a space here? In GitHub the comments no longer line up in their indentation.