Heed::SimpleAtomPhotoAbsCS Class Reference

#include <PhotoAbsCS.h>

Inheritance diagram for Heed::SimpleAtomPhotoAbsCS:

Public Member Functions
virtual double	get_threshold (int nshell) const
virtual double	get_ACS (double energy) const
virtual double	get_integral_ACS (double energy1, double energy2) const
virtual double	get_ACS (int nshell, double energy) const
virtual double	get_integral_ACS (int nshell, double energy1, double energy2y) const
virtual double	get_ICS (double energy) const
virtual double	get_integral_ICS (double energy1, double energy2) const
virtual double	get_ICS (int nshell, double energy) const
virtual double	get_integral_ICS (int nshell, double energy1, double energy2) const
virtual int	get_main_shell_number (int nshell) const
virtual void	print (std::ostream &file, int l) const
	macro_copy_total (SimpleAtomPhotoAbsCS)
	SimpleAtomPhotoAbsCS (void)
	SimpleAtomPhotoAbsCS (int fZ, const String &ffile_name)
	SimpleAtomPhotoAbsCS (int fZ, const PhotoAbsCS &fasc)
virtual	~SimpleAtomPhotoAbsCS ()
Public Member Functions inherited from Heed::AtomPhotoAbsCS
int	get_Z () const
int	get_qshell () const
virtual double	get_threshold (int nshell) const =0
virtual double	get_I_min (void) const
virtual double	get_ACS (double energy) const =0
virtual double	get_integral_ACS (double energy1, double energy2) const =0
virtual double	get_ACS (int nshell, double energy) const =0
virtual double	get_integral_ACS (int nshell, double energy1, double energy2) const =0
virtual double	get_ICS (double energy) const =0
virtual double	get_TICS (double energy, double factual_minimal_threshold) const
virtual double	get_integral_ICS (double energy1, double energy2) const =0
virtual double	get_integral_TICS (double energy1, double energy2, double factual_minimal_threshold) const
virtual double	get_ICS (int nshell, double energy) const =0
virtual double	get_TICS (int nshell, double energy, double factual_minimal_threshold) const
virtual double	get_integral_ICS (int nshell, double energy1, double energy2) const =0
virtual double	get_integral_TICS (int nshell, double energy1, double energy2, double factual_minimal_threshold) const
virtual void	get_escape_particles (int nshell, double energy, DynLinArr< double > &el_energy, DynLinArr< double > &ph_energy) const
virtual int	get_main_shell_number (int nshell) const =0
virtual void	remove_shell (int nshell)
virtual void	restore_shell (int nshell)
virtual void	print (std::ostream &file, int l) const
	macro_copy_total_zero (AtomPhotoAbsCS)
	AtomPhotoAbsCS (void)
AtomicSecondaryProducts *	get_asp (int nshell)

Protected Attributes
String	file_name
DynLinArr< ActivePtr< PhotoAbsCS > >	acs
Protected Attributes inherited from Heed::AtomPhotoAbsCS
String	name
int	Z
int	qshell
DynLinArr< int >	s_ignore_shell
DynLinArr< AtomicSecondaryProducts >	asp

Detailed Description

Definition at line 401 of file PhotoAbsCS.h.

Constructor & Destructor Documentation

SimpleAtomPhotoAbsCS() [1/3]

Heed::SimpleAtomPhotoAbsCS::SimpleAtomPhotoAbsCS

(

void

)

Definition at line 1345 of file PhotoAbsCS.cpp.

: AtomPhotoAbsCS() {}

SimpleAtomPhotoAbsCS() [2/3]

Heed::SimpleAtomPhotoAbsCS::SimpleAtomPhotoAbsCS	(	int	fZ,
		const String &	ffile_name
	)

