Heed::EnTransfCS Class Reference

#include <EnTransfCS.h>

Inheritance diagram for Heed::EnTransfCS:

Public Member Functions
	EnTransfCS ()=default
	Default constructor.
	EnTransfCS (double fparticle_mass, double fgamma_1, bool fs_primary_electron, HeedMatterDef *fhmd, long fparticle_charge=1)
	Constructor.
void	print (std::ostream &file, int l) const
EnTransfCS *	copy () const

Public Attributes
double	particle_mass = 0.
	Particle mass [MeV].
long	particle_charge = 0
	Charge in units of electron charge (used square, sign does not matter).
double	gamma_1 = 0.
	Lorentz factor - 1 (the best dimensionless measurement of speed).
double	max_etransf = 0.
	Max. energy transfer [MeV].
bool	s_simple_form = true
bool	s_primary_electron = false
	Flag indicating whether the primary particle is an electron.
HeedMatterDef *	hmd = nullptr
std::vector< double >	log1C
	common first log without cs
std::vector< double >	log2C
	common second log without cs
std::vector< double >	chereC
	Cherenkov's radiation.
std::vector< double >	chereCangle
	angle of Cherenkov's radiation
std::vector< double >	Rruth
	term called R in my paper
std::vector< double >	addaC
	Sum of (ionization) differential cross-section terms.
double	quanC = 0.
	Integrated (ionization) cross-section.
double	meanC = 0.
double	meanC1 = 0.
std::vector< std::vector< std::vector< double > > >	cher
std::vector< std::vector< std::vector< double > > >	fruth
	Rutherford term.
std::vector< std::vector< std::vector< double > > >	adda
	Sum.
std::vector< std::vector< std::vector< double > > >	fadda
	Integral, normalised to unity.
std::vector< std::vector< double > >	quan
	Number of collisions / cm, for each atom and shell.
std::vector< std::vector< double > >	mean
	First moment, for each atom and shell.
std::vector< double >	length_y0

Detailed Description

The PAI cross section of energy transfers from charged particle to media. The particle has fixed parameters (energy, speed, etc.), which are not affected by energy transfers, since they are considered too small compared with the particle energy.

2003, I. Smirnov

Definition at line 17 of file EnTransfCS.h.

Constructor & Destructor Documentation

EnTransfCS() [1/2]

Heed::EnTransfCS::EnTransfCS ( )

default

Default constructor.

Referenced by copy().

EnTransfCS() [2/2]

Heed::EnTransfCS::EnTransfCS	(	double	fparticle_mass,
		double	fgamma_1,
		bool	fs_primary_electron,
		HeedMatterDef *	fhmd,
		long	fparticle_charge = 1
	)

Constructor.

