GENIE-MC · jtenavidal · Oct 20, 2023 · Oct 17, 2025 · Oct 17, 2025 · Oct 17, 2025
diff --git a/config/BY00DISNuclearModel.xml b/config/BY00DISNuclearModel.xml
@@ -0,0 +1,17 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+
+<alg_conf>
+
+<!--
+Configuration for the BY DIS Nuclear model 2000. 
+In this version, parameters are hardcoded. This model is kept for historical reasons.                                                                        
+This model is not advised, as a new version of this parameterization is available in BY21DISNuclearModel 
+-->
+
+ <param_set name="Default"> 
+
+
+ </param_set>
+
+</alg_conf>
+
diff --git a/config/BY21DISNuclearModel.xml b/config/BY21DISNuclearModel.xml
@@ -0,0 +1,78 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+
+<alg_conf>
+
+<!--
+Configuration for the BY DIS Nuclear model arxiv.org/pdf/2108.09240
+
+Configurable Parameters:
+...............................................................................................................................................
+Name                            Type    Optional   Comment                     
+...............................................................................................................................................
+BY21-NuclModel-2H-Scale         double  No         Overall scaling factor applied to deuterium to free nucleon correction
+BY21-NuclModel-2H-f0            double  No         Constant multiplying x_F^0 contribution of the 2H to free nucleon parameterization 
+BY21-NuclModel-2H-f1            double  No         Constant multiplying x_F^1 contribution of the 2H to free nucleon parameterization 
+BY21-NuclModel-2H-f2            double  No         Constant multiplying x_F^2 contribution of the 2H to free nucleon parameterization 
+BY21-NuclModel-2H-f3            double  No         Constant multiplying x_F^3 contribution of the 2H to free nucleon parameterization 
+BY21-NuclModel-2H-f4            double  No         Constant multiplying x_F^4 contribution of the 2H to free nucleon parameterization 
+BY21-NuclModel-2H-f5            double  No         Constant multiplying x_F^5 contribution of the 2H to free nucleon parameterization 
+
+BY21-NuclModel-Fe-f0            double  No         Constant used in the Isoscalar Iron to Deuterium parameterization
+BY21-NuclModel-Fe-f1            double  No         Constant used in the Isoscalar Iron to Deuterium parameterization
+BY21-NuclModel-Fe-f2            double  No         Constant used in the Isoscalar Iron to Deuterium parameterization
+BY21-NuclModel-Fe-f3            double  No         Constant used in the Isoscalar Iron to Deuterium parameterization
+BY21-NuclModel-Fe-f4            double  No         Constant used in the Isoscalar Iron to Deuterium parameterization
+
+BY21-NuclModel-C-f0            double  No         Constant used in the Isoscalar Carbon to Iron parameterization
+BY21-NuclModel-C-f1            double  No         Constant used in the Isoscalar Carbon to Iron parameterization
+BY21-NuclModel-C-f2            double  No         Constant used in the Isoscalar Carbon to Iron parameterization
+BY21-NuclModel-C-f3            double  No         Constant used in the Isoscalar Carbon to Iron parameterization
+BY21-NuclModel-C-f4            double  No         Constant used in the Isoscalar Carbon to Iron parameterization
+BY21-NuclModel-C-f5            double  No         Constant used in the Isoscalar Carbon to Iron parameterization
+
+BY21-NuclModel-Au-f0            double  No         Constant used in the Isoscalar Au,Pb  to Iron parameterization
+BY21-NuclModel-Au-f1            double  No         Constant used in the Isoscalar Au,Pb  to Iron parameterization
+BY21-NuclModel-Au-f2            double  No         Constant used in the Isoscalar Au,Pb  to Iron parameterization
+BY21-NuclModel-Au-f3            double  No         Constant used in the Isoscalar Au,Pb  to Iron parameterization
+BY21-NuclModel-Au-f4            double  No         Constant used in the Isoscalar Au,Pb  to Iron parameterization
+BY21-NuclModel-Au-f5            double  No         Constant used in the Isoscalar Au,Pb  to Iron parameterization
+BY21-NuclModel-Au-f6            double  No         Constant used in the Isoscalar Au,Pb  to Iron parameterization
+
+
+-->
+
+ <param_set name="Default"> 
+   <param type="double" name="BY21-NuclModel-2H-Scale">  0.985   </param>
+   <param type="double" name="BY21-NuclModel-2H-f0">     1       </param>
+   <param type="double" name="BY21-NuclModel-2H-f1">     0.422   </param>
+   <param type="double" name="BY21-NuclModel-2H-f2">     -2.745  </param>
+   <param type="double" name="BY21-NuclModel-2H-f3">     7.570   </param>
+   <param type="double" name="BY21-NuclModel-2H-f4">     -10.335 </param>
+   <param type="double" name="BY21-NuclModel-2H-f5">     5.422   </param>
+
+   <param type="double" name="BY21-NuclModel-Fe-f0">     1.096   </param>
+   <param type="double" name="BY21-NuclModel-Fe-f1">    -0.38    </param>
+   <param type="double" name="BY21-NuclModel-Fe-f2">    -0.30    </param>
+   <param type="double" name="BY21-NuclModel-Fe-f3">    -23.0    </param>
+   <param type="double" name="BY21-NuclModel-Fe-f4">     8.0     </param>
+
+   <param type="double" name="BY21-NuclModel-C-f0">     0.919   </param>
+   <param type="double" name="BY21-NuclModel-C-f1">     1.844   </param>
+   <param type="double" name="BY21-NuclModel-C-f2">    -12.73   </param>
+   <param type="double" name="BY21-NuclModel-C-f3">     36.89   </param>
+   <param type="double" name="BY21-NuclModel-C-f4">    -46.77   </param>
+   <param type="double" name="BY21-NuclModel-C-f5">     21.22   </param>
+
+   <param type="double" name="BY21-NuclModel-Au-f0">     0.932   </param>
+   <param type="double" name="BY21-NuclModel-Au-f1">     2.461   </param>
+   <param type="double" name="BY21-NuclModel-Au-f2">    -24.23   </param>
+   <param type="double" name="BY21-NuclModel-Au-f3">    101.03   </param>
+   <param type="double" name="BY21-NuclModel-Au-f4">   -203.47   </param>
+   <param type="double" name="BY21-NuclModel-Au-f5">    193.85   </param>
+   <param type="double" name="BY21-NuclModel-Au-f6">    -69.82   </param>
+
+
+ </param_set>
+
+</alg_conf>
+
diff --git a/config/BYPDF.xml b/config/BYPDF.xml
@@ -14,6 +14,8 @@ BY-X1                      double  No    Bodek-Yang param X1
 BY-X2                      double  No    Bodek-Yang param X2                    
 PDF-Q2min                  double  Yes   Min Q2 for PDF evaluation              from "Uncorr-PDF-Set" register
 Uncorr-PDF-Set             alg     No    Uncorrected PDF model 
