Heed::SimpleTablePhotoAbsCS Class Reference

#include <PhotoAbsCS.h>

Inheritance diagram for Heed::SimpleTablePhotoAbsCS:

Public Member Functions
	SimpleTablePhotoAbsCS ()
	Default constructor.
	SimpleTablePhotoAbsCS (const std::string &fname, int fZ, double fthreshold, const std::string &ffile_name)
	Constructor for reading table from file.
	SimpleTablePhotoAbsCS (const std::string &fname, int fZ, double fthreshold, const std::vector< double > &fener, const std::vector< double > &fcs)
	Constructor from given energy and cross-section tables.
	SimpleTablePhotoAbsCS (const std::string &fname, int fZ, double fthreshold, int l, double E0, double yw, double ya, double P, double sigma)
	SimpleTablePhotoAbsCS (const SimpleTablePhotoAbsCS &total, const SimpleTablePhotoAbsCS &part, double emax_repl)
	Replace part of the cross-section table.
virtual	~SimpleTablePhotoAbsCS ()
	Destructor.
void	remove_leading_zeros ()
	Remove points with zero cross-section from the table.
void	remove_leading_tiny (double level)
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].
const std::vector< double > &	get_arr_ener () const
const std::vector< double > &	get_arr_CS () const
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 SimpleTablePhotoAbsCS *	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

Typically this is for reading Experimental CS, for example of argon, as if there is one shell. File is two-column table, the first row is energy in eV, the second one is CS in Mbarn. The points are understood as local points. The interpolation between them are either straight or by power function. the choice is determined by function sign_nonlinear_interpolation (see above). If the first point is not zero cross section, and the threshold is to the left, the straight extrapolation is performed from the two edge points. The extrapolation to the right is performed to the end of the energy mesh by power function with power -2.75. The minimal number of points is 2, as in PointCoorMesh from tline.h The zero number of points is allowed as well. Then the cross section is assumed to be zero.

Definition at line 160 of file PhotoAbsCS.h.

Constructor & Destructor Documentation

SimpleTablePhotoAbsCS() [1/5]

Heed::SimpleTablePhotoAbsCS::SimpleTablePhotoAbsCS

(

)

Default constructor.

Referenced by copy().

SimpleTablePhotoAbsCS() [2/5]

Heed::SimpleTablePhotoAbsCS::SimpleTablePhotoAbsCS	(	const std::string &	fname,
		int	fZ,
		double	fthreshold,
		const std::string &	ffile_name
	)

Constructor for reading table from file.

Definition at line 194 of file PhotoAbsCS.cpp.

    : PhotoAbsCS(fname, fZ, fthreshold), file_name(ffile_name) {
  mfunnamep("SimpleTablePhotoAbsCS::SimpleTablePhotoAbsCS(...)");
  std::ifstream file(file_name.c_str());
  if (!file) {
    funnw.ehdr(mcerr);
    mcerr << "cannot open file " << file_name << std::endl;
    spexit(mcerr);
  }
  ener.reserve(20);
  cs.reserve(20);
  do {
    // Read the energy.
    double x = 0.;
    file >> x;
    if (!file.good()) break;
    // Make sure it is non-negative and in ascending order.
    check_econd11(x, < 0.0, mcerr);
    if (!ener.empty()) check_econd12(x, <, ener.back(), mcerr);
    // Read the cross-section.
    double y = 0.;
    file >> y;
    if (!file.good()) break;
    check_econd11(y, < 0.0, mcerr);
    // Add the point to the table.
    ener.push_back(x * 1.e-6);
    cs.push_back(y);
  } while (1);
}

SimpleTablePhotoAbsCS() [3/5]

Heed::SimpleTablePhotoAbsCS::SimpleTablePhotoAbsCS	(	const std::string &	fname,
		int	fZ,
		double	fthreshold,
		const std::vector< double > &	fener,
		const std::vector< double > &	fcs
	)

Constructor from given energy and cross-section tables.

Definition at line 226 of file PhotoAbsCS.cpp.

    : PhotoAbsCS(fname, fZ, fthreshold),
      file_name("none"),
      ener(fener),
      cs(fcs) {
  mfunname("SimpleTablePhotoAbsCS::SimpleTablePhotoAbsCS(...)");
  check_econd12(ener.size(), !=, cs.size(), mcerr);
}

SimpleTablePhotoAbsCS() [4/5]

Heed::SimpleTablePhotoAbsCS::SimpleTablePhotoAbsCS	(	const std::string &	fname,
		int	fZ,
		double	fthreshold,
		int	l,
		double	E0,
		double	yw,
		double	ya,
		double	P,
		double	sigma
	)

Constructor from fit parameters. Fit formula from Band-Band-Trzaskovskaya et al. It is difficult to integrate those formulas analytically, so I create numerical array and treat it as input data.

