Heed::PairProd Class Reference

#include <PairProd.h>

Inheritance diagram for Heed::PairProd:

Public Member Functions
	PairProd ()
	PairProd (const std::string &file_name, double fwa, double ffactorFano=0.19)
double	get_eloss () const
	Calculate energy loss (in eV).
double	get_eloss (const double ecur) const
	Calculate energy loss (in eV) assuming V = wa / 2.
virtual void	print (std::ostream &file, int l) const
Public Member Functions inherited from Heed::RegPassivePtr
	RegPassivePtr (void)
	RegPassivePtr (char fs_ban_del, char fs_ban_sub, char fs_ban_cop=0)
	RegPassivePtr (const RegPassivePtr &f)
RegPassivePtr &	operator= (const RegPassivePtr &f)
CountPP_ns::CountPassivePtr *	book (void) const
void	clear_pointers (void) const
virtual RegPassivePtr *	copy () const
virtual	~RegPassivePtr ()
virtual void	print (std::ostream &file, int l=1) const
void	set_s_ban_del (char fs_ban_del)
char	get_s_ban_del (void) const
void	set_s_ban_sub (char fs_ban_sub)
char	get_s_ban_sub (void) const
void	set_s_ban_cop (char fs_ban_cop)
char	get_s_ban_cop (void) const
void	set_s_allow_del_at_zero_count (char fs_allow_del_at_zero_count)
char	get_s_allow_del_at_zero_count (void) const
long	get_total_number_of_references (void) const

Additional Inherited Members
Static Public Member Functions inherited from Heed::RegPassivePtr
static void	set_s_ban_del_ignore (char fs_ban_del_ignore)
static char	get_s_ban_del_ignore (void)
static void	set_s_print_adr_cpp (char fs_print_adr_cpp)
static char	get_s_print_adr_cpp (void)

Detailed Description

Production of electron-ion pairs by delta-electrons. All energies are assumed in eV. But it seems that it is not strict requirement. If to change this in file "file_name" and in interface, the energies could be in any units.

2003, I. Smirnov

Definition at line 17 of file PairProd.h.

Constructor & Destructor Documentation

PairProd() [1/2]

Heed::PairProd::PairProd ( )

inline

Definition at line 19 of file PairProd.h.

{}

PairProd() [2/2]

Heed::PairProd::PairProd	(	const std::string &	file_name,
		double	fwa,
		double	ffactorFano = 0.19
	)

Definition at line 14 of file PairProd.cpp.

    : wa(fwa), factorFano(ffactorFano) {
  mfunnamep(
      "PairProd::PairProd(const std::string& file_name, double fwa, double "
      "ffactorFano)");
 
  std::ifstream file(file_name.c_str());
  if (!file) {
    funnw.ehdr(mcerr);
    mcerr << "cannot open file " << file_name << std::endl;
    spexit(mcerr);
  }
  long q;
  file >> wa_table >> I_table >> J_table >> factorFano_table >> q;
  if (!file.good()) {
    funnw.ehdr(mcerr);
    mcerr << "error at reading file" << std::endl;
    spexit(mcerr);
  }
  std::vector<double> xx(q);
  std::vector<double> yy(q);
  for (long n = 0; n < q; n++) file >> xx[n] >> yy[n];
  pran = PointsRan(xx, yy, I_table, J_table);
  k = sqrt(factorFano * wa * wa / (factorFano_table * wa_table * wa_table));
  s = wa - k * wa_table;
  // s = wa - wa / factorFano * ( 2.0 * factorFano - factorFano_table);
  // k = wa / (factorFano * wa_table) * ( 2.0 * factorFano - factorFano_table);
}

Member Function Documentation

get_eloss() [1/2]

double Heed::PairProd::get_eloss

(

)

const

Calculate energy loss (in eV).

Definition at line 43 of file PairProd.cpp.

                                 {
  mfunname("double PairProd::get_eloss() const");
  return k * pran.ran(SRANLUX()) + s;
}

get_eloss() [2/2]

double Heed::PairProd::get_eloss

(

const double

ecur

)

const

Calculate energy loss (in eV) assuming V = wa / 2.

Definition at line 59 of file PairProd.cpp.

                                                   {
  mfunname("double PairProd::get_eloss(const double ecur) const");
  const double e_loss = k * pran.ran(SRANLUX()) + s;
  const double V_ratio = 0.5;
  const double v = V_ratio * wa / e_cur; 
  // const double c = 1. / (1. - v);
  const double c = v < 1. ? 1. / (1. - v * v) : DBL_MAX;
  return e_loss * c;
}

print()

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

virtual

Reimplemented from Heed::RegPassivePtr.

Definition at line 71 of file PairProd.cpp.

                                                  {
  if (l <= 0) return;
  Ifile << "PairProd:\n";
  indn.n += 2;
  Ifile << "wa=" << wa << " factorFano=" << factorFano << '\n';
  Ifile << "wa_table=" << wa_table << " factorFano_table=" << factorFano_table
        << '\n';
  Ifile << " I_table=" << I_table << " J_table=" << J_table << " k=" << k
        << " s=" << s << '\n';
  pran.print(file);
  indn.n -= 2;
}

---

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

• PairProd.h
• PairProd.cpp