Definition at line 112 of file EnTransfCS.cpp.

    : particle_mass(fparticle_mass),
      particle_charge(fparticle_charge),
      gamma_1(fgamma_1),
      s_primary_electron(fs_primary_electron),
      hmd(fhmd) {
  mfunnamep("EnTransfCS::EnTransfCS(...)");
 
  const double beta = lorbeta(fgamma_1);
  const double beta2 = beta * beta;
  const double beta12 = 1.0 - beta2;
  const double gamma = fgamma_1 + 1.;
  // Particle kinetic energy.
  const double tkin = particle_mass * gamma_1;
  const double ener = particle_mass * gamma;
  // Calculate the max. energy transfer.
  if (s_primary_electron) {
    max_etransf = 0.5 * tkin;
  } else {
    double rm2 = particle_mass * particle_mass;
    double rme = electron_mass_c2;
    if (beta12 > 1.0e-10) {
      max_etransf =
          2.0 * rm2 * electron_mass_c2 * beta2 /
          ((rm2 + rme * rme + 2.0 * rme * gamma * particle_mass) * beta12);
      if (max_etransf > tkin) max_etransf = tkin;
    } else {
      max_etransf = tkin;
    }
  }
  const long qe = hmd->energy_mesh->get_q();
  log1C.assign(qe, 0.0);
  log2C.assign(qe, 0.0);
  chereC.assign(qe, 0.0);
  chereCangle.assign(qe, 0.0);
  Rruth.assign(qe, 0.0);
  addaC.assign(qe, 0.0);
#ifndef EXCLUDE_A_VALUES
  addaC_a.assign(qe, 0.0);
#endif
 
  const long qa = hmd->matter->qatom();
  cher.resize(qa);
  fruth.resize(qa);
  adda.resize(qa);
  fadda.resize(qa);
  quan.resize(qa);
  mean.resize(qa);
#ifndef EXCLUDE_A_VALUES
  cher_a.resize(qa);
  adda_a.resize(qa);
  fadda_a.resize(qa);
  quan_a.resize(qa);
  mean_a.resize(qa);
#endif
 
  for (long na = 0; na < qa; na++) {
    const long qs = hmd->apacs[na]->get_qshell();
    cher[na].assign(qs, std::vector<double>(qe, 0.));
    fruth[na].assign(qs, std::vector<double>(qe, 0.));
    adda[na].assign(qs, std::vector<double>(qe, 0.));
    fadda[na].assign(qs, std::vector<double>(qe, 0.));
    quan[na].assign(qs, 0.);
    mean[na].assign(qs, 0.);
#ifndef EXCLUDE_A_VALUES
    cher_a[na].assign(qs, std::vector<double>(qe, 0.));
    adda_a[na].assign(qs, std::vector<double>(qe, 0.));
    fadda_a[na].assign(qs, std::vector<double>(qe, 0.));
    quan_a[na].assign(qs, 0.);
    mean_a[na].assign(qs, 0.);
#endif
  }
 
  const long q2 = particle_charge * particle_charge;
  double coefpa = fine_structure_const * q2 / (beta2 * CLHEP::pi);
  const double coefCh = coefpa / hmd->eldens;
  for (long ne = 0; ne < qe; ne++) {
    const double eps1 = hmd->epsi1[ne];
    const double eps2 = hmd->epsi2[ne];
    const double a0 = 1. + eps1;
    const double a1 = -eps1 + a0 * beta12;
    const double a2 = beta2 * eps2;
    log1C[ne] = log(1. / sqrt(a1 * a1 + a2 * a2));
 
    const double ec = hmd->energy_mesh->get_ec(ne);
    const double a3 = 2. * electron_mass_c2 * beta2 / ec;
    log2C[ne] = a3 > 0. ? log(a3) : 0.;
 
    double a4 = atan(a2 / a1);
    if (a1 < 0) a4 += CLHEP::pi;
    chereCangle[ne] = a4;
 
    const double a5 = eps2 * eps2;
    const double a6 = (-a0 * a1 + beta2 * a5) / (a0 * a0 + a5);
    chereC[ne] = coefCh * a6 * a4;
 
    if (s_simple_form) {
      Rruth[ne] = 1. / (ec * ec);
      if (!s_primary_electron) {
        Rruth[ne] *= (1. - beta2 * ec / max_etransf);
      }
    } else {
      if (!s_primary_electron) {
        Rruth[ne] = 1. / (ec * ec) * (1. - beta2 * ec / max_etransf +
                                      ec * ec / (2. * ener * ener));
      } else {
        const double delta = ec / particle_mass;
        const double pg2 = gamma * gamma;
        const double dgd = delta * (gamma_1 - delta);
        Rruth[ne] = beta2 / (particle_mass * particle_mass) * 1.0 /
                    (pg2 - 1.0) * (gamma_1 * gamma_1 * pg2 / (dgd * dgd) -
                                   (2.0 * pg2 + 2.0 * gamma - 1.0) / dgd + 1.0);
      }
    }
  }
 
  double Z_mean = hmd->matter->Z_mean();
  const double ethr = hmd->min_ioniz_pot;
  for (long na = 0; na < qa; na++) {
    auto pacs = hmd->apacs[na];
    const double awq = hmd->matter->weight_quan(na);
    const long qs = pacs->get_qshell();
    for (long ns = 0; ns < qs; ns++) {
      for (long ne = 0; ne < qe; ne++) {
        double e1 = hmd->energy_mesh->get_e(ne);
        double e2 = hmd->energy_mesh->get_e(ne + 1);
        double ics = 0.;
        if (hmd->s_use_mixture_thresholds == 1) {
          ics = pacs->get_integral_TICS(ns, e1, e2, ethr) / (e2 - e1);
        } else {
          ics = pacs->get_integral_ICS(ns, e1, e2) / (e2 - e1);
        }
        check_econd11a(ics, < 0,
                       "na=" << na << " ns=" << ns << " ne=" << ne << '\n',
                       mcerr);
        const double tacs = hmd->ACS[ne];
        if (tacs <= 0.0) continue;
        cher[na][ns][ne] = chereC[ne] * awq * ics / tacs;
#ifndef EXCLUDE_A_VALUES
        double acs = pacs->get_integral_ACS(ns, e1, e2) / (e2 - e1);
        cher_a[na][ns][ne] = chereC[ne] * awq * acs / tacs;
#endif
      }
      // Calculate the integral.
      const double cR = C1_MEV2_MBN * awq * coefpa / Z_mean;
      double s = 0.;
      for (long ne = 0; ne < qe; ne++) {
        const double e1 = hmd->energy_mesh->get_e(ne);
        const double ec = hmd->energy_mesh->get_ec(ne);
        const double e2 = hmd->energy_mesh->get_e(ne + 1);
        const double r = pacs->get_integral_ACS(ns, e1, e2) * cR;
        // Here it must be ACS to satisfy sum rule for Rutherford
        check_econd11a(r, < 0.0, "na=" << na << " ns=" << ns << " ne=" << ne,
                       mcerr);
        if (ec > ethr && ec < max_etransf) {
          fruth[na][ns][ne] = (s + 0.5 * r) * Rruth[ne];
          check_econd11a(fruth[na][ns][ne], < 0,
                         "na=" << na << " ns=" << ns << " na=" << na, mcerr);
        } 
        s += r;
      }
    }
  }
  for (long ne = 0; ne < qe; ++ne) {
    double s = 0.0;
#ifndef EXCLUDE_A_VALUES
    double s_a = 0.0;
#endif
    double e1 = hmd->energy_mesh->get_e(ne);
    double ec = hmd->energy_mesh->get_ec(ne);
    double e2 = hmd->energy_mesh->get_e(ne + 1);
    const double eps1 = hmd->epsi1[ne];
    const double eps2 = hmd->epsi2[ne];
    const double eps11 = 1. + eps1;
    const double sqepsi = eps11 * eps11 + eps2 * eps2;
    const double cL = C1_MEV2_MBN * coefpa * (log1C[ne] + log2C[ne]) / 
                      (ec * Z_mean * sqepsi); 
    for (long na = 0; na < qa; na++) {
      double awq = hmd->matter->weight_quan(na);
      auto pacs = hmd->apacs[na];
      const long qs = pacs->get_qshell();
      for (long ns = 0; ns < qs; ns++) {
        double ics = 0.;
        if (hmd->s_use_mixture_thresholds == 1) {
          ics = pacs->get_integral_TICS(ns, e1, e2, ethr) / (e2 - e1);
        } else {
          ics = pacs->get_integral_ICS(ns, e1, e2) / (e2 - e1);
        }
        double r = std::max(cL * awq * ics + fruth[na][ns][ne], 0.);
#ifndef EXCLUDE_A_VALUES
        double acs = pacs->get_integral_ACS(ns, e1, e2) / (e2 - e1);
        double r_a = std::max(cL * awq * acs + fruth[na][ns][ne], 0.);
#endif
        if (ec > ethr) {
          r += cher[na][ns][ne];
          if (r < 0.0) {
            funnw.whdr(mcout);
            mcout << "negative adda\n";
            mcout << "na=" << na << " ns=" << ns << " ne=" << ne
                  << ": " << r << '\n';
            r = 0.;
          }
        }
#ifndef EXCLUDE_A_VALUES
        r_a += cher[na][ns][ne];
        check_econd11a(r_a, < 0,
                       "na=" << na << " ns=" << ns << " na=" << na, mcerr);
#endif
        adda[na][ns][ne] = r;
        s += r;
#ifndef EXCLUDE_A_VALUES
        adda_a[na][ns][ne] = r_a;
        s_a += r_a;
#endif
      }
    }
    addaC[ne] = s;
#ifndef EXCLUDE_A_VALUES
    addaC_a[ne] = s_a;
#endif
  }
 
  const double* aetemp = hmd->energy_mesh->get_ae();
  PointCoorMesh<double, const double*> pcm_e(qe + 1, &(aetemp));
  double emin = hmd->energy_mesh->get_emin();
  double emax = hmd->energy_mesh->get_emax();
 
  const double rho = hmd->xeldens;
  quanC = integrate(pcm_e, addaC, emin, emax, 0) * rho;
  meanC = integrate(pcm_e, addaC, emin, emax, 1) * rho;
 
#ifndef EXCLUDE_A_VALUES
  quanC_a = integrate(pcm_e, addaC_a, emin, emax, 0) * rho;
  meanC_a = integrate(pcm_e, addaC_a, emin, emax, 1) * rho;
#endif
  meanC1 = meanC;
  const double coef = fine_structure_const * fine_structure_const * q2 * twopi /
                      (electron_mass_c2 * beta2) * rho;
  if (s_simple_form) {
    if (!s_primary_electron) {
      if (max_etransf > hmd->energy_mesh->get_e(qe)) {
        double e1 = hmd->energy_mesh->get_e(qe);
        double e2 = max_etransf;
        meanC1 += coef * (log(e2 / e1) - beta2 / max_etransf * (e2 - e1));
      }
    } else {
      if (max_etransf > hmd->energy_mesh->get_e(qe)) {
        double e1 = hmd->energy_mesh->get_e(qe);
        double e2 = max_etransf;
        meanC1 += coef * log(e2 / e1);
      }
    }
  } else {
    if (!s_primary_electron) {
      if (max_etransf > hmd->energy_mesh->get_e(qe)) {
        double e1 = hmd->energy_mesh->get_e(qe);
        double e2 = max_etransf;
        meanC1 += coef *
                  (log(e2 / e1) - beta2 / max_etransf * (e2 - e1) +
                   (e2 * e2 - e1 * e1) / (4.0 * ener * ener));
      }
#ifndef EXCLUDE_A_VALUES
      meanC1_a = meanC_a;
      if (max_etransf > hmd->energy_mesh->get_e(qe)) {
        double e1 = hmd->energy_mesh->get_e(qe);
        double e2 = max_etransf;
        meanC1_a += coef * (log(e2 / e1) - beta2 / max_etransf * (e2 - e1) +
                            (e2 * e2 - e1 * e1) /
                                (4.0 * ener * ener));
      }
#endif
    }
  }
 
  for (long na = 0; na < qa; na++) {
    const long qs = hmd->apacs[na]->get_qshell();
    for (long ns = 0; ns < qs; ns++) {
      quan[na][ns] = integrate(pcm_e, adda[na][ns], emin, emax, 0) * rho;
      mean[na][ns] = integrate(pcm_e, adda[na][ns], emin, emax, 1) * rho;
#ifndef EXCLUDE_A_VALUES
      quan_a[na][ns] = integrate(pcm_e, adda_a[na][ns], emin, emax, 0) * rho;
      mean_a[na][ns] = integrate(pcm_e, adda_a[na][ns], emin, emax, 1) * rho;
#endif
    }
  }
 
  for (long na = 0; na < qa; na++) {
    const long qs = hmd->apacs[na]->get_qshell();
    for (long ns = 0; ns < qs; ns++) {
      if (quan[na][ns] > 0.0) {
        const double s = cdf(pcm_e, adda[na][ns], fadda[na][ns]);
        if (fabs(s * rho - quan[na][ns]) > 1.e-10) {
          std::cerr << "Heed::EnTransfCS: Integrals differ (warning).\n";
        }
      }
#ifndef EXCLUDE_A_VALUES
      if (quan_a[na][ns] > 0.0) {
        const double s = cdf(pcm_e, adda_a[na][ns], fadda_a[na][ns]);
        if (fabs(s * rho - quan_a[na][ns]) > 1.e-10) {
          std::cerr << "Heed::EnTransfCS: Integrals differ (warning).\n";
        }
      }
#endif
    }
  }
 
  length_y0.resize(qe, 0.);
  for (long ne = 0; ne < qe; ne++) {
    const double k0 = hmd->energy_mesh->get_ec(ne) / (hbarc / cm);
    const double det_value = 1.0 / (gamma * gamma) - hmd->epsi1[ne] * beta2;
    length_y0[ne] = det_value > 0. ? beta / k0 * 1.0 / sqrt(det_value) : 0.;
  }
 
  /*
  std::ofstream dcsfile;
  dcsfile.open("dcs.txt", std::ios::out);
  dcsfile << "# energy [MeV] vs. diff. cs per electron [Mbarn / MeV]\n";
  for (int i = 0; i < qe; ++i) {
    double sumR = 0.;
    double sumC = 0.;
    double sumL = 0.; 
    for (long na = 0; na < qa; ++na) {
      const long qs = hmd->apacs[na]->get_qshell();
      for (long ns = 0; ns < qs; ++ns) {
        sumR += fruth[na][ns][i];
        sumC += cher[na][ns][i];
        sumL += adda[na][ns][i] - cher[na][ns][i] - fruth[na][ns][i];
      }
    }
    const double f1 = log1C[i] / (log1C[i] + log2C[i]);
    const double f2 = 1. - f1;
    sumR /= C1_MEV2_MBN;
    sumC /= C1_MEV2_MBN;
    sumL /= C1_MEV2_MBN;
    const double sumL1 = f1 * sumL;
    const double sumL2 = f2 * sumL;
    dcsfile << hmd->energy_mesh->get_ec(i) << "  " << addaC[i] / C1_MEV2_MBN
            << "  " << sumL1 << "  " << "  " << sumC << "  " << sumL2 << "  " << sumR 
            << "\n";
  }
  dcsfile.close();
  //*/
 
  log1C.clear();
  log2C.clear();
  chereC.clear();
  chereCangle.clear();
  Rruth.clear();
 
  addaC.clear();
  cher.clear();
  fruth.clear();
  adda.clear();
  mean.clear();
#ifndef EXCLUDE_A_VALUES
  addaC_a.clear();
  cher_a.clear();
  adda_a.clear();
  fadda_a.clear();
  mean_a.clear();
#endif
}

