#include <GSLXSecFunc.h>

Inheritance diagram for genie::utils::gsl::dXSec_dElep_AR:

Collaboration diagram for genie::utils::gsl::dXSec_dElep_AR:

Public Member Functions
	dXSec_dElep_AR (const XSecAlgorithmI m, const Interaction i, string gsl_nd_integrator_type, double gsl_relative_tolerance, unsigned int max_n_calls)

	dXSec_dElep_AR ()

	~dXSec_dElep_AR ()

dXSec_dElep_AR *	Clone (void) const

double	DoEval (double xin) const

Private Attributes
const XSecAlgorithmI *	fModel

const Interaction *	fInteraction

const genie::utils::gsl::d3Xsec_dOmegaldThetapi *	func

ROOT::Math::IntegratorMultiDim	integrator

double	kine_min [3]

double	kine_max [3]

string	fGSLIntegratorType

double	fGSLRelTol

unsigned int	fGSLMaxCalls

Detailed Description

.....................................................................................

genie::utils::gsl::dXSec_dElep_AR A 1-D cross section function: dxsec/dElep Used for Alvarez-Ruso coherent.

Definition at line 422 of file GSLXSecFunc.h.

Constructor & Destructor Documentation

genie::utils::gsl::dXSec_dElep_AR::dXSec_dElep_AR	(	const XSecAlgorithmI *	m,
		const Interaction *	i,
		string	gsl_nd_integrator_type,
		double	gsl_relative_tolerance,
		unsigned int	max_n_calls
	)

Definition at line 1001 of file GSLXSecFunc.cxx.

References fInteraction, fModel, func, integrator, genie::controls::kASmallNum, kine_max, kine_min, and genie::constants::kPi.

                               :
 ROOT::Math::IBaseFunctionOneDim(),
 fModel(m),
 fInteraction(i),
 integrator(utils::gsl::IntegrationNDimTypeFromString(gsl_nd_integrator_type),1, gsl_relative_tolerance, max_n_calls),
 fGSLIntegratorType(gsl_nd_integrator_type),
 fGSLRelTol(gsl_relative_tolerance),
 fGSLMaxCalls(max_n_calls)
 {
   func = new utils::gsl::d3Xsec_dOmegaldThetapi(fModel, fInteraction);
 
   integrator.SetRelTolerance(gsl_relative_tolerance);
   integrator.SetFunction(*func);
 
   kine_min[0] = kine_min[1] = kine_min[2] = controls::kASmallNum;
 
   kine_max[0] = kine_max[2] = constants::kPi-controls::kASmallNum;
   kine_max[1] = 2 * constants::kPi-controls::kASmallNum;
 }

genie::utils::gsl::dXSec_dElep_AR::dXSec_dElep_AR ( )

inline

Definition at line 428 of file GSLXSecFunc.h.

428 {};

genie::utils::gsl::dXSec_dElep_AR::~dXSec_dElep_AR ( )

Definition at line 1023 of file GSLXSecFunc.cxx.

References func().

 {
   delete func;
 }

Member Function Documentation

genie::utils::gsl::dXSec_dElep_AR * genie::utils::gsl::dXSec_dElep_AR::Clone ( void ) const

Definition at line 1036 of file GSLXSecFunc.cxx.

 {
   return
     new genie::utils::gsl::dXSec_dElep_AR(fModel,fInteraction, fGSLIntegratorType, fGSLRelTol, fGSLMaxCalls);
 }

double genie::utils::gsl::dXSec_dElep_AR::DoEval ( double xin ) const

Definition at line 1027 of file GSLXSecFunc.cxx.

References func(), LOG, and pINFO.

 {
   double Elep = xin;
   func->SetE_lep(Elep);
   double xsec = integrator.Integral(&kine_min[0], &kine_max[0]) ;
   LOG("GSLXSecFunc",pINFO) << "dXSec_dElep_AR >> "<<func->NDim()<<"d integral done. (Elep = " <<Elep<< " , dxsec/dElep = "<<xsec << ")";
   return xsec;
 }