Project

General

Profile

R-2_11_2 variants

As in the table on the GENIE_Configuration_Files page, there are a variety of configurations for R-2_11_2.

Base R-2_11_2 vs. DefaultPlusMECWithNCHadWfix

--- ./genie_xsec/v2_11_2/NULL/DefaultPlusMECWithNCHadWFix/data/UserPhysicsOptions.original.tidyReduced.xml    2016-10-17 16:12:06.112907095 -0500
+++ ./genie_xsec/v2_11_2/NULL/DefaultPlusMECWithNCHadWFix/data/UserPhysicsOptions.tidyReduced.xml    2016-10-17 16:12:06.177726583 -0500
@@ -186,8 +186,8 @@
     <param type="double" name="PYTHIA-GaussianPt2">0.44</param>
     <param type="double" name="PYTHIA-NonGaussianPt2Tail">0.01</param>
     <param type="double" name="PYTHIA-RemainingEnergyCutoff">0.20</param>
-    <param type="double" name="KNO2PYTHIA-Wmin">2000.30</param>
-    <param type="double" name="KNO2PYTHIA-Wmax">3000.00</param>
+    <param type="double" name="KNO2PYTHIA-Wmin">2.30</param>
+    <param type="double" name="KNO2PYTHIA-Wmax">3.00</param>
     <param type="double" name="DIS-HMultWgt-vp-CC-m2">0.100</param>
     <param type="double" name="DIS-HMultWgt-vp-CC-m3">1.000</param>
     <param type="double" name="DIS-HMultWgt-vp-NC-m2">0.100</param>
@@ -231,18 +231,18 @@
     <param type="bool" name="INUKE-DoCompoundNucleus">true</param>
     <param type="bool" name="INUKE-XsecNNCorr">false</param>
     <param type="bool" name="DecayParticleWithCode=111">false</param>
-    <param type="bool" name="DecayParticleWithCode=411">false</param>
-    <param type="bool" name="DecayParticleWithCode=-411">false</param>
-    <param type="bool" name="DecayParticleWithCode=421">false</param>
-    <param type="bool" name="DecayParticleWithCode=-421">false</param>
-    <param type="bool" name="DecayParticleWithCode=431">false</param>
-    <param type="bool" name="DecayParticleWithCode=-431">false</param>
-    <param type="bool" name="DecayParticleWithCode=4122">false</param>
-    <param type="bool" name="DecayParticleWithCode=4112">false</param>
-    <param type="bool" name="DecayParticleWithCode=4212">false</param>
-    <param type="bool" name="DecayParticleWithCode=4222">false</param>
-    <param type="bool" name="DecayParticleWithCode=15">false</param>
-    <param type="bool" name="DecayParticleWithCode=-15">false</param>
+    <param type="bool" name="DecayParticleWithCode=411">true</param>
+    <param type="bool" name="DecayParticleWithCode=-411">true</param>
+    <param type="bool" name="DecayParticleWithCode=421">true</param>
+    <param type="bool" name="DecayParticleWithCode=-421">true</param>
+    <param type="bool" name="DecayParticleWithCode=431">true</param>
+    <param type="bool" name="DecayParticleWithCode=-431">true</param>
+    <param type="bool" name="DecayParticleWithCode=4122">true</param>
+    <param type="bool" name="DecayParticleWithCode=4112">true</param>
+    <param type="bool" name="DecayParticleWithCode=4212">true</param>
+    <param type="bool" name="DecayParticleWithCode=4222">true</param>
+    <param type="bool" name="DecayParticleWithCode=15">true</param>
+    <param type="bool" name="DecayParticleWithCode=-15">true</param>
     <param type="bool" name="DecayParticleWithCode=221">true</param>
     <param type="bool" name="DecayParticleWithCode=331">true</param>
     <param type="bool" name="DecayParticleWithCode=113">true</param>

The changes to KNO2PYTHIA-Wmin and KNO2PYTHIA-Wmax were necessitated by a mistake that ended up in the official R-2_11_2 and R-2_11_0 tags. This affect how the hadronic recoil is hadronized by incorrectly pushing the transition between using the KNO prescription and the Pythia6 scheme to values of W much higher than normally found. The transition region of 2.3 to 3.0 GeV is the traditional (and "official" intended) range.

