Heed::PhenoPhotoAbsCS Class Reference

Simple phenomenological CS for any shell (analytic formula). More...

#include <PhotoAbsCS.h>

Inheritance diagram for Heed::PhenoPhotoAbsCS:

Public Member Functions
	PhenoPhotoAbsCS ()
	Default constructor.
	PhenoPhotoAbsCS (const std::string &fname, int fZ, double fthreshold, double fpower=2.75)
virtual	~PhenoPhotoAbsCS ()
	Destructor.
virtual double	get_CS (double energy) const
	Retrieve cross-section [Mb] at a given energy [MeV].
virtual double	get_integral_CS (double energy1, double energy2) const
	Retrieve integral cross-section [Mb * MeV] in a given interval [MeV].
virtual void	scale (double fact)
	Multiply by some factor. Can be useful for debugging and other purposes.
virtual void	print (std::ostream &file, int l) const
virtual PhenoPhotoAbsCS *	copy () const
Public Member Functions inherited from Heed::PhotoAbsCS
	PhotoAbsCS ()
	Default constructor.
	PhotoAbsCS (const std::string &fname, int fZ, double fthreshold)
	Constructor.
virtual	~PhotoAbsCS ()
	Destructor.
const std::string &	get_name () const
	Name of this shell or atom.
int	get_number () const
	Number of this shell.
int	get_Z () const
double	get_threshold () const
	Return the threshold energy.
virtual double	get_CS (double energy) const =0
	Retrieve cross-section [Mb] at a given energy [MeV].
virtual double	get_integral_CS (double energy1, double energy2) const =0
	Retrieve integral cross-section [Mb * MeV] in a given interval [MeV].
virtual void	scale (double fact)=0
	Multiply by some factor. Can be useful for debugging and other purposes.
virtual void	print (std::ostream &file, int l) const
virtual PhotoAbsCS *	copy () const =0

Additional Inherited Members
Protected Attributes inherited from Heed::PhotoAbsCS
std::string	name
int	number
int	Z
double	threshold

Detailed Description

Simple phenomenological CS for any shell (analytic formula).

Definition at line 211 of file PhotoAbsCS.h.

Constructor & Destructor Documentation

PhenoPhotoAbsCS() [1/2]

Heed::PhenoPhotoAbsCS::PhenoPhotoAbsCS

(

)

Default constructor.

Definition at line 425 of file PhotoAbsCS.cpp.

: PhotoAbsCS("none", 0, 0.0), power(0.0) {}

Referenced by copy().

PhenoPhotoAbsCS() [2/2]

Heed::PhenoPhotoAbsCS::PhenoPhotoAbsCS	(	const std::string &	fname,
		int	fZ,
		double	fthreshold,
		double	fpower = 2.75
	)

Constructor

Parameters

fname	name of the shell or atom
fZ	number of electrons
fthreshold	threshold level
fpower	positive number in

Definition at line 427 of file PhotoAbsCS.cpp.

    : PhotoAbsCS(fname, fZ, fthreshold), power(fpower) {
  mfunname("PhenoPhotoAbsCS::PhenoPhotoAbsCS(...)");
  check_econd11a(power, <= 2, " value not allowed, integral would be infinite",
                 mcerr);
  const double a = power - 1.;
  factor = pow(threshold, a) * Thomas_sum_rule_const_Mb * Z * a;
}

~PhenoPhotoAbsCS()

virtual Heed::PhenoPhotoAbsCS::~PhenoPhotoAbsCS ( )

inlinevirtual

Destructor.

Definition at line 224 of file PhotoAbsCS.h.

{}

Member Function Documentation

copy()

virtual PhenoPhotoAbsCS * Heed::PhenoPhotoAbsCS::copy ( ) const

inlinevirtual

Implements Heed::PhotoAbsCS.

Definition at line 229 of file PhotoAbsCS.h.

{ return new PhenoPhotoAbsCS(*this); }

get_CS()

double Heed::PhenoPhotoAbsCS::get_CS ( double  energy ) const

virtual

Retrieve cross-section [Mb] at a given energy [MeV].

Implements Heed::PhotoAbsCS.

Definition at line 437 of file PhotoAbsCS.cpp.

                                                  {
  if (energy < threshold || energy == DBL_MAX) return 0.0;
  return factor * pow(energy, -power);
}

get_integral_CS()

double Heed::PhenoPhotoAbsCS::get_integral_CS ( double  energy1, double  energy2  ) const

virtual

Retrieve integral cross-section [Mb * MeV] in a given interval [MeV].

Implements Heed::PhotoAbsCS.

Definition at line 442 of file PhotoAbsCS.cpp.

                                                                            {
  if (energy2 < threshold) return 0.0;
  if (energy1 < threshold) energy1 = threshold;
  const double a = power - 1.;
  double s;
  if (energy2 == DBL_MAX) {
    s = factor / a * (1. / pow(energy1, a));
  } else {
    s = factor / a * (1. / pow(energy1, a) - 1. / pow(energy2, a));
  }
  return s;
}

print()

void Heed::PhenoPhotoAbsCS::print ( std::ostream &  file, int  l  ) const

virtual

Reimplemented from Heed::PhotoAbsCS.

Definition at line 460 of file PhotoAbsCS.cpp.

                                                         {
  if (l <= 0) return;
  Ifile << "PhenoPhotoAbsCS: name=" << name << " Z = " << Z << std::endl;
  Ifile << " threshold = " << threshold << " power=" << power
        << " factor=" << factor << std::endl;
}

scale()

void Heed::PhenoPhotoAbsCS::scale ( double  fact )

virtual

Multiply by some factor. Can be useful for debugging and other purposes.

Implements Heed::PhotoAbsCS.

Definition at line 455 of file PhotoAbsCS.cpp.

                                       {
  mfunnamep("void PhenoPhotoAbsCS::scale(double fact)");
  factor *= fact;
}

---

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

• PhotoAbsCS.h
• PhotoAbsCS.cpp