Heed::ExAtomPhotoAbsCS Class Reference

#include <PhotoAbsCS.h>

Inheritance diagram for Heed::ExAtomPhotoAbsCS:

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 energy2) 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
void	replace_shells_by_overage (double fwidth, double fstep, long fmax_q_step)
virtual void	print (std::ostream &file, int l) const
	macro_copy_total (ExAtomPhotoAbsCS)
	ExAtomPhotoAbsCS (void)
	ExAtomPhotoAbsCS (int fZ, const String &fthreshold_file_name, const String &fsimple_table_file_name, const String &fname="none", double fminimal_threshold=0.0)
	ExAtomPhotoAbsCS (int fZ, const String &fname, const String &fBT_file_name, int id, double fminimal_threshold=0.0)
	ExAtomPhotoAbsCS (int fZ, const String &fname, const String &fFitBT_file_name, int id, int s_no_scale, double fminimal_threshold=0.0)
	ExAtomPhotoAbsCS (int fZ, const String &fname, const String &fFitBT_file_name, const String &fsimple_table_file_name, double emax_repl, int id, double fminimal_threshold=0.0)
virtual	~ExAtomPhotoAbsCS ()
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	threshold_file_name
String	simple_table_file_name
String	BT_file_name
DynLinArr< ActivePtr< PhotoAbsCS > >	acs
double	integ_abs_before_corr
double	integ_abs_after_corr
double	integ_ioniz_after_corr
double	height_of_excitation
double	exener [2]
double	minimal_threshold
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 445 of file PhotoAbsCS.h.

Constructor & Destructor Documentation

ExAtomPhotoAbsCS() [1/5]

Heed::ExAtomPhotoAbsCS::ExAtomPhotoAbsCS ( void  )

inline

Definition at line 466 of file PhotoAbsCS.h.

: AtomPhotoAbsCS() {}

Referenced by ExAtomPhotoAbsCS().

ExAtomPhotoAbsCS() [2/5]

Heed::ExAtomPhotoAbsCS::ExAtomPhotoAbsCS	(	int	fZ,
		const String &	fthreshold_file_name,
		const String &	fsimple_table_file_name,
		const String &	fname = "none",
		double	fminimal_threshold = 0.0
	)

Definition at line 1575 of file PhotoAbsCS.cpp.

    : threshold_file_name(fthreshold_file_name),
      simple_table_file_name(fsimple_table_file_name),
      BT_file_name("none"),
      minimal_threshold(fminimal_threshold) {
  mfunnamep("SimpleAtomPhotoAbsCS::SimpleAtomPhotoAbsCS(int fZ, const String& "
            "fname, , const String& fthreshold_file_name, const String& "
            "fsimple_table_file_name, double fminimal_threshold)");
  //Imcout<<"ExAtomPhotoAbsCS::ExAtomPhotoAbsCS is run for fZ="<<fZ<<std::endl;
  check_econd11(fZ, < 1, mcerr);
#ifdef USE_STLSTRING
  std::ifstream threshold_file(threshold_file_name.c_str());
#else
  std::ifstream threshold_file(threshold_file_name);
#endif
  if (!threshold_file) {
    funnw.ehdr(mcerr);
    mcerr << "cannot open file " << threshold_file_name << std::endl;
    spexit(mcerr);
  }
  DynLinArr<double> thr(0, 0.0);
  DynLinArr<int> Zshell(0, 0);
  DynLinArr<double> fl(0, 0.0);
  DynLinArr<String> shell_name(0);
  while (findmark(threshold_file, "#") == 1) {
    threshold_file >> Z;
    if (Z == fZ) {
      threshold_file >> qshell;
      check_econd21(qshell, < 1 ||, > 10000, mcerr);
      s_ignore_shell.put_qel(qshell, 0);
      //Iprintn(mcout, qshell);
      thr = DynLinArr<double>(qshell, 0.0);
      Zshell = DynLinArr<int>(qshell, 0);
      fl = DynLinArr<double>(qshell, 0.0);
      shell_name = DynLinArr<String>(qshell);
      acs = DynLinArr<ActivePtr<PhotoAbsCS> >(qshell);
      asp = DynLinArr<AtomicSecondaryProducts>(qshell);
      String temp_name;
      threshold_file >> temp_name;
      if (fname == "none")
        name = temp_name;
      else
        name = fname;
      int nshell;
      int sZshell = 0;
      for (nshell = 0; nshell < qshell; nshell++) {
        threshold_file >> thr[nshell];
        check_econd11(thr[nshell], <= 0.0, mcerr);
        thr[nshell] *= 1.0e-6;
        threshold_file >> Zshell[nshell];
        check_econd11(Zshell[nshell], <= 0, mcerr);
        sZshell += Zshell[nshell];
        threshold_file >> fl[nshell];
        findmark(threshold_file, "!");
        threshold_file >> shell_name[nshell];
        //PhotoAbsCS* ap = new PhenoPhotoAbsCS(shell_name, Zshell, thr);
        //acs[nshell].pass( ap );
      }
      check_econd12(sZshell, !=, Z, mcerr);
      int n_min = 0;
      double st = DBL_MAX;
      for (nshell = 0; nshell < qshell;
           nshell++) {  // currently the minimal shell is the last,
                        // but to avoid this assumption
                        // we check all
        if (thr[nshell] < st) {
          n_min = nshell;
          st = thr[nshell];
        }
      }
      for (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] = thr[nshell] - thr[n_min];
          asp[nshell].add_channel(fl[nshell], felectron_energy, fphoton_energy);
        }
      }
      goto mark1;
    }
  }
  funnw.ehdr(mcerr);
  mcerr << "there is no element Z=" << fZ << " in file " << threshold_file_name
        << std::endl;
  spexit(mcerr);
