TIonizationLoss Class Reference

Ionization energy losses. More...

#include <galaxy.h>

Inheritance diagram for TIonizationLoss:

List of all members.

Public Member Functions
	TIonizationLoss (TGrid *, vector< TGas * >, TGas *, double, double)
	Initialize ionization losses.
virtual	~TIonizationLoss ()

---

Detailed Description

Ionization energy losses.

Author:

Luca Maccione luca.maccione@desy.de

Definition at line 215 of file galaxy.h.

---

Constructor & Destructor Documentation

TIonizationLoss::TIonizationLoss	(	TGrid *	coord,
		vector< TGas * >	gas,
		TGas *	totalgas,
		double	A,
		double	Z
	)

Initialize ionization losses.

Parameters:

	coord	geometry
	gas	ISM gas components
	totalgas	Total ISM gas
	A	Nucleus mass
	Z	Nucleus charge

Definition at line 24 of file galaxy.cc.

References TEnergyLoss::dimE, TEnergyLoss::dimr, TEnergyLoss::dimz, TEnergyLoss::dpdt, EH, EHe, EiH, EiHe, TGrid::GetBeta(), TGrid::GetBetaEl(), TGrid::GetGamma(), TGrid::GetGammaEl(), TGas::GetGas(), He_abundance, TEnergyLoss::index(), TGrid::index(), Mele, mp, Myr, PIR02MC2C, and ZHe.

                                                                                                    : TEnergyLoss(coord) {

  const double factor = 1e-9*Myr;

  if (A != 0) { // hadrons

    vector<double> bet = coord->GetBeta();
    vector<double> gamma = coord->GetGamma();

    const double MA = A*mp*1.e3;  // nucleus mass MeV
    
    for (int k = 0; k < dimE; ++k) {
      double gammak = gamma[k];
      double betak = bet[k];
      double qmax = 2.*Mele*(gammak*gammak-1.)/(1. +2.*gammak*Mele/MA);
      double bh   = log(2.*Mele*(gammak*gammak-1.)*qmax /(EH*EH))  -2.*betak*betak;
      double bhe  = log(2.*Mele*(gammak*gammak-1.)*qmax /(EHe*EHe))-2.*betak*betak;  

      for (int i = 0; i < dimr; ++i) {
        for (int j = 0; j < dimz; ++j) {
          
          int ind = coord->index(i,j);
          double nh = totalgas->GetGas(ind) - gas[2]->GetGas(ind);
          dpdt[index(i,j,k)] = 2.*PIR02MC2C*fabs(Z)*fabs(Z)/betak*nh*(bh +He_abundance*bhe) * factor / betak;
          
        }
      }
    }
  }
  else { // leptons

    vector<double> bet = coord->GetBetaEl();
    vector<double> gamma = coord->GetGammaEl();
    
    for (int k = 0; k < dimE; ++k) {
      double gammak = gamma[k];
      double betk = bet[k];

      for (int i = 0; i < dimr; ++i) {
        for (int j = 0; j < dimz; ++j) {
          
          int ind = coord->index(i,j);
          double nh = totalgas->GetGas(ind) - gas[2]->GetGas(ind);
          dpdt[index(i,j,k)] = 2.*PIR02MC2C/betk*(  nh*(1.0+ZHe*He_abundance)*( log(gammak-1.)- M_LN2 +1./8.
                                                                                +2.*log(gammak*betk*Mele) )-nh*2.*log(EiH)-ZHe*nh*He_abundance*2.*log(EiHe)  ) * factor / betk;
        }
      }
    }
  }
}

virtual TIonizationLoss::~TIonizationLoss

(

)

[inline, virtual]

Destructor.

Definition at line 223 of file galaxy.h.

---

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

• galaxy.h
• galaxy.cc