Integrand for the calculation of the k^th nuclear density moment: {0}^{} {A}(r) r^k d^{3}r where {A}(r) is the nuclear density for a nucleus with atomic mass number A. More...

#include <NuclDensityMomentIntegrand.h>

Inheritance diagram for genie::utils::gsl::wrap::NuclDensityMomentIntegrand:

[legend]

Collaboration diagram for genie::utils::gsl::wrap::NuclDensityMomentIntegrand:

[legend]

Public Member Functions
	NuclDensityMomentIntegrand (int A, int k)

	~NuclDensityMomentIntegrand ()

unsigned int	NDim (void) const

double	DoEval (double xin) const

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

Private Attributes
int	fA

int	fK

Detailed Description

Integrand for the calculation of the k^th nuclear density moment: {0}^{} {A}(r) r^k d^{3}r where {A}(r) is the nuclear density for a nucleus with atomic mass number A.

Author: Costas Andreopoulos <c.andreopoulos cern.ch> University of Liverpool

Created:: July 12, 2019

License:: Copyright (c) 2003-2024, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 30 of file NuclDensityMomentIntegrand.h.

Constructor & Destructor Documentation

utils::gsl::wrap::NuclDensityMomentIntegrand::NuclDensityMomentIntegrand	(	int	A,
		int	k
	)

Definition at line 20 of file NuclDensityMomentIntegrand.cxx.

References genie::units::A, fA, fK, LOG, pDEBUG, and pWARN.

               :
 ROOT::Math::IBaseFunctionOneDim()
 {
   fA  = A;
   fK  = k;
 
   if(fK < 0) {
      LOG("Nuclear", pDEBUG)
         << "Sure you want to calculate an inverse nuclear density moment ("
         << "E[r^{" << fK << "}]) ?";
   }
   if(fA <= 1) {
      LOG("Nuclear", pWARN)
         << "The atomic mass number A should be >1 (input value was: "
         << fA << ")";
   }
 }

utils::gsl::wrap::NuclDensityMomentIntegrand::~NuclDensityMomentIntegrand ( )

Definition at line 39 of file NuclDensityMomentIntegrand.cxx.

 {
 
 }

Member Function Documentation

ROOT::Math::IBaseFunctionOneDim * utils::gsl::wrap::NuclDensityMomentIntegrand::Clone ( void ) const

Definition at line 61 of file NuclDensityMomentIntegrand.cxx.

 {
   return new utils::gsl::wrap::NuclDensityMomentIntegrand(fA, fK);
 }

double utils::gsl::wrap::NuclDensityMomentIntegrand::DoEval ( double xin ) const

Definition at line 49 of file NuclDensityMomentIntegrand.cxx.

References genie::utils::nuclear::Density().

 {
   if(fA <= 1) return 0;
 
   double rho  = utils::nuclear::Density(r,fA);
   double rkth = TMath::Power(r,fK);
 
   double integrand = rho * rkth;
   return integrand;
 }