<?xml version="1.0" encoding="ISO-8859-1"?>

<global_param_list>

<!--
***************************************************************************************************
This file defines defaults for the most important GENIE physics user configuration parameters. 

Each GENIE algorithm would pick the value of a configuration parameter from here if it is not 
set at its own configuration set. Typically, in the "Default" configuration set of each algorithm, 
these params are left undefined so that they can be set from a single place (here).

Multilple named lists of user physics options can be defined. 

Use the "Default" one to get the best GENIE physics & tuning.
***************************************************************************************************
-->

 <param_set name="Default"> 

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  CKM quark mixing parameters 
  Use values listed at the July 2006 particle physics booklet
  -->
  <param type="double" name="CKM-Vud">  0.97377  </param>
  <param type="double" name="CKM-Vus">  0.2257   </param>
  <param type="double" name="CKM-Vcd">  0.230    </param>  <!-- 0.2239 -->
  <param type="double" name="CKM-Vcs">  0.957    </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Cabbibo and Weinberg angles 
  -->
  <param type="double" name="CabbiboAngle">  0.22853207 </param> <!-- 0.226893 -->  
  <param type="double" name="WeinbergAngle"> 0.49744211 </param> <!-- 0.501568 --> 

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Constants in Feynman-Kislinger-Ravndal model for baryon excitation used in the Rein-Seghal 
  resonance model
  -->
  <param type="double" name="RS-Omega"> 1.05  </param>
  <param type="double" name="RS-Zeta">  0.762 </param> <!-- 0.75 -->

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Constants in the Bodek-Yang DIS structure function model and in the Bodek-Yang GRV-LO-9* PDF 
  corrections
  -->
  <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-X0">   -0.00817 </param>
  <param type="double" name="BY-X1">    0.0506  </param>
  <param type="double" name="BY-X2">    0.0798  </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Minimum Q2 for PDF evaluation - PDF evolution frozen for Q2<Q2min
  -->
  <param type="double" name="PDF-Q2min"> 0.800 </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) 
  -->
  <param type="bool" name="DISSF-IncludeR">       true </param>
  <param type="bool" name="DISSF-IncludeNuclMod"> true </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  -->
  <param type="double" name="DIS-XSecScale"> 1.032 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  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>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Axial and vector masses for elastic scattering  
  -->
  <param type="double" name="EL-Ma">  0.990  </param>
  <param type="double" name="EL-Mv">  0.840  </param> 

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Axial and vector masses for quasi-elastic scattering  
  -->
  <param type="double" name="QEL-Ma"> 0.990 </param>
  <param type="double" name="QEL-Mv"> 0.840 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Axial and vector masses for resonance production  
  -->
  <param type="double" name="RES-Ma"> 1.120 </param>
  <param type="double" name="RES-Mv"> 0.840 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Axial mass for coherent scattering 
  -->
  <param type="double" name="COH-Ma"> 1.000 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Value of axial form factor at Q2=0  
  -->
  <param type="double" name="QEL-FA0"> -1.2670 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameter eta controlling the strange axial form factor 
  -->
  <param type="double" name="EL-Axial-Eta"> 0.12 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Nuclear scale parameter in coherent scattering (in fm)
  -->
  <param type="double" name="COH-Ro"> 1.000 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Real over imaginary fwd scattering amplitude in coherent interactions 
  -->
  <param type="double" name="COH-ReImAmpl"> 0.300 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Use a modification to Adler's PCAC theorem that includes the final state lepton mass. 
  See Rein-Seghal, hep-ph/0606185. Affects the COH CC cross section.
  -->
  <param type="bool" name="COH-UseModifiedPCAC"> true </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameters for neutrino-electron elastic scattering 
  -->

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameters for anomaly mediated neutrino-gamma interactions
  -->
 <param type="double" name="AMNuGamma-Gw"> 20.0 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  BBA2003 elastic form factor parameters
  - BBA03-Gen-* are parameters for the Krutov Gen parameterization used in BBA2003
  - BBA03-Gep/Gmp/Gmn-* are coefficients for the BBA2003 inverse polynomial fit for Gep,Gmp,Gmn
  - BBA03-Q2Max is a Q2 value above which Gep/Gmp is assummed const
  H.Budd, NuINT-02 proceedings
  -->
  <param type="double" name="BBA03-Gen-a">    0.942       </param>
  <param type="double" name="BBA03-Gen-b">    4.610       </param>
  <param type="double" name="BBA03-Gep-a2">   3.253       </param>
  <param type="double" name="BBA03-Gep-a4">   1.422       </param>
  <param type="double" name="BBA03-Gep-a6">   0.08582     </param>
  <param type="double" name="BBA03-Gep-a8">   0.3318      </param>
  <param type="double" name="BBA03-Gep-a10"> -0.09371     </param>
  <param type="double" name="BBA03-Gep-a12">  0.01076     </param>
  <param type="double" name="BBA03-Gmp-a2">   3.104       </param>
  <param type="double" name="BBA03-Gmp-a4">   1.428       </param>
  <param type="double" name="BBA03-Gmp-a6">   0.1112      </param>                               
  <param type="double" name="BBA03-Gmp-a8">  -0.006981    </param>
  <param type="double" name="BBA03-Gmp-a10">  0.0003705   </param>
  <param type="double" name="BBA03-Gmp-a12"> -0.7063E-5   </param>
  <param type="double" name="BBA03-Gmn-a2">   3.043       </param>
  <param type="double" name="BBA03-Gmn-a4">   0.8548      </param>
  <param type="double" name="BBA03-Gmn-a6">   0.6806      </param>
  <param type="double" name="BBA03-Gmn-a8">  -0.1287      </param>
  <param type="double" name="BBA03-Gmn-a10">  0.008912    </param>
  <param type="double" name="BBA03-Gmn-a12">  0.          </param>
  <param type="double" name="BBA03-Q2Max">    6.0         </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  BBA2005 elastic form factor parameters
  R.Bradford, A.Bodek, H.Budd and J.Arrington, hep-ex/0602017
  -->
  <param type="double" name="BBA05-Gep-a0">     1.      </param>
  <param type="double" name="BBA05-Gep-a1">    -0.0578  </param>  <!--  +/- 0.1660 -->
  <param type="double" name="BBA05-Gep-a2">     0.      </param>
  <param type="double" name="BBA05-Gep-b1">    11.100   </param>  <!--  +/- 0.217  -->
  <param type="double" name="BBA05-Gep-b2">    13.60    </param>  <!--  +/- 1.39   -->
  <param type="double" name="BBA05-Gep-b3">    33.00    </param>  <!--  +/- 8.95   -->
  <param type="double" name="BBA05-Gep-b4">     0.      </param> 
  <param type="double" name="BBA05-Gmp-a0">     1.      </param> 
  <param type="double" name="BBA05-Gmp-a1">     0.1500  </param>  <!--  +/- 0.0312 -->
  <param type="double" name="BBA05-Gmp-a2">     0.      </param>
  <param type="double" name="BBA05-Gmp-b1">    11.100   </param>  <!--  +/- 0.103  -->
  <param type="double" name="BBA05-Gmp-b2">    19.600   </param>  <!--  +/- 0.281  -->
  <param type="double" name="BBA05-Gmp-b3">     7.540   </param>  <!--  +/- 0.967  -->
  <param type="double" name="BBA05-Gmp-b4">     0.      </param>
  <param type="double" name="BBA05-Gen-a0">     0.      </param>
  <param type="double" name="BBA05-Gen-a1">     1.250   </param>  <!--  +/- 0.368  -->
  <param type="double" name="BBA05-Gen-a2">     1.30    </param>  <!--  +/- 1.99   -->
  <param type="double" name="BBA05-Gen-b1">    -9.86    </param>  <!--  +/- 6.46   -->
  <param type="double" name="BBA05-Gen-b2">   305.0     </param>  <!--  +/- 28.6   -->
  <param type="double" name="BBA05-Gen-b3">  -758.0     </param>  <!--  +/- 77.5   -->
  <param type="double" name="BBA05-Gen-b4">   802.0     </param>  <!--  +/- 156.   -->
  <param type="double" name="BBA05-Gmn-a0">     1.      </param>
  <param type="double" name="BBA05-Gmn-a1">     1.810   </param>  <!--  +/- 0.402  -->
  <param type="double" name="BBA05-Gmn-a2">     0.      </param>
  <param type="double" name="BBA05-Gmn-b1">    14.100   </param>  <!--  +/- 0.597  -->
  <param type="double" name="BBA05-Gmn-b2">    20.70    </param>  <!--  +/- 2.55   -->
  <param type="double" name="BBA05-Gmn-b3">    68.7     </param>  <!--  +/- 14.1   -->
  <param type="double" name="BBA05-Gmn-b4">     0.      </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Charm mass 
  -->
  <param type="double" name="Charm-Mass"> 1.430 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  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),P13(1720),F15(1680),P31(1910),
         P33(1920),F35(1905),F37(1950),P11(1710) </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Nuclear model selection. 
  Options:
   - genie::FGMBodekRitchie/Default
   - genie::SpectralFunc1d/Default
  The 'NuclearModel' option defines the default basic model which should work for _any_ nuclei 
  (typically a Fermi Gas model with the Bodek-Ritchie NN corellatin tail). Refinements for specific 
  nuclei are possible, by specifying the 'NuclearModel@Pdg=10LZZZAAAI' option. 
  Currently the same nuclear model is forced for all isotopes.
  -->
  <param type="alg"  name="NuclearModel">                 genie::FGMBodekRitchie/Default </param>