mark1:
  // Here it reads the PACS as an one shell curve:
  SimpleTablePhotoAbsCS stpacs(name, Z, 0.0, fsimple_table_file_name);
  const DynLinArr<double>& ener = stpacs.get_arr_ener();
  const DynLinArr<double>& CS = stpacs.get_arr_CS();
  //Iprintn(mcout, ener.get_qel());
  DynLinArr<double> left_CS = CS;             // used in sequencial algorithm
  DynLinArr<DynLinArr<double> > SCS(qshell);  // here cs is saved
  int n;
  for (n = 0; n < qshell; n++)
    SCS[n] = DynLinArr<double>(ener.get_qel(), 0.0);
  /*
  if(ener[0] < thr[qshell-1])
  {
    funnw.ehdr(cerr);
    mcerr<<"ener[0] < thr[qshell-1]\n";
    mcerr<<"This possibility is not bad, but it is not implemented \n"
         <<"in the program yet, sorry\n";
    spexit(mcerr);
  }
  */
  int nct = qshell - 1;  //"current" threshold index
  long nce = 0;          // "current" energy index
                         // Let uss ignore values below the lowest threshold
  // It is not clear whether it is right, perhaps this is one of possible ways
  if (ener[0] < thr[qshell - 1]) {
    long ne;
    for (ne = 0; ne < ener.get_qel() && (ener[ne] < thr[qshell - 1] ||
                                         (ener[ne] >= thr[qshell - 1] &&
                                          ne > 1 && CS[ne - 1] <= CS[ne - 2]));
         ne++) {
      if (ne > 0) left_CS[ne - 1] = 0.0;
      nce = ne;
    }
  }
  //Iprintn(mcout, nce);
  //mcout<<"ener[nce]="<<ener[nce]<<" CS[nce]="<<CS[nce]
  //     <<"  CS[nce+1]="<<CS[nce+1]<<std::endl;
  int s_more;
  int nt2 = 0;  // < nt1
  int s_spes = 0;
  do  // Actually this is a loop by the group of thresholds
      {
    //mcout<<"nct="<<nct<<" nce="<<nce<<std::endl;
    // Find all thresholds which are less then the current energy
    int nt;
    s_more = 0;  // sign that there are thresholds more than
                 // the current energy
    for (nt = nct; nt >= 0; nt--) {
      if (s_spes == 0) {
        if (thr[nt] > ener[nce]) {
          s_more = 1;
          break;
        }
      } else {
        if (thr[nt] > ener[nce + 1]) {
          s_more = 1;
          break;
        }
      }
    }
    // nt is now index of the next threshold or -1, if the thresholds are
    // made up.
    int nt1 = nct;
    int nce_next = ener.get_qel();
    nt2 = nt + 1;
    //mcout<<"nt="<<nt<<" nt1="<<nt1<<" nt2="<<nt2<<" s_more="<<s_more<<'\n';
    if (s_more == 1) {
      //if(nt >= 0)  // so if there are other larger thresholds,
      //{        // we should check how far we can pass at this step
      long ne;
      for (ne = nce; ne < ener.get_qel();
           ne++) {  // finding energy larger than the next threshold
        if (thr[nt] <= ener[ne]) {
          nce_next = ne;
          s_spes = 0;
          break;
        } else {  // At the following condition energy could be less then
                  // threshold,
          // but this point will anyway be assotiated with the next shell
          // corresponding to this threshold.
          // This is related to not precise measurement of cross section
          // and not precise knowledge of shell energies.
          // Occurence of this condition is marked by s_spes = 1.
          // At the next passing of this loop the thresholds are compared with
          // the next energy.
          if (ne > 1 && ne < ener.get_qel() - 1 && ne > nce + 2 &&
              thr[nt] <= ener[ne + 1] &&
              (thr[nt] - ener[ne]) / (ener[ne + 1] - ener[ne]) < 0.1 &&
              CS[ne] > CS[ne - 1]) {
            //mcout<<"special condition is satisf.\n";
            nce_next = ne;
            s_spes = 1;
            break;
          }
        }
      }
      if (ne == ener.get_qel())  // threshold is larger then energy mesh
        s_more = 0;              // to finish the loop
    }
    //Iprintn(mcout, nce_next);
    //Iprintn(mcout, ener[nce_next-1]);
    int qt = nt1 - nt2 + 1;   // quantity of the new thresholds
                              //Iprintn(mcout, qt);
    DynLinArr<double> w(qt);  // weights accouring to charges
    int s = 0;                // sum of Z
    for (nt = 0; nt < qt; nt++) {
      int nshell = nct - nt;
      s += Zshell[nshell];
    }
    for (nt = 0; nt < qt; nt++) {
      int nshell = nct - nt;
      w[nt] = double(Zshell[nshell]) / s;
    }
    double save_left_CS = left_CS[nce_next - 1];
    long ne;
    for (ne = 0; ne < nce_next; ne++) {
      for (nt = 0; nt < qt; nt++) {
        int nshell = nct - nt;
        SCS[nshell][ne] = left_CS[ne] * w[nt];
      }
      left_CS[ne] = 0.0;
    }
    for (ne = nce_next; ne < ener.get_qel(); ne++) {
      double extrap_CS =
          save_left_CS * pow(ener[nce_next - 1], 2.75) / pow(ener[ne], 2.75);
      if (extrap_CS > left_CS[ne]) extrap_CS = left_CS[ne];
      for (nt = 0; nt < qt; nt++) {
        int nshell = nct - nt;
        SCS[nshell][ne] = extrap_CS * w[nt];
      }
      left_CS[ne] -= extrap_CS;
    }
    nce = nce_next;
    nct = nt2 - 1;
  } while (s_more != 0);
  // now nt2 will be index of last filled shell
  // Now to fill the shells which are absent in the input table.
  // They will be initialized phenomenologically, basing on the sum rule:
  int ns;
  //Iprintn(mcout, nt2);
  for (ns = 0; ns < nt2; ns++) {
    acs[ns].pass(new PhenoPhotoAbsCS(shell_name[ns], Zshell[ns], thr[ns]));
  }
  // Initialization of input shells:
  for (ns = nt2; ns < qshell; ns++) {
    //acs[ns].pass( new SimpleTablePhotoAbsCS(shell_name[ns], Zshell[ns],
    //                                thr[ns],
    //                                ener, SCS[ns]) );
    SimpleTablePhotoAbsCS* adr = new SimpleTablePhotoAbsCS(
        shell_name[ns], Zshell[ns], thr[ns], ener, SCS[ns]);
    adr->remove_leading_zeros();
    acs[ns].pass(adr);
  }
  height_of_excitation = 0.0;
  exener[0] = exener[1] = 0.0;
  double integ = get_integral_ACS(0.0, DBL_MAX);
  //Iprintn(mcout, integ);
  integ_abs_before_corr = integ;
  double pred_integ = Thomas_sum_rule_const_Mb * Z;
  //Iprintn(mcout, pred_integ);
  if (pred_integ > integ) {
    if (s_add_excitations_to_normalize == 1) {
      // add excitation
      exener[0] =
          low_boundary_of_excitations * acs[qshell - 1]->get_threshold();
      exener[1] = 1.0 * acs[qshell - 1]->get_threshold();
      height_of_excitation = (pred_integ - integ) / (exener[1] - exener[0]);
      if (minimal_threshold > 0.0) {
        if (minimal_threshold > acs[qshell - 1]->get_threshold()) {
          // currently the minimal shell is the last one
          exener[0] += minimal_threshold - acs[qshell - 1]->get_threshold();
          exener[1] += minimal_threshold - acs[qshell - 1]->get_threshold();
        }
      }
    }
    /*
    else
    {
      height_of_excitation = 0.0;
      exener[0] = exener[1] = 0.0;
    }
    */
  } else if (pred_integ < integ) {
    if (s_scale_to_normalize_if_more == 1) {
      const double fact = pred_integ / integ;
      for (int nshell = 0; nshell < qshell; ++nshell) {
        acs[nshell]->scale(fact);
      }
    }
  }
 
  integ_abs_after_corr = get_integral_ACS(0.0, DBL_MAX);
  integ_ioniz_after_corr = get_integral_ICS(0.0, DBL_MAX);
}

