GENIE-MC · mroda88 · May 26, 2020 · May 29, 2020 · Jul 8, 2020 · Jul 14, 2020
diff --git a/config/AR23_20i/CommonParam.xml b/config/AR23_20i/CommonParam.xml
@@ -291,6 +291,7 @@ Or changing the name of this parameter set
       P31(1910),P33(1920),F35(1905),F37(1950),P11(1710)
     </param>
 
+    <param type="double" name="N-Delta-f_Pi">  0.0924 </param>
   </param_set>
 
   <param_set name="Lepton">

diff --git a/config/AR23_20i/ModelConfiguration.xml b/config/AR23_20i/ModelConfiguration.xml
@@ -120,6 +120,9 @@ STFC, Rutherford Appleton Laboratory
   <param type="alg" name="XSecModel@genie::EventGenerator/COH-NC-PION">       genie::AlvarezRusoCOHPiPXSec/Fast       </param>
   -->
 
+ <!-- Alvarez-Ruso et al COH gamma -->
+  <param type="alg" name="XSecModel@genie::EventGenerator/COH-NC-GAMMA">   genie::AlvarezRusoSalaCOHGammaPXSec/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 --git a/config/AR23_20i/TuneGeneratorList.xml b/config/AR23_20i/TuneGeneratorList.xml
@@ -47,7 +47,7 @@ Generator-%d                alg     No
 -->
 
   <param_set name="Default"> 
-     <param type="int" name="NGenerators">   18                                    </param>
+     <param type="int" name="NGenerators">   19                                    </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>
@@ -66,6 +66,7 @@ Generator-%d                alg     No
      <param type="alg" name="Generator-15">  genie::EventGenerator/QEL-CC-LAMBDA   </param>
      <param type="alg" name="Generator-16">  genie::EventGenerator/MEC-CC          </param>
      <param type="alg" name="Generator-17">  genie::EventGenerator/MEC-NC          </param>
+     <param type="alg" name="Generator-18">  genie::EventGenerator/COH-NC-GAMMA    </param>
   </param_set>
 
 </tune_generator_list>

diff --git a/config/AlvarezRusoCOHPiPXSec.xml b/config/AlvarezRusoCOHPiPXSec.xml
@@ -39,13 +39,13 @@ CoHAR_mA	 double   yes       Coherent axial mass        1.28
 <alg_conf>
 
   <param_set name="Default"> 
-      <param type="alg" name="XSec-Integrator"> genie::COHXSecAR/Default </param>
+      <param type="alg" name="XSec-Integrator"> genie::COHXSecAR/Pion </param>
       <param type="bool" name="UseLookupTable"> false </param>
 
    </param_set>
 
   <param_set name="Fast"> 
-      <param type="alg" name="XSec-Integrator"> genie::COHXSecAR/Fast </param>
+      <param type="alg" name="XSec-Integrator"> genie::COHXSecAR/Fast_Pion </param>
       <param type="bool" name="UseLookupTable"> false </param>
   </param_set>
 