Member Function Documentation

copy()

EnTransfCS * Heed::EnTransfCS::copy ( ) const

inline

Definition at line 26 of file EnTransfCS.h.

{ return new EnTransfCS(*this); }

print()

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

Definition at line 477 of file EnTransfCS.cpp.

                                                    {
  if (l <= 0) return;
  Ifile << "EnTransfCS(l=" << l << "):\n";
  indn.n += 2;
  Ifile << "particle_mass=" << particle_mass
        << "particle_ener=" << particle_mass * (gamma_1 + 1.)
        << " particle_charge=" << particle_charge << std::endl;
  Ifile << "max_etransf=" << max_etransf << std::endl;
  Ifile << "s_primary_electron=" << s_primary_electron << std::endl;
  Ifile << "hmd:\n";
  hmd->print(file, 1);
#ifndef EXCLUDE_A_VALUES
  Ifile << "quanC=" << quanC << " quanC_a=" << quanC_a << '\n';
  Ifile << "meanC=" << meanC << " meanC_a=" << meanC_a << '\n';
  Ifile << "meanC1=" << meanC1 << " meanC1_a=" << meanC1_a << '\n';
#else
  Ifile << "quanC=" << quanC << '\n';
  Ifile << "meanC=" << meanC << '\n';
  Ifile << "meanC1=" << meanC1 << '\n';
#endif
  if (l > 2) {
    long qe = hmd->energy_mesh->get_q();
    long ne;
    if (l > 4) {
      Ifile << "       enerc,      log1C,      log2C,      chereC,     addaC, "
               "chereCangle   Rruth   length_y0\n";
      for (ne = 0; ne < qe; ne++) {
        Ifile << std::setw(12) << hmd->energy_mesh->get_ec(ne) << std::setw(12)
              << log1C[ne] << std::setw(12) << log2C[ne] << std::setw(12)
              << chereC[ne] << std::setw(12) << addaC[ne] << std::setw(12)
              << chereCangle[ne] << std::setw(12) << Rruth[ne]
              << std::setw(12) << length_y0[ne] << '\n';
      }
    }
    if (l > 3) {
      long qa = hmd->matter->qatom();
      long na;
      Iprintn(file, hmd->matter->qatom());
      for (na = 0; na < qa; na++) {
        Iprintn(file, na);
        long qs = hmd->apacs[na]->get_qshell();
        long ns;
        Iprintn(file, hmd->apacs[na]->get_qshell());
        for (ns = 0; ns < qs; ns++) {
          Iprintn(file, ns);
          Ifile << "quan      =" << std::setw(13) << quan[na][ns]
                << " mean  =" << std::setw(13) << mean[na][ns] << '\n';
#ifndef EXCLUDE_A_VALUES
          Ifile << "quan_a    =" << std::setw(13) << quan_a[na][ns]
                << " mean_a=" << std::setw(13) << mean_a[na][ns] << '\n';
#endif
          if (l > 5) {
            Ifile << "   enerc,        cher,       cher_a,     fruth,   adda, "
                     "  adda_a,  fadda,   fadda_a\n";
            for (ne = 0; ne < qe; ne++) {
              Ifile << std::setw(12)
                    << hmd->energy_mesh->get_ec(ne)
                    << std::setw(12) << cher[na][ns][ne]
#ifndef EXCLUDE_A_VALUES
                    << std::setw(12) << cher_a[na][ns][ne]
#endif
                    << std::setw(12) << fruth[na][ns][ne] << std::setw(12)
                    << adda[na][ns][ne]
#ifndef EXCLUDE_A_VALUES
                    << std::setw(12) << adda_a[na][ns][ne]
#endif
                    << std::setw(12) << fadda[na][ns][ne]
#ifndef EXCLUDE_A_VALUES
                    << std::setw(12) << fadda_a[na][ns][ne]
#endif
                    << '\n';
            }
          }
        }
      }
    }
  }
  indn.n -= 2;
}