ExAtomPhotoAbsCS() [3/5]

Heed::ExAtomPhotoAbsCS::ExAtomPhotoAbsCS	(	int	fZ,
		const String &	fname,
		const String &	fBT_file_name,
		int	id,
		double	fminimal_threshold = 0.0
	)

Definition at line 1860 of file PhotoAbsCS.cpp.

    : threshold_file_name("none"),
      simple_table_file_name("none"),
      BT_file_name(fBT_file_name),
      minimal_threshold(fminimal_threshold) {
  mfunnamep("ExAtomPhotoAbsCS::ExAtomPhotoAbsCS(int fZ, const String& fname, "
            "const String& fBT_file_name, int id, double fminimal_threshold)");
  check_econd11(fZ, < 1, mcerr);
  check_econd21(id, < 1 ||, > 2, mcerr);
 
  name = fname;
#ifdef USE_STLSTRING
  std::ifstream BT_file(BT_file_name.c_str());
#else
  std::ifstream BT_file(BT_file_name);
#endif
  if (!BT_file) {
    funnw.ehdr(mcerr);
    mcerr << "cannot open file " << BT_file_name << std::endl;
    spexit(mcerr);
  }
  DynLinArr<double> thresh(0, 0.0);
  DynLinArr<double> fl(0, 0.0);
  Z = fZ;
  int i = findmark(BT_file, "NUCLEAR CHARGE =");
  check_econd11a(i, != 1, "wrong file format", mcerr);
  int Z_from_file;
  BT_file >> Z_from_file;
  check_econd12(Z_from_file, !=, Z, mcerr);
  qshell = 0;
  while ((i = findmark(BT_file, "Z =")) == 1) {
    BT_file >> i;
    check_econd11(i, != Z, mcerr);
    String shellname;
    BT_file >> shellname;
    //Iprintn(mcout, shellname);
    i = findmark(BT_file, "$");
    check_econd11(i, != 1, mcerr);
    long qen;
    BT_file >> qen;
    check_econd11(qen, <= 0, mcerr);
    DynLinArr<double> fener(qen, 0.0);
    DynLinArr<double> fcs(qen, 0.0);
    double thr = 0.0;
    BT_file >> thr;
    check_econd11(thr, <= 0, mcerr);
    thr *= 1.0e-3;  // pass from keV to MeV
    if (id == 2) {
      thresh.increment();
      fl.increment();
      thresh[qshell] = thr;
      BT_file >> fl[qshell];
      check_econd21(fl[qshell], < 0.0 ||, > 1.0, mcerr);
      //Iprintn(mcout, fl[qshell]);
    }
    long nen;
    for (nen = 0; nen < qen; nen++) {
      BT_file >> fener[nen] >> fcs[nen];
      check_econd11(fener[nen], <= 0.0, mcerr);
      check_econd11(fcs[nen], < 0.0, mcerr);
      fener[nen] *= 1.0e-3;  // pass from keV to MeV
    }
    qshell++;
    acs.put_qel(qshell);
    acs[qshell - 1]
        .pass(new SimpleTablePhotoAbsCS(shellname, 0,  // unknown here
                                        thr, fener, fcs));
  }
  if (id == 2) {
    // a copy of similar thing from subroutine above
    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 (thresh[nshell] < st) {
        n_min = nshell;
        st = thresh[nshell];
      }
    }
    asp = DynLinArr<AtomicSecondaryProducts>(qshell);
    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] = thresh[nshell] - thresh[n_min];
        asp[nshell].add_channel(fl[nshell], felectron_energy, fphoton_energy);
      }
    }
  }
 
  check_econd11(qshell, <= 0, mcerr);
  s_ignore_shell.put_qel(qshell, 0);
  height_of_excitation = 0.0;
  exener[0] = exener[1] = 0.0;
  double integ = get_integral_ACS(0.0, DBL_MAX);
  //Iprintn(mcout, integ);
  integ_abs_before_corr = integ;
  double pred_integ = Thomas_sum_rule_const_Mb * Z;
  //Iprintn(mcout, pred_integ);
  if (pred_integ > integ) {
    if (s_add_excitations_to_normalize == 1) {
      // add excitation
      exener[0] =
          low_boundary_of_excitations * acs[qshell - 1]->get_threshold();
      exener[1] = 1.0 * acs[qshell - 1]->get_threshold();
      height_of_excitation = (pred_integ - integ) / (exener[1] - exener[0]);
      if (minimal_threshold > 0.0) {
        if (minimal_threshold > acs[qshell - 1]->get_threshold()) {
          // currently the minimal shell is the last one
          exener[0] += minimal_threshold - acs[qshell - 1]->get_threshold();
          exener[1] += minimal_threshold - acs[qshell - 1]->get_threshold();
        }
      }
    }
  } else {
    if (s_scale_to_normalize_if_more == 1) {
      const double fact = pred_integ / integ;
      for (int nshell = 0; nshell < qshell; ++nshell) {
        acs[nshell]->scale(fact);
      }
    }
  }
  integ_abs_after_corr = get_integral_ACS(0.0, DBL_MAX);
  integ_ioniz_after_corr = get_integral_ICS(0.0, DBL_MAX);
}

