Abstract base class for media. More...

#include <Medium.hh>

Inheritance diagram for Garfield::Medium:

Public Member Functions
	Medium ()

virtual	~Medium ()

int	GetId () const

const std::string &	GetName () const

virtual bool	IsGas () const

virtual bool	IsSemiconductor () const

void	SetTemperature (const double t)

double	GetTemperature () const

void	SetPressure (const double p)

double	GetPressure () const

void	SetDielectricConstant (const double eps)

double	GetDielectricConstant () const

unsigned int	GetNumberOfComponents () const

virtual void	GetComponent (const unsigned int i, std::string &label, double &f)

virtual void	SetAtomicNumber (const double z)

virtual double	GetAtomicNumber () const

virtual void	SetAtomicWeight (const double a)

virtual double	GetAtomicWeight () const

virtual void	SetNumberDensity (const double n)

virtual double	GetNumberDensity () const

virtual void	SetMassDensity (const double rho)

virtual double	GetMassDensity () const

virtual void	EnableDrift ()

void	DisableDrift ()

virtual void	EnablePrimaryIonisation ()

void	DisablePrimaryIonisation ()

bool	IsDriftable () const

bool	IsMicroscopic () const

bool	IsIonisable () const

void	SetW (const double w)

double	GetW ()

void	SetFanoFactor (const double f)

double	GetFanoFactor ()

virtual bool	ElectronVelocity (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &vx, double &vy, double &vz)

virtual bool	ElectronDiffusion (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &dl, double &dt)

virtual bool	ElectronDiffusion (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double cov[3][3])

virtual bool	ElectronTownsend (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &alpha)

virtual bool	ElectronAttachment (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &eta)

virtual bool	ElectronLorentzAngle (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &lor)

virtual double	GetElectronEnergy (const double px, const double py, const double pz, double &vx, double &vy, double &vz, const int band=0)

virtual void	GetElectronMomentum (const double e, double &px, double &py, double &pz, int &band)

virtual double	GetElectronNullCollisionRate (const int band=0)

virtual double	GetElectronCollisionRate (const double e, const int band=0)

virtual bool	GetElectronCollision (const double e, int &type, int &level, double &e1, double &dx, double &dy, double &dz, int &nion, int &ndxc, int &band)

virtual unsigned int	GetNumberOfIonisationProducts () const

virtual bool	GetIonisationProduct (const unsigned int i, int &type, double &energy) const

virtual unsigned int	GetNumberOfDeexcitationProducts () const

virtual bool	GetDeexcitationProduct (const unsigned int i, double &t, double &s, int &type, double &energy) const

virtual bool	HoleVelocity (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &vx, double &vy, double &vz)

virtual bool	HoleDiffusion (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &dl, double &dt)

virtual bool	HoleDiffusion (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double cov[3][3])

virtual bool	HoleTownsend (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &alpha)

virtual bool	HoleAttachment (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &eta)

virtual bool	IonVelocity (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &vx, double &vy, double &vz)

virtual bool	IonDiffusion (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &dl, double &dt)

virtual bool	IonDissociation (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, double &diss)

void	SetFieldGrid (double emin, double emax, int ne, bool logE, double bmin=0., double bmax=0., int nb=1, double amin=0., double amax=0., int na=1)

void	SetFieldGrid (const std::vector< double > &efields, const std::vector< double > &bfields, const std::vector< double > &angles)

void	GetFieldGrid (std::vector< double > &efields, std::vector< double > &bfields, std::vector< double > &angles)

bool	GetElectronVelocityE (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &v)

bool	GetElectronVelocityExB (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &v)

bool	GetElectronVelocityB (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &v)

bool	GetElectronLongitudinalDiffusion (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &dl)

bool	GetElectronTransverseDiffusion (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &dt)

bool	GetElectronTownsend (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &alpha)

bool	GetElectronAttachment (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &eta)

bool	GetElectronLorentzAngle (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &lor)

bool	GetHoleVelocityE (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &v)

bool	GetHoleVelocityExB (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &v)

bool	GetHoleVelocityB (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &v)

bool	GetHoleLongitudinalDiffusion (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &dl)

bool	GetHoleTransverseDiffusion (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &dt)

bool	GetHoleTownsend (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &alpha)

bool	GetHoleAttachment (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &eta)

bool	GetIonMobility (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &mu)

bool	GetIonLongitudinalDiffusion (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &dl)

bool	GetIonTransverseDiffusion (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &dt)

bool	GetIonDissociation (const unsigned int ie, const unsigned int ib, const unsigned int ia, double &diss)

void	ResetElectronVelocity ()

void	ResetElectronDiffusion ()

void	ResetElectronTownsend ()

void	ResetElectronAttachment ()

void	ResetElectronLorentzAngle ()

void	ResetHoleVelocity ()

void	ResetHoleDiffusion ()

void	ResetHoleTownsend ()

void	ResetHoleAttachment ()

void	ResetIonMobility ()

void	ResetIonDiffusion ()

void	ResetIonDissociation ()

bool	SetIonMobility (const unsigned int ie, const unsigned int ib, const unsigned int ia, const double mu)

bool	SetIonMobility (const std::vector< double > &fields, const std::vector< double > &mobilities)

void	SetExtrapolationMethodVelocity (const std::string &extrLow, const std::string &extrHigh)

void	SetExtrapolationMethodDiffusion (const std::string &extrLow, const std::string &extrHigh)

void	SetExtrapolationMethodTownsend (const std::string &extrLow, const std::string &extrHigh)

void	SetExtrapolationMethodAttachment (const std::string &extrLow, const std::string &extrHigh)

void	SetExtrapolationMethodIonMobility (const std::string &extrLow, const std::string &extrHigh)

void	SetExtrapolationMethodIonDissociation (const std::string &extrLow, const std::string &extrHigh)

void	SetInterpolationMethodVelocity (const unsigned int intrp)

void	SetInterpolationMethodDiffusion (const unsigned int intrp)

void	SetInterpolationMethodTownsend (const unsigned int intrp)

void	SetInterpolationMethodAttachment (const unsigned int intrp)

void	SetInterpolationMethodIonMobility (const unsigned int intrp)

void	SetInterpolationMethodIonDissociation (const unsigned int intrp)

virtual double	ScaleElectricField (const double e) const

virtual double	UnScaleElectricField (const double e) const

virtual double	ScaleVelocity (const double v) const

virtual double	ScaleDiffusion (const double d) const

virtual double	ScaleDiffusionTensor (const double d) const

virtual double	ScaleTownsend (const double alpha) const

virtual double	ScaleAttachment (const double eta) const

virtual double	ScaleLorentzAngle (const double lor) const

virtual double	ScaleDissociation (const double diss) const

virtual bool	GetOpticalDataRange (double &emin, double &emax, const unsigned int i=0)

virtual bool	GetDielectricFunction (const double e, double &eps1, double &eps2, const unsigned int i=0)

virtual bool	GetPhotoAbsorptionCrossSection (const double e, double &sigma, const unsigned int i=0)

virtual double	GetPhotonCollisionRate (const double e)

virtual bool	GetPhotonCollision (const double e, int &type, int &level, double &e1, double &ctheta, int &nsec, double &esec)

void	EnableDebugging ()

void	DisableDebugging ()

Protected Member Functions
double	GetAngle (const double ex, const double ey, const double ez, const double bx, const double by, const double bz, const double e, const double b) const

double	Interpolate1D (const double e, const std::vector< double > &table, const std::vector< double > &fields, const unsigned int intpMeth, const int jExtr, const int iExtr)

bool	GetExtrapolationIndex (std::string extrStr, unsigned int &extrNb)

void	CloneTable (std::vector< std::vector< std::vector< double > > > &tab, const std::vector< double > &efields, const std::vector< double > &bfields, const std::vector< double > &angles, const unsigned int intp, const unsigned int extrLow, const unsigned int extrHigh, const double init, const std::string &label)

void	CloneTensor (std::vector< std::vector< std::vector< std::vector< double > > > > &tab, const unsigned int n, const std::vector< double > &efields, const std::vector< double > &bfields, const std::vector< double > &angles, const unsigned int intp, const unsigned int extrLow, const unsigned int extrHigh, const double init, const std::string &label)

void	InitParamArrays (const unsigned int eRes, const unsigned int bRes, const unsigned int aRes, std::vector< std::vector< std::vector< double > > > &tab, const double val)

void	InitParamTensor (const unsigned int eRes, const unsigned int bRes, const unsigned int aRes, const unsigned int tRes, std::vector< std::vector< std::vector< std::vector< double > > > > &tab, const double val)

Protected Attributes
std::string	m_className

int	m_id

std::string	m_name

double	m_temperature

double	m_pressure

double	m_epsilon

unsigned int	m_nComponents

double	m_z

double	m_a

double	m_density

bool	m_driftable

bool	m_microscopic

bool	m_ionisable

double	m_w

double	m_fano

bool	m_isChanged

bool	m_debug

std::vector< double >	m_eFields

std::vector< double >	m_bFields

std::vector< double >	m_bAngles

bool	m_map2d

bool	m_hasElectronVelocityE

bool	m_hasElectronVelocityB

bool	m_hasElectronVelocityExB

bool	m_hasElectronDiffLong

bool	m_hasElectronDiffTrans

bool	m_hasElectronDiffTens

bool	m_hasElectronAttachment

bool	m_hasElectronLorentzAngle

std::vector< std::vector< std::vector< double > > >	tabElectronVelocityE

std::vector< std::vector< std::vector< double > > >	tabElectronVelocityExB

std::vector< std::vector< std::vector< double > > >	tabElectronVelocityB

std::vector< std::vector< std::vector< double > > >	tabElectronDiffLong

std::vector< std::vector< std::vector< double > > >	tabElectronDiffTrans

std::vector< std::vector< std::vector< double > > >	tabElectronTownsend

std::vector< std::vector< std::vector< double > > >	tabElectronAttachment

std::vector< std::vector< std::vector< double > > >	tabElectronLorentzAngle

std::vector< std::vector< std::vector< std::vector< double > > > >	tabElectronDiffTens

bool	m_hasHoleVelocityE

bool	m_hasHoleVelocityB

bool	m_hasHoleVelocityExB

bool	m_hasHoleDiffLong

bool	m_hasHoleDiffTrans

bool	m_hasHoleDiffTens

bool	m_hasHoleTownsend

bool	m_hasHoleAttachment

std::vector< std::vector< std::vector< double > > >	tabHoleVelocityE

std::vector< std::vector< std::vector< double > > >	tabHoleVelocityExB

std::vector< std::vector< std::vector< double > > >	tabHoleVelocityB

std::vector< std::vector< std::vector< double > > >	tabHoleDiffLong

std::vector< std::vector< std::vector< double > > >	tabHoleDiffTrans

std::vector< std::vector< std::vector< double > > >	tabHoleTownsend

std::vector< std::vector< std::vector< double > > >	tabHoleAttachment

std::vector< std::vector< std::vector< std::vector< double > > > >	tabHoleDiffTens

bool	m_hasIonMobility

bool	m_hasIonDiffLong

bool	m_hasIonDiffTrans

bool	m_hasIonDissociation

std::vector< std::vector< std::vector< double > > >	tabIonMobility

std::vector< std::vector< std::vector< double > > >	tabIonDiffLong

std::vector< std::vector< std::vector< double > > >	tabIonDiffTrans

std::vector< std::vector< std::vector< double > > >	tabIonDissociation

int	thrElectronTownsend

int	thrElectronAttachment

int	thrHoleTownsend

int	thrHoleAttachment

int	thrIonDissociation

unsigned int	m_extrLowVelocity

unsigned int	m_extrHighVelocity

unsigned int	m_extrLowDiffusion

unsigned int	m_extrHighDiffusion

unsigned int	m_extrLowTownsend

unsigned int	m_extrHighTownsend

unsigned int	m_extrLowAttachment

unsigned int	m_extrHighAttachment

unsigned int	m_extrLowLorentzAngle

unsigned int	m_extrHighLorentzAngle

unsigned int	m_extrLowMobility

unsigned int	m_extrHighMobility

unsigned int	m_extrLowDissociation

unsigned int	m_extrHighDissociation

unsigned int	m_intpVelocity

unsigned int	m_intpDiffusion

unsigned int	m_intpTownsend

unsigned int	m_intpAttachment

unsigned int	m_intpLorentzAngle

unsigned int	m_intpMobility

unsigned int	m_intpDissociation

Static Protected Attributes
static int	m_idCounter = -1

Detailed Description

Abstract base class for media.

Definition at line 11 of file Medium.hh.

Constructor & Destructor Documentation

◆ Medium()

Garfield::Medium::Medium ( )

Definition at line 24 of file Medium.cc.

    : m_className("Medium"),
      m_id(++m_idCounter),
      m_name(""),
      m_temperature(293.15),
      m_pressure(760.),
      m_epsilon(1.),
      m_nComponents(1),
      m_z(1.),
      m_a(0.),
      m_density(0.),
      m_driftable(false),
      m_microscopic(false),
      m_ionisable(false),
      m_w(0.),
      m_fano(0.),
      m_isChanged(true),
      m_debug(false),
      m_map2d(false) {
 
  // Initialise the transport tables.
  m_bFields.assign(1, 0.);
  m_bAngles.assign(1, 0.);
 
  m_hasElectronVelocityE = false;
  m_hasElectronVelocityB = false;
  m_hasElectronVelocityExB = false;
  m_hasElectronDiffLong = false;
  m_hasElectronDiffTrans = false;
  m_hasElectronAttachment = false;
  m_hasElectronLorentzAngle = false;
  m_hasElectronDiffTens = false;
 
  m_hasHoleVelocityE = false;
  m_hasHoleVelocityB = false;
  m_hasHoleVelocityExB = false;
  m_hasHoleDiffLong = false;
  m_hasHoleDiffTrans = false;
  m_hasHoleTownsend = false;
  m_hasHoleAttachment = false;
  m_hasHoleDiffTens = false;
 
  m_hasIonMobility = false;
  m_hasIonDiffLong = false;
  m_hasIonDiffTrans = false;
  m_hasIonDissociation = false;
 
  m_extrLowVelocity = 0;
  m_extrHighVelocity = 1;
  m_extrLowDiffusion = 0;
  m_extrHighDiffusion = 1;
  m_extrLowTownsend = 0;
  m_extrHighTownsend = 1;
  m_extrLowAttachment = 0;
  m_extrHighAttachment = 1;
  m_extrLowLorentzAngle = 0;
  m_extrHighLorentzAngle = 1;
  m_extrLowMobility = 0;
  m_extrHighMobility = 1;
  m_extrLowDissociation = 0;
  m_extrHighDissociation = 1;
 
  m_intpVelocity = 2;
  m_intpDiffusion = 2;
  m_intpTownsend = 2;
  m_intpAttachment = 2;
  m_intpLorentzAngle = 2;
  m_intpMobility = 2;
  m_intpDissociation = 2;
 
  thrElectronTownsend = thrElectronAttachment = 0;
  thrHoleTownsend = thrHoleAttachment = 0;
  thrIonDissociation = 0;
 
  // Set the default grid.
  SetFieldGrid(100., 100000., 20, true, 0., 0., 1, 0., 0., 1);
}

◆ ~Medium()

Garfield::Medium::~Medium ( )

virtual

Definition at line 102 of file Medium.cc.

102{}

Member Function Documentation

◆ CloneTable()

void Garfield::Medium::CloneTable	(	std::vector< std::vector< std::vector< double > > > &	tab,
		const std::vector< double > &	efields,
		const std::vector< double > &	bfields,
		const std::vector< double > &	angles,
		const unsigned int	intp,
		const unsigned int	extrLow,
		const unsigned int	extrHigh,
		const double	init,
		const std::string &	label
	)

protected

Definition at line 2137 of file Medium.cc.

                                                                   {
 
  if (m_debug) {
    std::cout << m_className << "::CloneTable:\n";
    std::cout << "    Copying values of " << label << " to new grid.\n";
  }
 
  if (tab.empty()) {
    if (m_debug) std::cout << "    Table is empty.\n";
    return;
  }
  // Get the dimensions of the new grid.
  const int nEfieldsNew = efields.size();
  const int nBfieldsNew = bfields.size();
  const int nAnglesNew = angles.size();
 
  // Create a temporary table to store the values at the new grid points.
  std::vector<std::vector<std::vector<double> > > tabClone;
  InitParamArrays(nEfieldsNew, nBfieldsNew, nAnglesNew, tabClone, init);
 
  // Fill the temporary table.
  for (int i = 0; i < nEfieldsNew; ++i) {
    for (int j = 0; j < nBfieldsNew; ++j) {
      for (int k = 0; k < nAnglesNew; ++k) {
        double val = 0.;
        if (m_map2d) {
          if (!Numerics::Boxin3(tab, m_bAngles, m_bFields, m_eFields, m_bAngles.size(),
                                m_bFields.size(), m_eFields.size(), angles[k], bfields[j],
                                efields[i], val, intp)) {
            std::cerr << m_className << "::SetFieldGrid:\n";
            std::cerr << "    Interpolation of " << label << " failed.\n";
            std::cerr << "    Cannot copy value to new grid at: \n";
            std::cerr << "      E = " << efields[i] << "\n";
            std::cerr << "      B = " << bfields[j] << "\n";
            std::cerr << "      angle: " << angles[k] << "\n";
          } else {
            tabClone[k][j][i] = val;
          }
        } else {
          val = Interpolate1D(efields[i], tab[0][0], m_eFields, intp, extrLow,
                              extrHigh);
          tabClone[k][j][i] = val;
        }
      }
    }
  }
  // Re-dimension the original table.
  InitParamArrays(nEfieldsNew, nBfieldsNew, nAnglesNew, tab, init);
  // Copy the values to the original table.
  for (int i = 0; i < nEfieldsNew; ++i) {
    for (int j = 0; j < nBfieldsNew; ++j) {
      for (int k = 0; k < nAnglesNew; ++k) {
        tab[k][j][i] = tabClone[k][j][i];
      }
    }
  }
  tabClone.clear();
}

Referenced by SetFieldGrid().

◆ CloneTensor()

void Garfield::Medium::CloneTensor	(	std::vector< std::vector< std::vector< std::vector< double > > > > &	tab,
		const unsigned int	n,
		const std::vector< double > &	efields,
		const std::vector< double > &	bfields,
		const std::vector< double > &	angles,
		const unsigned int	intp,
		const unsigned int	extrLow,
		const unsigned int	extrHigh,
		const double	init,
		const std::string &	label
	)

protected