+BY-UpScale                 double  Yes   Scaling the up distribution. Default 0.
+By-DownScale               double  Yes   Scaling the down distribution. Default 0.
 -->
 
   <param_set name="Default"> 
@@ -26,9 +28,7 @@ Uncorr-PDF-Set             alg     No    Uncorrected PDF model
 
   </param_set>
 
-
   <param_set name="Optional"> 
-
     <param type="alg" name="Uncorr-PDF-Set">  genie::LHAPDF5/GRVLO  </param>
 
   </param_set>

diff --git a/config/BYStrucFunc.xml b/config/BYStrucFunc.xml
@@ -39,30 +39,32 @@ WeinbergAngle              double  No
   <param_set name="Default"> 
     <param type="string" name="CommonParam"> CKM,Masses,WeakInt </param>
 
-     <param type="alg" name="PDF-Set"> genie::BYPDF/Default </param>
+    <param type="alg" name="PDF-Set"> genie::BYPDF/Default </param>
 
-     <param type="double" name="BY-A">     0.538   </param>
-     <param type="double" name="BY-B">     0.305   </param>
-     <param type="double" name="BY-CsU">   0.363   </param>
-     <param type="double" name="BY-CsD">   0.621   </param>
-     <param type="double" name="BY-Cv1U">  0.291   </param>
-     <param type="double" name="BY-Cv2U">  0.189   </param>
-     <param type="double" name="BY-Cv1D">  0.202   </param>
-     <param type="double" name="BY-Cv2D">  0.255   </param>
+    <param type="double" name="BY-A">     0.538   </param>
+    <param type="double" name="BY-B">     0.305   </param>
+    <param type="double" name="BY-CsU">   0.363   </param>
+    <param type="double" name="BY-CsD">   0.621   </param>
+    <param type="double" name="BY-Cv1U">  0.291   </param>
+    <param type="double" name="BY-Cv2U">  0.189   </param>
+    <param type="double" name="BY-Cv1D">  0.202   </param>
+    <param type="double" name="BY-Cv2D">  0.255   </param>
 