ExAtomPhotoAbsCS() [4/5]

Heed::ExAtomPhotoAbsCS::ExAtomPhotoAbsCS	(	int	fZ,
		const String &	fname,
		const String &	fFitBT_file_name,
		int	id,
		int	s_no_scale,
		double	fminimal_threshold = 0.0
	)

Definition at line 1992 of file PhotoAbsCS.cpp.

    : threshold_file_name("none"),
      simple_table_file_name("none"),
      BT_file_name(fFitBT_file_name),
      minimal_threshold(fminimal_threshold) {
  mfunnamep(
      "ExAtomPhotoAbsCS::ExAtomPhotoAbsCS(int fZ, const String& fname, const "
      "String& fFitBT_file_name, int id, int id1, double fminimal_threshold)");
  check_econd11(fZ, < 1, mcerr);
  check_econd21(id, < 1 ||, > 2, mcerr);
  Z = fZ;
  name = fname;
#ifdef USE_STLSTRING
  std::ifstream BT_file(fFitBT_file_name.c_str());
#else
  std::ifstream BT_file(fFitBT_file_name);
#endif
  if (!BT_file) {
    funnw.ehdr(mcerr);
    mcerr << "cannot open file " << BT_file_name << std::endl;
    spexit(mcerr);
  }
  DynLinArr<double> thresh(0, 0.0);
  DynLinArr<double> fl(0, 0.0);
 
  int i = 0;
  while ((i = findmark(BT_file, "$")) == 1) {
    long iZ;
    BT_file >> iZ;
    //Iprintn(mcout, iZ);
    if (iZ == Z) {
      BT_file >> qshell;
      //Iprintn(mcout, qshell);
      check_econd11(qshell, <= 0, mcerr);
      check_econd11(qshell, > 1000, mcerr);
      acs.put_qel(qshell);
      if (id == 2) {
        thresh.put_qel(qshell);
        fl.put_qel(qshell);
      }
      for (int nshell = 0; nshell < qshell; ++nshell) {
#ifdef READ_FILE_WITH_PRINCIPAL_NUMBERS
        int n_princ = 0;
#endif
        int l;
        double threshold;
        double E0;
        double yw;
        double ya;
        double P;
        double sigma;
        if (BT_file.eof()) {
          mcerr << "unexpected end of file " << BT_file_name << '\n';
          spexit(mcerr);
        }
        if (!BT_file.good()) {
          mcerr << "bad format of file " << BT_file_name << '\n';
          spexit(mcerr);
        }
#ifdef READ_FILE_WITH_PRINCIPAL_NUMBERS
        BT_file >> n_princ;
        if (!BT_file.good()) {
          mcerr << "bad format of file " << BT_file_name << '\n';
          spexit(mcerr);
        }
        check_econd21(n_princ, < 0 ||, > 10, mcerr);
#endif
        BT_file >> l >> threshold >> E0 >> sigma >> ya >> P >> yw;
        check_econd11(l, < 0, mcerr);
        check_econd11(l, > 20, mcerr);
        threshold *= 1.0e-6;
        E0 *= 1.0e-6;
 
        check_econd11a(threshold, <= 2.0e-6,
                       "n_princ=" << n_princ << " l=" << l << '\n', mcerr);
        check_econd11(E0, <= 0, mcerr);
        double flu = 0.0;
        if (id == 2) {
          if (BT_file.eof()) {
            mcerr << "unexpected end of file " << BT_file_name << '\n';
            spexit(mcerr);
          }
          if (!BT_file.good()) {
            mcerr << "bad format of file " << BT_file_name << '\n';
            spexit(mcerr);
          }
          BT_file >> flu;
          check_econd11(flu, < 0.0, mcerr);
          check_econd11(flu, > 1.0, mcerr);
          thresh[nshell] = threshold;
          fl[nshell] = flu;
        }
#ifdef READ_FILE_WITH_PRINCIPAL_NUMBERS
        String shellname(long_to_String(n_princ) +  // necessary
                         // for generation escape products
                         String(" shell number ") + long_to_String(nshell));
#else
        String shellname(String("shell number ") + long_to_String(nshell));
#endif
        acs[nshell].pass(
            new SimpleTablePhotoAbsCS(shellname, 0,  // unknown here
                                      threshold, l, E0, yw, ya, P, sigma));
        //Iprintn(mcout, nshell);
        //Iprint3n(mcout, l, threshold, E0);
        //Iprint4n(mcout, yw, ya, P, sigma);
        //acs[nshell]->print(mcout, 5);
      }
      goto mark1;
    }
  }
  funnw.ehdr(mcerr);
  mcerr << "there is no element Z=" << fZ << " in file " << fFitBT_file_name
        << std::endl;
  spexit(mcerr);
mark1:
  if (id == 2) {
    // a copy of similar thing from subroutine above
    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 (thresh[nshell] < st) {
        n_min = nshell;
        st = thresh[nshell];
      }
    }
    asp = DynLinArr<AtomicSecondaryProducts>(qshell);
    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] = thresh[nshell] - thresh[n_min];
        asp[nshell].add_channel(fl[nshell], felectron_energy, fphoton_energy);
      }
    }
  }
 
  check_econd11(qshell, <= 0, mcerr);
  s_ignore_shell.put_qel(qshell, 0);
  height_of_excitation = 0.0;
  exener[0] = exener[1] = 0.0;
  double integ = get_integral_ACS(0.0, DBL_MAX);
  //Iprintn(mcout, integ);
  integ_abs_before_corr = integ;
  double pred_integ = Thomas_sum_rule_const_Mb * Z;
  //Iprintn(mcout, pred_integ);
  if (pred_integ > integ) {
    if (s_add_excitations_to_normalize == 1) {
      // add excitation
      exener[0] =
          low_boundary_of_excitations * acs[qshell - 1]->get_threshold();
      exener[1] = 1.0 * acs[qshell - 1]->get_threshold();
      height_of_excitation = (pred_integ - integ) / (exener[1] - exener[0]);
      if (minimal_threshold > 0.0) {
        if (minimal_threshold > acs[qshell - 1]->get_threshold()) {
          // currently the minimal shell is the last one
          exener[0] += minimal_threshold - acs[qshell - 1]->get_threshold();
          exener[1] += minimal_threshold - acs[qshell - 1]->get_threshold();
        }
      }
    }
  } else {
    if (s_scale_to_normalize_if_more == 1 && s_no_scale == 0) {
      const double fact = pred_integ / integ;
      for (int nshell = 0; nshell < qshell; ++nshell) {
        acs[nshell]->scale(fact);
      }
    }
  }
  integ_abs_after_corr = get_integral_ACS(0.0, DBL_MAX);
  integ_ioniz_after_corr = get_integral_ICS(0.0, DBL_MAX);
}