Definition at line 2203 of file Medium.cc.

                            {
 
  // If the table does not exist, do nothing.
  if (tab.empty()) return;
 
  // Get the dimensions of the new grid.
  const unsigned int nEfieldsNew = efields.size();
  const unsigned int nBfieldsNew = bfields.size();
  const unsigned int nAnglesNew = angles.size();
 
  // Create a temporary table to store the values at the new grid points.
  std::vector<std::vector<std::vector<std::vector<double> > > > tabClone;
  tabClone.clear();
  InitParamTensor(nEfieldsNew, nBfieldsNew, nAnglesNew, n, tabClone, init);
 
  // Fill the temporary table.
  for (unsigned int l = 0; l < n; ++l) {
    for (unsigned int i = 0; i < nEfieldsNew; ++i) {
      for (unsigned int j = 0; j < nBfieldsNew; ++j) {
        for (unsigned int k = 0; k < nAnglesNew; ++k) {
          double val = 0.;
          if (m_map2d) {
            if (!Numerics::Boxin3(tab[l], m_bAngles, m_bFields, m_eFields, m_bAngles.size(),
                                  m_bFields.size(), m_eFields.size(), angles[k], bfields[j],
                                  efields[i], val, intp)) {
              std::cerr << m_className << "::SetFieldGrid:\n";
              std::cerr << "    Interpolation of " << label << " failed.\n";
              std::cerr << "    Cannot copy value to new grid at: \n";
              std::cerr << "      Index: " << l << "\n";
              std::cerr << "      E = " << efields[i] << "\n";
              std::cerr << "      B = " << bfields[j] << "\n";
              std::cerr << "      angle: " << angles[k] << "\n";
            } else {
              tabClone[l][k][j][i] = val;
            }
          } else {
            val = Interpolate1D(efields[i], tab[l][0][0], m_eFields, intp,
                                extrLow, extrHigh);
            tabClone[l][k][j][i] = val;
          }
        }
      }
    }
  }
  // Re-dimension the original table.
  InitParamTensor(nEfieldsNew, nBfieldsNew, nAnglesNew, n, tab, 0.);
  // Copy the values to the original table.
  for (unsigned int l = 0; l < n; ++l) {
    for (unsigned int i = 0; i < nEfieldsNew; ++i) {
      for (unsigned int j = 0; j < nBfieldsNew; ++j) {
        for (unsigned int k = 0; k < nAnglesNew; ++k) {
          tab[l][k][j][i] = tabClone[l][k][j][i];
        }
      }
    }
  }
}

Referenced by SetFieldGrid().

◆ DisableDebugging()

void Garfield::Medium::DisableDebugging ( )

inline

Definition at line 291 of file Medium.hh.

291{ m_debug = false; }

Referenced by GarfieldPhysics::InitializePhysics().

◆ DisableDrift()

void Garfield::Medium::DisableDrift ( )

inline

Definition at line 53 of file Medium.hh.

53{ m_driftable = false; }

◆ DisablePrimaryIonisation()

void Garfield::Medium::DisablePrimaryIonisation ( )

inline

Definition at line 55 of file Medium.hh.

55{ m_ionisable = false; }

◆ ElectronAttachment()

bool Garfield::Medium::ElectronAttachment	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	eta
	)

virtual

Reimplemented in Garfield::MediumCdTe, Garfield::MediumGaAs, and Garfield::MediumSilicon.

Definition at line 544 of file Medium.cc.

                                                                               {
 
  eta = 0.;
  if (!m_hasElectronAttachment) return false;
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    // Interpolate.
    if (e0 < m_eFields[thrElectronAttachment]) {
      if (!Numerics::Boxin3(tabElectronAttachment, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, eta,
                            1)) {
        eta = -30.;
      }
    } else {
      if (!Numerics::Boxin3(tabElectronAttachment, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, eta,
                            m_intpAttachment)) {
        eta = -30.;
      }
    }
  } else {
    // Interpolate.
    if (e0 < m_eFields[thrElectronAttachment]) {
      eta = Interpolate1D(e0, tabElectronAttachment[0][0], m_eFields, 1,
                          m_extrLowAttachment, m_extrHighAttachment);
    } else {
      eta =
          Interpolate1D(e0, tabElectronAttachment[0][0], m_eFields,
                        m_intpAttachment, 
                        m_extrLowAttachment, m_extrHighAttachment);
    }
  }
 
  if (eta < -20.) {
    eta = 0.;
  } else {
    eta = exp(eta);
  }
 
  // Apply scaling.
  eta = ScaleAttachment(eta);
  return true;
}

Referenced by Garfield::MediumCdTe::ElectronAttachment(), Garfield::MediumGaAs::ElectronAttachment(), Garfield::MediumSilicon::ElectronAttachment(), and Garfield::ViewMedium::EvaluateFunction().

◆ ElectronDiffusion() [1/2]

bool Garfield::Medium::ElectronDiffusion	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	dl,
		double &	dt
	)

virtual

Definition at line 362 of file Medium.cc.

                                           {
 
  dl = dt = 0.;
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    // Interpolate.
    if (m_hasElectronDiffLong) {
      if (!Numerics::Boxin3(tabElectronDiffLong, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, dl,
                            m_intpDiffusion)) {
        dl = 0.;
      }
    }
    if (m_hasElectronDiffTrans) {
      if (!Numerics::Boxin3(tabElectronDiffTrans, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, dt,
                            m_intpDiffusion)) {
        dt = 0.;
      }
    }
  } else {
    if (m_hasElectronDiffLong) {
      dl = Interpolate1D(e0, tabElectronDiffLong[0][0], m_eFields, 
                         m_intpDiffusion,
                         m_extrLowDiffusion, m_extrHighDiffusion);
    }
    if (m_hasElectronDiffTrans) {
      dt = Interpolate1D(e0, tabElectronDiffTrans[0][0], m_eFields, 
                         m_intpDiffusion,
                         m_extrLowDiffusion, m_extrHighDiffusion);
    }
  }
 
  // If no data available, calculate
  // the diffusion coefficients using the Einstein relation
  if (!m_hasElectronDiffLong || !m_hasElectronDiffTrans) {
    const double d = sqrt(2. * BoltzmannConstant * m_temperature / e);
    if (!m_hasElectronDiffLong) dl = d;
    if (!m_hasElectronDiffTrans) dt = d;
  }
  // Verify values and apply scaling.
  if (dl < 0.) dl = 0.;
  if (dt < 0.) dt = 0.;
  dl = ScaleDiffusion(dl);
  dt = ScaleDiffusion(dt);
 
  return true;
}

Referenced by Garfield::ViewMedium::EvaluateFunction().

◆ ElectronDiffusion() [2/2]

bool Garfield::Medium::ElectronDiffusion	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double	cov[3][3]
	)

virtual

Definition at line 423 of file Medium.cc.

                                                 {
 
  // Initialise the tensor.
  cov[0][0] = cov[0][1] = cov[0][2] = 0.;
  cov[1][0] = cov[1][1] = cov[1][2] = 0.;
  cov[2][0] = cov[2][1] = cov[2][2] = 0.;
 
  if (!m_hasElectronDiffTens) return false;
 
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    // Interpolate.
    double diff = 0.;
    for (int l = 0; l < 6; ++l) {
      if (!Numerics::Boxin3(tabElectronDiffTens[l], m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, diff,
                            m_intpDiffusion)) {
        diff = 0.;
      }
      // Apply scaling.
      diff = ScaleDiffusionTensor(diff);
      if (l < 3) {
        cov[l][l] = diff;
      } else if (l == 3) {
        cov[0][1] = cov[1][0] = diff;
      } else if (l == 4) {
        cov[0][2] = cov[2][0] = diff;
      } else if (l == 5) {
        cov[1][2] = cov[2][1] = diff;
      }
    }
  } else {
    // Interpolate.
    for (int l = 0; l < 6; ++l) {
      double diff =
          Interpolate1D(e0, tabElectronDiffTens[l][0][0], m_eFields,
                        m_intpDiffusion, 
                        m_extrLowDiffusion, m_extrHighDiffusion);
      // Apply scaling.
      diff = ScaleDiffusionTensor(diff);
      if (l < 3) {
        cov[l][l] = diff;
      } else if (l == 3) {
        cov[0][1] = cov[1][0] = diff;
      } else if (l == 4) {
        cov[0][2] = cov[2][0] = diff;
      } else if (l == 5) {
        cov[1][2] = cov[2][1] = diff;
      }
    }
  }
 
  return true;
}

◆ ElectronLorentzAngle()

bool Garfield::Medium::ElectronLorentzAngle	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	lor
	)

virtual

Definition at line 599 of file Medium.cc.

                                               {
 
  lor = 0.;
  if (!m_hasElectronLorentzAngle) return false;
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    // Interpolate.
    if (!Numerics::Boxin3(tabElectronLorentzAngle, m_bAngles, m_bFields, m_eFields,
                          m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, lor,
                          m_intpLorentzAngle)) {
      lor = 0.;
    }
  } else {
    // Interpolate.
    lor =
        Interpolate1D(e0, tabElectronLorentzAngle[0][0], m_eFields,
                      m_intpLorentzAngle, 
                      m_extrLowLorentzAngle, m_extrHighLorentzAngle);
  }
  // Apply scaling.
  lor = ScaleLorentzAngle(lor);
  return true;
}

Referenced by Garfield::ViewMedium::EvaluateFunction().

◆ ElectronTownsend()

bool Garfield::Medium::ElectronTownsend	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	alpha
	)

virtual

Reimplemented in Garfield::MediumCdTe, Garfield::MediumGaAs, and Garfield::MediumSilicon.

Definition at line 490 of file Medium.cc.

                                             {
 
  alpha = 0.;
  if (tabElectronTownsend.empty()) return false;
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    // Interpolate.
    if (e0 < m_eFields[thrElectronTownsend]) {
      if (!Numerics::Boxin3(tabElectronTownsend, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, alpha,
                            1)) {
        alpha = -30.;
      }
    } else {
      if (!Numerics::Boxin3(tabElectronTownsend, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, alpha,
                            m_intpTownsend)) {
        alpha = -30.;
      }
    }
  } else {
    // Interpolate.
    if (e0 < m_eFields[thrElectronTownsend]) {
      alpha = Interpolate1D(e0, tabElectronTownsend[0][0], m_eFields, 1,
                            m_extrLowTownsend, m_extrHighTownsend);
    } else {
      alpha = Interpolate1D(e0, tabElectronTownsend[0][0], m_eFields,
                            m_intpTownsend, 
                            m_extrLowTownsend, m_extrHighTownsend);
    }
  }
 
  if (alpha < -20.) {
    alpha = 0.;
  } else {
    alpha = exp(alpha);
  }
 
  // Apply scaling.
  alpha = ScaleTownsend(alpha);
  return true;
}

Referenced by Garfield::MediumCdTe::ElectronTownsend(), Garfield::MediumGaAs::ElectronTownsend(), Garfield::MediumSilicon::ElectronTownsend(), and Garfield::ViewMedium::EvaluateFunction().

◆ ElectronVelocity()

bool Garfield::Medium::ElectronVelocity	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	vx,
		double &	vy,
		double &	vz
	)

virtual

Reimplemented in Garfield::MediumCdTe, Garfield::MediumGaAs, and Garfield::MediumSilicon.

Definition at line 204 of file Medium.cc.

                                                                  {
 
  vx = vy = vz = 0.;
  // Make sure there is at least a table of velocities along E.
  if (!m_hasElectronVelocityE) return false;
 
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return false;
 
  // Compute the magnitude of the magnetic field.
  const double b = sqrt(bx * bx + by * by + bz * bz);
 
  // Compute the angle between B field and E field.
  const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
  if (b < Small) {
    // No magnetic field.
 
    // Calculate the velocity along E.
    double ve = 0.;
    if (m_map2d) {
      if (!Numerics::Boxin3(tabElectronVelocityE, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, ve,
                            m_intpVelocity)) {
        std::cerr << m_className << "::ElectronVelocity:\n";
        std::cerr << "    Interpolation of velocity along E failed.\n";
        return false;
      }
    } else {
      ve = Interpolate1D(e0, tabElectronVelocityE[0][0], m_eFields, m_intpVelocity,
                         m_extrLowVelocity, m_extrHighVelocity);
    }
    const double q = -1.;
    const double mu = q * ve / e;
    vx = mu * ex;
    vy = mu * ey;
    vz = mu * ez;
 
  } else if (m_hasElectronVelocityB && m_hasElectronVelocityExB) {
    // Magnetic field, velocities along ExB and Bt available
 
    // Compute unit vectors along E, E x B and Bt.
    double ue[3] = {ex / e, ey / e, ez / e};
    double uexb[3] = {ey * bz - ez * by, ez * bx - ex * bz, ex * by - ey * bx};
    const double exb =
        sqrt(uexb[0] * uexb[0] + uexb[1] * uexb[1] + uexb[2] * uexb[2]);
    if (exb > 0.) {
      uexb[0] /= exb;
      uexb[1] /= exb;
      uexb[2] /= exb;
    } else {
      uexb[0] = ue[0];
      uexb[1] = ue[1];
      uexb[2] = ue[2];
    }
 
    double ubt[3] = {uexb[1] * ez - uexb[2] * ey, uexb[2] * ex - uexb[0] * ez,
                     uexb[0] * ey - uexb[1] * ex};
    const double bt = sqrt(ubt[0] * ubt[0] + ubt[1] * ubt[1] + ubt[2] * ubt[2]);
 
    if (bt > 0.) {
      ubt[0] /= bt;
      ubt[1] /= bt;
      ubt[2] /= bt;
    } else {
      ubt[0] = ue[0];
      ubt[1] = ue[1];
      ubt[2] = ue[2];
    }
 
    if (m_debug) {
      std::cout << std::setprecision(5);
      std::cout << m_className << "::ElectronVelocity:\n";
      std::cout << "    unit vector along E:     (" << ue[0] << ", " << ue[1]
                << ", " << ue[2] << ")\n";
      std::cout << "    unit vector along E x B: (" << uexb[0] << ", "
                << uexb[1] << ", " << uexb[2] << ")\n";
      std::cout << "    unit vector along Bt:    (" << ubt[0] << ", " << ubt[1]
                << ", " << ubt[2] << ")\n";
    }
 
    // Calculate the velocities in all directions.
    double ve = 0., vbt = 0., vexb = 0.;
    if (m_map2d) {
      if (!Numerics::Boxin3(tabElectronVelocityE, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, ve,
                            m_intpVelocity)) {
        std::cerr << m_className << "::ElectronVelocity:\n";
        std::cerr << "    Interpolation of velocity along E failed.\n";
        return false;
      }
      if (!Numerics::Boxin3(tabElectronVelocityExB, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, vexb,
                            m_intpVelocity)) {
        std::cerr << m_className << "::ElectronVelocity:\n";
        std::cerr << "    Interpolation of velocity along ExB failed.\n";
        return false;
      }
      if (!Numerics::Boxin3(tabElectronVelocityB, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, vbt,
                            m_intpVelocity)) {
        std::cerr << m_className << "::ElectronVelocity:\n";
        std::cerr << "    Interpolation of velocity along Bt failed.\n";
        return false;
      }
    } else {
      ve = Interpolate1D(e0, tabElectronVelocityE[0][0], m_eFields, 
                         m_intpVelocity,
                         m_extrLowVelocity, m_extrHighVelocity);
      vbt = Interpolate1D(e0, tabElectronVelocityB[0][0], m_eFields, 
                          m_intpVelocity,
                          m_extrLowVelocity, m_extrHighVelocity);
      vexb = Interpolate1D(e0, tabElectronVelocityExB[0][0], m_eFields,
                           m_intpVelocity, 
                           m_extrLowVelocity, m_extrHighVelocity);
    }
    const double q = -1.;
    if (ex * bx + ey * by + ez * bz > 0.) vbt = fabs(vbt);
     else vbt = -fabs(vbt);
    vx = q * (ve * ue[0] + q * q * vbt * ubt[0] + q * vexb * uexb[0]);
    vy = q * (ve * ue[1] + q * q * vbt * ubt[1] + q * vexb * uexb[1]);
    vz = q * (ve * ue[2] + q * q * vbt * ubt[2] + q * vexb * uexb[2]);
 
  } else {
    // Magnetic field, velocities along ExB, Bt not available
 
    // Calculate the velocity along E.
    double ve = 0.;
    if (m_map2d) {
      if (!Numerics::Boxin3(tabElectronVelocityE, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, ve,
                            m_intpVelocity)) {
        std::cerr << m_className << "::ElectronVelocity:\n";
        std::cerr << "    Interpolation of velocity along E failed.\n";
        return false;
      }
    } else {
      ve = Interpolate1D(e0, tabElectronVelocityE[0][0], m_eFields, 
                         m_intpVelocity,
                         m_extrLowVelocity, m_extrHighVelocity);
    }
 
    const double q = -1.;
    const double mu = q * ve / e;
    const double eb = bx * ex + by * ey + bz * ez;
    const double nom = 1. + pow(mu * b, 2);
    vx = mu * (ex + mu * (ey * bz - ez * by) + mu * mu * bx * eb) / nom;
    vy = mu * (ey + mu * (ez * bx - ex * bz) + mu * mu * by * eb) / nom;
    vz = mu * (ez + mu * (ex * by - ey * bx) + mu * mu * bz * eb) / nom;
  }
 
  return true;
}

Referenced by Garfield::MediumCdTe::ElectronVelocity(), Garfield::MediumGaAs::ElectronVelocity(), Garfield::MediumSilicon::ElectronVelocity(), and Garfield::ViewMedium::EvaluateFunction().

◆ EnableDebugging()

void Garfield::Medium::EnableDebugging ( )

inline

Definition at line 290 of file Medium.hh.

290{ m_debug = true; }

Referenced by GarfieldPhysics::InitializePhysics().

◆ EnableDrift()

virtual void Garfield::Medium::EnableDrift ( )

inlinevirtual

Reimplemented in Garfield::MediumConductor, and Garfield::MediumPlastic.

Definition at line 52 of file Medium.hh.

52{ m_driftable = true; }

Referenced by Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumGaAs::MediumGaAs(), Garfield::MediumGas::MediumGas(), Garfield::MediumMagboltz::MediumMagboltz(), and Garfield::MediumSilicon::MediumSilicon().

◆ EnablePrimaryIonisation()

virtual void Garfield::Medium::EnablePrimaryIonisation ( )

inlinevirtual

Definition at line 54 of file Medium.hh.

54{ m_ionisable = true; }

Referenced by Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumGaAs::MediumGaAs(), Garfield::MediumGas::MediumGas(), Garfield::MediumMagboltz::MediumMagboltz(), and Garfield::MediumSilicon::MediumSilicon().

◆ GetAngle()

double Garfield::Medium::GetAngle	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		const double	e,
		const double	b
	)		const

protected