Definition at line 1347 of file PhotoAbsCS.cpp.

    : file_name(ffile_name) {
  mfunnamep("SimpleAtomPhotoAbsCS::SimpleAtomPhotoAbsCS(int fZ, const String& "
            "ffile_name)");
  check_econd11(fZ, < 1, mcerr);
#ifdef USE_STLSTRING
  std::ifstream file(file_name.c_str());
#else
  std::ifstream file(file_name);
#endif
  if (!file) {
    funnw.ehdr(mcerr);
    mcerr << "cannot open file " << file_name << std::endl;
    spexit(mcerr);
  }
  while (findmark(file, "#") == 1) {
    file >> Z;
    if (Z == fZ) {
      file >> qshell;
      check_econd21(qshell, < 1 ||, > 10000, mcerr);
      s_ignore_shell.put_qel(qshell, 0);
      file >> name;
      acs = DynLinArr<ActivePtr<PhotoAbsCS> >(qshell);
      asp = DynLinArr<AtomicSecondaryProducts>(qshell);
      DynLinArr<double> fl(qshell);
      int sZshell = 0;
      for (int nshell = 0; nshell < qshell; ++nshell) {
        double thr = 0.0;
        int Zshell = 0;
        //double fl;
        String shell_name;
        file >> thr;
        check_econd11(thr, <= 0.0, mcerr);
        file >> Zshell;
        check_econd11(Zshell, <= 0, mcerr);
        sZshell += Zshell;
        file >> fl[nshell];
        findmark(file, "!");
        file >> shell_name;
        //PhotoAbsCS* ap = new PhenoPhotoAbsCS(shell_name, Zshell, thr);
        //acs[nshell].pass( ap );
        acs[nshell].pass(new PhenoPhotoAbsCS(shell_name, Zshell, thr * 1.0e-6));
      }
      check_econd12(sZshell, !=, Z, mcerr);
 
      int n_min = 0;
      double st = DBL_MAX;
      for (int nshell = 0; nshell < qshell; ++nshell) {
        // currently the minimal shell is the last,
        // but to avoid this assumption we check all
        if (get_threshold(nshell) < st) n_min = nshell;
      }
      for (int nshell = 0; nshell < qshell; ++nshell) {
        if (fl[nshell] > 0) {
          check_econd12(nshell, ==, n_min, mcerr);
          DynLinArr<double> felectron_energy(0);
          DynLinArr<double> fphoton_energy(1);
          fphoton_energy[0] = get_threshold(nshell) - get_threshold(n_min);
          asp[nshell].add_channel(fl[nshell], felectron_energy, fphoton_energy);
        }
      }
      return;
    }
  }
  funnw.ehdr(mcerr);
  mcerr << "there is no element Z=" << fZ << " in file " << file_name
        << std::endl;
  spexit(mcerr);
}

SimpleAtomPhotoAbsCS() [3/3]

Heed::SimpleAtomPhotoAbsCS::SimpleAtomPhotoAbsCS	(	int	fZ,
		const PhotoAbsCS &	fasc
	)

Definition at line 1417 of file PhotoAbsCS.cpp.

                                                                         {
  mfunname("SimpleAtomPhotoAbsCS::SimpleAtomPhotoAbsCS(int fZ, const "
           "PhotoAbsCS& facs)");
  check_econd11(fZ, <= 0, mcerr);
  check_econd12(fZ, !=, facs.get_Z(), mcerr);
  Z = fZ;
  qshell = 1;
  s_ignore_shell.put_qel(qshell, 0);
  name = facs.get_name();
  acs.put_qel(1);
  acs[0].put(&facs);
}

~SimpleAtomPhotoAbsCS()

virtual Heed::SimpleAtomPhotoAbsCS::~SimpleAtomPhotoAbsCS ( )

inlinevirtual

Definition at line 426 of file PhotoAbsCS.h.

{}

Member Function Documentation

get_ACS() [1/2]