ExAtomPhotoAbsCS() [5/5]

Heed::ExAtomPhotoAbsCS::ExAtomPhotoAbsCS	(	int	fZ,
		const String &	fname,
		const String &	fFitBT_file_name,
		const String &	fsimple_table_file_name,
		double	emax_repl,
		int	id,
		double	fminimal_threshold = 0.0
	)

Definition at line 2171 of file PhotoAbsCS.cpp.

                               {
  mfunname("ExAtomPhotoAbsCS::ExAtomPhotoAbsCS(...)");
  Z = fZ;
  name = fname;
  int s_no_scale = 1;
  *this =
      ExAtomPhotoAbsCS(fZ, fname,  // just name of this atom
                       fFitBT_file_name, id, s_no_scale, fminimal_threshold);
 
  height_of_excitation = 0.0;
  exener[0] = exener[1] = 0.0;
 
  double thrmin = DBL_MAX;
  long nsmin = -1;
  // look for minimal shell (usually the last)
  for (long ns = 0; ns < qshell; ++ns) {
    //Iprint2n(mcout, ns, acs[ns]->get_threshold());
    if (thrmin > acs[ns]->get_threshold()) {
      nsmin = ns;
      thrmin = acs[ns]->get_threshold();
    }
  }
  //Iprint3n(mcout, nsmin, acs[nsmin]->get_threshold(), thrmin);
  check_econd11(nsmin, < 0, mcerr);
  check_econd11(nsmin, != qshell - 1, mcerr);  // now it has to be by this way
  ActivePtr<PhotoAbsCS> facs = acs[nsmin];     // copying the valence shell
  PhotoAbsCS* apacs = facs.get();
  //PhotoAbsCS* apacs = acs[nsmin]->getver();
  SimpleTablePhotoAbsCS* first_shell =
      dynamic_cast<SimpleTablePhotoAbsCS*>(apacs);
  //mcout<<"first_shell:\n";
  //first_shell->print(mcout, 2);
  // Strange why the following does not work (because of something being not
  // implemented in commpiler)? - may be due to an error before, which is now
  // corrected.
  // assumes the first is the last
  // SimpleTablePhotoAbsCS* first_shell =
  // dynamic_cast<SimpleTablePhotoAbsCS*>(&(acs[nsmin]->getver()));
 
  check_econd11(first_shell, == NULL, mcerr);
  /*
#ifdef USE_STLSTRING
  std::ifstream threshold_file(threshold_file_name.c_str());
#else
  std::ifstream threshold_file(threshold_file_name);
#endif
  if( !threshold_file )
  {
    funnw.ehdr(mcerr);
    mcerr<<"cannot open file "<<threshold_file_name<<std::endl;
    spexit(mcerr);
  }
  DynLinArr< double > thr(0, 0.0);
  DynLinArr< int > Zshell(0, 0);
  DynLinArr< double > fl(0, 0.0);
  DynLinArr< String > shell_name(0);
  int n_min=0;
  double thr_min = DBL_MAX;
  while (findmark(threshold_file, "#") == 1) {
    threshold_file >> Z;
    if (Z == fZ) {
      threshold_file >> qshell;
      check_econd21( qshell , < 1 || , > 10000 , mcerr);
      s_ignore_shell.put_qel(qshell, 0);
      //Iprintn(mcout, qshell);
      thr = DynLinArr< double >(qshell, 0.0);
      Zshell = DynLinArr< int >(qshell, 0);
      fl = DynLinArr< double >(qshell, 0.0);
      shell_name = DynLinArr< String >(qshell);
      String temp_name;
      threshold_file >> temp_name;
      if(fname == "none") name = temp_name;
      else name = fname;
      int sZshell = 0;
      for (int nshell = 0; nshell < qshell; nshell++) {
        threshold_file >> thr[nshell];
        check_econd11(thr[nshell] , <= 0.0, mcerr);
        thr[nshell] *= 1.0e-6;
        threshold_file >> Zshell[nshell];
        check_econd11(Zshell[nshell] , <= 0, mcerr);
        sZshell += Zshell[nshell];
        threshold_file >> fl[nshell];
        findmark(threshold_file, "!");
        threshold_file >> shell_name[nshell];
        //PhotoAbsCS* ap = new PhenoPhotoAbsCS(shell_name, Zshell, thr);
        //acs[nshell].pass( ap );
      }
      check_econd12( sZshell , != , Z , mcerr);
      for(nshell=0; nshell<qshell; nshell++) {
        // currently the minimal shell is the last,
        // but to avoid this assumption we check all
        if (thr[nshell] < thr_min) {
          n_min = nshell;
          thr_min = thr[nshell];
        }
      }
      goto mark1;
    }
  }
  funnw.ehdr(mcerr);
  mcerr<<"there is no element Z="<<fZ<<" in file
"<<threshold_file_name<<std::endl;
  spexit(mcerr);
  mark1:
  */
 
  SimpleTablePhotoAbsCS stpacs(name, Z, 0.0, fsimple_table_file_name);
  stpacs.remove_leading_tiny(1.0e-10);
 
  // Merging shells:
  acs[nsmin].pass(new SimpleTablePhotoAbsCS(*first_shell, stpacs, emax_repl));
 
  s_ignore_shell.put_qel(qshell, 0);
  height_of_excitation = 0.0;
  exener[0] = exener[1] = 0.0;
  double integ = get_integral_ACS(0.0, DBL_MAX);
  //Iprintn(mcout, integ);
  integ_abs_before_corr = integ;
  double pred_integ = Thomas_sum_rule_const_Mb * Z;
  //Iprintn(mcout, pred_integ);
  if (pred_integ > integ) {
    if (s_add_excitations_to_normalize == 1) {
      // add excitation
      exener[0] =
          low_boundary_of_excitations * acs[qshell - 1]->get_threshold();
      exener[1] = 1.0 * acs[qshell - 1]->get_threshold();
      height_of_excitation = (pred_integ - integ) / (exener[1] - exener[0]);
      if (minimal_threshold > 0.0) {
        if (minimal_threshold > acs[qshell - 1]->get_threshold()) {
          // currently the minimal shell is the last one
          exener[0] += minimal_threshold - acs[qshell - 1]->get_threshold();
          exener[1] += minimal_threshold - acs[qshell - 1]->get_threshold();
        }
      }
    }
  } else {
    if (s_scale_to_normalize_if_more == 1) {
      const double fact = pred_integ / integ;
      for (int nshell = 0; nshell < qshell; ++nshell) {
        acs[nshell]->scale(fact);
      }
    }
  }
  integ_abs_after_corr = get_integral_ACS(0.0, DBL_MAX);
  integ_ioniz_after_corr = get_integral_ICS(0.0, DBL_MAX);
}