Definition at line 2503 of file Medium.cc.

                                                              {
 
  // Ion
  /*
    if (e * b > 0.) {
      const double eb = fabs(ex * bx + ey * by + ez * bz);
      if (eb > 0.2 * e * b) {
        ebang = asin(std::min(
            1., sqrt(pow(ex * by - ey * bx, 2) + pow(ex * bz - ez * bx, 2) +
                     pow(ez * by - ey * bz, 2)) /
                    (e * b)));
      } else {
        ebang = acos(std::min(1., eb / (e * b)));
      }
    } else {
      ebang = m_bAngles[0];
    }
  */
  // Hole
  /*
    if (e * b > 0.) {
      const double eb = fabs(ex * bx + ey * by + ez * bz);
      if (eb > 0.2 * e * b) {
        ebang = asin(std::min(
            1., sqrt(pow(ex * by - ey * bx, 2) + pow(ex * bz - ez * bx, 2) +
                     pow(ez * by - ey * bz, 2)) /
                    (e * b)));
      } else {
        ebang = acos(std::min(1., eb / (e * b)));
      }
    } else {
      ebang = m_bAngles[0];
    }
  */
  /*
    if (e * b > 0.) {
      const double eb = fabs(ex * bx + ey * by + ez * bz);
      if (eb > 0.2 * e * b) {
        ebang = asin(std::min(
            1., sqrt(pow(ex * by - ey * bx, 2) + pow(ex * bz - ez * bx, 2) +
                     pow(ez * by - ey * bz, 2)) /
                    (e * b)));
      } else {
        ebang = acos(std::min(1., eb / (e * b)));
      }
    } else {
      ebang = m_bAngles[0];
    }
  */
  if (e * b <= 0.) return m_bAngles[0];
  const double eb = fabs(ex * bx + ey * by + ez * bz);
  if (eb > 0.2 * e * b) {
    const double ebxy = ex * by - ey * bx;
    const double ebxz = ex * bz - ez * bx;
    const double ebzy = ez * by - ey * bz;
    return asin(std::min(1., sqrt(ebxy * ebxy + ebxz * ebxz + ebzy * ebzy) / (e * b)));
  }
  return acos(std::min(1., eb / (e * b)));
} 

Referenced by ElectronAttachment(), ElectronDiffusion(), ElectronLorentzAngle(), ElectronTownsend(), ElectronVelocity(), HoleAttachment(), HoleDiffusion(), HoleTownsend(), HoleVelocity(), IonDiffusion(), IonDissociation(), and IonVelocity().

◆ GetAtomicNumber()

virtual double Garfield::Medium::GetAtomicNumber ( ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 42 of file Medium.hh.

42{ return m_z; }

◆ GetAtomicWeight()

virtual double Garfield::Medium::GetAtomicWeight ( ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 44 of file Medium.hh.

44{ return m_a; }

◆ GetComponent()

void Garfield::Medium::GetComponent	(	const unsigned int	i,
		std::string &	label,
		double &	f
	)

virtual

Reimplemented in Garfield::MediumCdTe, Garfield::MediumGaAs, and Garfield::MediumGas.

Definition at line 142 of file Medium.cc.

                                                       {
 
  if (i >= m_nComponents) {
    std::cerr << m_className << "::GetComponent:\n";
    std::cerr << "    Index out of range.\n";
  }
 
  label = m_name;
  f = 1.;
}

◆ GetDeexcitationProduct()

bool Garfield::Medium::GetDeexcitationProduct	(	const unsigned int	i,
		double &	t,
		double &	s,
		int &	type,
		double &	energy
	)		const

virtual

Reimplemented in Garfield::MediumMagboltz.

Definition at line 694 of file Medium.cc.

                                                                     {
 
  if (m_debug) PrintNotImplemented(m_className, "GetDeexcitationProduct");
  t = s = energy = 0.;
  type = 0;
  return false;
}

◆ GetDielectricConstant()

double Garfield::Medium::GetDielectricConstant ( ) const

inline

Definition at line 34 of file Medium.hh.

34{ return m_epsilon; }

◆ GetDielectricFunction()

bool Garfield::Medium::GetDielectricFunction	(	const double	e,
		double &	eps1,
		double &	eps2,
		const unsigned int	i = `0`
	)

virtual

Reimplemented in Garfield::MediumCdTe, Garfield::MediumGaAs, and Garfield::MediumSilicon.

Definition at line 1243 of file Medium.cc.

                                                         {
 
  if (i >= m_nComponents) {
    std::cerr << m_className << "::GetDielectricFunction:\n";
    std::cerr << "    Component " << i << " does not exist.\n";
    return false;
  }
 
  if (e < 0.) {
    std::cerr << m_className << "::GetDielectricFunction:\n";
    std::cerr << "    Energy must be > 0.\n";
    return false;
  }
 
  if (m_debug) PrintNotImplemented(m_className, "GetDielectricFunction");
  eps1 = 1.;
  eps2 = 0.;
  return false;
}

◆ GetElectronAttachment()

bool Garfield::Medium::GetElectronAttachment	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	eta
	)

Definition at line 1826 of file Medium.cc.

                                                                       {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetElectronAttachment:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    eta = 0.;
    return false;
  }
  if (!m_hasElectronAttachment) {
    if (m_debug) {
      std::cerr << m_className << "::GetElectronAttachment:\n";
      std::cerr << "    Data not available.\n";
    }
    eta = 0.;
    return false;
  }
 
  eta = tabElectronAttachment[ia][ib][ie];
  return true;
}

◆ GetElectronCollision()

bool Garfield::Medium::GetElectronCollision	(	const double	e,
		int &	type,
		int &	level,
		double &	e1,
		double &	dx,
		double &	dy,
		double &	dz,
		int &	nion,
		int &	ndxc,
		int &	band
	)

virtual

Reimplemented in Garfield::MediumMagboltz, and Garfield::MediumSilicon.

Definition at line 672 of file Medium.cc.

                                                                               {
 
  type = level = -1;
  e1 = e;
  nion = ndxc = band = 0;
  RndmDirection(dx, dy, dz);
 
  if (m_debug) PrintNotImplemented(m_className, "GetElectronCollision");
  return false;
}

◆ GetElectronCollisionRate()

double Garfield::Medium::GetElectronCollisionRate	(	const double	e,
		const int	band = `0`
	)

virtual

Reimplemented in Garfield::MediumMagboltz, and Garfield::MediumSilicon.

Definition at line 665 of file Medium.cc.

                                                     {
 
  if (m_debug) PrintNotImplemented(m_className, "GetElectronCollisionRate");
  return 0.;
}

◆ GetElectronEnergy()

double Garfield::Medium::GetElectronEnergy	(	const double	px,
		const double	py,
		const double	pz,
		double &	vx,
		double &	vy,
		double &	vz,
		const int	band = `0`
	)

virtual

Reimplemented in Garfield::MediumSilicon.

Definition at line 635 of file Medium.cc.

                                                             {
 
  if (band > 0) {
    std::cerr << m_className << "::GetElectronEnergy:\n";
    std::cerr << "    Unknown band index.\n";
  }
 
  vx = SpeedOfLight * px / ElectronMass;
  vy = SpeedOfLight * py / ElectronMass;
  vz = SpeedOfLight * pz / ElectronMass;
 
  return 0.5 * (px * px + py * py + pz * pz) / ElectronMass;
}

◆ GetElectronLongitudinalDiffusion()

bool Garfield::Medium::GetElectronLongitudinalDiffusion	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	dl
	)

Definition at line 1752 of file Medium.cc.

                                                          {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetElectronLongitudinalDiffusion:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    dl = 0.;
    return false;
  }
  if (!m_hasElectronDiffLong) {
    if (m_debug) {
      std::cerr << m_className << "::GetElectronLongitudinalDiffusion:\n";
      std::cerr << "    Data not available.\n";
    }
    dl = 0.;
    return false;
  }
 
  dl = tabElectronDiffLong[ia][ib][ie];
  return true;
}

◆ GetElectronLorentzAngle()

bool Garfield::Medium::GetElectronLorentzAngle	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	lor
	)

Definition at line 1850 of file Medium.cc.

                                                                         {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetElectronLorentzAngle:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    lor = 0.;
    return false;
  }
  if (!m_hasElectronLorentzAngle) {
    if (m_debug) {
      std::cerr << m_className << "::GetElectronLorentzAngle:\n";
      std::cerr << "    Data not available.\n";
    }
    lor = 0.;
    return false;
  }
 
  lor = tabElectronLorentzAngle[ia][ib][ie];
  return true;
}

◆ GetElectronMomentum()

void Garfield::Medium::GetElectronMomentum	(	const double	e,
		double &	px,
		double &	py,
		double &	pz,
		int &	band
	)

virtual

Reimplemented in Garfield::MediumSilicon.

Definition at line 651 of file Medium.cc.

                                                        {
 
  const double p = sqrt(2. * ElectronMass * e) / SpeedOfLight;
  RndmDirection(px, py, pz, p);
  band = -1;
}

◆ GetElectronNullCollisionRate()

double Garfield::Medium::GetElectronNullCollisionRate ( const int band = 0 )

virtual

Reimplemented in Garfield::MediumMagboltz, and Garfield::MediumSilicon.

Definition at line 659 of file Medium.cc.

                                                       {
 
  if (m_debug) PrintNotImplemented(m_className, "GetElectronNullCollisionRate");
  return 0.;
}

◆ GetElectronTownsend()

bool Garfield::Medium::GetElectronTownsend	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	alpha
	)

Definition at line 1802 of file Medium.cc.

                                                                       {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetElectronTownsend:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    alpha = 0.;
    return false;
  }
  if (tabElectronTownsend.empty()) {
    if (m_debug) {
      std::cerr << m_className << "::GetElectronTownsend:\n";
      std::cerr << "    Data not available.\n";
    }
    alpha = 0.;
    return false;
  }
 
  alpha = tabElectronTownsend[ia][ib][ie];
  return true;
}

◆ GetElectronTransverseDiffusion()

bool Garfield::Medium::GetElectronTransverseDiffusion	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	dt
	)

Definition at line 1777 of file Medium.cc.

                                                        {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetElectronTransverseDiffusion:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    dt = 0.;
    return false;
  }
  if (!m_hasElectronDiffTrans) {
    if (m_debug) {
      std::cerr << m_className << "::GetElectronTransverseDiffusion:\n";
      std::cerr << "    Data not available.\n";
    }
    dt = 0.;
    return false;
  }
 
  dt = tabElectronDiffTrans[ia][ib][ie];
  return true;
}

◆ GetElectronVelocityB()

bool Garfield::Medium::GetElectronVelocityB	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	v
	)

Definition at line 1728 of file Medium.cc.

                                                                    {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetElectronVelocityB:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    v = 0.;
    return false;
  }
  if (!m_hasElectronVelocityB) {
    if (m_debug) {
      std::cerr << m_className << "::GetElectronVelocityB:\n";
      std::cerr << "    Data not available.\n";
    }
    v = 0.;
    return false;
  }
 
  v = tabElectronVelocityB[ia][ib][ie];
  return true;
}

◆ GetElectronVelocityE()

bool Garfield::Medium::GetElectronVelocityE	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	v
	)

Definition at line 1680 of file Medium.cc.

                                                                    {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetElectronVelocityE:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    v = 0.;
    return false;
  }
  if (!m_hasElectronVelocityE) {
    if (m_debug) {
      std::cerr << m_className << "::GetElectronVelocityE:\n";
      std::cerr << "    Data not available.\n";
    }
    v = 0.;
    return false;
  }
 
  v = tabElectronVelocityE[ia][ib][ie];
  return true;
}

◆ GetElectronVelocityExB()

bool Garfield::Medium::GetElectronVelocityExB	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	v
	)

Definition at line 1704 of file Medium.cc.

                                                                      {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetElectronVelocityExB:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    v = 0.;
    return false;
  }
  if (!m_hasElectronVelocityExB) {
    if (m_debug) {
      std::cerr << m_className << "::GetElectronVelocityExB:\n";
      std::cerr << "    Data not available.\n";
    }
    v = 0.;
    return false;
  }
 
  v = tabElectronVelocityExB[ia][ib][ie];
  return true;
}

◆ GetExtrapolationIndex()

bool Garfield::Medium::GetExtrapolationIndex	(	std::string	extrStr,
		unsigned int &	extrNb
	)

protected

Definition at line 2453 of file Medium.cc.

                                                                          {
 
  // Convert to upper-case
  for (unsigned int i = 0; i < extrStr.length(); ++i) {
    extrStr[i] = toupper(extrStr[i]);
  }
 
  if (extrStr == "CONST" || extrStr == "CONSTANT") {
    extrNb = 0;
  } else if (extrStr == "LIN" || extrStr == "LINEAR") {
    extrNb = 1;
  } else if (extrStr == "EXP" || extrStr == "EXPONENTIAL") {
    extrNb = 2;
  } else {
    return false;
  }
 
  return true;
}

Referenced by SetExtrapolationMethodAttachment(), SetExtrapolationMethodDiffusion(), Garfield::MediumGas::SetExtrapolationMethodExcitationRates(), SetExtrapolationMethodIonDissociation(), Garfield::MediumGas::SetExtrapolationMethodIonisationRates(), SetExtrapolationMethodIonMobility(), SetExtrapolationMethodTownsend(), and SetExtrapolationMethodVelocity().

◆ GetFanoFactor()

double Garfield::Medium::GetFanoFactor ( )

inline

Definition at line 65 of file Medium.hh.

65{ return m_fano; }

◆ GetFieldGrid()

void Garfield::Medium::GetFieldGrid	(	std::vector< double > &	efields,
		std::vector< double > &	bfields,
		std::vector< double > &	angles
	)

Definition at line 1671 of file Medium.cc.

                                                     {
 
  efields = m_eFields;
  bfields = m_bFields;
  angles = m_bAngles;
}

◆ GetHoleAttachment()

bool Garfield::Medium::GetHoleAttachment	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	eta
	)

Definition at line 2018 of file Medium.cc.

                                                                   {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetHoleAttachment:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    eta = 0.;
    return false;
  }
  if (!m_hasHoleAttachment) {
    if (m_debug) {
      std::cerr << m_className << "::GetHoleAttachment:\n";
      std::cerr << "    Data not available.\n";
    }
    eta = 0.;
    return false;
  }
 
  eta = tabHoleAttachment[ia][ib][ie];
  return true;
}

◆ GetHoleLongitudinalDiffusion()

bool Garfield::Medium::GetHoleLongitudinalDiffusion	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	dl
	)

Definition at line 1946 of file Medium.cc.

                                                                             {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetHoleLongitudinalDiffusion:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    dl = 0.;
    return false;
  }
  if (!m_hasHoleDiffLong) {
    if (m_debug) {
      std::cerr << m_className << "::GetHoleLongitudinalDiffusion:\n";
      std::cerr << "    Data not available.\n";
    }
    dl = 0.;
    return false;
  }
 
  dl = tabHoleDiffLong[ia][ib][ie];
  return true;
}

◆ GetHoleTownsend()

bool Garfield::Medium::GetHoleTownsend	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	alpha
	)

Definition at line 1994 of file Medium.cc.

                                                                   {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetHoleTownsend:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    alpha = 0.;
    return false;
  }
  if (!m_hasHoleTownsend) {
    if (m_debug) {
      std::cerr << m_className << "::GetHoleTownsend:\n";
      std::cerr << "    Data not available.\n";
    }
    alpha = 0.;
    return false;
  }
 
  alpha = tabHoleTownsend[ia][ib][ie];
  return true;
}

◆ GetHoleTransverseDiffusion()

bool Garfield::Medium::GetHoleTransverseDiffusion	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	dt
	)

Definition at line 1970 of file Medium.cc.

                                                                           {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetHoleTransverseDiffusion:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    dt = 0.;
    return false;
  }
  if (!m_hasHoleDiffTrans) {
    if (m_debug) {
      std::cerr << m_className << "::GetHoleTransverseDiffusion:\n";
      std::cerr << "    Data not available.\n";
    }
    dt = 0.;
    return false;
  }
 
  dt = tabHoleDiffTrans[ia][ib][ie];
  return true;
}

◆ GetHoleVelocityB()

bool Garfield::Medium::GetHoleVelocityB	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	v
	)

Definition at line 1922 of file Medium.cc.

                                                                {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetHoleVelocityB:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    v = 0.;
    return false;
  }
  if (!m_hasHoleVelocityB) {
    if (m_debug) {
      std::cerr << m_className << "::GetHoleVelocityB:\n";
      std::cerr << "    Data not available.\n";
    }
    v = 0.;
    return false;
  }
 
  v = tabHoleVelocityB[ia][ib][ie];
  return true;
}

◆ GetHoleVelocityE()

bool Garfield::Medium::GetHoleVelocityE	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	v
	)

Definition at line 1874 of file Medium.cc.

                                                                {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetHoleVelocityE:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    v = 0.;
    return false;
  }
  if (!m_hasHoleVelocityE) {
    if (m_debug) {
      std::cerr << m_className << "::GetHoleVelocityE:\n";
      std::cerr << "    Data not available.\n";
    }
    v = 0.;
    return false;
  }
 
  v = tabHoleVelocityE[ia][ib][ie];
  return true;
}

◆ GetHoleVelocityExB()

bool Garfield::Medium::GetHoleVelocityExB	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	v
	)

Definition at line 1898 of file Medium.cc.

                                                                  {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetHoleVelocityExB:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    v = 0.;
    return false;
  }
  if (!m_hasHoleVelocityExB) {
    if (m_debug) {
      std::cerr << m_className << "::GetHoleVelocityExB:\n";
      std::cerr << "    Data not available.\n";
    }
    v = 0.;
    return false;
  }
 
  v = tabHoleVelocityExB[ia][ib][ie];
  return true;
}

◆ GetId()

int Garfield::Medium::GetId ( ) const

inline

Definition at line 20 of file Medium.hh.

20{ return m_id; }

Referenced by Garfield::GeometrySimple::AddSolid(), and Garfield::ViewGeometry::Plot().

◆ GetIonDissociation()

bool Garfield::Medium::GetIonDissociation	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	diss
	)

Definition at line 2113 of file Medium.cc.

                                                                     {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetIonDissociation:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    diss = 0.;
    return false;
  }
  if (!m_hasIonDissociation) {
    if (m_debug) {
      std::cerr << m_className << "::GetIonDissociation:\n";
      std::cerr << "    Data not available.\n";
    }
    diss = 0.;
    return false;
  }
 
  diss = tabIonDissociation[ia][ib][ie];
  return true;
}

◆ GetIonisationProduct()

bool Garfield::Medium::GetIonisationProduct	(	const unsigned int	i,
		int &	type,
		double &	energy
	)		const

virtual

Reimplemented in Garfield::MediumMagboltz, and Garfield::MediumSilicon.