This second block of changes is in order to force GENIE (via Pythia6) to do the decays of charmed hadrons and taus. We do this because Geant4 does not properly generate secondaries for charm decay (the energy goes "poof"). The same used to be true for taus as well so historically we've also forced this; this is no longer absolutely necessary but we retained it in order to maintain consistency.

In order to highlight the important differences the remaining differences are relative to DefaultPlusMECWithNC (since there are no generated splines for the base configuration).

The DefaultPlusMECWithNC set also has a modified EventGeneratorListAssembler.xml file:

$ diff -u GXSPLINES-v2_11_2-DefaultPlusMECWithNC/cfg/EventGeneratorListAssembler.*tidyReduced.xml
--- GXSPLINES-v2_11_2-DefaultPlusMECWithNC/cfg/EventGeneratorListAssembler.original.tidyReduced.xml    2016-10-07 15:22:41.766621260 -0500
+++ GXSPLINES-v2_11_2-DefaultPlusMECWithNC/cfg/EventGeneratorListAssembler.tidyReduced.xml    2016-10-07 15:22:41.758540017 -0500
@@ -177,6 +177,24 @@
     <param type="alg" name="Generator-6">genie::EventGenerator/MEC-CC</param>
   </param_set>
   <param_set name="Default+MEC">
+    <param type="int" name="NGenerators">15</param>
+    <param type="alg" name="Generator-0">genie::EventGenerator/QEL-CC</param>
+    <param type="alg" name="Generator-1">genie::EventGenerator/QEL-NC</param>
+    <param type="alg" name="Generator-2">genie::EventGenerator/RES-CC</param>
+    <param type="alg" name="Generator-3">genie::EventGenerator/RES-NC</param>
+    <param type="alg" name="Generator-4">genie::EventGenerator/DIS-CC</param>
+    <param type="alg" name="Generator-5">genie::EventGenerator/DIS-NC</param>
+    <param type="alg" name="Generator-6">genie::EventGenerator/COH-CC</param>
+    <param type="alg" name="Generator-7">genie::EventGenerator/COH-NC</param>
+    <param type="alg" name="Generator-8">genie::EventGenerator/DIS-CC-CHARM</param>
+    <param type="alg" name="Generator-9">genie::EventGenerator/QEL-CC-CHARM</param>
+    <param type="alg" name="Generator-10">genie::EventGenerator/NUE-EL</param>
+    <param type="alg" name="Generator-11">genie::EventGenerator/IMD</param>
+    <param type="alg" name="Generator-12">genie::EventGenerator/IMD-ANH</param>
+    <param type="alg" name="Generator-13">genie::EventGenerator/MEC-CC</param>
+    <param type="alg" name="Generator-14">genie::EventGenerator/MEC-NC</param>
+  </param_set>
+  <param_set name="Default+CCMEC">
     <param type="int" name="NGenerators">14</param>
     <param type="alg" name="Generator-0">genie::EventGenerator/QEL-CC</param>
     <param type="alg" name="Generator-1">genie::EventGenerator/QEL-NC</param>

This replaces Default+MEC with a set that include both MEC-CC and MEC-NC, and adds in an explicit Default+CCMEC set to allow generation similar to what was in R-2_10_X.

DefaultPlusMECWithNCHadWFix vs. EffSFTEMHadWFix

--- ./genie_xsec/v2_11_2/NULL/EffSFTEMHadWFix/data/UserPhysicsOptions.original.tidyReduced.xml    2016-10-17 16:20:39.198029289 -0500
+++ ./genie_xsec/v2_11_2/NULL/EffSFTEMHadWFix/data/UserPhysicsOptions.tidyReduced.xml    2016-10-17 16:20:39.405183442 -0500
@@ -137,9 +137,9 @@
     <param type="double" name="BBA05-Gmn-b4">0.</param>
     <param type="double" name="Charm-Mass">1.430</param>
     <param type="string" name="ResonanceNameList">P33(1232),S11(1535),D13(1520),S11(1650),D13(1700),D15(1675), S31(1620),D33(1700),P11(1440),P13(1720),F15(1680),P31(1910), P33(1920),F35(1905),F37(1950),P11(1710)</param>