double Heed::SimpleAtomPhotoAbsCS::get_ACS ( double  energy ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 1436 of file PhotoAbsCS.cpp.

                                                        {
  mfunname("double SimpleAtomPhotoAbsCS::get_ACS(double energy) const");
  double s = 0.0;
  for (int n = 0; n < qshell; ++n) {
    if (s_ignore_shell[n] == 0) {
      s += acs[n]->get_CS(energy);
    }
  }
  return s;
}

get_ACS() [2/2]

double Heed::SimpleAtomPhotoAbsCS::get_ACS ( int  nshell, double  energy  ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 1468 of file PhotoAbsCS.cpp.

                                                                    {
  mfunname("double SimpleAtomPhotoAbsCS::get_ACS(int nshell, double energy)");
  check_econd21(nshell, < 0 ||, > qshell, mcerr);
  if (s_ignore_shell[nshell] == 0) {
    return acs[nshell]->get_CS(energy);
  }
  return 0.0;
}

get_ICS() [1/2]

double Heed::SimpleAtomPhotoAbsCS::get_ICS ( double  energy ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 1488 of file PhotoAbsCS.cpp.

                                                        {
  mfunname("double SimpleAtomPhotoAbsCS::get_ICS(double energy) const");
  double s = 0.0;
  for (int n = 0; n < qshell; ++n) {
    if (s_ignore_shell[n] == 0) {
      s += acs[n]->get_CS(energy);
    }
  }
  return s;
}

get_ICS() [2/2]

double Heed::SimpleAtomPhotoAbsCS::get_ICS ( int  nshell, double  energy  ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 1512 of file PhotoAbsCS.cpp.

                                                                    {
  mfunname("double SimpleAtomPhotoAbsCS::get_ICS(int nshell, double energy)");
  check_econd21(nshell, < 0 ||, > qshell, mcerr);
  if (s_ignore_shell[nshell] == 0) {
    return acs[nshell]->get_CS(energy);
  }
  return 0.0;
}

get_integral_ACS() [1/2]

double Heed::SimpleAtomPhotoAbsCS::get_integral_ACS ( double  energy1, double  energy2  ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 1446 of file PhotoAbsCS.cpp.

                                                                    {
  mfunnamep(
      "double SimpleAtomPhotoAbsCS::get_integral_ACS(double energy)const");
  double s = 0.0;
  for (int n = 0; n < qshell; ++n) {
    if (s_ignore_shell[n] == 0) {
      double t;
      s += (t = acs[n]->get_integral_CS(energy1, energy2));
      if (t < 0) {
        funnw.ehdr(mcout);
        mcout << "t < 0\n";
        Iprintn(mcout, t);
        print(mcout, 4);
        spexit(mcout);
      }
    }
    //check_econd11a(t , < 0 , "n="<<n<<'\n',  mcerr);
  }
  return s;
}

get_integral_ACS() [2/2]

double Heed::SimpleAtomPhotoAbsCS::get_integral_ACS ( int  nshell, double  energy1, double  energy2y  ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 1477 of file PhotoAbsCS.cpp.

                                                                    {
  mfunname("double SimpleAtomPhotoAbsCS::get_integral_ACS(int nshell, double "
           "energy1, double energy2)");
  check_econd21(nshell, < 0 ||, > qshell, mcerr);
  if (s_ignore_shell[nshell] == 0) {
    return acs[nshell]->get_integral_CS(energy1, energy2);
  }
  return 0.0;
}

get_integral_ICS() [1/2]

double Heed::SimpleAtomPhotoAbsCS::get_integral_ICS ( double  energy1, double  energy2  ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 1499 of file PhotoAbsCS.cpp.

                                                                    {
  mfunname("double SimpleAtomPhotoAbsCS::get_integral_ICS(double energy1, "
           "double energy2) const");
  double s = 0.0;
  for (int n = 0; n < qshell; ++n) {
    if (s_ignore_shell[n] == 0) {
      s += acs[n]->get_integral_CS(energy1, energy2);
    }
  }
  return s;
}

get_integral_ICS() [2/2]

double Heed::SimpleAtomPhotoAbsCS::get_integral_ICS ( int  nshell, double  energy1, double  energy2  ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 1521 of file PhotoAbsCS.cpp.

                                                                    {
  mfunname("double SimpleAtomPhotoAbsCS::get_integral_ICS(int nshell, double "
           "energy1, double energy2)");
  check_econd21(nshell, < 0 ||, > qshell, mcerr);
  if (s_ignore_shell[nshell] == 0) {
    return acs[nshell]->get_integral_CS(energy1, energy2);
  }
  return 0.0;
}

get_main_shell_number()

int Heed::SimpleAtomPhotoAbsCS::get_main_shell_number ( int  nshell ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 1550 of file PhotoAbsCS.cpp.

                                                                {
  mfunname("int SimpleAtomPhotoAbsCS::get_main_shell_number(int nshell) const");
  String shell_name = acs[nshell]->get_name();
#ifdef STRSTREAM_AVAILABLE
#ifdef USE_STLSTRING
  istrstream sfile(shell_name.c_str());
#else
  istrstream sfile(shell_name);
#endif
#else
  std::istringstream sfile(shell_name.c_str());
#endif
  int i = -100;
  sfile >> i;
  //Iprintn(mcout, i);
  //check_econd(i <  1  || i > 50 , "i="<<i<<" shell_name="<<shell_name<<'\n' ,
  //                 mcerr);
  if (i < 1 || i > 50) {
    i = -1;
  }
  return i;
}

get_threshold()

double Heed::SimpleAtomPhotoAbsCS::get_threshold ( int  nshell ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 1430 of file PhotoAbsCS.cpp.

                                                           {
  mfunname("double SimpleAtomPhotoAbsCS::get_threshold(int nshell) const");
  check_econd21(nshell, < 0 ||, > qshell, mcerr);
  return acs[nshell]->get_threshold();
}

Referenced by SimpleAtomPhotoAbsCS().

macro_copy_total()

Heed::SimpleAtomPhotoAbsCS::macro_copy_total

(

SimpleAtomPhotoAbsCS

)

print()

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

virtual

Reimplemented from Heed::AtomPhotoAbsCS.

Definition at line 1532 of file PhotoAbsCS.cpp.

                                                              {
  if (l > 0) {
    Ifile << "SimpleAtomPhotoAbsCS(l=" << l << "): name=" << name
          << " Z = " << Z << " qshell = " << qshell
          << " file_name=" << file_name << std::endl;
    l--;
    if (l > 0) {
      indn.n += 2;
      for (int n = 0; n < qshell; ++n) {
        Ifile << "nshell=" << n << std::endl;
        acs[n].print(file, l);
      }
      AtomPhotoAbsCS::print(file, l);
      indn.n -= 2;
    }
  }
}

Referenced by get_integral_ACS().

Member Data Documentation

acs

DynLinArr<ActivePtr<PhotoAbsCS> > Heed::SimpleAtomPhotoAbsCS::acs

protected

Definition at line 430 of file PhotoAbsCS.h.

Referenced by get_ACS(), get_ICS(), get_integral_ACS(), get_integral_ICS(), get_main_shell_number(), get_threshold(), print(), and SimpleAtomPhotoAbsCS().

file_name

String Heed::SimpleAtomPhotoAbsCS::file_name

protected

Definition at line 429 of file PhotoAbsCS.h.

Referenced by print(), and SimpleAtomPhotoAbsCS().

---

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

• PhotoAbsCS.h
• PhotoAbsCS.cpp