genie::utils::gsl::d3Xsec_dTldTkdCosThetal Class Reference

#include <AlamSimoAtharVacasSKXSec.h>

Inheritance diagram for genie::utils::gsl::d3Xsec_dTldTkdCosThetal:

Collaboration diagram for genie::utils::gsl::d3Xsec_dTldTkdCosThetal:

Public Member Functions
	d3Xsec_dTldTkdCosThetal (const XSecAlgorithmI *m, const Interaction *i)
	~d3Xsec_dTldTkdCosThetal ()
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

Detailed Description

Definition at line 64 of file AlamSimoAtharVacasSKXSec.h.

Constructor & Destructor Documentation

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

Definition at line 171 of file AlamSimoAtharVacasSKXSec.cxx.

                                                      :
ROOT::Math::IBaseFunctionMultiDim(),
fModel(m),
fInteraction(i)
{

}

genie::utils::gsl::d3Xsec_dTldTkdCosThetal::~d3Xsec_dTldTkdCosThetal

(

)

Definition at line 180 of file AlamSimoAtharVacasSKXSec.cxx.

{

}

Member Function Documentation

ROOT::Math::IBaseFunctionMultiDim * genie::utils::gsl::d3Xsec_dTldTkdCosThetal::Clone

(

void

)

const

Definition at line 237 of file AlamSimoAtharVacasSKXSec.cxx.

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

double genie::utils::gsl::d3Xsec_dTldTkdCosThetal::DoEval

(

const double *

xin

)

const

Definition at line 192 of file AlamSimoAtharVacasSKXSec.cxx.

References genie::units::cm2, genie::utils::mec::J(), genie::kKVctl, genie::kKVphikq, genie::kKVTk, genie::kKVTl, genie::constants::kPi, genie::kRfLab, LOG, and pDEBUG.

{
// inputs:
//    Tl [GeV]
//    Tk [GeV]
//    cosine theta l
//    * calculate phi_kq based on neutrino energy -- this is for the integral only
// outputs:
//   differential cross section [10^-38 cm^2]
//

  double Enu = fInteraction->InitState().ProbeE(kRfLab);

  double phikq = 0.5*constants::kPi;
  if( Enu > 3.0 ) phikq = 0.55*constants::kPi;
  else if( Enu > 1.0 ) phikq = constants::kPi*(0.5 + 0.025*(Enu-1.0));

  Kinematics * kinematics = fInteraction->KinePtr();

  double T_l         = xin[0];
  double T_k         = xin[1];
  //double cos_theta_l = xin[2];
  double log_oneminuscostheta = xin[2];
  double cos_theta_l = 1.0 - TMath::Exp(log_oneminuscostheta);
  double J = 1.0 - cos_theta_l; // Jacobian for transformation

  kinematics->SetKV(kKVTl, T_l);
  kinematics->SetKV(kKVTk, T_k);
  kinematics->SetKV(kKVctl, cos_theta_l);
  kinematics->SetKV(kKVphikq, phikq);

  double xsec = fModel->XSec(fInteraction);
  LOG( "GXSecFunc", pDEBUG )
     << "t_l = " << T_l << " t_k = " << T_k
     << " costhetal = " << cos_theta_l << " phikq = " << phikq
     << " enu = " << Enu << " Xsec = " << xsec;

  // integrate out phi_kq by multiplying by 2pi

  xsec *= 2.0 * genie::constants::kPi * J;

  return xsec/(1E-38 * units::cm2);
}

unsigned int genie::utils::gsl::d3Xsec_dTldTkdCosThetal::NDim

(

void

)

const

Definition at line 185 of file AlamSimoAtharVacasSKXSec.cxx.

{
  // phi_kq is important for kinematics generation
  // But dependence is weak so we will not use it in the integration
  return 3;
}

Member Data Documentation

const Interaction* genie::utils::gsl::d3Xsec_dTldTkdCosThetal::fInteraction

private

Definition at line 75 of file AlamSimoAtharVacasSKXSec.h.

const XSecAlgorithmI* genie::utils::gsl::d3Xsec_dTldTkdCosThetal::fModel

private

Definition at line 74 of file AlamSimoAtharVacasSKXSec.h.

---

The documentation for this class was generated from the following files:

• AlamSimoAtharVacasSKXSec.h
• AlamSimoAtharVacasSKXSec.cxx