-    <param type="alg" name="NuclearModel">genie::FGMBodekRitchie/Default</param>
+    <param type="alg" name="NuclearModel">genie::EffectiveSF/Default</param>
     <param type="bool" name="IsotopesUseSameNuclearModel">true</param>
-    <param type="bool" name="UseElFFTransverseEnhancement">false</param>
+    <param type="bool" name="UseElFFTransverseEnhancement">true</param>
     <param type="alg" name="TransverseEnhancement">genie::TransverseEnhancementFFModel/Default</param>
     <param type="double" name="RFG-MomentumCutOff">0.5000</param>
     <param type="double" name="RFG-NucRemovalE@Pdg=1000030060">0.0170</param>
@@ -318,9 +318,6 @@
     <param type="alg" name="XSecModel@genie::EventGenerator/AM-NUGAMMA">genie::H3AMNuGammaPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/DFR-CC">genie::ReinDFRPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/DFR-NC">genie::ReinDFRPXSec/Default</param>
-    <param type="alg" name="XSecModel@genie::EventGenerator/MEC-CC">genie::EmpiricalMECPXSec2015/Default</param>
-    <param type="alg" name="XSecModel@genie::EventGenerator/MEC-NC">genie::EmpiricalMECPXSec2015/Default</param>
-    <param type="alg" name="XSecModel@genie::EventGenerator/MEC-EM">genie::EmpiricalMECPXSec2015/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/GLRES">genie::GLRESPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/NucleonDecay">genie::DummyPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/NeutronOsc">genie::NOscDummyPXSec/Default</param>

DefaultPlusMECWithNCHadWfix vs. LocalFGNievesQEAndMECDiNucOffHadWFix

--- ./genie_xsec/v2_11_2/NULL/LocalFGNievesQEAndMECHadWFix/data/UserPhysicsOptions.original.tidyReduced.xml    2016-10-17 15:48:49.054296004 -0500
+++ ./genie_xsec/v2_11_2/NULL/LocalFGNievesQEAndMECHadWFix/data/UserPhysicsOptions.tidyReduced.xml    2016-10-17 15:49:04.402934007 -0500
@@ -34,7 +34,7 @@
     <param type="double" name="EL-Mv">0.840</param>
     <param type="double" name="QEL-Ma">0.990</param>
     <param type="double" name="QEL-Mv">0.840</param>
-    <param type="bool" name="SetDiNucleonCode">true</param>
+    <param type="bool" name="SetDiNucleonCode">false</param>
     <param type="double" name="EmpiricalMEC-Mq2d">0.4</param>
     <param type="double" name="EmpiricalMEC-Mass">2.1</param>
     <param type="double" name="EmpiricalMEC-Width">0.05</param>
@@ -137,7 +137,7 @@
     <param type="double" name="BBA05-Gmn-b4">0.</param>
     <param type="double" name="Charm-Mass">1.430</param>
     <param type="string" name="ResonanceNameList">P33(1232),S11(1535),D13(1520),S11(1650),D13(1700),D15(1675), S31(1620),D33(1700),P11(1440),P13(1720),F15(1680),P31(1910), P33(1920),F35(1905),F37(1950),P11(1710)</param>
-    <param type="alg" name="NuclearModel">genie::FGMBodekRitchie/Default</param>
+    <param type="alg" name="NuclearModel">genie::LocalFGM/Default</param>
     <param type="bool" name="IsotopesUseSameNuclearModel">true</param>
     <param type="bool" name="UseElFFTransverseEnhancement">false</param>
     <param type="alg" name="TransverseEnhancement">genie::TransverseEnhancementFFModel/Default</param>
@@ -297,7 +297,7 @@
     <param type="bool" name="DecayParticleWithCode=42112">true</param>
     <param type="bool" name="DecayParticleWithCode=42212">true</param>
     <param type="bool" name="DecayParticleWithCode=-24">false</param>