<!--
  <param type="alg" name="NuclearModel@Pdg=1000060120">   genie::SpectralFunc1d/Default  </param>
  <param type="alg" name="NuclearModel@Pdg=1000260560">   genie::SpectralFunc1d/Default  </param>
-->
  <param type="bool" name="IsotopesUseSameNuclearModel">  true                           </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  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>
  <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=1000260560">  0.0360 </param> <!-- Fe56   -->
  <param type="double" name="RFG-NucRemovalE@Pdg=1000280580">  0.0360 </param> <!-- Ni58   -->
  <param type="double" name="RFG-NucRemovalE@Pdg=1000822080">  0.0440 </param> <!-- Pb208  -->
  <param type="bool"   name="RFG-IsotopesHaveSameRemovalE">    true  </param> 

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameters related with GENIE's 1-d/effective spectral function implementation  
  -->
  <param type="bool"   name="SF1d-UseRFGRemovalE">       false </param> 
  <param type="bool"   name="SF1d-UseRFGMomentumCutoff"> true  </param> 

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameters related with GENIE's full spectral function implementation  
  -->

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  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>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Options relevant for the simulation of the 2nd nucleon when the scattering if off a correlated
  NN pair
  -->
  <param type="bool"   name="CorrelNN-Enable">            false </param>
  <param type="double" name="CorrelNN-MomentumThreshold"> 0.3   </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameters controlling corrections to final state primary lepton simulation
  -->
  <param type="bool" name="FSL-ApplyCoulombCorrection"> false </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameters controlling the hadronic multiplicities in the KNO model
  -->
  <param type="double" name="KNO-Alpha-vp">  0.40 </param> 
  <param type="double" name="KNO-Alpha-vn"> -0.20 </param> 
  <param type="double" name="KNO-Alpha-vbp"> 0.02 </param> 
  <param type="double" name="KNO-Alpha-vbn"> 0.80 </param> 
  <param type="double" name="KNO-Beta-vp">   1.42 </param> 
  <param type="double" name="KNO-Beta-vn">   1.42 </param> 
  <param type="double" name="KNO-Beta-vbp">  1.28 </param> 
  <param type="double" name="KNO-Beta-vbn">  0.95 </param> 

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Levy function (KNO parameterization) parameter c at kno(z) = 2*exp(-c)*pow(c,cz+1)/Gamma(cz+1)
  v+p    : 7.93 +/- 0.34 source: Tingjun's fit
  v+n    : 5.22 +/- 0.15 source: Tingjun's fit
  vbar+p : same as vn
  vbar+n : same as vp
  -->
  <param type="double" name="KNO-LevyC-vp">  7.93 </param> 
  <param type="double" name="KNO-LevyC-vn">  5.22 </param> 
  <param type="double" name="KNO-LevyC-vbp"> 5.22 </param> 
  <param type="double" name="KNO-LevyC-vbn"> 7.93 </param> 

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  KNO hadronization model probabilities for producing hadron pairs
  -->
  <param type="double" name="KNO-ProbPi0Pi0">        0.30 </param>
  <param type="double" name="KNO-ProbPiplusPiminus"> 0.60 </param>
  <param type="double" name="KNO-ProbKplusKminus">   0.05 </param>
  <param type="double" name="KNO-ProbK0K0bar">       0.05 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameters controlling whether to reweight the KNO phase space decay reweighting and the actual
  reweighting function.
  See: A.B.Clegg, A.Donnachie, A Description of Jet Structure by pT-limited Phase Space.
  -->
  <param type="bool"   name="KNO-PhaseSpDec-Reweight">     true  </param>
  <param type="double" name="KNO-PhaseSpDec-ReweightParm"> 3.5   </param> 

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameters controlling whether to use the baryon xF and pT2 pdfs in the KNO hadronization.
  The option to use isotropic, non-reweighted 2-body phase space decays is used for compatibility
  with neugen/daikon.
  -->
  <param type="bool" name="KNO-UseBaryonPdfs-xFpT2">  true  </param>
  <param type="bool" name="KNO-UseIsotropic2BodyDec"> true  </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameters controlling the PYTHIA hadronization model:
  - ssbar:   [PYTHIA def:0.30    ] [NUX tune:0.21    ]
  - pt2:     [PYTHIA def:0.36 GeV] [NUX tune:0.44 GeV]
  - pt2tail: [PYTHIA def:0.01    ] [NUX tune:0.01    ]
  - cutoff:  [PYTHIA def:0.80 GeV] [NUX tune:0.20 GeV]
  -->
  <param type="double" name="PYTHIA-SSBarSuppression">      0.21 </param>
  <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>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  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>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  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.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>
  <param type="double" name="DIS-HMultWgt-vp-NC-m3">  1.000  </param>
  <param type="double" name="DIS-HMultWgt-vn-CC-m2">  0.300  </param>
  <param type="double" name="DIS-HMultWgt-vn-CC-m3">  1.000  </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-vbp-CC-m2"> 0.300  </param>
  <param type="double" name="DIS-HMultWgt-vbp-CC-m3"> 1.000  </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-vbn-CC-m2"> 0.100  </param>
  <param type="double" name="DIS-HMultWgt-vbn-CC-m3"> 1.000  </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>

 <!-- 	
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Parameters relevant to formation zone simulation
  - ct0 is the formation time times the speed of light (given in fm)
  - KPt2 is the parameter multiplying pT2 in formation zone calc.
  -->
  <param type="double" name="FZONE-ct0">    0.342   </param>
  <param type="double" name="FZONE-KPt2">   0.0     </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  RES/DIS joining parameter Wcut 
  DRJoinScheme:
   - ...
   - ...
  -->
  <param type="bool"   name="UseDRJoinScheme"> true </param>
  <param type="double" name="Wcut">            1.7  </param>


 <!-- 
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  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>

 <!-- 
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  NUCL-VtxGenerationMethod defines the intranuclear vertex generation method:
   - 0 : generated uniformly in volume within the nucleus 
   - 1 : generated according to a realistic nuclear density profile
  -->
  <param type="int" name="NUCL-VtxGenerationMethod"> 1 </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Intranuclear rescattering 
  Use the HadronTransportEnable option to toggle intranuclear rescattering on/off. 
  Possible model options include INTRANUKE, FLUKA and GIBUU
  -->
  <param type="bool" name="HadronTransp-Enable"> true </param>
  <param type="alg"  name="HadronTransp-Model">  genie::Intranuke/Default  </param>

 <!-- 	
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  INTRANUKE-specific parameters:
  - Mode options are: hA, hN
  - NucRemovalE is the binding E to subtract from cascade nucleons (in GeV)
  -->
  <param type="string" name="INUKE-Mode">            hA    </param>
  <param type="double" name="INUKE-NucRemovalE">     0.008 </param>
  <param type="double" name="INUKE-HadStep">         0.05  </param>
  <param type="double" name="INUKE-DelRPion">        0.5   </param>
  <param type="double" name="INUKE-DelRNucleon">     1.0   </param>

 <!-- 	
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Particle decayer parameters
  -->
  <param type="bool" name="ForceCharmedHadronDecay"> false </param>
  <param type="bool" name="InhibitPi0Decay">         true  </param>
  

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  Specify which cross section model is to be used by each GENIE event generation thread.
  The parameter name is build as: "XSecModel@[name of thread]"
  -->
  <param type="alg" name="XSecModel@genie::EventGenerator/QEL-CC">       genie::QELPXSec/CC-Default            </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/QEL-NC">       genie::NuNucElasticPXSec/Default      </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/RES-CC">       genie::ReinSeghalRESPXSec/Default     </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/RES-NC">       genie::ReinSeghalRESPXSec/Default     </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/DIS-CC">       genie::DISPartonModelPXSec/CC-Default </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/DIS-NC">       genie::DISPartonModelPXSec/NC-Default </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/COHPI-CC">     genie::ReinSeghalCOHPiPXSec/Default   </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/COHPI-NC">     genie::ReinSeghalCOHPiPXSec/Default   </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/IMD">          genie::BardinIMDRadCorPXSec/Default   </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/NUE-EL">       genie::NuElectronPXSec/Default        </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/QEL-CC-CHARM"> genie::KovalenkoQELCharmPXSec/Default </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/DIS-CC-CHARM"> genie::AivazisCharmPXSecLO/CC-Default </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/AM-NUGAMMA">   genie::H3AMNuGammaPXSec/Default       </param>
  <param type="alg" name="XSecModel@genie::EventGenerator/MEC">          genie::MECPXSec/Default               </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  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
  -->
  <param type="alg" name="ElasticFormFactorsModel"> genie::BBA05ELFormFactorsModel/Default </param>

 <!-- 
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  GENIE validity range
  -->
  <param type="double" name="GVLD-Emin">    0.010  </param>
  <param type="double" name="GVLD-Emax">  300.000  </param>

 </param_set>

</global_param_list>