Definition at line 685 of file Medium.cc.

                                                                   {
 
  if (m_debug) PrintNotImplemented(m_className, "GetIonisationProduct");
  type = 0;
  energy = 0.;
  return false;
}

◆ GetIonLongitudinalDiffusion()

bool Garfield::Medium::GetIonLongitudinalDiffusion	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	dl
	)

Definition at line 2065 of file Medium.cc.

                                                                            {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetIonLongitudinalDiffusion:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    dl = 0.;
    return false;
  }
  if (!m_hasIonDiffLong) {
    if (m_debug) {
      std::cerr << m_className << "::GetIonLongitudinalDiffusion:\n";
      std::cerr << "    Data not available.\n";
    }
    dl = 0.;
    return false;
  }
 
  dl = tabIonDiffLong[ia][ib][ie];
  return true;
}

◆ GetIonMobility()

bool Garfield::Medium::GetIonMobility	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	mu
	)

Definition at line 2042 of file Medium.cc.

                                                               {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetIonMobility:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    mu = 0.;
    return false;
  }
  if (!m_hasIonMobility) {
    if (m_debug) {
      std::cerr << m_className << "::GetIonMobility:\n";
      std::cerr << "    Data not available.\n";
    }
    mu = 0.;
    return false;
  }
 
  mu = tabIonMobility[ia][ib][ie];
  return true;
}

◆ GetIonTransverseDiffusion()

bool Garfield::Medium::GetIonTransverseDiffusion	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		double &	dt
	)

Definition at line 2089 of file Medium.cc.

                                                                          {
 
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || ia >= m_bAngles.size()) {
    std::cerr << m_className << "::GetIonTransverseDiffusion:\n";
    std::cerr << "     Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    dt = 0.;
    return false;
  }
  if (!m_hasIonDiffTrans) {
    if (m_debug) {
      std::cerr << m_className << "::GetIonTransverseDiffusion:\n";
      std::cerr << "    Data not available.\n";
    }
    dt = 0.;
    return false;
  }
 
  dt = tabIonDiffTrans[ia][ib][ie];
  return true;
}

◆ GetMassDensity()

double Garfield::Medium::GetMassDensity ( ) const

virtual

Reimplemented in Garfield::MediumGas.

Definition at line 137 of file Medium.cc.

                                    {
 
  return m_density * AtomicMassUnit * m_a;
}

Referenced by Garfield::TrackHeed::NewTrack(), Garfield::GeometryRoot::SetMedium(), Garfield::TrackHeed::TransportDeltaElectron(), and Garfield::TrackHeed::TransportPhoton().

◆ GetName()

const std::string & Garfield::Medium::GetName ( ) const

inline

Definition at line 22 of file Medium.hh.

22{ return m_name; }

Referenced by Garfield::TrackBichsel::GetCluster(), Garfield::TrackElectron::GetCluster(), Garfield::TrackPAI::GetCluster(), Garfield::TrackBichsel::NewTrack(), Garfield::TrackElectron::NewTrack(), Garfield::TrackHeed::NewTrack(), Garfield::TrackPAI::NewTrack(), Garfield::GeometryRoot::SetMedium(), Garfield::ComponentFieldMap::SetMedium(), Garfield::TrackHeed::TransportDeltaElectron(), and Garfield::TrackHeed::TransportPhoton().

◆ GetNumberDensity()

virtual double Garfield::Medium::GetNumberDensity ( ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 47 of file Medium.hh.

47{ return m_density; }

Referenced by Garfield::TrackElectron::GetCluster(), Garfield::TrackPAI::GetCluster(), and Garfield::TrackPAI::NewTrack().

◆ GetNumberOfComponents()

unsigned int Garfield::Medium::GetNumberOfComponents ( ) const

inline

Definition at line 37 of file Medium.hh.

37{ return m_nComponents; }

◆ GetNumberOfDeexcitationProducts()

virtual unsigned int Garfield::Medium::GetNumberOfDeexcitationProducts ( ) const

inlinevirtual

Reimplemented in Garfield::MediumMagboltz.

Definition at line 123 of file Medium.hh.

123{ return 0; }

◆ GetNumberOfIonisationProducts()

virtual unsigned int Garfield::Medium::GetNumberOfIonisationProducts ( ) const

inlinevirtual

Reimplemented in Garfield::MediumMagboltz, and Garfield::MediumSilicon.

Definition at line 119 of file Medium.hh.

119{ return 0; }

◆ GetOpticalDataRange()

bool Garfield::Medium::GetOpticalDataRange	(	double &	emin,
		double &	emax,
		const unsigned int	i = `0`
	)

virtual

Reimplemented in Garfield::MediumCdTe, Garfield::MediumGaAs, and Garfield::MediumSilicon.

Definition at line 1229 of file Medium.cc.

                                                       {
 
  if (i >= m_nComponents) {
    std::cerr << m_className << "::GetOpticalDataRange:\n";
    std::cerr << "    Component " << i << " does not exist.\n";
    return false;
  }
 
  if (m_debug) PrintNotImplemented(m_className, "GetOpticalDataRange");
  emin = emax = 0.;
  return false;
}

◆ GetPhotoAbsorptionCrossSection()

bool Garfield::Medium::GetPhotoAbsorptionCrossSection	(	const double	e,
		double &	sigma,
		const unsigned int	i = `0`
	)

virtual

Definition at line 1264 of file Medium.cc.

                                                                  {
 
  if (i >= m_nComponents) {
    std::cerr << m_className << "::GetPhotoAbsorptionCrossSection:\n";
    std::cerr << "    Component " << i << " does not exist.\n";
    return false;
  }
 
  if (e < 0.) {
    std::cerr << m_className << "::GetPhotoAbsorptionCrossSection:\n";
    std::cerr << "    Energy must be > 0.\n";
    return false;
  }
 
  if (m_debug) {
    PrintNotImplemented(m_className, "GetPhotoAbsorptionCrossSection");
  }
  sigma = 0.;
  return false;
}

Referenced by GetPhotonCollisionRate().

◆ GetPhotonCollision()

bool Garfield::Medium::GetPhotonCollision	(	const double	e,
		int &	type,
		int &	level,
		double &	e1,
		double &	ctheta,
		int &	nsec,
		double &	esec
	)

virtual

Reimplemented in Garfield::MediumMagboltz.

Definition at line 1294 of file Medium.cc.

                                              {
 
  type = level = -1;
  e1 = e;
  ctheta = 1.;
  nsec = 0;
  esec = 0.;
  return false;
}

◆ GetPhotonCollisionRate()

double Garfield::Medium::GetPhotonCollisionRate ( const double e )

virtual

Reimplemented in Garfield::MediumMagboltz.

Definition at line 1286 of file Medium.cc.

                                                    {
 
  double sigma = 0.;
  if (!GetPhotoAbsorptionCrossSection(e, sigma)) return 0.;
 
  return sigma * m_density * SpeedOfLight;
}

◆ GetPressure()

double Garfield::Medium::GetPressure ( ) const

inline

Definition at line 31 of file Medium.hh.

31{ return m_pressure; }

◆ GetTemperature()

double Garfield::Medium::GetTemperature ( ) const

inline

Definition at line 28 of file Medium.hh.

28{ return m_temperature; }

◆ GetW()

double Garfield::Medium::GetW ( )

inline

Definition at line 63 of file Medium.hh.

63{ return m_w; }

◆ HoleAttachment()

bool Garfield::Medium::HoleAttachment	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	eta
	)

virtual

Reimplemented in Garfield::MediumCdTe, Garfield::MediumGaAs, and Garfield::MediumSilicon.

Definition at line 1023 of file Medium.cc.

                                         {
 
  eta = 0.;
  if (!m_hasHoleAttachment) return false;
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    // Interpolate.
    if (e0 < m_eFields[thrHoleAttachment]) {
      if (!Numerics::Boxin3(tabHoleAttachment, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, eta,
                            1)) {
        eta = -30.;
      }
    } else {
      if (!Numerics::Boxin3(tabHoleAttachment, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, eta,
                            m_intpAttachment)) {
        eta = -30.;
      }
    }
  } else {
    // Interpolate.
    if (e0 < m_eFields[thrHoleAttachment]) {
      eta = Interpolate1D(e0, tabHoleAttachment[0][0], m_eFields, 1,
                          m_extrLowAttachment, m_extrHighAttachment);
    } else {
      eta = Interpolate1D(e0, tabHoleAttachment[0][0], m_eFields, 
                          m_intpAttachment,
                          m_extrLowAttachment, m_extrHighAttachment);
    }
  }
 
  if (eta < -20.) {
    eta = 0.;
  } else {
    eta = exp(eta);
  }
 
  // Apply scaling.
  eta = ScaleAttachment(eta);
  return true;
}

Referenced by Garfield::ViewMedium::EvaluateFunction(), Garfield::MediumCdTe::HoleAttachment(), Garfield::MediumGaAs::HoleAttachment(), and Garfield::MediumSilicon::HoleAttachment().

◆ HoleDiffusion() [1/2]

bool Garfield::Medium::HoleDiffusion	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	dl,
		double &	dt
	)

virtual

Definition at line 846 of file Medium.cc.

                                                   {
 
  dl = dt = 0.;
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    // Interpolate.
    if (m_hasHoleDiffLong) {
      if (!Numerics::Boxin3(tabHoleDiffLong, m_bAngles, m_bFields, m_eFields, m_bAngles.size(),
                            m_bFields.size(), m_eFields.size(), ebang, b, e0, dl,
                            m_intpDiffusion)) {
        dl = 0.;
      }
    }
    if (m_hasHoleDiffTrans) {
      if (!Numerics::Boxin3(tabHoleDiffTrans, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, dt,
                            m_intpDiffusion)) {
        dt = 0.;
      }
    }
  } else {
    if (m_hasHoleDiffLong) {
      dl = Interpolate1D(e0, tabHoleDiffLong[0][0], m_eFields, m_intpDiffusion,
                         m_extrLowDiffusion, m_extrHighDiffusion);
    }
    if (m_hasHoleDiffTrans) {
      dt = Interpolate1D(e0, tabHoleDiffTrans[0][0], m_eFields, m_intpDiffusion,
                         m_extrLowDiffusion, m_extrHighDiffusion);
    }
  }
 
  // If no data available, calculate
  // the diffusion coefficients using the Einstein relation
  if (!m_hasHoleDiffLong) {
    dl = sqrt(2. * BoltzmannConstant * m_temperature / e);
  }
  if (!m_hasHoleDiffTrans) {
    dt = sqrt(2. * BoltzmannConstant * m_temperature / e);
  }
 
  // Verify values and apply scaling.
  if (dl < 0.) dl = 0.;
  if (dt < 0.) dt = 0.;
  dl = ScaleDiffusion(dl);
  dt = ScaleDiffusion(dt);
 
  return true;
}

Referenced by Garfield::ViewMedium::EvaluateFunction().

◆ HoleDiffusion() [2/2]

bool Garfield::Medium::HoleDiffusion	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double	cov[3][3]
	)

virtual

Definition at line 906 of file Medium.cc.

                                             {
 
  // Initialise the tensor.
  cov[0][0] = cov[0][1] = cov[0][2] = 0.;
  cov[1][0] = cov[1][1] = cov[1][2] = 0.;
  cov[2][0] = cov[2][1] = cov[2][2] = 0.;
 
  if (!m_hasHoleDiffTens) return false;
 
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    // Interpolate.
    double diff = 0.;
    for (int l = 0; l < 6; ++l) {
      if (!Numerics::Boxin3(tabHoleDiffTens[l], m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, diff,
                            m_intpDiffusion)) {
        diff = 0.;
      }
      // Apply scaling.
      diff = ScaleDiffusionTensor(diff);
      if (l < 3) {
        cov[l][l] = diff;
      } else if (l == 3) {
        cov[0][1] = cov[1][0] = diff;
      } else if (l == 4) {
        cov[0][2] = cov[2][0] = diff;
      } else if (l == 5) {
        cov[1][2] = cov[2][1] = diff;
      }
    }
  } else {
    // Interpolate.
    for (int l = 0; l < 6; ++l) {
      double diff =
          Interpolate1D(e0, tabHoleDiffTens[l][0][0], m_eFields, m_intpDiffusion,
                        m_extrLowDiffusion, m_extrHighDiffusion);
      // Apply scaling.
      diff = ScaleDiffusionTensor(diff);
      if (l < 3) {
        cov[l][l] = diff;
      } else if (l == 3) {
        cov[0][1] = cov[1][0] = diff;
      } else if (l == 4) {
        cov[0][2] = cov[2][0] = diff;
      } else if (l == 5) {
        cov[1][2] = cov[2][1] = diff;
      }
    }
  }
 
  return true;
}

◆ HoleTownsend()

bool Garfield::Medium::HoleTownsend	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	alpha
	)

virtual

Reimplemented in Garfield::MediumCdTe, Garfield::MediumGaAs, and Garfield::MediumSilicon.

Definition at line 971 of file Medium.cc.

                                         {
 
  alpha = 0.;
  if (!m_hasHoleTownsend) return false;
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    // Interpolate.
    if (e0 < m_eFields[thrHoleTownsend]) {
      if (!Numerics::Boxin3(tabHoleTownsend, m_bAngles, m_bFields, m_eFields, m_bAngles.size(),
                            m_bFields.size(), m_eFields.size(), ebang, b, e0, alpha, 1)) {
        alpha = -30.;
      }
    } else {
      if (!Numerics::Boxin3(tabHoleTownsend, m_bAngles, m_bFields, m_eFields, m_bAngles.size(),
                            m_bFields.size(), m_eFields.size(), ebang, b, e0, alpha,
                            m_intpTownsend)) {
        alpha = -30.;
      }
    }
  } else {
    // Interpolate.
    if (e0 < m_eFields[thrHoleTownsend]) {
      alpha = Interpolate1D(e0, tabHoleTownsend[0][0], m_eFields, 1,
                            m_extrLowTownsend, m_extrHighTownsend);
    } else {
      alpha = Interpolate1D(e0, tabHoleTownsend[0][0], m_eFields, m_intpTownsend,
                            m_extrLowTownsend, m_extrHighTownsend);
    }
  }
 
  if (alpha < -20.) {
    alpha = 0.;
  } else {
    alpha = exp(alpha);
  }
 
  // Apply scaling.
  alpha = ScaleTownsend(alpha);
  return true;
}

Referenced by Garfield::ViewMedium::EvaluateFunction(), Garfield::MediumCdTe::HoleTownsend(), Garfield::MediumGaAs::HoleTownsend(), and Garfield::MediumSilicon::HoleTownsend().

◆ HoleVelocity()

bool Garfield::Medium::HoleVelocity	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	vx,
		double &	vy,
		double &	vz
	)

virtual

Reimplemented in Garfield::MediumCdTe, Garfield::MediumGaAs, and Garfield::MediumSilicon.

Definition at line 704 of file Medium.cc.

                                                              {
 
  vx = vy = vz = 0.;
  // Make sure there is at least a table of velocities along E.
  if (!m_hasHoleVelocityE) return false;
 
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  // Compute the magnitude of the magnetic field.
  const double b = sqrt(bx * bx + by * by + bz * bz);
 
  // Compute the angle between B field and E field.
  const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
  if (b < Small) {
    // No magnetic field.
    // Calculate the velocity along E.
    double ve = 0.;
    if (m_map2d) {
      if (!Numerics::Boxin3(tabHoleVelocityE, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, ve,
                            m_intpVelocity)) {
        std::cerr << m_className << "::HoleVelocity:\n";
        std::cerr << "    Interpolation of velocity along E failed.\n";
        return false;
      }
    } else {
      ve = Interpolate1D(e0, tabHoleVelocityE[0][0], m_eFields, m_intpVelocity,
                         m_extrLowVelocity, m_extrHighVelocity);
    }
    const double q = 1.;
    const double mu = q * ve / e;
    vx = mu * ex;
    vy = mu * ey;
    vz = mu * ez;
 
  } else if (m_hasHoleVelocityB && m_hasHoleVelocityExB) {
    // Magnetic field, velocities along ExB and Bt available
 
    // Compute unit vectors along E, E x B and Bt.
    double ue[3] = {ex / e, ey / e, ez / e};
    double uexb[3] = {ey * bz - ez * by, ez * bx - ex * bz, ex * by - ey * bx};
    const double exb =
        sqrt(uexb[0] * uexb[0] + uexb[1] * uexb[1] + uexb[2] * uexb[2]);
    if (exb > 0.) {
      uexb[0] /= exb;
      uexb[1] /= exb;
      uexb[2] /= exb;
    } else {
      uexb[0] = ue[0];
      uexb[1] = ue[1];
      uexb[2] = ue[2];
    }
 
    double ubt[3] = {uexb[1] * ez - uexb[2] * ey, uexb[2] * ex - uexb[0] * ez,
                     uexb[0] * ey - uexb[1] * ex};
    const double bt = sqrt(ubt[0] * ubt[0] + ubt[1] * ubt[1] + ubt[2] * ubt[2]);
 
    if (bt > 0.) {
      ubt[0] /= bt;
      ubt[1] /= bt;
      ubt[2] /= bt;
    } else {
      ubt[0] = ue[0];
      ubt[1] = ue[1];
      ubt[2] = ue[2];
    }
 
    // Calculate the velocities in all directions.
    double ve = 0., vbt = 0., vexb = 0.;
    if (m_map2d) {
      if (!Numerics::Boxin3(tabHoleVelocityE, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, ve,
                            m_intpVelocity)) {
        std::cerr << m_className << "::HoleVelocity:\n";
        std::cerr << "    Interpolation of velocity along E failed.\n";
        return false;
      }
      if (!Numerics::Boxin3(tabHoleVelocityExB, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, vexb,
                            m_intpVelocity)) {
        std::cerr << m_className << "::HoleVelocity:\n";
        std::cerr << "    Interpolation of velocity along ExB failed.\n";
        return false;
      }
      if (!Numerics::Boxin3(tabHoleVelocityB, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, vbt,
                            m_intpVelocity)) {
        std::cerr << m_className << "::HoleVelocity:\n";
        std::cerr << "    Interpolation of velocity along Bt failed.\n";
        return false;
      }
    } else {
      ve = Interpolate1D(e0, tabHoleVelocityE[0][0], m_eFields, m_intpVelocity,
                         m_extrLowVelocity, m_extrHighVelocity);
      vbt = Interpolate1D(e0, tabHoleVelocityB[0][0], m_eFields, m_intpVelocity,
                          m_extrLowVelocity, m_extrHighVelocity);
      vexb = Interpolate1D(e0, tabHoleVelocityExB[0][0], m_eFields, m_intpVelocity,
                           m_extrLowVelocity, m_extrHighVelocity);
    }
    const double q = 1.;
    if (ex * bx + ey * by + ez * bz > 0.) vbt = fabs(vbt);
     else vbt = -fabs(vbt);
    vx = q * (ve * ue[0] + q * q * vbt * ubt[0] + q * vexb * uexb[0]);
    vy = q * (ve * ue[1] + q * q * vbt * ubt[1] + q * vexb * uexb[1]);
    vz = q * (ve * ue[2] + q * q * vbt * ubt[2] + q * vexb * uexb[2]);
 
  } else {
    // Magnetic field, velocities along ExB, Bt not available
 
    // Calculate the velocity along E.
    double ve = 0.;
    if (m_map2d) {
      if (!Numerics::Boxin3(tabHoleVelocityE, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, ve,
                            m_intpVelocity)) {
        std::cerr << m_className << "::HoleVelocity:\n";
        std::cerr << "    Interpolation of velocity along E failed.\n";
        return false;
      }
    } else {
      ve = Interpolate1D(e0, tabHoleVelocityE[0][0], m_eFields, m_intpVelocity,
                         m_extrLowVelocity, m_extrHighVelocity);
    }
 
    const double q = 1.;
    const double mu = q * ve / e;
    const double eb = bx * ex + by * ey + bz * ez;
    const double nom = 1. + pow(mu * b, 2);
    vx = mu * (ex + mu * (ey * bz - ez * by) + mu * mu * bx * eb) / nom;
    vy = mu * (ey + mu * (ez * bx - ex * bz) + mu * mu * by * eb) / nom;
    vz = mu * (ez + mu * (ex * by - ey * bx) + mu * mu * bz * eb) / nom;
  }
 
  return true;
}