-    <param type="alg" name="XSecModel@genie::EventGenerator/QEL-CC">genie::LwlynSmithQELCCPXSec/Default</param>
+    <param type="alg" name="XSecModel@genie::EventGenerator/QEL-CC">genie::NievesQELCCPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/QEL-NC">genie::AhrensNCELPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/QEL-EM">genie::RosenbluthPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/RES-CC">genie::ReinSehgalRESPXSec/Default</param>
@@ -318,7 +318,7 @@
     <param type="alg" name="XSecModel@genie::EventGenerator/AM-NUGAMMA">genie::H3AMNuGammaPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/DFR-CC">genie::ReinDFRPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/DFR-NC">genie::ReinDFRPXSec/Default</param>
-    <param type="alg" name="XSecModel@genie::EventGenerator/MEC-CC">genie::EmpiricalMECPXSec2015/Default</param>
+    <param type="alg" name="XSecModel@genie::EventGenerator/MEC-CC">genie::NievesSimoVacasMECPXSec2016/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/MEC-NC">genie::EmpiricalMECPXSec2015/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/MEC-EM">genie::EmpiricalMECPXSec2015/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/GLRES">genie::GLRESPXSec/Default</param>
--- ./genie_xsec/v2_11_2/NULL/LocalFGNievesQEAndMECHadWFix/data/EventGenerator.original.tidyReduced.xml    2016-10-17 15:42:04.664127081 -0500
+++ ./genie_xsec/v2_11_2/NULL/LocalFGNievesQEAndMECHadWFix/data/EventGenerator.tidyReduced.xml    2016-10-17 15:49:04.376635885 -0500
@@ -20,19 +20,16 @@
   </param_set>
   <param_set name="QEL-CC">
     <param type="string" name="VldContext"></param>
-    <param type="int" name="NModules">12</param>
+    <param type="int" name="NModules">9</param>
     <param type="alg" name="Module-0">genie::InitialStateAppender/Default</param>
     <param type="alg" name="Module-1">genie::VertexGenerator/Default</param>
-    <param type="alg" name="Module-2">genie::FermiMover/Default</param>
-    <param type="alg" name="Module-3">genie::QELKinematicsGenerator/CC-Default</param>
-    <param type="alg" name="Module-4">genie::QELPrimaryLeptonGenerator/Default</param>
-    <param type="alg" name="Module-5">genie::QELHadronicSystemGenerator/Default</param>
-    <param type="alg" name="Module-6">genie::PauliBlocker/Default</param>
-    <param type="alg" name="Module-7">genie::UnstableParticleDecayer/BeforeHadronTransport</param>
-    <param type="alg" name="Module-8">genie::NucDeExcitationSim/Default</param>
-    <param type="alg" name="Module-9">genie::HadronTransporter/Default</param>
-    <param type="alg" name="Module-10">genie::NucBindEnergyAggregator/Default</param>
-    <param type="alg" name="Module-11">genie::UnstableParticleDecayer/AfterHadronTransport</param>
+    <param type="alg" name="Module-2">genie::QELEventGenerator/Default</param>
+    <param type="alg" name="Module-3">genie::PauliBlocker/Default</param>
+    <param type="alg" name="Module-4">genie::UnstableParticleDecayer/BeforeHadronTransport</param>
+    <param type="alg" name="Module-5">genie::NucDeExcitationSim/Default</param>
+    <param type="alg" name="Module-6">genie::HadronTransporter/Default</param>
+    <param type="alg" name="Module-7">genie::NucBindEnergyAggregator/Default</param>
+    <param type="alg" name="Module-8">genie::UnstableParticleDecayer/AfterHadronTransport</param>
     <param type="alg" name="ILstGen">genie::QELInteractionListGenerator/CC-Default</param>
   </param_set>
   <param_set name="QEL-NC">

DefaultPlusMECWithNC vs. AltPion