-     <!-- 
-	  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-	  Switches controlling the DIS structure function calculation. Those switches specify whether to:
-	  - include an R (~FL) factor into the calculation
-	  - include a nuclear factor (accounting for shadowing/anti-shadowing)
-	  - include corrections for calculating relation between 2xF1 and F2
-	  - Select value for Q2 cutoff in relation between 2xF1 and F2
-     -->
-     <param type="bool"   name="IncludeR">           true  </param>
-     <param type="bool"   name="IncludeNuclMod">     true  </param>
-     <param type="bool"   name="Use2016Corrections"> false </param>
-     <param type="double" name="LowQ2CutoffF1F2">    0.8   </param>
-
+    <!-- 
+	 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+	 Switches controlling the DIS structure function calculation. Those switches specify whether to:
+	 - include an R (~FL) factor into the calculation
+	 - include a nuclear factor (accounting for shadowing/anti-shadowing)
+	 - include corrections for calculating relation between 2xF1 and F2
+	 - Select value for Q2 cutoff in relation between 2xF1 and F2
+    -->
+    <param type="bool"   name="IncludeR">           true  </param>
+    <param type="bool"   name="IncludeNuclMod">     true  </param>
+    <param type="bool"   name="Use2016Corrections"> false </param>
+    <param type="double" name="LowQ2CutoffF1F2">    0.8   </param>
+
+    <!-- Specify DIS nuclear correction to use-->
+    <param type="alg" name="DISNuclModel">  genie::BYDISNuclearModel/Default </param>
   </param_set>
 
 </alg_conf>