Referenced by Garfield::ViewMedium::EvaluateFunction(), Garfield::MediumCdTe::HoleVelocity(), Garfield::MediumGaAs::HoleVelocity(), and Garfield::MediumSilicon::HoleVelocity().

◆ InitParamArrays()

void Garfield::Medium::InitParamArrays	(	const unsigned int	eRes,
		const unsigned int	bRes,
		const unsigned int	aRes,
		std::vector< std::vector< std::vector< double > > > &	tab,
		const double	val
	)

protected

Definition at line 2640 of file Medium.cc.

                                                                      {
 
  if (eRes == 0 || bRes == 0 || aRes == 0) {
    std::cerr << m_className << "::InitParamArrays:\n";
    std::cerr << "    Invalid grid.\n";
    return;
  }
 
  tab.assign(aRes, std::vector<std::vector<double> >(bRes, std::vector<double>(eRes, val))); 
  /*
  tab.resize(aRes);
  for (unsigned int i = 0; i < aRes; ++i) {
    tab[i].resize(bRes);
    for (unsigned int j = 0; j < bRes; ++j) {
      tab[i][j].assign(eRes, val);
    }
  }
  */
}

Referenced by CloneTable(), Garfield::MediumMagboltz::GenerateGasTable(), Garfield::MediumGas::LoadGasFile(), and SetIonMobility().

◆ InitParamTensor()

void Garfield::Medium::InitParamTensor	(	const unsigned int	eRes,
		const unsigned int	bRes,
		const unsigned int	aRes,
		const unsigned int	tRes,
		std::vector< std::vector< std::vector< std::vector< double > > > > &	tab,
		const double	val
	)

protected

Definition at line 2663 of file Medium.cc.

                      {
 
  if (eRes == 0 || bRes == 0 || aRes == 0 || tRes == 0) {
    std::cerr << m_className << "::InitParamArrays: Invalid grid.\n";
    return;
  }
 
  tab.resize(tRes);
  for (unsigned int l = 0; l < tRes; ++l) {
    tab[l].resize(aRes);
    for (unsigned int i = 0; i < aRes; ++i) {
      tab[l][i].resize(bRes);
      for (unsigned int j = 0; j < bRes; ++j) {
        tab[l][i][j].assign(eRes, val);
      }
    }
  }
}

Referenced by CloneTensor(), and Garfield::MediumGas::LoadGasFile().

◆ Interpolate1D()

double Garfield::Medium::Interpolate1D	(	const double	e,
		const std::vector< double > &	table,
		const std::vector< double > &	fields,
		const unsigned int	intpMeth,
		const int	jExtr,
		const int	iExtr
	)

protected

Definition at line 2565 of file Medium.cc.

                                                 {
 
  // This function is a generalized version of the Fortran functions
  // GASVEL, GASVT1, GASVT2, GASLOR, GASMOB, GASDFT, and GASDFL
  // for the case of a 1D table. All variables are generic.
 
  const int nSizeTable = fields.size();
 
  if (e < 0. || nSizeTable < 1) return 0.;
 
  double result = 0.;
 
  if (nSizeTable == 1) {
    // Only one point
    result = table[0];
  } else if (e < fields[0]) {
    // Extrapolation towards small fields
    if (fields[0] >= fields[1]) {
      if (m_debug) {
        std::cerr << m_className << "::Interpolate1D:\n";
        std::cerr << "    First two field values coincide.\n";
        std::cerr << "    No extrapolation to lower fields.\n";
      }
      result = table[0];
    } else if (extrLow == 1) {
      // Linear extrapolation
      const double extr4 = (table[1] - table[0]) / (fields[1] - fields[0]);
      const double extr3 = table[0] - extr4 * fields[0];
      result = extr3 + extr4 * e;
    } else if (extrLow == 2) {
      // Logarithmic extrapolation
      const double extr4 = log(table[1] / table[0]) / (fields[1] - fields[0]);
      const double extr3 = log(table[0] - extr4 * fields[0]);
      result = std::exp(std::min(50., extr3 + extr4 * e));
    } else {
      result = table[0];
    }
  } else if (e > fields[nSizeTable - 1]) {
    // Extrapolation towards large fields
    if (fields[nSizeTable - 1] <= fields[nSizeTable - 2]) {
      if (m_debug) {
        std::cerr << m_className << "::Interpolate1D:\n";
        std::cerr << "    Last two field values coincide.\n";
        std::cerr << "    No extrapolation to higher fields.\n";
      }
      result = table[nSizeTable - 1];
    } else if (extrHigh == 1) {
      // Linear extrapolation
      const double extr2 = (table[nSizeTable - 1] - table[nSizeTable - 2]) /
                           (fields[nSizeTable - 1] - fields[nSizeTable - 2]);
      const double extr1 =
          table[nSizeTable - 1] - extr2 * fields[nSizeTable - 1];
      result = extr1 + extr2 * e;
    } else if (extrHigh == 2) {
      // Logarithmic extrapolation
      const double extr2 = log(table[nSizeTable - 1] / table[nSizeTable - 2]) /
                           (fields[nSizeTable - 1] - fields[nSizeTable - 2]);
      const double extr1 =
          log(table[nSizeTable - 1]) - extr2 * fields[nSizeTable - 1];
      result = exp(std::min(50., extr1 + extr2 * e));
    } else {
      result = table[nSizeTable - 1];
    }
  } else {
    // Intermediate points, spline interpolation (not implemented).
    // Intermediate points, Newtonian interpolation
    result = Numerics::Divdif(table, fields, nSizeTable, e, intpMeth);
  }
 
  return result;
}

Referenced by CloneTable(), CloneTensor(), ElectronAttachment(), ElectronDiffusion(), ElectronLorentzAngle(), ElectronTownsend(), ElectronVelocity(), HoleAttachment(), HoleDiffusion(), HoleTownsend(), HoleVelocity(), IonDiffusion(), IonDissociation(), IonVelocity(), and SetIonMobility().

◆ IonDiffusion()

bool Garfield::Medium::IonDiffusion	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	dl,
		double &	dt
	)

virtual

Definition at line 1122 of file Medium.cc.

                                                  {
 
  dl = dt = 0.;
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    // Interpolate.
    if (m_hasIonDiffLong) {
      if (!Numerics::Boxin3(tabIonDiffLong, m_bAngles, m_bFields, m_eFields, m_bAngles.size(),
                            m_bFields.size(), m_eFields.size(), ebang, b, e0, dl,
                            m_intpDiffusion)) {
        dl = 0.;
      }
    }
    if (m_hasIonDiffTrans) {
      if (!Numerics::Boxin3(tabIonDiffTrans, m_bAngles, m_bFields, m_eFields, m_bAngles.size(),
                            m_bFields.size(), m_eFields.size(), ebang, b, e0, dt,
                            m_intpDiffusion)) {
        dt = 0.;
      }
    }
  } else {
    if (m_hasIonDiffLong) {
      dl = Interpolate1D(e0, tabIonDiffLong[0][0], m_eFields, m_intpDiffusion,
                         m_extrLowDiffusion, m_extrHighDiffusion);
    }
    if (m_hasIonDiffTrans) {
      dt = Interpolate1D(e0, tabIonDiffTrans[0][0], m_eFields, m_intpDiffusion,
                         m_extrLowDiffusion, m_extrHighDiffusion);
    }
  }
 
  // If no data available, calculate
  // the diffusion coefficients using the Einstein relation
  if (!m_hasIonDiffLong) {
    dl = sqrt(2. * BoltzmannConstant * m_temperature / e);
  }
  if (!m_hasIonDiffTrans) {
    dt = sqrt(2. * BoltzmannConstant * m_temperature / e);
  }
 
  return true;
}

Referenced by Garfield::ViewMedium::EvaluateFunction().

◆ IonDissociation()

bool Garfield::Medium::IonDissociation	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	diss
	)

virtual

Definition at line 1176 of file Medium.cc.

                                           {
 
  diss = 0.;
  if (!m_hasIonDissociation) return false;
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
 
  if (m_map2d) {
    // Compute the magnitude of the magnetic field.
    const double b = sqrt(bx * bx + by * by + bz * bz);
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
    // Interpolate.
    if (e0 < m_eFields[thrIonDissociation]) {
      if (!Numerics::Boxin3(tabIonDissociation, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, diss,
                            1)) {
        diss = -30.;
      }
    } else {
      if (!Numerics::Boxin3(tabIonDissociation, m_bAngles, m_bFields, m_eFields,
                            m_bAngles.size(), m_bFields.size(), m_eFields.size(), ebang, b, e0, diss,
                            m_intpDissociation)) {
        diss = -30.;
      }
    }
  } else {
    // Interpolate.
    if (e0 < m_eFields[thrIonDissociation]) {
      diss = Interpolate1D(e0, tabIonDissociation[0][0], m_eFields, 1,
                           m_extrLowDissociation, m_extrHighDissociation);
    } else {
      diss = Interpolate1D(e0, tabHoleTownsend[0][0], m_eFields, 
                           m_intpDissociation,
                           m_extrLowDissociation, m_extrHighDissociation);
    }
  }
 
  if (diss < -20.) {
    diss = 0.;
  } else {
    diss = exp(diss);
  }
 
  // Apply scaling.
  diss = ScaleDissociation(diss);
  return true;
}

◆ IonVelocity()

bool Garfield::Medium::IonVelocity	(	const double	ex,
		const double	ey,
		const double	ez,
		const double	bx,
		const double	by,
		const double	bz,
		double &	vx,
		double &	vy,
		double &	vz
	)

virtual

Definition at line 1077 of file Medium.cc.

                                                             {
 
  vx = vy = vz = 0.;
  if (!m_hasIonMobility) return false;
  // Compute the magnitude of the electric field.
  const double e = sqrt(ex * ex + ey * ey + ez * ez);
  const double e0 = ScaleElectricField(e);
  if (e < Small || e0 < Small) return true;
  // Compute the magnitude of the electric field.
  const double b = sqrt(bx * bx + by * by + bz * bz);
 
  double mu = 0.;
  if (m_map2d) {
    // Compute the angle between B field and E field.
    const double ebang = GetAngle(ex, ey, ez, bx, by, bz, e, b);
 
    if (!Numerics::Boxin3(tabIonMobility, m_bAngles, m_bFields, m_eFields, m_bAngles.size(),
                          m_bFields.size(), m_eFields.size(), ebang, b, e0, mu, 
                          m_intpMobility)) {
      mu = 0.;
    }
  } else {
    mu = Interpolate1D(e0, tabIonMobility[0][0], m_eFields, m_intpMobility,
                       m_extrLowMobility, m_extrHighMobility);
  }
 
  const double q = 1.;
  mu *= q;
  if (b < Small) {
    vx = mu * ex;
    vy = mu * ey;
    vz = mu * ez;
  } else {
    const double eb = bx * ex + by * ey + bz * ez;
    const double nom = 1. + pow(mu * b, 2);
    vx = mu * (ex + mu * (ey * bz - ez * by) + mu * mu * bx * eb) / nom;
    vy = mu * (ey + mu * (ez * bx - ex * bz) + mu * mu * by * eb) / nom;
    vz = mu * (ez + mu * (ex * by - ey * bx) + mu * mu * bz * eb) / nom;
  }
 
  return true;
}

Referenced by Garfield::ViewMedium::EvaluateFunction().

◆ IsDriftable()

bool Garfield::Medium::IsDriftable ( ) const

inline

Definition at line 57 of file Medium.hh.

57{ return m_driftable; }

Referenced by Garfield::ComponentAnsys121::ElectricField(), Garfield::ComponentAnsys123::ElectricField(), Garfield::ComponentComsol::ElectricField(), Garfield::ComponentElmer::ElectricField(), Garfield::ComponentAnalyticField::ElectricField(), Garfield::ComponentConstant::ElectricField(), and Garfield::ComponentUser::ElectricField().

◆ IsGas()