~ExAtomPhotoAbsCS()

virtual Heed::ExAtomPhotoAbsCS::~ExAtomPhotoAbsCS ( )

inlinevirtual

Definition at line 511 of file PhotoAbsCS.h.

{}

Member Function Documentation

get_ACS() [1/2]

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

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 2400 of file PhotoAbsCS.cpp.

                                                    {
  mfunname("double ExAtomPhotoAbsCS::get_ACS(double energy) const");
  double s = 0.0;
  for (int n = 0; n < qshell; ++n) {
    if (s_ignore_shell[n] == 0) {
      double shift = 0.0;
      const double t = acs[n]->get_threshold();
      if (minimal_threshold > 0.0) {
        if (t < minimal_threshold) shift = minimal_threshold - t;
      }
      s += acs[n]->get_CS(energy - shift);
    }
  }
  if (energy >= exener[0] && energy <= exener[1]) {
    s += height_of_excitation;
  }
  return s;
}

get_ACS() [2/2]

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

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 2448 of file PhotoAbsCS.cpp.

                                                                {
  mfunname("double ExAtomPhotoAbsCS::get_ACS(int nshell, double energy)");
  check_econd21(nshell, < 0 ||, > qshell, mcerr);
  if (s_ignore_shell[nshell] == 0) {
    double shift = 0.0;
    const double t = acs[nshell]->get_threshold();
    if (minimal_threshold > 0.0) {
      if (t < minimal_threshold) shift = minimal_threshold - t;
    }
    double s = acs[nshell]->get_CS(energy - shift);
    if (nshell == qshell - 1 && energy >= exener[0] && energy <= exener[1]) {
      s += height_of_excitation;
    }
    return s;
  }
  return 0.0;
}