$ diff -u GXSPLINES-v2_11_2-AltPion/cfg/User*.tidyReduced.xml
--- GXSPLINES-v2_11_2-AltPion/cfg/UserPhysicsOptions.original.tidyReduced.xml    2016-10-07 15:34:04.953026968 -0500
+++ GXSPLINES-v2_11_2-AltPion/cfg/UserPhysicsOptions.tidyReduced.xml    2016-10-07 15:34:04.923719175 -0500
@@ -300,14 +300,14 @@
     <param type="alg" name="XSecModel@genie::EventGenerator/QEL-CC">genie::LwlynSmithQELCCPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/QEL-NC">genie::AhrensNCELPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/QEL-EM">genie::RosenbluthPXSec/Default</param>
-    <param type="alg" name="XSecModel@genie::EventGenerator/RES-CC">genie::ReinSehgalRESPXSec/Default</param>
-    <param type="alg" name="XSecModel@genie::EventGenerator/RES-NC">genie::ReinSehgalRESPXSec/Default</param>
-    <param type="alg" name="XSecModel@genie::EventGenerator/RES-EM">genie::ReinSehgalRESPXSec/Default</param>
+    <param type="alg" name="XSecModel@genie::EventGenerator/RES-CC">genie::BergerSehgalRESPXSec2014/Default</param>
+    <param type="alg" name="XSecModel@genie::EventGenerator/RES-NC">genie::BergerSehgalRESPXSec2014/Default</param>
+    <param type="alg" name="XSecModel@genie::EventGenerator/RES-EM">genie::BergerSehgalRESPXSec2014/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/DIS-CC">genie::QPMDISPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/DIS-NC">genie::QPMDISPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/DIS-EM">genie::QPMDISPXSec/Default</param>
-    <param type="alg" name="XSecModel@genie::EventGenerator/COH-CC">genie::ReinSehgalCOHPiPXSec/Default</param>
-    <param type="alg" name="XSecModel@genie::EventGenerator/COH-NC">genie::ReinSehgalCOHPiPXSec/Default</param>
+    <param type="alg" name="XSecModel@genie::EventGenerator/COH-CC">genie::BergerSehgalCOHPiPXSec2015/Default</param>
+    <param type="alg" name="XSecModel@genie::EventGenerator/COH-NC">genie::BergerSehgalCOHPiPXSec2015/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/IMD">genie::BardinIMDRadCorPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/IMD-ANH">genie::IMDAnnihilationPXSec/Default</param>
     <param type="alg" name="XSecModel@genie::EventGenerator/NUE-EL">genie::NuElectronPXSec/Default</param>
$ diff -u GXSPLINES-v2_11_2-AltPion/cfg/Event*.tidyReduced.xml
--- GXSPLINES-v2_11_2-AltPion/cfg/EventGeneratorListAssembler.original.tidyReduced.xml    2016-10-07 15:34:04.945929653 -0500
+++ GXSPLINES-v2_11_2-AltPion/cfg/EventGeneratorListAssembler.tidyReduced.xml    2016-10-07 15:34:04.938884195 -0500
@@ -1,7 +1,7 @@
 <?xml version="1.0" encoding="iso-8859-1"?>
 <alg_conf>
   <param_set name="Default">
-    <param type="int" name="NGenerators">13</param>
+    <param type="int" name="NGenerators">15</param>
     <param type="alg" name="Generator-0">genie::EventGenerator/QEL-CC</param>
     <param type="alg" name="Generator-1">genie::EventGenerator/QEL-NC</param>
     <param type="alg" name="Generator-2">genie::EventGenerator/RES-CC</param>
@@ -15,6 +15,8 @@
     <param type="alg" name="Generator-10">genie::EventGenerator/NUE-EL</param>
     <param type="alg" name="Generator-11">genie::EventGenerator/IMD</param>
     <param type="alg" name="Generator-12">genie::EventGenerator/IMD-ANH</param>
+    <param type="alg" name="Generator-13">genie::EventGenerator/DFR-CC</param>
+    <param type="alg" name="Generator-14">genie::EventGenerator/DFR-NC</param>
   </param_set>
   <param_set name="AnomalyMediatedNuGamma">
     <param type="int" name="NGenerators">1</param>