virtual bool Garfield::Medium::IsGas ( ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 23 of file Medium.hh.

23{ return false; }

Referenced by Garfield::TrackElectron::NewTrack(), and Garfield::ViewGeometry::Plot().

◆ IsIonisable()

bool Garfield::Medium::IsIonisable ( ) const

inline

Definition at line 59 of file Medium.hh.

59{ return m_ionisable; }

Referenced by Garfield::TrackBichsel::GetCluster(), Garfield::TrackElectron::GetCluster(), Garfield::TrackPAI::GetCluster(), Heed::HeedFieldMap::inside(), Garfield::TrackBichsel::NewTrack(), Garfield::TrackElectron::NewTrack(), Garfield::TrackHeed::NewTrack(), Garfield::TrackPAI::NewTrack(), Garfield::TrackSrim::NewTrack(), Garfield::TrackHeed::TransportDeltaElectron(), and Garfield::TrackHeed::TransportPhoton().

◆ IsMicroscopic()

bool Garfield::Medium::IsMicroscopic ( ) const

inline

Definition at line 58 of file Medium.hh.

58{ return m_microscopic; }

◆ IsSemiconductor()

virtual bool Garfield::Medium::IsSemiconductor ( ) const

inlinevirtual

Reimplemented in Garfield::MediumCdTe, Garfield::MediumGaAs, and Garfield::MediumSilicon.

Definition at line 24 of file Medium.hh.

24{ return false; }

Referenced by Garfield::ViewGeometry::Plot().

◆ ResetElectronAttachment()

void Garfield::Medium::ResetElectronAttachment ( )

Definition at line 1331 of file Medium.cc.

                                     {
 
  tabElectronAttachment.clear();
  m_hasElectronAttachment = false;
}

◆ ResetElectronDiffusion()

void Garfield::Medium::ResetElectronDiffusion ( )

Definition at line 1316 of file Medium.cc.

                                    {
 
  tabElectronDiffLong.clear();
  tabElectronDiffTrans.clear();
  tabElectronDiffTens.clear();
  m_hasElectronDiffLong = false;
  m_hasElectronDiffTrans = false;
  m_hasElectronDiffTens = false;
}

◆ ResetElectronLorentzAngle()

void Garfield::Medium::ResetElectronLorentzAngle ( )

Definition at line 1337 of file Medium.cc.

                                       {
 
  tabElectronLorentzAngle.clear();
  m_hasElectronLorentzAngle = false;
}

◆ ResetElectronTownsend()

void Garfield::Medium::ResetElectronTownsend ( )

Definition at line 1326 of file Medium.cc.

                                   {
 
  tabElectronTownsend.clear();
}

◆ ResetElectronVelocity()

void Garfield::Medium::ResetElectronVelocity ( )

Definition at line 1306 of file Medium.cc.

                                   {
 
  tabElectronVelocityE.clear();
  tabElectronVelocityB.clear();
  tabElectronVelocityExB.clear();
  m_hasElectronVelocityE = false;
  m_hasElectronVelocityB = false;
  m_hasElectronVelocityExB = false;
}

◆ ResetHoleAttachment()

void Garfield::Medium::ResetHoleAttachment ( )

Definition at line 1369 of file Medium.cc.

                                 {
 
  tabHoleAttachment.clear();
  m_hasHoleAttachment = false;
}

◆ ResetHoleDiffusion()

void Garfield::Medium::ResetHoleDiffusion ( )

Definition at line 1353 of file Medium.cc.

                                {
 
  tabHoleDiffLong.clear();
  tabHoleDiffTrans.clear();
  tabHoleDiffTens.clear();
  m_hasHoleDiffLong = false;
  m_hasHoleDiffTrans = false;
  m_hasHoleDiffTens = false;
}

◆ ResetHoleTownsend()

void Garfield::Medium::ResetHoleTownsend ( )

Definition at line 1363 of file Medium.cc.

                               {
 
  tabHoleTownsend.clear();
  m_hasHoleTownsend = false;
}

◆ ResetHoleVelocity()

void Garfield::Medium::ResetHoleVelocity ( )

Definition at line 1343 of file Medium.cc.

                               {
 
  tabHoleVelocityE.clear();
  tabHoleVelocityB.clear();
  tabHoleVelocityExB.clear();
  m_hasHoleVelocityE = false;
  m_hasHoleVelocityB = false;
  m_hasHoleVelocityExB = false;
}

◆ ResetIonDiffusion()

void Garfield::Medium::ResetIonDiffusion ( )

Definition at line 1381 of file Medium.cc.

                               {
 
  tabIonDiffLong.clear();
  tabIonDiffTrans.clear();
  m_hasIonDiffLong = false;
  m_hasIonDiffTrans = false;
}

◆ ResetIonDissociation()

void Garfield::Medium::ResetIonDissociation ( )

Definition at line 1389 of file Medium.cc.

                                  {
 
  tabIonDissociation.clear();
  m_hasIonDissociation = false;
}

◆ ResetIonMobility()

void Garfield::Medium::ResetIonMobility ( )

Definition at line 1375 of file Medium.cc.

                              {
 
  tabIonMobility.clear();
  m_hasIonMobility = false;
}

Referenced by SetIonMobility().

◆ ScaleAttachment()

virtual double Garfield::Medium::ScaleAttachment ( const double eta ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 270 of file Medium.hh.

270{ return eta; }

Referenced by ElectronAttachment(), and HoleAttachment().

◆ ScaleDiffusion()

virtual double Garfield::Medium::ScaleDiffusion ( const double d ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 267 of file Medium.hh.

267{ return d; }

Referenced by ElectronDiffusion(), and HoleDiffusion().

◆ ScaleDiffusionTensor()

virtual double Garfield::Medium::ScaleDiffusionTensor ( const double d ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 268 of file Medium.hh.

268{ return d; }

Referenced by ElectronDiffusion(), and HoleDiffusion().

◆ ScaleDissociation()

virtual double Garfield::Medium::ScaleDissociation ( const double diss ) const

inlinevirtual

Definition at line 272 of file Medium.hh.

272{ return diss; }

Referenced by IonDissociation().

◆ ScaleElectricField()

virtual double Garfield::Medium::ScaleElectricField ( const double e ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 264 of file Medium.hh.

264{ return e; }

Referenced by ElectronAttachment(), ElectronDiffusion(), ElectronLorentzAngle(), ElectronTownsend(), ElectronVelocity(), HoleAttachment(), HoleDiffusion(), HoleTownsend(), HoleVelocity(), IonDiffusion(), IonDissociation(), and IonVelocity().

◆ ScaleLorentzAngle()

virtual double Garfield::Medium::ScaleLorentzAngle ( const double lor ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 271 of file Medium.hh.

271{ return lor; }

Referenced by ElectronLorentzAngle().

◆ ScaleTownsend()

virtual double Garfield::Medium::ScaleTownsend ( const double alpha ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 269 of file Medium.hh.

269{ return alpha; }

Referenced by ElectronTownsend(), and HoleTownsend().

◆ ScaleVelocity()

virtual double Garfield::Medium::ScaleVelocity ( const double v ) const

inlinevirtual

Definition at line 266 of file Medium.hh.

266{ return v; }

◆ SetAtomicNumber()

void Garfield::Medium::SetAtomicNumber ( const double z )

virtual

Reimplemented in Garfield::MediumGas.

Definition at line 154 of file Medium.cc.

                                           {
 
  if (z < 1.) {
    std::cerr << m_className << "::SetAtomicNumber:\n";
    std::cerr << "    Atomic number must be >= 1.\n";
    return;
  }
  m_z = z;
  m_isChanged = true;
}

Referenced by Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumGaAs::MediumGaAs(), and Garfield::MediumSilicon::MediumSilicon().

◆ SetAtomicWeight()

void Garfield::Medium::SetAtomicWeight ( const double a )

virtual

Reimplemented in Garfield::MediumGas.

Definition at line 165 of file Medium.cc.

                                           {
 
  if (a <= 0.) {
    std::cerr << m_className << "::SetAtomicWeight:\n";
    std::cerr << "    Atomic weight must be greater than zero.\n";
    return;
  }
  m_a = a;
  m_isChanged = true;
}

Referenced by Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumGaAs::MediumGaAs(), and Garfield::MediumSilicon::MediumSilicon().

◆ SetDielectricConstant()

void Garfield::Medium::SetDielectricConstant ( const double eps )

Definition at line 126 of file Medium.cc.

                                                   {
 
  if (eps < 1.) {
    std::cerr << m_className << "::SetDielectricConstant:\n";
    std::cerr << "    Dielectric constant must be >= 1.\n";
    return;
  }
  m_epsilon = eps;
  m_isChanged = true;
}

Referenced by Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumGaAs::MediumGaAs(), and Garfield::MediumSilicon::MediumSilicon().

◆ SetExtrapolationMethodAttachment()

void Garfield::Medium::SetExtrapolationMethodAttachment	(	const std::string &	extrLow,
		const std::string &	extrHigh
	)

Definition at line 2397 of file Medium.cc.

                                                                         {
 
  unsigned int iExtr = 0;
  if (GetExtrapolationIndex(extrLow, iExtr)) {
    m_extrLowAttachment = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodAttachment:\n";
    std::cerr << "    Unknown extrapolation method (" << extrLow << ")\n";
  }
 
  if (GetExtrapolationIndex(extrHigh, iExtr)) {
    m_extrHighAttachment = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodAttachment:\n";
    std::cerr << "    Unknown extrapolation method (" << extrHigh << ")\n";
  }
}

◆ SetExtrapolationMethodDiffusion()

void Garfield::Medium::SetExtrapolationMethodDiffusion	(	const std::string &	extrLow,
		const std::string &	extrHigh
	)

Definition at line 2359 of file Medium.cc.

                                                                        {
 
  unsigned int iExtr = 0;
  if (GetExtrapolationIndex(extrLow, iExtr)) {
    m_extrLowDiffusion = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodDiffusion:\n";
    std::cerr << "    Unknown extrapolation method (" << extrLow << ")\n";
  }
 
  if (GetExtrapolationIndex(extrHigh, iExtr)) {
    m_extrHighDiffusion = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodDiffusion:\n";
    std::cerr << "    Unknown extrapolation method (" << extrHigh << ")\n";
  }
}

◆ SetExtrapolationMethodIonDissociation()

void Garfield::Medium::SetExtrapolationMethodIonDissociation	(	const std::string &	extrLow,
		const std::string &	extrHigh
	)

Definition at line 2434 of file Medium.cc.

                                                                              {
 
  unsigned int iExtr = 0;
  if (GetExtrapolationIndex(extrLow, iExtr)) {
    m_extrLowDissociation = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodIonDissociation:\n";
    std::cerr << "    Unknown extrapolation method (" << extrLow << ")\n";
  }
 
  if (GetExtrapolationIndex(extrHigh, iExtr)) {
    m_extrHighDissociation = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodIonDissociation:\n";
    std::cerr << "    Unknown extrapolation method (" << extrHigh << ")\n";
  }
}

◆ SetExtrapolationMethodIonMobility()

void Garfield::Medium::SetExtrapolationMethodIonMobility	(	const std::string &	extrLow,
		const std::string &	extrHigh
	)

Definition at line 2416 of file Medium.cc.

                                                                          {
 
  unsigned int iExtr = 0;
  if (GetExtrapolationIndex(extrLow, iExtr)) {
    m_extrLowMobility = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodIonMobility:\n";
    std::cerr << "    Unknown extrapolation method (" << extrLow << ")\n";
  }
  if (GetExtrapolationIndex(extrHigh, iExtr)) {
    m_extrHighMobility = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodIonMobility:\n";
    std::cerr << "    Unknown extrapolation method (" << extrHigh << ")\n";
  }
}

◆ SetExtrapolationMethodTownsend()

void Garfield::Medium::SetExtrapolationMethodTownsend	(	const std::string &	extrLow,
		const std::string &	extrHigh
	)

Definition at line 2378 of file Medium.cc.

                                                                       {
 
  unsigned int iExtr = 0;
  if (GetExtrapolationIndex(extrLow, iExtr)) {
    m_extrLowTownsend = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodTownsend:\n";
    std::cerr << "    Unknown extrapolation method (" << extrLow << ")\n";
  }
 
  if (GetExtrapolationIndex(extrHigh, iExtr)) {
    m_extrHighTownsend = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodTownsend:\n";
    std::cerr << "    Unknown extrapolation method (" << extrHigh << ")\n";
  }
}

◆ SetExtrapolationMethodVelocity()

void Garfield::Medium::SetExtrapolationMethodVelocity	(	const std::string &	extrLow,
		const std::string &	extrHigh
	)

Definition at line 2340 of file Medium.cc.

                                                                       {
 
  unsigned int iExtr = 0;
  if (GetExtrapolationIndex(extrLow, iExtr)) {
    m_extrLowVelocity = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodVelocity:\n";
    std::cerr << "    Unknown extrapolation method (" << extrLow << ")\n";
  }
 
  if (GetExtrapolationIndex(extrHigh, iExtr)) {
    m_extrHighVelocity = iExtr;
  } else {
    std::cerr << m_className << "::SetExtrapolationMethodVelocity:\n";
    std::cerr << "    Unknown extrapolation method (" << extrHigh << ")\n";
  }
}

◆ SetFanoFactor()

void Garfield::Medium::SetFanoFactor ( const double f )

inline

Definition at line 64 of file Medium.hh.

64{ m_fano = f; }

◆ SetFieldGrid() [1/2]

void Garfield::Medium::SetFieldGrid	(	const std::vector< double > &	efields,
		const std::vector< double > &	bfields,
		const std::vector< double > &	angles
	)

Definition at line 1506 of file Medium.cc.

                                                           {
 
  if (efields.empty()) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Number of E-fields must be > 0.\n";
    return;
  }
  if (bfields.empty()) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Number of B-fields must be > 0.\n";
    return;
  }
  if (angles.empty()) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Number of angles must be > 0.\n";
    return;
  }
 
  // Make sure the values are not negative.
  if (efields[0] < 0.) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    E-fields must be >= 0.\n";
  }
  if (bfields[0] < 0.) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    B-fields must be >= 0.\n";
  }
  if (angles[0] < 0.) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Angles must be >= 0.\n";
  }
 
  const unsigned int nEfieldsNew = efields.size();
  const unsigned int nBfieldsNew = bfields.size();
  const unsigned int nAnglesNew = angles.size();
  // Make sure the values are in strictly monotonic, ascending order.
  if (nEfieldsNew > 1) {
    for (unsigned int i = 1; i < nEfieldsNew; ++i) {
      if (efields[i] <= efields[i - 1]) {
        std::cerr << m_className << "::SetFieldGrid:\n";
        std::cerr << "    E-fields are not in ascending order.\n";
        return;
      }
    }
  }
  if (nBfieldsNew > 1) {
    for (unsigned int i = 1; i < nBfieldsNew; ++i) {
      if (bfields[i] <= bfields[i - 1]) {
        std::cerr << m_className << "::SetFieldGrid:\n";
        std::cerr << "    B-fields are not in ascending order.\n";
        return;
      }
    }
  }
  if (nAnglesNew > 1) {
    for (unsigned int i = 1; i < nAnglesNew; ++i) {
      if (angles[i] <= angles[i - 1]) {
        std::cerr << m_className << "::SetFieldGrid:\n";
        std::cerr << "    Angles are not in ascending order.\n";
        return;
      }
    }
  }
 
  if (m_debug) {
    std::cout << m_className << "::SetFieldGrid:\n";
    std::cout << "    E-fields:\n";
    for (unsigned int i = 0; i < nEfieldsNew; ++i) {
      std::cout << "      " << efields[i] << "\n";
    }
    std::cout << "    B-fields:\n";
    for (unsigned int i = 0; i < nBfieldsNew; ++i) {
      std::cout << "      " << bfields[i] << "\n";
    }
    std::cout << "    Angles:\n";
    for (unsigned int i = 0; i < nAnglesNew; ++i) {
      std::cout << "      " << angles[i] << "\n";
    }
  }
 
  // Clone the existing tables.
  // Electrons
  CloneTable(tabElectronVelocityE, efields, bfields, angles, m_intpVelocity,
             m_extrLowVelocity, m_extrHighVelocity, 0.,
             "electron velocity along E");
  CloneTable(tabElectronVelocityB, efields, bfields, angles, m_intpVelocity,
             m_extrLowVelocity, m_extrHighVelocity, 0.,
             "electron velocity along Bt");
  CloneTable(tabElectronVelocityExB, efields, bfields, angles, m_intpVelocity,
             m_extrLowVelocity, m_extrHighVelocity, 0.,
             "electron velocity along ExB");
  CloneTable(tabElectronDiffLong, efields, bfields, angles, m_intpDiffusion,
             m_extrLowDiffusion, m_extrHighDiffusion, 0.,
             "electron longitudinal diffusion");
  CloneTable(tabElectronDiffTrans, efields, bfields, angles, m_intpDiffusion,
             m_extrLowDiffusion, m_extrHighDiffusion, 0.,
             "electron transverse diffusion");
  CloneTable(tabElectronTownsend, efields, bfields, angles, m_intpTownsend,
             m_extrLowTownsend, m_extrHighTownsend, -30.,
             "electron Townsend coefficient");
  CloneTable(tabElectronAttachment, efields, bfields, angles, m_intpAttachment,
             m_extrLowAttachment, m_extrHighAttachment, -30.,
             "electron attachment coefficient");
  CloneTable(tabElectronLorentzAngle, efields, bfields, angles, m_intpLorentzAngle,
             m_extrLowLorentzAngle, m_extrHighLorentzAngle, 0.,
             "electron attachment coefficient");
  if (m_hasElectronDiffTens) {
    CloneTensor(tabElectronDiffTens, 6, efields, bfields, angles, m_intpDiffusion,
                m_extrLowDiffusion, m_extrHighDiffusion, 0.,
                "electron diffusion tensor");
  }
 
  // Holes
  CloneTable(tabHoleVelocityE, efields, bfields, angles, m_intpVelocity,
             m_extrLowVelocity, m_extrHighVelocity, 0., 
             "hole velocity along E");
  CloneTable(tabHoleVelocityB, efields, bfields, angles, m_intpVelocity,
             m_extrLowVelocity, m_extrHighVelocity, 0., 
             "hole velocity along Bt");
  CloneTable(tabHoleVelocityExB, efields, bfields, angles, m_intpVelocity,
             m_extrLowVelocity, m_extrHighVelocity, 0.,
             "hole velocity along ExB");
  CloneTable(tabHoleDiffLong, efields, bfields, angles, m_intpDiffusion,
             m_extrLowDiffusion, m_extrHighDiffusion, 0.,
             "hole longitudinal diffusion");
  CloneTable(tabHoleDiffTrans, efields, bfields, angles, m_intpDiffusion,
             m_extrLowDiffusion, m_extrHighDiffusion, 0.,
             "hole transverse diffusion");
  CloneTable(tabHoleTownsend, efields, bfields, angles, m_intpTownsend,
             m_extrLowTownsend, m_extrHighTownsend, -30.,
             "hole Townsend coefficient");
  CloneTable(tabHoleAttachment, efields, bfields, angles, m_intpAttachment,
             m_extrLowAttachment, m_extrHighAttachment, -30.,
             "hole attachment coefficient");
  if (m_hasHoleDiffTens) {
    CloneTensor(tabHoleDiffTens, 6, efields, bfields, angles, m_intpDiffusion,
                m_extrLowDiffusion, m_extrHighDiffusion, 0.,
                "hole diffusion tensor");
  }
 
  // Ions
  CloneTable(tabIonMobility, efields, bfields, angles, m_intpMobility,
             m_extrLowMobility, m_extrHighMobility, 0., 
             "ion mobility");
  CloneTable(tabIonDiffLong, efields, bfields, angles, m_intpDiffusion,
             m_extrLowDiffusion, m_extrHighDiffusion, 0.,
             "ion longitudinal diffusion");
  CloneTable(tabIonDiffTrans, efields, bfields, angles, m_intpDiffusion,
             m_extrLowDiffusion, m_extrHighDiffusion, 0.,
             "ion transverse diffusion");
  CloneTable(tabIonDissociation, efields, bfields, angles, m_intpDissociation,
             m_extrLowDissociation, m_extrHighDissociation, -30.,
             "ion dissociation");
 
  if (nBfieldsNew > 1 || nAnglesNew > 1) m_map2d = true;
  m_eFields.resize(nEfieldsNew);
  m_bFields.resize(nBfieldsNew);
  m_bAngles.resize(nAnglesNew);
  m_eFields = efields;
  m_bFields = bfields;
  m_bAngles = angles;
}

◆ SetFieldGrid() [2/2]

void Garfield::Medium::SetFieldGrid	(	double	emin,
		double	emax,
		int	ne,
		bool	logE,
		double	bmin = `0.`,
		double	bmax = `0.`,
		int	nb = `1`,
		double	amin = `0.`,
		double	amax = `0.`,
		int	na = `1`
	)

Definition at line 1395 of file Medium.cc.

                                               {
 
  // Check if the requested E-field range makes sense.
  if (ne <= 0) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Number of E-fields must be > 0.\n";
    return;
  }
 
  if (emin < 0. || emax < 0.) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Electric fields must be positive.\n";
    return;
  }
 
  if (emax < emin) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Swapping min./max. E-field.\n";
    const double etemp = emin;
    emin = emax;
    emax = etemp;
  }
 
  double estep = 0.;
  if (logE) {
    // Logarithmic scale
    if (emin < Small) {
      std::cerr << m_className << "::SetFieldGrid:\n";
      std::cerr << "    Min. E-field must be non-zero for log. scale.\n";
      return;
    }
    if (ne == 1) {
      std::cerr << m_className << "::SetFieldGrid:\n";
      std::cerr << "    Number of E-fields must be > 1 for log. scale.\n";
      return;
    }
    estep = pow(emax / emin, 1. / (ne - 1.));
  } else {
    // Linear scale
    if (ne > 1) estep = (emax - emin) / (ne - 1.);
  }
 
  // Check if the requested B-field range makes sense.
  if (nb <= 0) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Number of B-fields must be > 0.\n";
    return;
  }
  if (bmax < 0. || bmin < 0.) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Magnetic fields must be positive.\n";
    return;
  }
  if (bmax < bmin) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Swapping min./max. B-field.\n";
    const double btemp = bmin;
    bmin = bmax;
    bmax = btemp;
  }
 
  double bstep = 0.;
  if (nb > 1) bstep = (bmax - bmin) / (nb - 1.);
 
  // Check if the requested angular range makes sense.
  if (na <= 0) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Number of angles must be > 0.\n";
    return;
  }
  if (amax < 0. || amin < 0.) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Angles must be positive.\n";
    return;
  }
  if (amax < amin) {
    std::cerr << m_className << "::SetFieldGrid:\n";
    std::cerr << "    Swapping min./max. angle.\n";
    const double atemp = amin;
    amin = amax;
    amax = atemp;
  }
  double astep = 0.;
  if (na > 1) astep = (amax - amin) / (na - 1.);
 
  // Setup the field grids.
  std::vector<double> eFieldsNew(ne);
  std::vector<double> bFieldsNew(nb);
  std::vector<double> bAnglesNew(na);
  for (int i = 0; i < ne; ++i) {
    if (logE) {
      eFieldsNew[i] = emin * pow(estep, i);
    } else {
      eFieldsNew[i] = emin + i * estep;
    }
  }
  for (int i = 0; i < nb; ++i) {
    bFieldsNew[i] = bmin + i * bstep;
  }
  if (na == 1 && nb == 1 && fabs(bmin) < 1.e-4) {
    bAnglesNew[0] = HalfPi;
  } else {
    for (int i = 0; i < na; ++i) {
      bAnglesNew[i] = amin + i * astep;
    }
  }
  SetFieldGrid(eFieldsNew, bFieldsNew, bAnglesNew);
}

Referenced by Medium(), and SetFieldGrid().

◆ SetInterpolationMethodAttachment()

void Garfield::Medium::SetInterpolationMethodAttachment ( const unsigned int intrp )

Definition at line 2488 of file Medium.cc.

                                                                      {
 
  if (intrp > 0) m_intpAttachment = intrp;
}