get_ICS() [1/2]

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

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 2335 of file PhotoAbsCS.cpp.

                                                    {
  mfunname("double ExAtomPhotoAbsCS::get_ACS(double energy) const");
  double s = 0.0;
  for (int n = 0; n < qshell; ++n) {
    if (s_ignore_shell[n] == 0) {
      double shift = 0.0;
      const double t = acs[n]->get_threshold();
      if (minimal_threshold > 0.0) {
        if (t < minimal_threshold) shift = minimal_threshold - t;
      }
      //Iprint3n(mcout, n, t, shift);
      s += acs[n]->get_CS(energy - shift);
      //Iprintn(mcout, s);
    }
  }
  return s;
}

get_ICS() [2/2]

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

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 2370 of file PhotoAbsCS.cpp.

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

get_integral_ACS() [1/2]

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

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 2419 of file PhotoAbsCS.cpp.

                                                                {
  mfunname("double ExAtomPhotoAbsCS::get_integral_ACS(double energy1, double "
           "energy2) const");
  double s = 0.0;
  for (int n = 0; n < qshell; ++n) {
    if (s_ignore_shell[n] == 0) {
      double shift = 0.0;
      const double t = acs[n]->get_threshold();
      if (minimal_threshold > 0.0) {
        if (t < minimal_threshold) shift = minimal_threshold - t;
      }
      s += acs[n]->get_integral_CS(energy1 - shift, energy2 - shift);
    }
  }
  //Imcout<<"energy1="<<setw(13)<<energy1
  //        <<"energy2="<<setw(13)<<energy2
  //        <<" s="<<s<<'\n';
  double b[2];
  b[0] = find_max(exener[0], energy1);
  b[1] = find_min(exener[1], energy2);
  //Iprint3n(mcout, b[0], b[1], height_of_excitation);
  if (b[1] >= b[0]) {
    s += height_of_excitation * (b[1] - b[0]);
  }
  //Iprintn(mcout, s);
  return s;
}

Referenced by ExAtomPhotoAbsCS(), and print().

get_integral_ACS() [2/2]

