Public Member Functions | Private Attributes | List of all members
genie::utils::gsl::wrap::NuclDensityMomentIntegrand Class Reference

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:

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 <constantinos.andreopoulos cern.ch> University of Liverpool & STFC Rutherford Appleton Laboratory

July 12, 2019

Copyright (c) 2003-2020, 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.

21  :
22 ROOT::Math::IBaseFunctionOneDim()
23 {
24  fA = A;
25  fK = k;
26 
27  if(fK < 0) {
28  LOG("Nuclear", pDEBUG)
29  << "Sure you want to calculate an inverse nuclear density moment ("
30  << "E[r^{" << fK << "}]) ?";
31  }
32  if(fA <= 1) {
33  LOG("Nuclear", pWARN)
34  << "The atomic mass number A should be >1 (input value was: "
35  << fA << ")";
36  }
37 }
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96
pWARN
#define pWARN
Definition: Messenger.h:60
A
#define A
Definition: memgrp.cpp:38
pDEBUG
#define pDEBUG
Definition: Messenger.h:63
utils::gsl::wrap::NuclDensityMomentIntegrand::~NuclDensityMomentIntegrand ( )

Definition at line 39 of file NuclDensityMomentIntegrand.cxx.

40 {
41 
42 }

Member Function Documentation

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

Definition at line 61 of file NuclDensityMomentIntegrand.cxx.

62 {
63  return new utils::gsl::wrap::NuclDensityMomentIntegrand(fA, fK);
64 }
genie::utils::gsl::wrap::NuclDensityMomentIntegrand
Integrand for the calculation of the k^th nuclear density moment: {0}^{} {A}(r) r^k d^{3}r where {A}(...
Definition: NuclDensityMomentIntegrand.h:30
double utils::gsl::wrap::NuclDensityMomentIntegrand::DoEval ( double  xin) const

Definition at line 49 of file NuclDensityMomentIntegrand.cxx.

50 {
51  if(fA <= 1) return 0;
52 
53  double rho = utils::nuclear::Density(r,fA);
54  double rkth = TMath::Power(r,fK);
55 
56  double integrand = rho * rkth;
57  return integrand;
58 }
genie::utils::nuclear::Density
double Density(double r, int A, double ring=0.)
Definition: NuclearUtils.cxx:396
plot_model.r
r
Definition: plot_model.py:57
unsigned int utils::gsl::wrap::NuclDensityMomentIntegrand::NDim ( void  ) const

Definition at line 44 of file NuclDensityMomentIntegrand.cxx.

45 {
46  return 1;
47 }

Member Data Documentation

int genie::utils::gsl::wrap::NuclDensityMomentIntegrand::fA
private

Definition at line 40 of file NuclDensityMomentIntegrand.h.

int genie::utils::gsl::wrap::NuclDensityMomentIntegrand::fK
private

Definition at line 41 of file NuclDensityMomentIntegrand.h.

The documentation for this class was generated from the following files:
Generated by   doxygen 1.8.11