Member Data Documentation

adda

std::vector<std::vector<std::vector<double> > > Heed::EnTransfCS::adda

Sum.

Definition at line 86 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

addaC

std::vector<double> Heed::EnTransfCS::addaC

Sum of (ionization) differential cross-section terms.

Definition at line 58 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

cher

std::vector<std::vector<std::vector<double> > > Heed::EnTransfCS::cher

In the following arrays there are three indices: atom number in the matter, shell number in atom, energy Fraction of Cherenkov term.

Definition at line 82 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

chereC

std::vector<double> Heed::EnTransfCS::chereC

Cherenkov's radiation.

Definition at line 53 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

chereCangle

std::vector<double> Heed::EnTransfCS::chereCangle

angle of Cherenkov's radiation

Definition at line 54 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

fadda

std::vector<std::vector<std::vector<double> > > Heed::EnTransfCS::fadda

Integral, normalised to unity.

Definition at line 88 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

fruth

std::vector<std::vector<std::vector<double> > > Heed::EnTransfCS::fruth

Rutherford term.

Definition at line 84 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

gamma_1

double Heed::EnTransfCS::gamma_1 = 0.

Lorentz factor - 1 (the best dimensionless measurement of speed).

Definition at line 34 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

hmd

HeedMatterDef* Heed::EnTransfCS::hmd = nullptr

