diff --git a/config/KNOTunedQPMDISPXSec.xml b/config/KNOTunedQPMDISPXSec.xml index 77cc15ad96..0098c0a2c6 100644 --- a/config/KNOTunedQPMDISPXSec.xml +++ b/config/KNOTunedQPMDISPXSec.xml @@ -12,6 +12,7 @@ DISModel alg no It must be QPMDISPXSec UseCache bool Yes Cache suppression factors? true Wcut double no CommonParam[NonResBackground] +Wcutmin double Yes 0 --> @@ -20,9 +21,10 @@ Wcut double no Co NonResBackground - genie::QPMDISPXSec/Default + genie::QPMDISPXSec/Default genie::AGKYLowW2019/Default 1. + 0. diff --git a/src/Physics/DeepInelastic/XSection/KNOTunedQPMDISPXSec.cxx b/src/Physics/DeepInelastic/XSection/KNOTunedQPMDISPXSec.cxx index f2f45a220a..f98594a2d0 100644 --- a/src/Physics/DeepInelastic/XSection/KNOTunedQPMDISPXSec.cxx +++ b/src/Physics/DeepInelastic/XSection/KNOTunedQPMDISPXSec.cxx @@ -103,7 +103,8 @@ double KNOTunedQPMDISPXSec::DISRESJoinSuppressionFactor( // double R=0, Ro=0; - const double Wmin = kNeutronMass + kPionMass + 1E-3; + double Wmin = kNeutronMass + kPionMass + 1E-3; + if( Wmin < fWcutmin ) Wmin = fWcutmin ; const InitialState & ist = in->InitState(); const ProcessInfo & pi = in->ProcInfo(); @@ -235,6 +236,7 @@ void KNOTunedQPMDISPXSec::LoadConfig(void) assert(fHadronizationModel); GetParam( "Wcut", fWcut ) ; + GetParam( "Wcutmin", fWcutmin ) ; GetParam( "NRB-EM-XSecScale", fNRBEMScale ); if ( fWcut <= 0. ) { diff --git a/src/Physics/DeepInelastic/XSection/KNOTunedQPMDISPXSec.h b/src/Physics/DeepInelastic/XSection/KNOTunedQPMDISPXSec.h index ad1cc99e31..8d51340c3d 100644 --- a/src/Physics/DeepInelastic/XSection/KNOTunedQPMDISPXSec.h +++ b/src/Physics/DeepInelastic/XSection/KNOTunedQPMDISPXSec.h @@ -56,6 +56,7 @@ class KNOTunedQPMDISPXSec : public XSecAlgorithmI { bool fUseCache; ///< cache reduction factors used in joining scheme double fWcut; ///< apply DIS/RES joining scheme < Wcut + double fWcutmin; ///< Limit to the Non-Resonant background coverage. Below < Wcutmin, only RES is considered double fNRBEMScale; ///< apply NRB EM Scale factor };