diff --git a/config/AlvarezRusoSalaCOHGammaPXSec.xml b/config/AlvarezRusoSalaCOHGammaPXSec.xml
@@ -0,0 +1,43 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+
+<!--
+
+Configuration sets for the AlvarezRusoSalaCOHGammaPXSec
+
+Configurable Parameters:
+.........................................................................................................................
+Name                          Type     Optional   Comment                                        Default
+.........................................................................................................................
+COH-HadronicCurrents          vec-alg  No         Vector of hadronic current sub algos                         -
+XSec-Integrator               alg      No         Cross section integrator sub algo                            - 
+COH-FormFactor                alg      yes        DeVries form facor sub algo                    The one of the integrator              
+AsymptoticMaxGammaEnergy      double   No         Configuration parameter for maximum W for      taken from subalgo
+                                                  the cross section. 
+-->
+
+
+
+
+<alg_conf>
+
+  <param_set name="Default">
+
+    <param type="vec-alg" name="COH-HadronicCurrents" delim=";">
+      genie::COHDeltaCurrent/Default
+    </param>
+
+    <param type="alg"  name="XSec-Integrator"> genie::COHXSecAR/Gamma </param>
+
+  </param_set>
+
+  <param_set name="Consistent">
+
+    <param type="vec-alg" name="COH-HadronicCurrents" delim=";">
+      genie::COHDeltaCurrent/FermiMomentumTable
+    </param>
+
+  </param_set>
+
+
+
+</alg_conf>
diff --git a/config/COHDeltaCurrent.xml b/config/COHDeltaCurrent.xml
@@ -0,0 +1,33 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+
+<!--
+
+Configuration sets for the COHDeltaCurrent
+
+Configurable Parameters:
+..........................................................................................................
+Name                          Type     Optional   Comment                                        Default
+..........................................................................................................
+COH-Delta-FormFactor          alg      No         Delta Transistion form facort sub algo           -
+COH-Delta-Corrections         alg      No         Delta Corrections in medium sub algo             -
+-->
+
+
+
+<alg_conf>
+
+  <param_set name="Default">
+
+    <param type="alg"  name="COH-Delta-FormFactor"> genie::DeltaTransitionFormFactor/Default </param>
+    <param type="alg"  name="COH-Delta-Corrections"> genie::DeltaInMediumCorrections/Default </param>
+
+  </param_set>
+
+  <param_set name="FermiMomentumTable">
+
+    <param type="alg"  name="COH-Delta-Corrections"> genie::DeltaInMediumCorrections/FermiMomentumTable </param>
+
+  </param_set>
+
+
+</alg_conf>
diff --git a/config/COHGammaIntegrationLimits.xml b/config/COHGammaIntegrationLimits.xml
@@ -0,0 +1,32 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+
+
+
+<!--
+
+Configuration sets for the COHGammaIntegrationLimits
+
+Configurable Parameters:
+.....................................................................................................................
+Name                          Type     Optional   Comment                               Default
+.....................................................................................................................
+COH-FormFactor                alg      No         form facor sub algo                                  -      
+AsymptoticMaxGammaEnergy      double   No         Energy cut off for gamma energy [GeV]
+                                                  in the limit for infinit mass nuclei
+MaxGammaTheta                 double   yes        Theta cut off for gamma [deg]         180-epsilon                    
+
+-->
+
+
+
+<alg_conf>
+
+  <param_set name="Default">
+    <param type="double" name="AsymptoticMaxGammaEnergy"> 0.9 </param>
+    <param type="double" name="MaxGammaTheta"> 150. </param>
+
+    <param type="alg"  name="COH-FormFactor"> genie::COHProtonFormFactorInterpolation/Default </param>
+
+  </param_set>
+
+</alg_conf>
diff --git a/config/COHGammaKinematicsGenerator.xml b/config/COHGammaKinematicsGenerator.xml
@@ -0,0 +1,54 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+
+<!--
+Configuration for the COHKinematicsGenerator EventRecordVisitorI
+
+Algorithm Configurable Parameters:
+.......................................................................................................................
+Name                     Type       Opt     Comment                                        Default
+.......................................................................................................................
+UniformOverPhaseSpace    bool       Yes     kinematics uniformly over allowd phase space   false
+                                            wgt = (phase_space_volume)*(diff_xsec)/(xsec)
+MaxXSec-SafetyFactor     double     Yes     multiplies max xsec in rejection method        1.6
+MaxXSec-DiffTolerance    double     Yes     max allowed 200*(xsec-xsecmax)/(xsec+xsecmax)  999999.00 (disable)
+                                            if xsec>xsecmax
+Cache-MinEnergy          double     Yes     minimum energy for which max xsec is cached    0.00
+MinimScanPoint           vec-int    yes     number of points in each phase space           all 1
+                                            coordinate to be used for a scan 
+                                            for Maximum XSec
+ScaleScanPoint           vec-double No      The number of points used for the Maximum 
+                                            Xsec scan increases with energy. 
+                                            These are the linear coefficients of the
+                                            scaling  
+
+
+IntegrationLimits        algo       No      Specificies integration limit algorithm                                  
+                                            It is required for COH gamma configurations 
+tPhaseSpace              bool       yes     Generate the events in the t phase space       false 
+                                            instead of theta_l 
+
+MaxIterations            int        yes     Maximum number of iterations to try            controls::kRjMaxIterations
+                                            acceptance / rejection method
+-->
+
+<alg_conf>
+
+  <param_set name="Default"> 
+
+    <param type="double" name="MaxXSec-SafetyFactor"> 1.5 </param>
+
+    <param type="double" name="MaxXSec-DiffTolerance"> 1. </param>
+
+    <param type="int" name="MaxIterations"> 100000 </param>
+
+    <param type="vec-int" name="MinimScanPoint" delim=";" > 11 ; 175 ; 11 ; 3 </param> 
+
+    <param type="vec-double" name="ScaleScanPoint" delim=";" > 0.4 ; 0.8 ; 0.4 ; 0.01 </param>
+
+    <param type="alg" name ="IntegrationLimits"  > genie::COHGammaIntegrationLimits/Default </param>
+
+    </param_set>
+
+
+</alg_conf>
+
diff --git a/config/COHKinematicsGenerator.xml b/config/COHKinematicsGenerator.xml
diff --git a/config/COHPionKinematicsGenerator.xml b/config/COHPionKinematicsGenerator.xml
@@ -0,0 +1,38 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+
+<!--
+Configuration for the COHPionKinematicsGenerator EventRecordVisitorI
+
+Algorithm Configurable Parameters:
+.......................................................................................................................
+Name                     Type    Opt     Comment                                        Default
+.......................................................................................................................
+UniformOverPhaseSpace    bool    Yes     kinematics uniformly over allowd phase space   false
+                                         wgt = (phase_space_volume)*(diff_xsec)/(xsec)
+MaxXSec-SafetyFactor     double  Yes     multiplies max xsec in rejection method        1.6
+MaxXSec-DiffTolerance    double  Yes     max allowed 200*(xsec-xsecmax)/(xsec+xsecmax)  999999.00 (disable)
+                                         if xsec>xsecmax
+Cache-MinEnergy          double  Yes     minimum energy for which max xsec is cached    0.00
+
+COH-Ro                   double  No      Nuclear size scale                             CommonParam[Coherent]
+COH-Q2-min               double  No      Minimum considered Q^2 for Berger-Sehgal       CommonParam[Coherent]
+                                         coherent reactions when estimating the max 
+                                         cross section. Units in GeV^2.
+COH-Q2-max               double  No      Maximum considered Q^2 for Berger-Sehgal       CommonParam[Coherent]
+                                         coherent reactions when estimating the max 
+                                         cross section. Units in GeV^2.
+COH-t-max                double  No      Maximum considered t for Berger-Sehgal         CommonParam[Coherent]
+                                         coherent reactions when estimating 
+                                         the max cross section.  Units in GeV^2.
+-->
+
+<alg_conf>
+
+  <param_set name="Default"> 
+
+    <param type="string" name="CommonParam"> Coherent </param>
+
+  </param_set>
+
+</alg_conf>
+
diff --git a/config/COHProtonFormFactorInterpolation.xml b/config/COHProtonFormFactorInterpolation.xml
@@ -0,0 +1,27 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+
+<!--
+
+Configuration sets for the COHProtonFormFactorInterpolation
+
+Configurable Parameters:
+.........................................................................................................................
+Name                          Type     Optional   Comment                                                    Default
+.........................................................................................................................
+BaseCOHFormFactor             alg      No         Base form factor algorithm 
+AllowExtrapolation            bool     yes        Allow map extrapolation at higher pdg than map maximum     false
+-->
+
+
+
+<alg_conf>
+
+  <param_set name="Default">
+    <param type="alg" name="BaseCOHFormFactor" >
+      genie::DeVriesFormFactorMap/Default
+    </param>
+
+    <param type="bool" name="AllowExtrapolation"> false </param>
+  </param_set>
+
+</alg_conf>