◆ SetInterpolationMethodDiffusion()

void Garfield::Medium::SetInterpolationMethodDiffusion ( const unsigned int intrp )

Definition at line 2478 of file Medium.cc.

                                                                     {
 
  if (intrp > 0) m_intpDiffusion = intrp;
}

◆ SetInterpolationMethodIonDissociation()

void Garfield::Medium::SetInterpolationMethodIonDissociation ( const unsigned int intrp )

Definition at line 2498 of file Medium.cc.

                                                                           {
 
  if (intrp > 0) m_intpDissociation = intrp;
}

◆ SetInterpolationMethodIonMobility()

void Garfield::Medium::SetInterpolationMethodIonMobility ( const unsigned int intrp )

Definition at line 2493 of file Medium.cc.

                                                                       {
 
  if (intrp > 0) m_intpMobility = intrp;
}

◆ SetInterpolationMethodTownsend()

void Garfield::Medium::SetInterpolationMethodTownsend ( const unsigned int intrp )

Definition at line 2483 of file Medium.cc.

                                                                    {
 
  if (intrp > 0) m_intpTownsend = intrp;
}

◆ SetInterpolationMethodVelocity()

void Garfield::Medium::SetInterpolationMethodVelocity ( const unsigned int intrp )

Definition at line 2473 of file Medium.cc.

                                                                    {
 
  if (intrp > 0) m_intpVelocity = intrp;
}

◆ SetIonMobility() [1/2]

bool Garfield::Medium::SetIonMobility	(	const std::vector< double > &	fields,
		const std::vector< double > &	mobilities
	)

Definition at line 2303 of file Medium.cc.

                                                                 {
 
  const int ne = efields.size();
  const int nm = mobilities.size();
  if (ne != nm) {
    std::cerr << m_className << "::SetIonMobility:\n";
    std::cerr << "    E-field and mobility arrays have different sizes.\n";
    return false;
  }
 
  ResetIonMobility();
  const unsigned int nEfields = m_eFields.size();
  const unsigned int nBfields = m_bFields.size();
  const unsigned int nAngles = m_bAngles.size();
  InitParamArrays(nEfields, nBfields, nAngles, tabIonMobility, 0.);
  for (unsigned int i = 0; i < nEfields; ++i) {
    const double e = m_eFields[i];
    const double mu = Interpolate1D(e, mobilities, efields, m_intpMobility,
                                    m_extrLowMobility, m_extrHighMobility);
    tabIonMobility[0][0][i] = mu;
  }
 
  if (m_map2d) {
    for (unsigned int i = 0; i < nAngles; ++i) {
      for (unsigned int j = 0; j< nBfields; ++j) {
        for (unsigned int k = 0; k < nEfields; ++k) {
          tabIonMobility[i][j][k] = tabIonMobility[0][0][k];
        }
      }
    }
  }
  m_hasIonMobility = true;
 
  return true;
}

◆ SetIonMobility() [2/2]

bool Garfield::Medium::SetIonMobility	(	const unsigned int	ie,
		const unsigned int	ib,
		const unsigned int	ia,
		const double	mu
	)

Definition at line 2269 of file Medium.cc.

                                                                    {
 
  // Check the index.
  if (ie >= m_eFields.size() || ib >= m_bFields.size() || 
      ia >= m_bAngles.size()) {
    std::cerr << m_className << "::SetIonMobility:\n";
    std::cerr << "    Index (" << ie << ", " << ib << ", " << ia
              << ") out of range.\n";
    return false;
  }
 
  if (!m_hasIonMobility) {
    std::cerr << m_className << "::SetIonMobility:\n";
    std::cerr << "    Ion mobility table not initialised.\n";
    return false;
  }
 
  if (mu == 0.) {
    std::cerr << m_className << "::SetIonMobility:\n";
    std::cerr << "    Zero value not permitted.\n";
    return false;
  }
 
  tabIonMobility[ia][ib][ie] = mu;
  if (m_debug) {
    std::cout << m_className << "::SetIonMobility:\n";
    std::cout << "   Ion mobility at E = " << m_eFields[ie]
              << " V/cm, B = " << m_bFields[ib] << " T, angle " << m_bAngles[ia]
              << " set to " << mu << " cm2/(V ns).\n";
  }
  return true;
}

Referenced by Garfield::MediumGas::LoadIonMobility().

◆ SetMassDensity()

void Garfield::Medium::SetMassDensity ( const double rho )

virtual

Reimplemented in Garfield::MediumGas.

Definition at line 187 of file Medium.cc.

                                            {
 
  if (rho <= 0.) {
    std::cerr << m_className << "::SetMassDensity:\n";
    std::cerr << "    Density [g/cm3] must be greater than zero.\n";
    return;
  }
 
  if (m_a <= 0.) {
    std::cerr << m_className << "::SetMassDensity:\n";
    std::cerr << "    Atomic weight is not defined.\n";
    return;
  }
  m_density = rho / (AtomicMassUnit * m_a);
  m_isChanged = true;
}

Referenced by Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumGaAs::MediumGaAs(), and Garfield::MediumSilicon::MediumSilicon().

◆ SetNumberDensity()

void Garfield::Medium::SetNumberDensity ( const double n )

virtual

Reimplemented in Garfield::MediumGas.

Definition at line 176 of file Medium.cc.

                                            {
 
  if (n <= 0.) {
    std::cerr << m_className << "::SetNumberDensity:\n";
    std::cerr << "    Density [cm-3] must be greater than zero.\n";
    return;
  }
  m_density = n;
  m_isChanged = true;
}

◆ SetPressure()

void Garfield::Medium::SetPressure ( const double p )

Definition at line 115 of file Medium.cc.

                                       {
 
  if (p <= 0.) {
    std::cerr << m_className << "::SetPressure:\n";
    std::cerr << "    Pressure [Torr] must be greater than zero.\n";
    return;
  }
  m_pressure = p;
  m_isChanged = true;
}

Referenced by GarfieldPhysics::InitializePhysics().

◆ SetTemperature()

void Garfield::Medium::SetTemperature ( const double t )

Definition at line 104 of file Medium.cc.

                                          {
 
  if (t <= 0.) {
    std::cerr << m_className << "::SetTemperature:\n";
    std::cerr << "    Temperature [K] must be greater than zero.\n";
    return;
  }
  m_temperature = t;
  m_isChanged = true;
}

Referenced by GarfieldPhysics::InitializePhysics(), Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumGaAs::MediumGaAs(), and Garfield::MediumSilicon::MediumSilicon().

◆ SetW()

void Garfield::Medium::SetW ( const double w )

inline

Definition at line 62 of file Medium.hh.

62{ m_w = w; }

◆ UnScaleElectricField()

virtual double Garfield::Medium::UnScaleElectricField ( const double e ) const

inlinevirtual

Reimplemented in Garfield::MediumGas.

Definition at line 265 of file Medium.hh.

265{ return e; }

Member Data Documentation

◆ m_a

double Garfield::Medium::m_a

protected

Definition at line 313 of file Medium.hh.

Referenced by GetAtomicWeight(), GetMassDensity(), SetAtomicWeight(), and SetMassDensity().

◆ m_bAngles

std::vector<double> Garfield::Medium::m_bAngles

protected

Definition at line 334 of file Medium.hh.

Referenced by CloneTable(), CloneTensor(), ElectronAttachment(), ElectronDiffusion(), ElectronLorentzAngle(), ElectronTownsend(), ElectronVelocity(), Garfield::MediumMagboltz::GenerateGasTable(), GetAngle(), GetElectronAttachment(), GetElectronLongitudinalDiffusion(), GetElectronLorentzAngle(), GetElectronTownsend(), GetElectronTransverseDiffusion(), GetElectronVelocityB(), GetElectronVelocityE(), GetElectronVelocityExB(), GetFieldGrid(), GetHoleAttachment(), GetHoleLongitudinalDiffusion(), GetHoleTownsend(), GetHoleTransverseDiffusion(), GetHoleVelocityB(), GetHoleVelocityE(), GetHoleVelocityExB(), GetIonDissociation(), GetIonLongitudinalDiffusion(), GetIonMobility(), GetIonTransverseDiffusion(), HoleAttachment(), HoleDiffusion(), HoleTownsend(), HoleVelocity(), IonDiffusion(), IonDissociation(), IonVelocity(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), SetIonMobility(), and Garfield::MediumGas::WriteGasFile().

◆ m_bFields

std::vector<double> Garfield::Medium::m_bFields

protected

Definition at line 333 of file Medium.hh.

