EnTransfCS.h

Go to the documentation of this file.

#ifndef ENTRANFCS_H
#define ENTRANFCS_H
 
#include "heed++/code/HeedMatterDef.h"
 
namespace Heed {
 
#define EXCLUDE_A_VALUES   // exclude absorption values
 
/// 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
 
class EnTransfCS {
 public:
  /// Default constructor
  EnTransfCS() = default;
  /// Constructor
  EnTransfCS(double fparticle_mass, double fgamma_1, bool fs_primary_electron,
             HeedMatterDef* fhmd, double fparticle_charge = 1.,
             const bool debug = false);
 
  void print(std::ostream& file, int l) const;
  EnTransfCS* copy() const { return new EnTransfCS(*this); }
 
  /// Flag indicating whether the calculation was successful.
  bool m_ok = true;
 
  /// Particle mass [MeV]
  double particle_mass = 0.;
  /// Charge in units of electron charge (used square, sign does not matter).
  double particle_charge = 1.;
 
  /// Lorentz factor - 1 (the best dimensionless measurement of speed).
  double gamma_1 = 0.;
 
  /// Max. energy transfer [MeV]
  double max_etransf = 0.;
  /// 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.
  bool s_simple_form = true;
  /// Flag indicating whether the primary particle is an electron.
  bool s_primary_electron = false;
 
  HeedMatterDef* hmd = nullptr;
 
  /// Integrated (ionization) cross-section
  double quanC = 0.;
  /// First moment (mean restricted energy loss) [MeV]
  double meanC = 0.;
  /// First moment with additional tail to max. kinematically allowed transfer,
  /// calculated only for heavy particles (integral for electrons non-trivial).
  double meanC1 = 0.;
 
#ifndef EXCLUDE_A_VALUES
  /// Integrated (absorption) cross-section
  double quanC_a = 0.;
  double meanC1_a = 0.;
  double meanC_a = 0.;
#endif
 
  /// Integral, normalised to unity for each atom, shell and energy.
  std::vector<std::vector<std::vector<double> > > fadda;
  /// Number of collisions / cm, for each atom and shell.
  std::vector<std::vector<double> > quan;
#ifndef EXCLUDE_A_VALUES
  // Integral (total absorption), normalised to unity.
  std::vector<std::vector<std::vector<double> > > fadda_a;
  /// Number of collisions / cm (total absorption), for each atom and shell.
  std::vector<std::vector<double> > quan_a;
#endif
 
  std::vector<double> length_y0;
 
  // Prefactor (without thickness dependence) of the Highland formula.
  double sigma_ms = 0.; 
};
}
 
#endif