Heed::AveragePhotoAbsCS Class Reference

Smoothed/smeared photoabsorption cross-section. More...

#include <PhotoAbsCS.h>

Inheritance diagram for Heed::AveragePhotoAbsCS:

Public Member Functions
	AveragePhotoAbsCS ()
	Default constructor.
	AveragePhotoAbsCS (PhotoAbsCS *apacs, double fwidth, double fstep, long fmax_q_step)
virtual	~AveragePhotoAbsCS ()
	Destructor.
double	get_CS (double energy) const override
	Retrieve cross-section [Mb] at a given energy [MeV].
double	get_integral_CS (double energy1, double energy2) const override
	Retrieve integral cross-section [Mb * MeV] in a given interval [MeV].
void	scale (double fact) override
	Multiply by some factor. Can be useful for debugging and other purposes.
void	print (std::ostream &file, int l) const override
AveragePhotoAbsCS *	copy () const override
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 = 0
int	Z
double	threshold

Detailed Description

Smoothed/smeared photoabsorption cross-section.

Definition at line 83 of file PhotoAbsCS.h.

Constructor & Destructor Documentation

AveragePhotoAbsCS() [1/2]

Heed::AveragePhotoAbsCS::AveragePhotoAbsCS ( )

inline

Default constructor.

Definition at line 94 of file PhotoAbsCS.h.

: PhotoAbsCS() {}

Referenced by copy().

AveragePhotoAbsCS() [2/2]

Heed::AveragePhotoAbsCS::AveragePhotoAbsCS	(	PhotoAbsCS *	apacs,
		double	fwidth,
		double	fstep,
		long	fmax_q_step
	)

Constructor.

Parameters

apacs	photoabsorption cross-section
fwidth	width [MeV] for smoothing
fstep	step size [MeV] for numerical integration
fmax_q_step	max number of integration steps

Definition at line 99 of file PhotoAbsCS.cpp.

    : width(fwidth),
      max_q_step(fmax_q_step),
      step(fstep) {
  mfunname("AveragePhotoAbsCS::AveragePhotoAbsCS(...)");
  check_econd11(apacs, == nullptr, mcerr);
  real_pacs.reset(apacs);
  // Check the parameters (step = 0.5 * width is bad but OK).
  if (fwidth > 0.0) check_econd11(fstep, >= 0.6 * fwidth, mcerr);
  // Copy the parameters of the "real" cross-section.
  name = real_pacs->get_name();
  number = real_pacs->get_number();
  Z = real_pacs->get_Z();
  threshold = real_pacs->get_threshold();
}

~AveragePhotoAbsCS()

virtual Heed::AveragePhotoAbsCS::~AveragePhotoAbsCS ( )

inlinevirtual

Destructor.

Definition at line 108 of file PhotoAbsCS.h.

{}

Member Function Documentation

copy()

AveragePhotoAbsCS * Heed::AveragePhotoAbsCS::copy ( ) const

inlineoverridevirtual

Implements Heed::PhotoAbsCS.

Definition at line 115 of file PhotoAbsCS.h.

                                           {
    return new AveragePhotoAbsCS(*this);
  }

get_CS()

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

overridevirtual

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

Implements Heed::PhotoAbsCS.

Definition at line 117 of file PhotoAbsCS.cpp.

                                                    {
  mfunname("double AveragePhotoAbsCS::get_CS(double energy) const");
  // In case of zero width, return the unmodified "real" cross-section.
  if (width == 0.0) return real_pacs->get_CS(energy);
  const double w2 = width * 0.5;
  const double e1 = std::max(energy - w2, 0.);
  return real_pacs->get_integral_CS(e1, energy + w2) / width;
}

Referenced by get_integral_CS().

get_integral_CS()

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

overridevirtual

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

Implements Heed::PhotoAbsCS.

Definition at line 126 of file PhotoAbsCS.cpp.

                                                                {
  mfunname("double AveragePhotoAbsCS::get_integral_CS(...) const");
  if (width == 0.0 || energy1 >= energy2) {
    // Return the integral of the unmodified "real" cross-section.
    return real_pacs->get_integral_CS(energy1, energy2);
  }
  long q = long((energy2 - energy1) / step);
  if (q > max_q_step) {
    return real_pacs->get_integral_CS(energy1, energy2);
  }
  q++;
  const double rstep = (energy2 - energy1) / q;
  double x0 = energy1 + 0.5 * rstep;
  double s = 0.0;
  for (long n = 0; n < q; n++) s += get_CS(x0 + rstep * n);
  return s * rstep;
}

print()

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

overridevirtual

Reimplemented from Heed::PhotoAbsCS.

Definition at line 150 of file PhotoAbsCS.cpp.

                                                           {
  mfunname("void PhotoAbsCS::print(std::ostream& file, int l) const");
  Ifile << "AveragePhotoAbsCS: width = " << width << " step=" << step
        << " max_q_step=" << max_q_step << '\n';
  indn.n += 2;
  real_pacs->print(file, l);
  indn.n -= 2;
}

scale()

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

overridevirtual

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

Implements Heed::PhotoAbsCS.

Definition at line 145 of file PhotoAbsCS.cpp.

                                         {
  mfunname("void AveragePhotoAbsCS::scale(double fact)");
  real_pacs->scale(fact);
}

---

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

• PhotoAbsCS.h
• PhotoAbsCS.cpp