Referenced by CloneTable(), CloneTensor(), ElectronAttachment(), ElectronDiffusion(), ElectronLorentzAngle(), ElectronTownsend(), ElectronVelocity(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronAttachment(), GetElectronLongitudinalDiffusion(), GetElectronLorentzAngle(), GetElectronTownsend(), GetElectronTransverseDiffusion(), GetElectronVelocityB(), GetElectronVelocityE(), GetElectronVelocityExB(), GetFieldGrid(), GetHoleAttachment(), GetHoleLongitudinalDiffusion(), GetHoleTownsend(), GetHoleTransverseDiffusion(), GetHoleVelocityB(), GetHoleVelocityE(), GetHoleVelocityExB(), GetIonDissociation(), GetIonLongitudinalDiffusion(), GetIonMobility(), GetIonTransverseDiffusion(), HoleAttachment(), HoleDiffusion(), HoleTownsend(), HoleVelocity(), IonDiffusion(), IonDissociation(), IonVelocity(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), SetIonMobility(), and Garfield::MediumGas::WriteGasFile().

◆ m_className

std::string Garfield::Medium::m_className

protected

Definition at line 294 of file Medium.hh.

◆ m_debug

bool Garfield::Medium::m_debug

protected

Definition at line 329 of file Medium.hh.

Referenced by CloneTable(), DisableDebugging(), ElectronVelocity(), EnableDebugging(), Garfield::MediumMagboltz::GenerateGasTable(), GetDeexcitationProduct(), GetDielectricFunction(), GetElectronAttachment(), GetElectronCollision(), Garfield::MediumMagboltz::GetElectronCollision(), Garfield::MediumMagboltz::GetElectronCollisionRate(), GetElectronCollisionRate(), GetElectronLongitudinalDiffusion(), GetElectronLorentzAngle(), Garfield::MediumSilicon::GetElectronMomentum(), Garfield::MediumMagboltz::GetElectronNullCollisionRate(), GetElectronNullCollisionRate(), GetElectronTownsend(), GetElectronTransverseDiffusion(), GetElectronVelocityB(), GetElectronVelocityE(), GetElectronVelocityExB(), GetHoleAttachment(), GetHoleLongitudinalDiffusion(), GetHoleTownsend(), GetHoleTransverseDiffusion(), GetHoleVelocityB(), GetHoleVelocityE(), GetHoleVelocityExB(), GetIonDissociation(), GetIonisationProduct(), GetIonLongitudinalDiffusion(), GetIonMobility(), GetIonTransverseDiffusion(), Garfield::MediumMagboltz::GetLevel(), GetOpticalDataRange(), Garfield::MediumCdTe::GetOpticalDataRange(), Garfield::MediumGaAs::GetOpticalDataRange(), Garfield::MediumSilicon::GetOpticalDataRange(), GetPhotoAbsorptionCrossSection(), Garfield::MediumSilicon::Initialise(), Garfield::MediumMagboltz::Initialise(), Interpolate1D(), Garfield::MediumGas::LoadGasFile(), Garfield::MediumGas::LoadIonMobility(), Garfield::MediumMagboltz::RunMagboltz(), SetFieldGrid(), SetIonMobility(), and Garfield::MediumGas::WriteGasFile().

◆ m_density

double Garfield::Medium::m_density

protected

Definition at line 315 of file Medium.hh.

Referenced by GetMassDensity(), GetNumberDensity(), GetPhotonCollisionRate(), SetMassDensity(), and SetNumberDensity().

◆ m_driftable

bool Garfield::Medium::m_driftable

protected

Definition at line 318 of file Medium.hh.

Referenced by DisableDrift(), EnableDrift(), and IsDriftable().

◆ m_eFields

std::vector<double> Garfield::Medium::m_eFields

protected

Definition at line 332 of file Medium.hh.

Referenced by CloneTable(), CloneTensor(), ElectronAttachment(), ElectronDiffusion(), ElectronLorentzAngle(), ElectronTownsend(), ElectronVelocity(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronAttachment(), GetElectronLongitudinalDiffusion(), GetElectronLorentzAngle(), GetElectronTownsend(), GetElectronTransverseDiffusion(), GetElectronVelocityB(), GetElectronVelocityE(), GetElectronVelocityExB(), GetFieldGrid(), GetHoleAttachment(), GetHoleLongitudinalDiffusion(), GetHoleTownsend(), GetHoleTransverseDiffusion(), GetHoleVelocityB(), GetHoleVelocityE(), GetHoleVelocityExB(), GetIonDissociation(), GetIonLongitudinalDiffusion(), GetIonMobility(), GetIonTransverseDiffusion(), HoleAttachment(), HoleDiffusion(), HoleTownsend(), HoleVelocity(), IonDiffusion(), IonDissociation(), IonVelocity(), Garfield::MediumGas::LoadGasFile(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), SetIonMobility(), and Garfield::MediumGas::WriteGasFile().

◆ m_epsilon

double Garfield::Medium::m_epsilon

protected

Definition at line 307 of file Medium.hh.

Referenced by GetDielectricConstant(), and SetDielectricConstant().

◆ m_extrHighAttachment

unsigned int Garfield::Medium::m_extrHighAttachment

protected

Definition at line 390 of file Medium.hh.

Referenced by ElectronAttachment(), HoleAttachment(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodAttachment(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrHighDiffusion

unsigned int Garfield::Medium::m_extrHighDiffusion

protected

Definition at line 388 of file Medium.hh.

Referenced by ElectronDiffusion(), HoleDiffusion(), IonDiffusion(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodDiffusion(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrHighDissociation

unsigned int Garfield::Medium::m_extrHighDissociation

protected

Definition at line 393 of file Medium.hh.

Referenced by IonDissociation(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodIonDissociation(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrHighLorentzAngle

unsigned int Garfield::Medium::m_extrHighLorentzAngle

protected

Definition at line 391 of file Medium.hh.

Referenced by ElectronLorentzAngle(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrHighMobility

unsigned int Garfield::Medium::m_extrHighMobility

protected

Definition at line 392 of file Medium.hh.

Referenced by IonVelocity(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodIonMobility(), SetFieldGrid(), SetIonMobility(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrHighTownsend

unsigned int Garfield::Medium::m_extrHighTownsend

protected

Definition at line 389 of file Medium.hh.

Referenced by ElectronTownsend(), HoleTownsend(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodTownsend(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrHighVelocity

unsigned int Garfield::Medium::m_extrHighVelocity

protected

Definition at line 387 of file Medium.hh.

Referenced by ElectronVelocity(), HoleVelocity(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodVelocity(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrLowAttachment

unsigned int Garfield::Medium::m_extrLowAttachment

protected

Definition at line 390 of file Medium.hh.

Referenced by ElectronAttachment(), HoleAttachment(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodAttachment(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrLowDiffusion

unsigned int Garfield::Medium::m_extrLowDiffusion

protected

Definition at line 388 of file Medium.hh.

Referenced by ElectronDiffusion(), HoleDiffusion(), IonDiffusion(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodDiffusion(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrLowDissociation

unsigned int Garfield::Medium::m_extrLowDissociation

protected

Definition at line 393 of file Medium.hh.

Referenced by IonDissociation(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodIonDissociation(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrLowLorentzAngle

unsigned int Garfield::Medium::m_extrLowLorentzAngle

protected

Definition at line 391 of file Medium.hh.

Referenced by ElectronLorentzAngle(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrLowMobility

unsigned int Garfield::Medium::m_extrLowMobility

protected

Definition at line 392 of file Medium.hh.

Referenced by IonVelocity(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodIonMobility(), SetFieldGrid(), SetIonMobility(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrLowTownsend

unsigned int Garfield::Medium::m_extrLowTownsend

protected

Definition at line 389 of file Medium.hh.

Referenced by ElectronTownsend(), HoleTownsend(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodTownsend(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_extrLowVelocity

unsigned int Garfield::Medium::m_extrLowVelocity

protected

Definition at line 387 of file Medium.hh.

Referenced by ElectronVelocity(), HoleVelocity(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetExtrapolationMethodVelocity(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_fano

double Garfield::Medium::m_fano

protected

Definition at line 323 of file Medium.hh.

Referenced by GetFanoFactor(), Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumGaAs::MediumGaAs(), Garfield::MediumSilicon::MediumSilicon(), and SetFanoFactor().

◆ m_hasElectronAttachment

bool Garfield::Medium::m_hasElectronAttachment

protected

Definition at line 341 of file Medium.hh.

Referenced by ElectronAttachment(), Garfield::MediumCdTe::ElectronAttachment(), Garfield::MediumGaAs::ElectronAttachment(), Garfield::MediumSilicon::ElectronAttachment(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronAttachment(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetElectronAttachment(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasElectronDiffLong

bool Garfield::Medium::m_hasElectronDiffLong

protected

Definition at line 340 of file Medium.hh.

Referenced by ElectronDiffusion(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronLongitudinalDiffusion(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetElectronDiffusion(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasElectronDiffTens

bool Garfield::Medium::m_hasElectronDiffTens

protected

Definition at line 340 of file Medium.hh.

Referenced by ElectronDiffusion(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetElectronDiffusion(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasElectronDiffTrans

bool Garfield::Medium::m_hasElectronDiffTrans

protected

Definition at line 340 of file Medium.hh.

Referenced by ElectronDiffusion(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronTransverseDiffusion(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetElectronDiffusion(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasElectronLorentzAngle

bool Garfield::Medium::m_hasElectronLorentzAngle

protected

Definition at line 342 of file Medium.hh.

Referenced by ElectronLorentzAngle(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronLorentzAngle(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetElectronLorentzAngle(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasElectronVelocityB

bool Garfield::Medium::m_hasElectronVelocityB

protected

Definition at line 339 of file Medium.hh.

Referenced by ElectronVelocity(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronVelocityB(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetElectronVelocity(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasElectronVelocityE

bool Garfield::Medium::m_hasElectronVelocityE

protected

Definition at line 339 of file Medium.hh.

Referenced by ElectronVelocity(), Garfield::MediumCdTe::ElectronVelocity(), Garfield::MediumGaAs::ElectronVelocity(), Garfield::MediumSilicon::ElectronVelocity(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronVelocityE(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetElectronVelocity(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasElectronVelocityExB

bool Garfield::Medium::m_hasElectronVelocityExB

protected

Definition at line 339 of file Medium.hh.

Referenced by ElectronVelocity(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronVelocityExB(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetElectronVelocity(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasHoleAttachment

bool Garfield::Medium::m_hasHoleAttachment

protected

Definition at line 358 of file Medium.hh.

Referenced by GetHoleAttachment(), HoleAttachment(), Garfield::MediumCdTe::HoleAttachment(), Garfield::MediumGaAs::HoleAttachment(), Garfield::MediumSilicon::HoleAttachment(), Medium(), and ResetHoleAttachment().

◆ m_hasHoleDiffLong

bool Garfield::Medium::m_hasHoleDiffLong

protected

Definition at line 357 of file Medium.hh.

Referenced by GetHoleLongitudinalDiffusion(), HoleDiffusion(), Medium(), and ResetHoleDiffusion().

◆ m_hasHoleDiffTens

bool Garfield::Medium::m_hasHoleDiffTens

protected

Definition at line 357 of file Medium.hh.

Referenced by HoleDiffusion(), Medium(), ResetHoleDiffusion(), and SetFieldGrid().

◆ m_hasHoleDiffTrans

bool Garfield::Medium::m_hasHoleDiffTrans

protected

Definition at line 357 of file Medium.hh.

Referenced by GetHoleTransverseDiffusion(), HoleDiffusion(), Medium(), and ResetHoleDiffusion().

◆ m_hasHoleTownsend

bool Garfield::Medium::m_hasHoleTownsend

protected

Definition at line 358 of file Medium.hh.

Referenced by GetHoleTownsend(), HoleTownsend(), Garfield::MediumCdTe::HoleTownsend(), Garfield::MediumGaAs::HoleTownsend(), Garfield::MediumSilicon::HoleTownsend(), Medium(), and ResetHoleTownsend().

◆ m_hasHoleVelocityB

bool Garfield::Medium::m_hasHoleVelocityB

protected

Definition at line 356 of file Medium.hh.

Referenced by GetHoleVelocityB(), HoleVelocity(), Medium(), and ResetHoleVelocity().

◆ m_hasHoleVelocityE

bool Garfield::Medium::m_hasHoleVelocityE

protected

Definition at line 356 of file Medium.hh.

Referenced by GetHoleVelocityE(), HoleVelocity(), Garfield::MediumCdTe::HoleVelocity(), Garfield::MediumGaAs::HoleVelocity(), Garfield::MediumSilicon::HoleVelocity(), Medium(), and ResetHoleVelocity().

◆ m_hasHoleVelocityExB

bool Garfield::Medium::m_hasHoleVelocityExB

protected

Definition at line 356 of file Medium.hh.

Referenced by GetHoleVelocityExB(), HoleVelocity(), Medium(), and ResetHoleVelocity().

◆ m_hasIonDiffLong

bool Garfield::Medium::m_hasIonDiffLong

protected

Definition at line 371 of file Medium.hh.

Referenced by Garfield::MediumMagboltz::GenerateGasTable(), GetIonLongitudinalDiffusion(), IonDiffusion(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetIonDiffusion(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasIonDiffTrans

bool Garfield::Medium::m_hasIonDiffTrans

protected

Definition at line 371 of file Medium.hh.

Referenced by Garfield::MediumMagboltz::GenerateGasTable(), GetIonTransverseDiffusion(), IonDiffusion(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetIonDiffusion(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasIonDissociation

bool Garfield::Medium::m_hasIonDissociation

protected

Definition at line 372 of file Medium.hh.

Referenced by Garfield::MediumMagboltz::GenerateGasTable(), GetIonDissociation(), IonDissociation(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetIonDissociation(), and Garfield::MediumGas::WriteGasFile().

◆ m_hasIonMobility

bool Garfield::Medium::m_hasIonMobility

protected

Definition at line 370 of file Medium.hh.

Referenced by Garfield::MediumMagboltz::GenerateGasTable(), GetIonMobility(), IonVelocity(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), ResetIonMobility(), SetIonMobility(), and Garfield::MediumGas::WriteGasFile().

◆ m_id

int Garfield::Medium::m_id

protected

Definition at line 299 of file Medium.hh.

Referenced by GetId().

◆ m_idCounter

int Garfield::Medium::m_idCounter = -1

staticprotected

Definition at line 296 of file Medium.hh.

◆ m_intpAttachment

unsigned int Garfield::Medium::m_intpAttachment

protected

Definition at line 399 of file Medium.hh.

Referenced by ElectronAttachment(), HoleAttachment(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), SetInterpolationMethodAttachment(), and Garfield::MediumGas::WriteGasFile().

◆ m_intpDiffusion

unsigned int Garfield::Medium::m_intpDiffusion

protected

Definition at line 397 of file Medium.hh.

Referenced by ElectronDiffusion(), HoleDiffusion(), IonDiffusion(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), SetInterpolationMethodDiffusion(), and Garfield::MediumGas::WriteGasFile().

◆ m_intpDissociation

unsigned int Garfield::Medium::m_intpDissociation

protected

Definition at line 402 of file Medium.hh.

Referenced by IonDissociation(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), SetInterpolationMethodIonDissociation(), and Garfield::MediumGas::WriteGasFile().

◆ m_intpLorentzAngle

unsigned int Garfield::Medium::m_intpLorentzAngle

protected

Definition at line 400 of file Medium.hh.

Referenced by ElectronLorentzAngle(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ m_intpMobility

unsigned int Garfield::Medium::m_intpMobility

protected

Definition at line 401 of file Medium.hh.

Referenced by IonVelocity(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), SetInterpolationMethodIonMobility(), SetIonMobility(), and Garfield::MediumGas::WriteGasFile().

◆ m_intpTownsend

unsigned int Garfield::Medium::m_intpTownsend

protected

Definition at line 398 of file Medium.hh.

Referenced by ElectronTownsend(), HoleTownsend(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), SetInterpolationMethodTownsend(), and Garfield::MediumGas::WriteGasFile().

◆ m_intpVelocity

unsigned int Garfield::Medium::m_intpVelocity

protected

Definition at line 396 of file Medium.hh.

Referenced by ElectronVelocity(), HoleVelocity(), Garfield::MediumGas::LoadGasFile(), Medium(), Garfield::MediumGas::PrintGas(), SetFieldGrid(), SetInterpolationMethodVelocity(), and Garfield::MediumGas::WriteGasFile().

◆ m_ionisable

bool Garfield::Medium::m_ionisable

protected

Definition at line 320 of file Medium.hh.

Referenced by DisablePrimaryIonisation(), EnablePrimaryIonisation(), and IsIonisable().

◆ m_isChanged

bool Garfield::Medium::m_isChanged

protected

Definition at line 326 of file Medium.hh.

◆ m_map2d

bool Garfield::Medium::m_map2d

protected

Definition at line 337 of file Medium.hh.

Referenced by CloneTable(), CloneTensor(), ElectronAttachment(), ElectronDiffusion(), ElectronLorentzAngle(), ElectronTownsend(), ElectronVelocity(), HoleAttachment(), HoleDiffusion(), HoleTownsend(), HoleVelocity(), IonDiffusion(), IonDissociation(), IonVelocity(), Garfield::MediumGas::LoadGasFile(), SetFieldGrid(), SetIonMobility(), and Garfield::MediumGas::WriteGasFile().

◆ m_microscopic

bool Garfield::Medium::m_microscopic

protected

Definition at line 319 of file Medium.hh.

Referenced by IsMicroscopic(), Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumGaAs::MediumGaAs(), Garfield::MediumMagboltz::MediumMagboltz(), and Garfield::MediumSilicon::MediumSilicon().

◆ m_name

std::string Garfield::Medium::m_name

protected

Definition at line 301 of file Medium.hh.

Referenced by GetComponent(), GetName(), Garfield::MediumGas::LoadGasFile(), Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumConductor::MediumConductor(), Garfield::MediumGaAs::MediumGaAs(), Garfield::MediumGas::MediumGas(), Garfield::MediumPlastic::MediumPlastic(), Garfield::MediumSilicon::MediumSilicon(), Garfield::MediumGas::PrintGas(), Garfield::MediumGas::SetComposition(), and Garfield::MediumGas::WriteGasFile().

◆ m_nComponents

unsigned int Garfield::Medium::m_nComponents

protected

◆ m_pressure

double Garfield::Medium::m_pressure

protected

Definition at line 305 of file Medium.hh.

◆ m_temperature

double Garfield::Medium::m_temperature

protected

Definition at line 303 of file Medium.hh.

Referenced by ElectronDiffusion(), Garfield::MediumMagboltz::GenerateGasTable(), Garfield::MediumGas::GetNumberDensity(), GetTemperature(), HoleDiffusion(), IonDiffusion(), Garfield::MediumGas::LoadGasFile(), Garfield::MediumGas::LoadIonMobility(), Garfield::MediumMagboltz::MediumMagboltz(), Garfield::MediumGas::PrintGas(), Garfield::MediumMagboltz::RunMagboltz(), and SetTemperature().

◆ m_w

double Garfield::Medium::m_w

protected

Definition at line 323 of file Medium.hh.

Referenced by GetW(), Garfield::MediumCdTe::MediumCdTe(), Garfield::MediumGaAs::MediumGaAs(), Garfield::MediumSilicon::MediumSilicon(), and SetW().

◆ m_z

double Garfield::Medium::m_z

protected

Definition at line 311 of file Medium.hh.

Referenced by GetAtomicNumber(), and SetAtomicNumber().

◆ tabElectronAttachment

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabElectronAttachment

protected

Definition at line 349 of file Medium.hh.

Referenced by ElectronAttachment(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronAttachment(), Garfield::MediumGas::LoadGasFile(), ResetElectronAttachment(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabElectronDiffLong

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabElectronDiffLong

protected

Definition at line 346 of file Medium.hh.

Referenced by ElectronDiffusion(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronLongitudinalDiffusion(), Garfield::MediumGas::LoadGasFile(), ResetElectronDiffusion(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabElectronDiffTens

std::vector<std::vector<std::vector<std::vector<double> > > > Garfield::Medium::tabElectronDiffTens

protected

Definition at line 353 of file Medium.hh.

Referenced by ElectronDiffusion(), Garfield::MediumGas::LoadGasFile(), ResetElectronDiffusion(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabElectronDiffTrans

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabElectronDiffTrans

protected

Definition at line 347 of file Medium.hh.

Referenced by ElectronDiffusion(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronTransverseDiffusion(), Garfield::MediumGas::LoadGasFile(), ResetElectronDiffusion(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabElectronLorentzAngle

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabElectronLorentzAngle

protected

Definition at line 350 of file Medium.hh.

Referenced by ElectronLorentzAngle(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronLorentzAngle(), Garfield::MediumGas::LoadGasFile(), ResetElectronLorentzAngle(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabElectronTownsend

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabElectronTownsend

protected

Definition at line 348 of file Medium.hh.

Referenced by ElectronTownsend(), Garfield::MediumCdTe::ElectronTownsend(), Garfield::MediumGaAs::ElectronTownsend(), Garfield::MediumSilicon::ElectronTownsend(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronTownsend(), Garfield::MediumGas::LoadGasFile(), Garfield::MediumGas::PrintGas(), ResetElectronTownsend(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabElectronVelocityB

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabElectronVelocityB

protected

Definition at line 345 of file Medium.hh.

Referenced by ElectronVelocity(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronVelocityB(), Garfield::MediumGas::LoadGasFile(), ResetElectronVelocity(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabElectronVelocityE

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabElectronVelocityE

protected

Definition at line 343 of file Medium.hh.

Referenced by ElectronVelocity(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronVelocityE(), Garfield::MediumGas::LoadGasFile(), ResetElectronVelocity(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabElectronVelocityExB

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabElectronVelocityExB

protected

Definition at line 344 of file Medium.hh.

Referenced by ElectronVelocity(), Garfield::MediumMagboltz::GenerateGasTable(), GetElectronVelocityExB(), Garfield::MediumGas::LoadGasFile(), ResetElectronVelocity(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabHoleAttachment

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabHoleAttachment

protected

Definition at line 365 of file Medium.hh.

Referenced by GetHoleAttachment(), HoleAttachment(), ResetHoleAttachment(), and SetFieldGrid().

◆ tabHoleDiffLong

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabHoleDiffLong

protected

Definition at line 362 of file Medium.hh.

Referenced by GetHoleLongitudinalDiffusion(), HoleDiffusion(), ResetHoleDiffusion(), and SetFieldGrid().

◆ tabHoleDiffTens

std::vector<std::vector<std::vector<std::vector<double> > > > Garfield::Medium::tabHoleDiffTens

protected

Definition at line 367 of file Medium.hh.

Referenced by HoleDiffusion(), ResetHoleDiffusion(), and SetFieldGrid().

◆ tabHoleDiffTrans

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabHoleDiffTrans

protected

Definition at line 363 of file Medium.hh.

Referenced by GetHoleTransverseDiffusion(), HoleDiffusion(), ResetHoleDiffusion(), and SetFieldGrid().

◆ tabHoleTownsend

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabHoleTownsend

protected

Definition at line 364 of file Medium.hh.

Referenced by GetHoleTownsend(), HoleTownsend(), IonDissociation(), ResetHoleTownsend(), and SetFieldGrid().

◆ tabHoleVelocityB

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabHoleVelocityB

protected

Definition at line 361 of file Medium.hh.

Referenced by GetHoleVelocityB(), HoleVelocity(), ResetHoleVelocity(), and SetFieldGrid().

◆ tabHoleVelocityE

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabHoleVelocityE

protected

Definition at line 359 of file Medium.hh.

Referenced by GetHoleVelocityE(), HoleVelocity(), ResetHoleVelocity(), and SetFieldGrid().

◆ tabHoleVelocityExB

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabHoleVelocityExB

protected

Definition at line 360 of file Medium.hh.

Referenced by GetHoleVelocityExB(), HoleVelocity(), ResetHoleVelocity(), and SetFieldGrid().

◆ tabIonDiffLong

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabIonDiffLong

protected

Definition at line 374 of file Medium.hh.

Referenced by GetIonLongitudinalDiffusion(), IonDiffusion(), Garfield::MediumGas::LoadGasFile(), ResetIonDiffusion(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabIonDiffTrans

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabIonDiffTrans

protected

Definition at line 375 of file Medium.hh.

Referenced by GetIonTransverseDiffusion(), IonDiffusion(), Garfield::MediumGas::LoadGasFile(), ResetIonDiffusion(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabIonDissociation

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabIonDissociation

protected

Definition at line 376 of file Medium.hh.

Referenced by GetIonDissociation(), IonDissociation(), Garfield::MediumGas::LoadGasFile(), ResetIonDissociation(), SetFieldGrid(), and Garfield::MediumGas::WriteGasFile().

◆ tabIonMobility

std::vector<std::vector<std::vector<double> > > Garfield::Medium::tabIonMobility

protected

Definition at line 373 of file Medium.hh.

Referenced by GetIonMobility(), IonVelocity(), Garfield::MediumGas::LoadGasFile(), ResetIonMobility(), SetFieldGrid(), SetIonMobility(), and Garfield::MediumGas::WriteGasFile().

◆ thrElectronAttachment

int Garfield::Medium::thrElectronAttachment

protected

Definition at line 380 of file Medium.hh.

Referenced by ElectronAttachment(), Garfield::MediumGas::LoadGasFile(), Medium(), and Garfield::MediumGas::WriteGasFile().

◆ thrElectronTownsend

int Garfield::Medium::thrElectronTownsend

protected

Definition at line 379 of file Medium.hh.

Referenced by ElectronTownsend(), Garfield::MediumGas::LoadGasFile(), Medium(), and Garfield::MediumGas::WriteGasFile().

◆ thrHoleAttachment

int Garfield::Medium::thrHoleAttachment

protected

Definition at line 383 of file Medium.hh.

Referenced by HoleAttachment(), and Medium().

◆ thrHoleTownsend

int Garfield::Medium::thrHoleTownsend

protected

Definition at line 382 of file Medium.hh.

Referenced by HoleTownsend(), and Medium().

◆ thrIonDissociation

int Garfield::Medium::thrIonDissociation

protected

Definition at line 384 of file Medium.hh.

Referenced by IonDissociation(), Garfield::MediumGas::LoadGasFile(), Medium(), and Garfield::MediumGas::WriteGasFile().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Medium()

◆ ~Medium()

Member Function Documentation

◆ CloneTable()

◆ CloneTensor()

◆ DisableDebugging()

◆ DisableDrift()

◆ DisablePrimaryIonisation()

◆ ElectronAttachment()

◆ ElectronDiffusion() [1/2]

◆ ElectronDiffusion() [2/2]

◆ ElectronLorentzAngle()

◆ ElectronTownsend()

◆ ElectronVelocity()

◆ EnableDebugging()

◆ EnableDrift()

◆ EnablePrimaryIonisation()

◆ GetAngle()

◆ GetAtomicNumber()

◆ GetAtomicWeight()

◆ GetComponent()

◆ GetDeexcitationProduct()

◆ GetDielectricConstant()

◆ GetDielectricFunction()

◆ GetElectronAttachment()

◆ GetElectronCollision()

◆ GetElectronCollisionRate()

◆ GetElectronEnergy()

◆ GetElectronLongitudinalDiffusion()

◆ GetElectronLorentzAngle()

◆ GetElectronMomentum()

◆ GetElectronNullCollisionRate()

◆ GetElectronTownsend()

◆ GetElectronTransverseDiffusion()

◆ GetElectronVelocityB()

◆ GetElectronVelocityE()

◆ GetElectronVelocityExB()

◆ GetExtrapolationIndex()

◆ GetFanoFactor()

◆ GetFieldGrid()

◆ GetHoleAttachment()

◆ GetHoleLongitudinalDiffusion()

◆ GetHoleTownsend()

◆ GetHoleTransverseDiffusion()

◆ GetHoleVelocityB()

◆ GetHoleVelocityE()

◆ GetHoleVelocityExB()

◆ GetId()

◆ GetIonDissociation()

◆ GetIonisationProduct()

◆ GetIonLongitudinalDiffusion()

◆ GetIonMobility()

◆ GetIonTransverseDiffusion()

◆ GetMassDensity()

◆ GetName()

◆ GetNumberDensity()

◆ GetNumberOfComponents()

◆ GetNumberOfDeexcitationProducts()

◆ GetNumberOfIonisationProducts()

◆ GetOpticalDataRange()

◆ GetPhotoAbsorptionCrossSection()

◆ GetPhotonCollision()

◆ GetPhotonCollisionRate()

◆ GetPressure()

◆ GetTemperature()

◆ GetW()

◆ HoleAttachment()

◆ HoleDiffusion() [1/2]

◆ HoleDiffusion() [2/2]

◆ HoleTownsend()

◆ HoleVelocity()

◆ InitParamArrays()

◆ InitParamTensor()

◆ Interpolate1D()