double Heed::ExAtomPhotoAbsCS::get_integral_ACS ( int  nshell, double  energy1, double  energy2  ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 2466 of file PhotoAbsCS.cpp.

                                                                {
  mfunname("double ExAtomPhotoAbsCS::get_integral_ACS(int nshell, double "
           "energy1, double energy2)");
  check_econd21(nshell, < 0 ||, > qshell, mcerr);
  if (s_ignore_shell[nshell] == 0) {
    double shift = 0.0;
    const double t = acs[nshell]->get_threshold();
    if (minimal_threshold > 0.0) {
      if (t < minimal_threshold) shift = minimal_threshold - t;
    }
    double s = acs[nshell]->get_integral_CS(energy1 - shift, energy2 - shift);
    //Imcout<<"energy1="<<setw(13)<<energy1<<" s="<<s<<'\n';
    if (nshell == qshell - 1) {
      double b[2];
      b[0] = find_max(exener[0], energy1);
      b[1] = find_min(exener[1], energy2);
      if (b[1] >= b[0]) {
        s += height_of_excitation * (b[1] - b[0]);
      }
    }
    return s;
  }
  return 0.0;
}

get_integral_ICS() [1/2]

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

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 2353 of file PhotoAbsCS.cpp.

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

Referenced by ExAtomPhotoAbsCS(), and print().

get_integral_ICS() [2/2]

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

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 2384 of file PhotoAbsCS.cpp.

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

get_main_shell_number()

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

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 2531 of file PhotoAbsCS.cpp.

                                                            {
  mfunname("int ExAtomPhotoAbsCS::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 = -1;
  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::ExAtomPhotoAbsCS::get_threshold ( int  nshell ) const

virtual

Implements Heed::AtomPhotoAbsCS.

Definition at line 2325 of file PhotoAbsCS.cpp.

                                                       {
  mfunname("double ExAtomPhotoAbsCS::get_threshold(int nshell) const");
  check_econd21(nshell, < 0 ||, > qshell, mcerr);
  double r = acs[nshell]->get_threshold();
  if (minimal_threshold > 0.0) {
    if (r < minimal_threshold) r = minimal_threshold;
  }
  return r;
}

Referenced by ExAtomPhotoAbsCS().

macro_copy_total()

Heed::ExAtomPhotoAbsCS::macro_copy_total

(

ExAtomPhotoAbsCS

)

print()

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

virtual

Reimplemented from Heed::AtomPhotoAbsCS.

Definition at line 2492 of file PhotoAbsCS.cpp.

                                                          {
  if (l > 0) {
    Ifile << "ExAtomPhotoAbsCS(l=" << l << "): name=" << name << " Z = " << Z
          << " qshell = " << qshell << std::endl;
    indn.n += 2;
    Ifile << "threshold_file_name=" << threshold_file_name << '\n';
    Ifile << "simple_table_file_name=" << simple_table_file_name << '\n';
    Ifile << "BT_file_name=" << BT_file_name << std::endl;
    Ifile << "Thomas_sum_rule_const_Mb * Z = " << Thomas_sum_rule_const_Mb* Z
          << '\n';
    Ifile << "integ_abs_before_corr        = " << integ_abs_before_corr << '\n';
    Ifile << "integ_abs_after_corr         = " << integ_abs_after_corr << '\n';
    Ifile << "integ_ioniz_after_corr       = " << integ_ioniz_after_corr
          << '\n';
    Ifile << "height_of_excitation=" << height_of_excitation
          << " exener=" << exener[0] << ' ' << exener[1] << '\n';
    Iprintn(file, minimal_threshold);
    Ifile << "integrals by shells:\n";
    Ifile << "nshell, int(acs), int(ics)\n";
    for (long n = 0; n < qshell; n++) {
      double ainteg = get_integral_ACS(n, 0.0, DBL_MAX);
      double iinteg = get_integral_ICS(n, 0.0, DBL_MAX);
      Ifile << n << "    " << ainteg << "    " << iinteg << '\n';
    }
 
    if (l > 1) {
      l--;
      indn.n += 2;
      for (long n = 0; n < qshell; ++n) {
        Ifile << "nshell=" << n << std::endl;
        acs[n].print(file, l);
      }
      AtomPhotoAbsCS::print(file, l);
      indn.n -= 2;
    }
    indn.n -= 2;
  }
}

replace_shells_by_overage()

void Heed::ExAtomPhotoAbsCS::replace_shells_by_overage	(	double	fwidth,
		double	fstep,
		long	fmax_q_step
	)

Definition at line 2554 of file PhotoAbsCS.cpp.

                                                                   {
  mfunname("void ExAtomPhotoAbsCS::replace_shells_by_overage(...)");
  for (long n = 0; n < qshell; n++) {
    //Iprintn(mcout, n);
    //mcout<<"----------------before replacement:\n";
    //acs[n]->print(mcout, 10);
    PhotoAbsCS* a =
        new OveragePhotoAbsCS(acs[n].getver(), fwidth, fstep, fmax_q_step);
    //mcout<<"----------------after replacement:\n";
    //acs[n]->print(mcout, 10);
    acs[n].pass(a);
  }
}

Member Data Documentation

acs

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

protected

Definition at line 517 of file PhotoAbsCS.h.

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

BT_file_name

String Heed::ExAtomPhotoAbsCS::BT_file_name

protected

Definition at line 516 of file PhotoAbsCS.h.

Referenced by ExAtomPhotoAbsCS(), and print().

exener

double Heed::ExAtomPhotoAbsCS::exener[2]

protected

Definition at line 529 of file PhotoAbsCS.h.

Referenced by ExAtomPhotoAbsCS(), get_ACS(), get_integral_ACS(), and print().

height_of_excitation

double Heed::ExAtomPhotoAbsCS::height_of_excitation

protected

Definition at line 528 of file PhotoAbsCS.h.

Referenced by ExAtomPhotoAbsCS(), get_ACS(), get_integral_ACS(), and print().

integ_abs_after_corr

double Heed::ExAtomPhotoAbsCS::integ_abs_after_corr

protected

Definition at line 526 of file PhotoAbsCS.h.

Referenced by ExAtomPhotoAbsCS(), and print().

integ_abs_before_corr

double Heed::ExAtomPhotoAbsCS::integ_abs_before_corr

protected

Definition at line 525 of file PhotoAbsCS.h.

Referenced by ExAtomPhotoAbsCS(), and print().

integ_ioniz_after_corr

double Heed::ExAtomPhotoAbsCS::integ_ioniz_after_corr

protected

Definition at line 527 of file PhotoAbsCS.h.

Referenced by ExAtomPhotoAbsCS(), and print().

minimal_threshold

double Heed::ExAtomPhotoAbsCS::minimal_threshold

protected

Definition at line 531 of file PhotoAbsCS.h.

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

simple_table_file_name

String Heed::ExAtomPhotoAbsCS::simple_table_file_name

protected

Definition at line 515 of file PhotoAbsCS.h.

Referenced by print().

threshold_file_name

String Heed::ExAtomPhotoAbsCS::threshold_file_name

protected

Definition at line 514 of file PhotoAbsCS.h.

Referenced by ExAtomPhotoAbsCS(), and print().

---

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

• PhotoAbsCS.h
• PhotoAbsCS.cpp