diff --git a/config/BYStrucFunc2021.xml b/config/BYStrucFunc2021.xml
@@ -0,0 +1,96 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+
+<!--
+Configuration for the BYStrucFunc DISStructureFuncModelI
+
+Algorithm Configurable Parameters:
+....................................................................................................
+Name                       Type    Opt   Comment                                Default
+....................................................................................................
+BY-A                       double  No    Bodek-Yang scaling var param A         
+By-B                       double  No    Bodek-Yang scaling var param B         
+BY-CsU                     double  No    Bodek-Yang u sea K factor param        
+BY-CsD                     double  No    Bodek-Yang d sea K factor param        
+BY-Cv1U                    double  No    Bodek-Yang u val K factor param        
+BY-Cv2U                    double  No    Bodek-Yang u val K factor param        
+BY-Cv1D                    double  No    Bodek-Yang d val K factor param        
+BY-Cv2D                    double  No    Bodek-Yang d val K factor param        
+BY-PsA                     double  No    Bodek-Yang P sea axial param
+BY-PvA                     double  No    Bodek-Yang P sea valance param
+BY-PsC                     double  No    Bodek-Yang C sea axial param
+BY-CaLW-nu                 double  No    Low-q0 axial neutrino parameter
+BY-CaLW-nubar              double  No    Low-q0 axial anti-neutrino parameter
+BY-H0                      double  No    Bodek-Yang high order QCD parameter
+BY-H1                      double  No    Bodek-Yang high order QCD parameter
+BY-H2                      double  No    Bodek-Yang high order QCD parameter
+BY-H3                      double  No    Bodek-Yang high order QCD parameter
+
++++ This model inherits from QPMDISStrucFunc and calls its configure
++++ The following parameters have to be the same as in QPMDISStrucFunc.xml
+
+PDF-Set                    alg     No    PDF model
+Charm-Mass                 double  No    charm mass                             CommonParam[Masses]
+Charm-Prod-Off             bool    Yes   charm production is turned off         false
+PDF-Q2min                  double  No    min Q2 for PDF evaluation              From "PDF-Set" 
+IncludeR                   bool    No    include R factor?              
+IncludeNuclMod             bool    No    include nuclear factor?            
+CKM-Vcd                    double  No                                           CommonParam[CKM]
+CKM-Vcs                    double  No                                           CommonParam[CKM]
+CKM-Vud                    double  No                                           CommonParam[CKM]
+CKM-Vus                    double  No                                           CommonParam[CKM]
+LowQ2CutoffF1F2            double  No    min for F1/F2 SF relation             
+WeinbergAngle              double  No                                           CommonParam[WeakInt]
+-->
+
+<alg_conf>
+
+  <param_set name="Default"> 
+    <param type="string" name="CommonParam"> CKM,Masses,WeakInt </param>
+    <param type="alg" name="PDF-Set"> genie::BYPDF/Default </param>
+    <param type="double" name="BY-A">     0.621   </param>
+    <param type="double" name="BY-B">     0.380   </param>
+    <param type="double" name="BY-CvLW">  0.218   </param>
+    <param type="double" name="BY-CsU">   0.369   </param>
+    <param type="double" name="BY-CsD">   0.561   </param>
+    <param type="double" name="BY-Cv1U">  0.417   </param>
+    <param type="double" name="BY-Cv2U">  0.264   </param>
+    <param type="double" name="BY-Cv1D">  0.341   </param>
+    <param type="double" name="BY-Cv2D">  0.323   </param>
+    <param type="double" name="BY-CsS">   0.561   </param>
+    <param type="double" name="BY-PsA">   0.55   </param>
+    <param type="double" name="BY-PvA">   0.30   </param>
+    <param type="double" name="BY-CsA">   0.75   </param>
+    <param type="double" name="BY-CaLW-nu"> 0.436 </param>
+    <param type="double" name="BY-CaLW-nubar"> 0.654 </param>
+    <param type="double" name="BY-H0">    0.914   </param>
+    <param type="double" name="BY-H1">    0.296   </param>
+    <param type="double" name="BY-H2">    -0.374   </param>
+    <param type="double" name="BY-H3">    0.165   </param>
+    <!-- 
+	 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+	 Switches controlling the DIS structure function calculation. Those switches specify whether to:
+	 - include an R (~FL) factor into the calculation
+	 - include a nuclear factor (accounting for shadowing/anti-shadowing)
+	 - include corrections for calculating relation between 2xF1 and F2
+	 - Select value for Q2 cutoff in relation between 2xF1 and F2
+    -->
+    <param type="bool"   name="IncludeR">           true  </param>
+    <param type="bool"   name="IncludeH">           true  </param>
+    <param type="bool"   name="IncludeNuclMod">     true  </param>
+    <param type="double" name="LowQ2CutoffF1F2">    0.8   </param>
+    <param type="bool"   name="Use2016Corrections"> true </param>
+    <!-- Specify DIS nuclear correction to use
+	 In this case we use the A-Dependent parameterization by J.Gomez et. al. 
+    -->
+    <param type="alg" name="DISNuclModel">  genie::JGDISNuclearModel/Default </param>
+  </param_set>
+
+  <param_set name="GRV98_2021">
+    <param type="alg" name="DISNuclModel">  genie::BY21DISNuclearModel/Default </param>
+  </param_set>
+
+  <param_set name="GRV98_2000">
+    <param type="alg" name="DISNuclModel">  genie::BY00DISNuclearModel/Default </param>
+  </param_set>
+</alg_conf>
+
diff --git a/config/BergerSehgalRESPXSec2014.xml b/config/BergerSehgalRESPXSec2014.xml
@@ -48,7 +48,6 @@ RES-Mb                      double  Yes   Mb parameter used for the GA calculati
   <param_set name="Default"> 
 
     <param type="string"  name="CommonParam"> WeakInt,NonResBackground,CKM,FermiGas </param>
-
 
      <!--
 	 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

diff --git a/config/DISHadronicSystemGenerator.xml b/config/DISHadronicSystemGenerator.xml
@@ -40,7 +40,7 @@ NUCL-NR             double  No
   </param_set>
 
   <param_set name="Charm">
