#include <ReinSehgalRESXSecWithCacheFast.h>

Inheritance diagram for genie::utils::gsl::d2XSecRESFast_dWQ2_E:

Collaboration diagram for genie::utils::gsl::d2XSecRESFast_dWQ2_E:

Public Member Functions
	d2XSecRESFast_dWQ2_E (const XSecAlgorithmI m, const Interaction i)

	~d2XSecRESFast_dWQ2_E ()

unsigned int	NDim (void) const

double	DoEval (const double *xin) const

ROOT::Math::IBaseFunctionMultiDim *	Clone (void) const

Private Attributes
const XSecAlgorithmI *	fModel

const Interaction *	fInteraction

double	fWmin

double	fWmax

bool	isfWcutLessfWmin

KPhaseSpace *	kps

Detailed Description

Definition at line 78 of file ReinSehgalRESXSecWithCacheFast.h.

Constructor & Destructor Documentation

genie::utils::gsl::d2XSecRESFast_dWQ2_E::d2XSecRESFast_dWQ2_E	(	const XSecAlgorithmI *	m,
		const Interaction *	i
	)

Definition at line 227 of file ReinSehgalRESXSecWithCacheFast.cxx.

References genie::Interaction::ExclTag(), fInteraction, fModel, fWmax, fWmin, genie::Registry::GetBool(), genie::Registry::GetBoolDef(), genie::Registry::GetDouble(), genie::Registry::GetDoubleDef(), isfWcutLessfWmin, kps, genie::utils::res::Mass(), genie::Range1D_t::max, genie::Range1D_t::min, genie::Interaction::PhaseSpacePtr(), genie::XclsTag::Resonance(), genie::utils::res::ResonanceIndex(), genie::utils::res::Width(), and genie::KPhaseSpace::WLim().

                                                       :
 ROOT::Math::IBaseFunctionMultiDim(),
 fModel(m),
 fInteraction(i)
 {
         kps = fInteraction->PhaseSpacePtr();
         Range1D_t Wl  = kps->WLim();
         fWmin = Wl.min;
         fWmax = Wl.max;
         Registry fConfig = (const_cast<XSecAlgorithmI *>(fModel))->GetConfig();
         bool fUsingDisResJoin = fConfig.GetBool("UseDRJoinScheme");
         double fWcut = 999999;
         if(fUsingDisResJoin)
         {
                 fWcut = fConfig.GetDouble("Wcut");
         }
         fWmax=TMath::Min(fWcut, fWmax);
         if (fWcut<fWmin)
                 isfWcutLessfWmin=true;
         else
                 isfWcutLessfWmin=false;
         bool fNormBW = fConfig.GetBoolDef("BreitWignerNorm", true);
         if (fNormBW)
         {
                 double fN2ResMaxNWidths = fConfig.GetDoubleDef("MaxNWidthForN2Res", 2.0);
                 double fN0ResMaxNWidths = fConfig.GetDoubleDef("MaxNWidthForN0Res", 6.0);
                 double fGNResMaxNWidths = fConfig.GetDoubleDef("MaxNWidthForGNRes", 4.0);
                 Resonance_t resonance = fInteraction->ExclTag().Resonance();
                 int    IR  = utils::res::ResonanceIndex    (resonance);
                 double MR  = utils::res::Mass              (resonance);
                 double WR  = utils::res::Width             (resonance);
                 if (IR==0)
                         fWcut = MR + fN0ResMaxNWidths * WR;
                 else if (IR==2)
                         fWcut = MR + fN2ResMaxNWidths * WR;
                 else
                         fWcut = MR + fGNResMaxNWidths * WR;
                 fWmax=TMath::Min(fWcut, fWmax);
                 if (fWcut<fWmin)
                         isfWcutLessfWmin=true;
         }
 
 }

genie::utils::gsl::d2XSecRESFast_dWQ2_E::~d2XSecRESFast_dWQ2_E ( )

Definition at line 271 of file ReinSehgalRESXSecWithCacheFast.cxx.

 {
 
 }

Member Function Documentation

ROOT::Math::IBaseFunctionMultiDim * genie::utils::gsl::d2XSecRESFast_dWQ2_E::Clone ( void ) const

Definition at line 301 of file ReinSehgalRESXSecWithCacheFast.cxx.

 {
   return
     new genie::utils::gsl::d2XSecRESFast_dWQ2_E(fModel,fInteraction);
 }

double genie::utils::gsl::d2XSecRESFast_dWQ2_E::DoEval ( const double * xin ) const

Definition at line 279 of file ReinSehgalRESXSecWithCacheFast.cxx.

References genie::units::cm2, genie::kPSWQ2fE, genie::Range1D_t::max, genie::Range1D_t::min, genie::utils::kinematics::Q2(), and genie::utils::kinematics::W().

 {
 // inputs:
 //    W
 //    Q2
 // outputs:
 //   differential cross section [10^-38 cm^2/GeV^3] for Resonance production
 //
 
   if (isfWcutLessfWmin)
         return 0;
   double W  = fWmin+(fWmax-fWmin)*xin[0];
   fInteraction->KinePtr()->SetW(W);
   Range1D_t Q2l = kps->Q2Lim_W();
   if (Q2l.min<0 || Q2l.max<0)
         return 0.0;
   double Q2 = Q2l.min+(Q2l.max-Q2l.min)*xin[1];
   fInteraction->KinePtr()->SetQ2(Q2);
   double xsec = fModel->XSec(fInteraction, kPSWQ2fE)*(fWmax-fWmin)*(Q2l.max-Q2l.min);
   return xsec/(1E-38 * units::cm2);
 }