Definition at line 47 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

length_y0

std::vector<double> Heed::EnTransfCS::length_y0

Definition at line 109 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

log1C

std::vector<double> Heed::EnTransfCS::log1C

common first log without cs

In the following arrays there is the only index: the energy. The meaning: the average value on the energy interval.

Definition at line 51 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

log2C

std::vector<double> Heed::EnTransfCS::log2C

common second log without cs

Definition at line 52 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

max_etransf

double Heed::EnTransfCS::max_etransf = 0.

Max. energy transfer [MeV].

Definition at line 37 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

mean

std::vector<std::vector<double> > Heed::EnTransfCS::mean

First moment, for each atom and shell.

Definition at line 101 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

meanC

double Heed::EnTransfCS::meanC = 0.

Definition at line 70 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

meanC1

double Heed::EnTransfCS::meanC1 = 0.

Definition at line 73 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

particle_charge

long Heed::EnTransfCS::particle_charge = 0

Charge in units of electron charge (used square, sign does not matter).

Definition at line 31 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

particle_mass

double Heed::EnTransfCS::particle_mass = 0.

Particle mass [MeV].

Definition at line 29 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

quan

std::vector<std::vector<double> > Heed::EnTransfCS::quan

Number of collisions / cm, for each atom and shell.

Definition at line 99 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

quanC

double Heed::EnTransfCS::quanC = 0.

Integrated (ionization) cross-section.

Definition at line 60 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

Rruth

std::vector<double> Heed::EnTransfCS::Rruth

term called R in my paper

Definition at line 55 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

s_primary_electron

bool Heed::EnTransfCS::s_primary_electron = false

Flag indicating whether the primary particle is an electron.

Definition at line 45 of file EnTransfCS.h.

Referenced by EnTransfCS(), and print().

s_simple_form

bool Heed::EnTransfCS::s_simple_form = true

Flag controlling the form of Rutherford scattering. For our purposes it is good to have simple form, so this variable is initialized to 1. Simple form means that there are two terms. The third term is assumed zero.

Definition at line 43 of file EnTransfCS.h.

Referenced by EnTransfCS().

---

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

• EnTransfCS.h
• EnTransfCS.cpp