-     <param type="alg"  name="Hadronizer">          genie::AGCharmPythia6Hadro2023/Default                     </param>
+     <param type="alg"  name="Hadronizer">          genie::AGCharmPythia8Hadro2023/Default                     </param>
      <param type="alg"  name="PreTransportDecayer"> genie::UnstableParticleDecayer/BeforeHadronTransport  </param>
      <param type="bool" name="FilterPreFragm">      false                                                 </param>
 

diff --git a/config/DISXSec.xml b/config/DISXSec.xml
@@ -2,6 +2,16 @@
 
 <!--
 Configuration for the DISXSec cross section algorithm
+
+Algorithm Configurable Parameters:
+..........................................................................................................................
+Name                        Type    Opt   Comment                                             Default
+..........................................................................................................................
+gsl-integration-type        string  Yes                                                       adaptative
+gsl-max-eval                int     Yes                                                       500000
+gsl-min-eval                int     Yes                                                       10000
+gsl-relative-tolerance      double  Yes                                                       1E-2
+DIS-NuclCorr                bool    Yes   Apply DIS Nuclear Correction                        false 
 -->
 
 <alg_conf>
@@ -15,5 +25,10 @@ Configuration for the DISXSec cross section algorithm
     <param type="double" name ="gsl-relative-tolerance">    0.0001  </param>
   </param_set>
 
+  <param_set name="GRV982010">
+    <param type="int"    name ="gsl-max-eval" >            50000000  </param>
+    <param type="double" name ="gsl-relative-tolerance">    0.00001  </param>
+    <param type="bool"   name="DIS-NuclCorr"> true </param>
+  </param_set>
 </alg_conf>
 
diff --git a/config/JGDISNuclearModel.xml b/config/JGDISNuclearModel.xml
@@ -0,0 +1,45 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+
+<alg_conf>
+
+<!--
+Configuration for the J. Gomez et. al. DIS Nuclear model journals.aps.org/prd/pdf/10.1103/PhysRevD.49.4348
+
+Configurable Parameters:
+.........................................................................................................
+Name                       Type    Optional   Comment                    
+.........................................................................................................
+JG-NuclModel-c0            double  No         
+JG-NuclModel-c1            double  No         
+JG-NuclModel-c2            double  No         
+
+JG-NuclModel-a0            double  No         
+JG-NuclModel-a1            double  No         
+JG-NuclModel-a2            double  No         
+JG-NuclModel-a3            double  No         
+JG-NuclModel-a4            double  No         
+JG-NuclModel-a5            double  No         
+JG-NuclModel-a6            double  No         
+JG-NuclModel-a7            double  No         
+JG-NuclModel-a8            double  No         
+-->
+
+ <param_set name="Default"> 
+   <param type="double" name="JG-NuclModel-c0">     0.017   </param>
+   <param type="double" name="JG-NuclModel-c1">     0.018   </param>
+   <param type="double" name="JG-NuclModel-c2">     0.005   </param>    
+
+   <param type="double" name="JG-NuclModel-a0">     -0.070  </param>
+   <param type="double" name="JG-NuclModel-a1">      2.189  </param>
+   <param type="double" name="JG-NuclModel-a2">    -24.667  </param>
+   <param type="double" name="JG-NuclModel-a3">    145.291  </param>
+   <param type="double" name="JG-NuclModel-a4">   -497.237  </param>
+   <param type="double" name="JG-NuclModel-a5">   1013.129  </param>
+   <param type="double" name="JG-NuclModel-a6">  -1208.393  </param>
+   <param type="double" name="JG-NuclModel-a7">    775.767  </param>
+   <param type="double" name="JG-NuclModel-a8">   -205.872  </param>
+
+ </param_set>
+
+</alg_conf>
+
-Original file line number
+Diff line change
@@ Expand Up / @@ -48,7 +48,6 @@ RES-Mb double Yes Mb parameter used for the GA calculati @@
       <param_set name="Default">
         <param type="string"  name="CommonParam"> WeakInt,NonResBackground,CKM,FermiGas </param>
          <!--
     	 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ Expand Down @@