Definition at line 238 of file PhotoAbsCS.cpp.

    : PhotoAbsCS(fname, fZ, fthreshold) {
  mfunname("SimpleTablePhotoAbsCS::SimpleTablePhotoAbsCS(...)");
  const long q = 1000;
  // Make a logarithmic energy mesh.
  ener.resize(q, 0.);
  const double emin = 2.e-6;
  const double emax = 2.e-1;
  const double rk = pow(emax / emin, (1.0 / double(q)));
  double e2 = emin;
  for (long n = 0; n < q; n++) {
    const double e1 = e2;
    e2 = e2 * rk;
    ener[n] = (e1 + e2) * 0.5;
  }
  cs.assign(q, 0.);
  for (long n = 0; n < q; n++) {
    double energy = ener[n];
    if (energy < threshold) continue;
    const double Q = 5.5 + l - 0.5 * P;
    const double y = energy / E0;
    const double Fpasc = ((y - 1) * (y - 1) + yw * yw) * pow(y, (-Q)) *
                         pow((1.0 + sqrt(y / ya)), (-P));
    cs[n] = Fpasc * sigma;
  }
  remove_leading_zeros();
}

SimpleTablePhotoAbsCS() [5/5]

Heed::SimpleTablePhotoAbsCS::SimpleTablePhotoAbsCS	(	const SimpleTablePhotoAbsCS &	total,
		const SimpleTablePhotoAbsCS &	part,
		double	emax_repl
	)

Replace part of the cross-section table.

Definition at line 269 of file PhotoAbsCS.cpp.

                                                               {
  mfunname(
      "SimpleTablePhotoAbsCS::SimpleTablePhotoAbsCS(const "
      "SimpleTablePhotoAbsCS& total,...)");
 
  *this = total;  // to assure that all is preserved
 
  long qe_i = total.ener.size();
 
  const std::vector<double>& ener_r = part.get_arr_ener();
  const std::vector<double>& cs_r = part.get_arr_CS();
  long qe_r = ener_r.size();
  std::vector<double> new_ener;
  std::vector<double> new_cs;
  // first write replacements
  for (long ne = 0; ne < qe_r; ne++) {
    if (ener_r[ne] >= total.get_threshold() && ener_r[ne] <= emax_repl) {
      new_ener.push_back(ener_r[ne]);
      new_cs.push_back(cs_r[ne]);
    }
  }
  for (long ne = 0; ne < qe_i; ne++) {
    if (ener[ne] >= total.get_threshold() && ener[ne] > emax_repl) {
      new_ener.push_back(total.ener[ne]);
      new_cs.push_back(total.cs[ne]);
    }
  }
  ener.swap(new_ener);
  cs.swap(new_cs);
}

~SimpleTablePhotoAbsCS()

virtual Heed::SimpleTablePhotoAbsCS::~SimpleTablePhotoAbsCS ( )

inlinevirtual

Destructor.

Definition at line 182 of file PhotoAbsCS.h.

{}

Member Function Documentation

copy()

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

inlinevirtual

Implements Heed::PhotoAbsCS.

Definition at line 197 of file PhotoAbsCS.h.

                                              {
    return new SimpleTablePhotoAbsCS(*this);
  }

get_arr_CS()

const std::vector< double > & Heed::SimpleTablePhotoAbsCS::get_arr_CS ( ) const

inline

Definition at line 194 of file PhotoAbsCS.h.

{ return cs; }

Referenced by Heed::ExAtomPhotoAbsCS::ExAtomPhotoAbsCS(), and SimpleTablePhotoAbsCS().

get_arr_ener()

const std::vector< double > & Heed::SimpleTablePhotoAbsCS::get_arr_ener ( ) const

inline

Definition at line 193 of file PhotoAbsCS.h.

{ return ener; }

Referenced by Heed::ExAtomPhotoAbsCS::ExAtomPhotoAbsCS(), and SimpleTablePhotoAbsCS().

get_CS()

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

virtual

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

Implements Heed::PhotoAbsCS.

Definition at line 337 of file PhotoAbsCS.cpp.

                                                        {
  mfunname("double SimpleTablePhotoAbsCS::get_CS(double energy) const");
  long q = ener.size();
  if (q == 0) return 0.0;
  check_econd11(q, == 1, mcerr);
  if (energy < threshold) return 0.;
  if (energy <= ener[q - 1]) {
    PointCoorMesh<double, const std::vector<double> > pcm(q, &ener);
    return t_value_generic_point_ar<
        double, std::vector<double>,
        PointCoorMesh<double, const std::vector<double> > >(
        pcm, cs, &my_val_fun, energy, 1, threshold, 0, DBL_MAX);
  } else {
    if (energy == DBL_MAX)
      return 0.0;
    else
      return cs[q - 1] * pow(energy, -2.75) / pow(ener[q - 1], -2.75);
  }
}

get_integral_CS()

