#include <SPPXSecWithCache.h>

Inheritance diagram for genie::utils::gsl::d3XSecMK_dWQ2CosTheta_E:

Collaboration diagram for genie::utils::gsl::d3XSecMK_dWQ2CosTheta_E:

Public Member Functions
	d3XSecMK_dWQ2CosTheta_E (const XSecAlgorithmI m, const Interaction i, double wcut)

	~d3XSecMK_dWQ2CosTheta_E ()

unsigned int	NDim (void) const

double	DoEval (const double *xin) const

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

Private Attributes
const XSecAlgorithmI *	fModel

Interaction *	fInteraction

Range1D_t	Wl

bool	isZero

KPhaseSpace *	kps

double	fWcut

Detailed Description

Definition at line 78 of file SPPXSecWithCache.h.

Constructor & Destructor Documentation

genie::utils::gsl::d3XSecMK_dWQ2CosTheta_E::d3XSecMK_dWQ2CosTheta_E	(	const XSecAlgorithmI *	m,
		const Interaction *	i,
		double	wcut
	)

Definition at line 218 of file SPPXSecWithCache.cxx.

References fInteraction, fWcut, isZero, genie::kISkipKinematicChk, genie::kISkipProcessChk, kps, genie::kRfHitNucRest, genie::Range1D_t::max, genie::Range1D_t::min, genie::Interaction::PhaseSpacePtr(), genie::InitialState::ProbeE(), Wl, and genie::KPhaseSpace::WLim_SPP_iso().

                                                                                                          :
   ROOT::Math::IBaseFunctionMultiDim(),
   fModel(m),
   fWcut(wcut)
 {
 
   isZero = false;
   fInteraction = const_cast<Interaction*>(interaction);
   fInteraction->SetBit(kISkipProcessChk);
   fInteraction->SetBit(kISkipKinematicChk);
   
   kps = fInteraction->PhaseSpacePtr();
   
   // Get kinematical parameters
   const InitialState & init_state = interaction -> InitState();
   double Enu = init_state.ProbeE(kRfHitNucRest);
 
 
   if (Enu < kps->Threshold_SPP_iso())
   {
     isZero = true;
     return;
   }
   
   Wl  = kps->WLim_SPP_iso();
   if (fWcut >= Wl.min)
     Wl.max = TMath::Min(fWcut,Wl.max);
   
   
 }

genie::utils::gsl::d3XSecMK_dWQ2CosTheta_E::~d3XSecMK_dWQ2CosTheta_E ( )

Definition at line 249 of file SPPXSecWithCache.cxx.

 {
 
 }

Member Function Documentation

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

Definition at line 279 of file SPPXSecWithCache.cxx.

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

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

Definition at line 257 of file SPPXSecWithCache.cxx.

References genie::units::cm2, genie::kKVctp, genie::kPSWQ2ctpfE, genie::Range1D_t::max, and genie::Range1D_t::min.

 {
   // outputs:
   //   differential cross section [1/GeV^3] for Resonance single pion production production
   //
 
   if (isZero) return 0.;
   
   double W2  = Wl.min*Wl.min + (Wl.max*Wl.max - Wl.min*Wl.min)*xin[0];
   fInteraction->KinePtr()->SetW(TMath::Sqrt(W2));
    
   Range1D_t Q2l = kps->Q2Lim_W_SPP_iso(); 
    
   double sqrt_Q2 = TMath::Sqrt(Q2l.min) + ( TMath::Sqrt(Q2l.max) - TMath::Sqrt(Q2l.min) )*xin[1];
   fInteraction->KinePtr()->SetQ2(sqrt_Q2*sqrt_Q2);
   
   fInteraction->KinePtr()->SetKV(kKVctp, -1. + 2.*xin[2]); //CosTheta
   
   double xsec = fModel->XSec(fInteraction, kPSWQ2ctpfE)*sqrt_Q2*(Wl.max*Wl.max - Wl.min*Wl.min)*(TMath::Sqrt(Q2l.max) - TMath::Sqrt(Q2l.min))*2/TMath::Sqrt(W2);
   return xsec/(1E-38 * units::cm2);
 }