double Heed::SimpleTablePhotoAbsCS::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 357 of file PhotoAbsCS.cpp.

                                                                    {
  mfunname("double SimpleTablePhotoAbsCS::get_integral_CS(...)");
 
  const long q = ener.size();
  if (q == 0) return 0.0;
  check_econd11(q, == 1, mcerr);
  if (energy2 < threshold) return 0.0;
  if (energy1 < threshold) energy1 = threshold;
  double s = 0.0;
  double energy21 = ener[q - 1];
  if (energy1 < energy21) {
    if (energy21 > energy2) energy21 = energy2;
    check_econd12(energy1, >, energy21, mcerr);
    PointCoorMesh<double, const std::vector<double> > pcm(q, &ener);
    s = t_integ_generic_point_ar<
        double, std::vector<double>,
        PointCoorMesh<double, const std::vector<double> > >(
        pcm, cs, &my_integr_fun, energy1, energy21, 1, threshold, 0, DBL_MAX);
  }
  // print(mcout, 3);
  // mcout << "energy1="<<energy1
  //      << " energy21="<<energy21
  //      << " ener[q-1]="<<ener[q-1]
  //      << " threshold="<<threshold
  //      << " s="<<s<<'\n';
  check_econd11(s, < 0.0, mcout);
  if (energy2 > ener[q - 1]) {
    // add tail
    if (energy2 == DBL_MAX) {
      if (energy1 < ener[q - 1]) energy1 = ener[q - 1];
      double c =
          cs[q - 1] / (1.75 * pow(ener[q - 1], -2.75)) * pow(energy1, -1.75);
      // check_econd11(c , < 0.0, mcout);
      s += c;
    } else {
      if (energy1 < ener[q - 1]) energy1 = ener[q - 1];
      double c = cs[q - 1] / (1.75 * pow(ener[q - 1], -2.75)) *
                 (pow(energy1, -1.75) - pow(energy2, -1.75));
      // check_econd11(c , < 0.0, mcout);
      s += c;
    }
  }
  return s;
}

print()

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

virtual

Reimplemented from Heed::PhotoAbsCS.

Definition at line 409 of file PhotoAbsCS.cpp.

                                                               {
  if (l <= 0) return;
  Ifile << "SimpleTablePhotoAbsCS: name=" << name << " Z = " << Z << "\n";
  Ifile << " threshold = " << threshold << " file_name=" << file_name << "\n";
  if (l > 1) {
    indn.n += 2;
    const long q = ener.size();
    for (long n = 0; n < q; ++n) {
      Ifile << "n=" << n << " ener=" << ener[n] << " cs=" << cs[n] << "\n";
    }
    indn.n -= 2;
  }
}

remove_leading_tiny()

void Heed::SimpleTablePhotoAbsCS::remove_leading_tiny

(

double

level

)

Remove points with cross section below a given level from the table. The function is designed for Henke tables, which are prepared for database with leading values like 1.0e-15 . Both functions allow to use the straight interpolation to threshold

Definition at line 319 of file PhotoAbsCS.cpp.

                                                            {
  const long q = ener.size();
  long ne = 0;
  for (ne = 0; ne < q; ne++) {
    if (cs[ne] > level) break;
  }
  if (ne <= 0) return;
  const long qn = q - ne;
  std::vector<double> enern(qn);
  std::vector<double> csn(qn);
  for (long nez = ne; nez < q; nez++) {
    enern[nez - ne] = ener[nez];
    csn[nez - ne] = cs[nez];
  }
  ener = enern;
  cs = csn;
}

Referenced by Heed::ExAtomPhotoAbsCS::ExAtomPhotoAbsCS().

remove_leading_zeros()

void Heed::SimpleTablePhotoAbsCS::remove_leading_zeros

(

)

Remove points with zero cross-section from the table.

Definition at line 302 of file PhotoAbsCS.cpp.

                                                 {
  const long q = ener.size();
  long ne = 0;
  for (ne = 0; ne < q; ne++) {
    if (cs[ne] > 0.0) break;
  }
  if (ne <= 0) return;
  const long qn = q - ne;
  std::vector<double> enern(qn);
  std::vector<double> csn(qn);
  for (long nez = ne; nez < q; nez++) {
    enern[nez - ne] = ener[nez];
    csn[nez - ne] = cs[nez];
  }
  ener = enern;
  cs = csn;
}

Referenced by Heed::ExAtomPhotoAbsCS::ExAtomPhotoAbsCS(), and SimpleTablePhotoAbsCS().

scale()

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

virtual

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

Implements Heed::PhotoAbsCS.

Definition at line 403 of file PhotoAbsCS.cpp.

                                             {
  mfunnamep("void SimpleTablePhotoAbsCS::scale(double fact)");
  const long q = ener.size();
  for (long n = 0; n < q; ++n) cs[n] *= fact;
}

---

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

• PhotoAbsCS.h
• PhotoAbsCS.cpp