constants.h File Reference

Definition of many constants and settings that are used in the code. More...

Go to the source code of this file.

Enumerations
enum	DMreaction { Decay, Annihilation }
	Chose whether to propagate products of DM decay or annihilation. More...
enum	DMspectrum { Delta, DarkSUSY }
	Chose whether to use DarkSUSY spectra or a Delta function. More...
enum	DMprofile {   ISO = 0, NFW = 1, Kra = 2, Moore = 3,   Einasto = 4 }
	Chose the spatial DM profile: Isothermal, NFW, Kra, Moore, Einasto. More...
enum	GasType { BronfFerr, NS, Galprop }
	Chose the gas distribution: Bronfmann + Ferriere, Nakanishi & Sophue or Galprop. More...
enum	SNRType { Galprop_, Ferriere, PointSource, Ring }
	Chose the SNR (CR source) distribution: Galprop, Ferriere, a Point Source at Galactic Center, a Ring (. More...
enum	DPerpType {   Constant, Exp, Expr, Radial,   ExpRadial }
	Spatial properties of the diffusion coefficient. More...
enum	SpallXSec { GalpropXSec, Webber03 }
	Choose the spallation cross section model. More...
enum	LeptonYields { GalpropTable, Kamae }
	Choose model for lepton yields in secondary electron and positron production. More...
Variables
const DMreaction	DMr = Annihilation
const DMspectrum	DMs = DarkSUSY
const double	rhos = 0.36
const double	EkDelta = 100.0
const double	small = 1e-45
const SpallXSec	spallationxsec = GalpropXSec
const LeptonYields	ly = GalpropTable
const bool	fullstore = true
const bool	partialstore = true
const std::string	ISRFfile = "data/MilkyWay_DR0.5_DZ0.1_DPHI10_RMAX20_ZMAX5_galprop_format.fits.gz"
const std::string	ElTablefile = "data/Electron_production.dat"
const std::string	PosTablefile = "data/Positron_production.dat"
const std::string	Webber03Data = "data/webber_xsec.dat"
const std::string	Webber03DataTotal = "data/webber_xsec_total.dat"
const std::string	BNLdata = "data/bnl_data.dat"
const double	C = 2.99792458e10
const double	h_planck = 6.6260755e-27
const double	year = 365.25*24.*3600.
const double	Myr = 1e6*year
const double	kpc = 3.08568e21
const double	km = 1e5
const double	Clight = C * Myr
const double	mp = .938
const double	Pi = 3.14159265358979312
const double	eV_to_erg = 1.60217733e-12
const double	erg_to_eV = 1./eV_to_erg
const double	u = 1./200.
const double	p = 3.0
const int	cross_section_option = 12
const double	t_half_limit = 1.e4
const int	Nrept = 20
const double	dtfactor = 0.25
const double	dtmin = 1e-6
const double	dtmax = 1e3
const bool	ADI = false
const bool	OpSplit = true
const double	tol = 3e-4
const double	alpha = 0.09
const double	rB = 10.
const double	zr = 2.0
const double	Bh = 5.
const double	D_ref_rig = 3.
const double	rd = 2.
const double	r0 = 2.
const double	dzzGal = 0.01
const int	diff_reacc = 1
const double	robs = 8.5
const double	zobs = 0.
const double	Rmin = 0.
const double	Rmax = 20.
const GasType	gas_model = Galprop
const SNRType	SNR_model = Ferriere
const double	ringmin = 3.
const double	ringmax = 4.
const double	He_abundance = 0.11
const double	ALPHAf = 1.0/137.035989561
const double	Mele = 0.5109990615
const double	MeleGeV = 1e-3*Mele
const double	Rele = 2.8179409238e-13
const double	H2PiC = 0.19732705359e-10
const double	PIR0H2C2 = Pi*Rele*H2PiC*H2PiC
const double	PIR02MC2C = Pi*Rele*Rele*C*(Mele*1.e6)
const double	PIR02C = Pi*Rele*Rele*C
const double	AFR02MC2C = ALPHAf*Rele*Rele*C*(Mele*1.e6)
const double	EiH = 13.6e-6
const double	EiHe = 24.59e-6
const double	TH = 62.8
const double	THe = 93.1
const double	ZHe = 2.
const double	EH = 19.e-6
const double	EHe = 44.e-6
const double	BK = 1.38066e-23/1.60218e-13
const double	Te = 1.0e4
const double	bet_e = sqrt(2.*BK*Te/Mele)
const double	xm = pow(pow(3.*sqrt(Pi)/4.,1./3.) *bet_e,3)
const double	MH = 1.67e-24
const double	MHe = 1.66e-24 *4.
const double	gam1 = 100.
const double	gam2 = 800.
const double	sigmaT = 6.6524e-25/(kpc*kpc)
const double	sp_ref_rig_norm = 1.e2
const double	sp_ref_rig_el = 33.20
const double	spect_norm = 5.e-2
const double	spect_norm_el = 4.2e-3
const double	Ekmax = 3.e3
const double	Ekmin = 0.01
const double	sp_ref_rig = 9.
const double	cutoff_rig = 5.e3
const double	spect_ind_el_low = 2.25
const double	spect_ind_el = 2.25
const double	sp_ref_rig_break_el = 4.
const double	ind_low_nucl = 2.75
const int	antiproton_cs = 2
const bool	scaling = false
const int	Zmax = 14
const int	Zmin = 1
const bool	prop_ap = false
const bool	prop_el = false
const bool	prop_extracomp = false
	Whether to propagate the extra component (electron/positron).
const double	ab_pr = 1.06e6
const double	ab_He3 = 9.033
const double	ab_He4 = 7.199e4
const double	ab_Li = 0.
const double	ab_Be9 = 0.
const double	ab_Be10 = 0.
const double	ab_B10 = 0.
const double	ab_B11 = 0.
const double	ab_C13 = 5.268e-07
const double	ab_C14 = 0.0
const double	ab_N15 = 5.961e-5
const double	ab_O = 3822.
const double	ab_O17 = 6.713e-7
const double	ab_O18 = 1.286
const double	ab_F19 = 2.664e-8
const double	ab_Ne20 = 312.5
const double	ab_Ne21 = 0.003556
const double	ab_Ne22 = 100.1
const double	ab_Na22 = 0.0
const double	ab_Na23 = 22.84
const double	ab_Mg24 = 658.1
const double	ab_Mg25 = 82.5
const double	ab_Mg26 = 104.7
const double	ab_Al26 = 0.0
const double	ab_Al27 = 76.42
const double	ab_Si28 = 725.7
const double	ab_Si29 = 35.02
const double	ab_Si30 = 24.68
const double	ab_Si32 = 24.68
const double	ab_P31 = 7.28
const double	ab_S32 = 94.40
const double	ab_S33 = 0.8
const double	ab_S34 = 4.45
const double	ab_S36 = 0.02
const double	ab_Cl35 = 2.11
const double	ab_Cl37 = 0.72
const double	ab_Ar36 = 11.17
const double	ab_Ar38 = 2.22
const double	ab_K39 = 4.47
const double	ab_Ca40 = 39.87
const double	ab_Ca42 = 0.28
const double	ab_Ca44 = 0.93
const double	ab_Ca48 = 0.09
const double	ab_Sc45 = 0.121
const double	ab_Ti46 = 0.17
const double	ab_Ti47 = 0.16
const double	ab_Ti48 = 1.64
const double	ab_Ti49 = 0.124
const double	ab_Ti50 = 0.12
const double	ab_V51 = 0.0
const double	ab_Cr50 = 0.72
const double	ab_Cr52 = 14.51
const double	ab_Cr53 = 1.69
const double	ab_Cr54 = 0.43
const double	ab_Mn55 = 16.15
const double	ab_Fe54 = 38.00
const double	ab_Fe56 = 620.2
const double	ab_Fe57 = 15.07
const double	ab_Fe58 = 2.31
const double	ab_Co59 = 1.24
const double	ab_Ni58 = 26.21
const double	ab_Ni60 = 10.44
const double	ab_Ni61 = 0.48
const double	ab_Ni62 = 1.51
const double	ab_Ni64 = 0.46
const double	spect_ind_pr = 2.75
const double	spect_ind_He = 2.75
const double	spect_ind_Be = 2.75
const double	spect_ind_B = 2.75
const double	spect_ind_C = 2.75
const double	spect_ind_N = 2.75
const double	spect_ind_O = 2.75
const double	spect_ind_F = 2.75
const double	spect_ind_Ne = 2.75
const double	spect_ind_Na = 2.75
const double	spect_ind_Mg = 2.75
const double	spect_ind_Al = 2.75
const double	spect_ind_Si = 2.75
const double	spect_ind_P = 2.75
const double	spect_ind_S = 2.75
const double	spect_ind_Cl = 2.75
const double	spect_ind_Ar = 2.75
const double	spect_ind_K = 2.75
const double	spect_ind_Ca = 2.75
const double	spect_ind_Sc = 2.75
const double	spect_ind_Ti = 2.75
const double	spect_ind_V = 2.75
const double	spect_ind_Cr = 2.75
const double	spect_ind_Mn = 2.75
const double	spect_ind_Fe = 2.75
const double	spect_ind_Co = 2.75
const double	spect_ind_Ni = 2.75

---

Detailed Description

Definition of many constants and settings that are used in the code.

Author:

Luca Maccione luca.maccione@desy.de

Definition in file constants.h.

---

Enumeration Type Documentation

enum DMprofile

Chose the spatial DM profile: Isothermal, NFW, Kra, Moore, Einasto.

See also:

include/input.h . For the DM density at solar position

rhos.

Enumerator:

ISO	Set Isothermal profile.
NFW	Set NFW profile.
Kra	Set Kra profile.
Moore	Set Moore profile.
Einasto	Einasto profile.

Definition at line 54 of file constants.h.

             {
  ISO=0, 
  NFW=1, 
  Kra=2, 
  Moore=3, 
  Einasto=4 
} DMprofile;

enum DMreaction

Chose whether to propagate products of DM decay or annihilation.

Enumerator:

Decay	Set DM decay.
Annihilation	SetDM annihilation.

Definition at line 36 of file constants.h.

             {
  Decay, 
  Annihilation 
} DMreaction;

enum DMspectrum

Chose whether to use DarkSUSY spectra or a Delta function.

Enumerator:

Delta	Set Delta function at kinetic energy See also: EkDelta.
DarkSUSY	Set DarkSUSY spectrum.

Definition at line 45 of file constants.h.

             {
  Delta, 
  DarkSUSY 
} DMspectrum;

enum DPerpType

Spatial properties of the diffusion coefficient.

Enumerator:

Constant	constant diffusion coefficient.
Exp	Exponentially growing with vertical scale. See also: zt .
Expr	Expr .
Radial	radial behavior according to SNR distribution.
ExpRadial	Radial behavior according to SNR distribution and exponential growth with vertical scale.

Definition at line 87 of file constants.h.

             {
  Constant, 
  Exp, 
  Expr, 
  Radial, 
  ExpRadial 
} DPerpType;

enum GasType

Chose the gas distribution: Bronfmann + Ferriere, Nakanishi & Sophue or Galprop.

Enumerator:

BronfFerr	Set Bronfmann + Ferriere distribution.
NS	Set Nakanishi & Sophue distribution.
Galprop	Set Galprop distribution.

Definition at line 66 of file constants.h.

             {
  BronfFerr, 
  NS, 
  Galprop 
} GasType; 

enum LeptonYields

Choose model for lepton yields in secondary electron and positron production.

Enumerator:

GalpropTable	Tabulated Galprop model.
Kamae	Kamae 2007 yields. See also: include/cparam.h .

Definition at line 108 of file constants.h.

             {
  GalpropTable, 
  Kamae 
} LeptonYields;

enum SNRType

Chose the SNR (CR source) distribution: Galprop, Ferriere, a Point Source at Galactic Center, a Ring (.

See also:

ringmin, ringmax).

Enumerator:

Galprop_	Set Galprop parameterized distribution.
Ferriere	Set Ferriere distribution.
PointSource	Set Point Source at GC.
Ring	Set Ring of sources: See also: ringmin ringmax .

Definition at line 76 of file constants.h.

             {
  Galprop_, 
  Ferriere, 
  PointSource, 
  Ring 
} SNRType; 

enum SpallXSec

Choose the spallation cross section model.

Enumerator:

GalpropXSec	Use Galprop cross sections.
Webber03	Use Webber 2003 cross sections where available.

Definition at line 99 of file constants.h.

             {
  GalpropXSec, 
  Webber03 
} SpallXSec;

---

Variable Documentation

const double ab_Al26 = 0.0

Definition at line 277 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Al27 = 76.42

Definition at line 278 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ar36 = 11.17

Definition at line 290 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ar38 = 2.22

Definition at line 291 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_B10 = 0.

Definition at line 258 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_B11 = 0.

Definition at line 259 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Be10 = 0.

Definition at line 257 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Be9 = 0.

Definition at line 256 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_C13 = 5.268e-07

Definition at line 261 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_C14 = 0.0

Definition at line 262 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ca40 = 39.87

Definition at line 293 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ca42 = 0.28

Definition at line 294 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ca44 = 0.93

Definition at line 295 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ca48 = 0.09

Definition at line 296 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Cl35 = 2.11

Definition at line 288 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Cl37 = 0.72

Definition at line 289 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Co59 = 1.24

Definition at line 313 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Cr50 = 0.72

Definition at line 304 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Cr52 = 14.51

Definition at line 305 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Cr53 = 1.69

Definition at line 306 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Cr54 = 0.43

Definition at line 307 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_F19 = 2.664e-8

Definition at line 268 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Fe54 = 38.00

Definition at line 309 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Fe56 = 620.2

Definition at line 310 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Fe57 = 15.07

Definition at line 311 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Fe58 = 2.31

Definition at line 312 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_He3 = 9.033

Definition at line 253 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_He4 = 7.199e4

Definition at line 254 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_K39 = 4.47

Definition at line 292 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Li = 0.

Definition at line 255 of file constants.h.

const double ab_Mg24 = 658.1

Definition at line 274 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Mg25 = 82.5

Definition at line 275 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Mg26 = 104.7

Definition at line 276 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Mn55 = 16.15

Definition at line 308 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_N15 = 5.961e-5

Definition at line 264 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Na22 = 0.0

Definition at line 272 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Na23 = 22.84

Definition at line 273 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ne20 = 312.5

Definition at line 269 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ne21 = 0.003556

Definition at line 270 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ne22 = 100.1

Definition at line 271 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ni58 = 26.21

Definition at line 314 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ni60 = 10.44

Definition at line 315 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ni61 = 0.48

Definition at line 316 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ni62 = 1.51

Definition at line 317 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ni64 = 0.46

Definition at line 318 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_O = 3822.

Definition at line 265 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_O17 = 6.713e-7

Definition at line 266 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_O18 = 1.286

Definition at line 267 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_P31 = 7.28

Definition at line 283 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_pr = 1.06e6

< Source abundances.

Definition at line 252 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_S32 = 94.40

Definition at line 284 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_S33 = 0.8

Definition at line 285 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_S34 = 4.45

Definition at line 286 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_S36 = 0.02

Definition at line 287 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Sc45 = 0.121

Definition at line 297 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Si28 = 725.7

Definition at line 279 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Si29 = 35.02

Definition at line 280 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Si30 = 24.68

Definition at line 281 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Si32 = 24.68

Definition at line 282 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ti46 = 0.17

Definition at line 298 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ti47 = 0.16

Definition at line 299 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ti48 = 1.64

Definition at line 300 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ti49 = 0.124

Definition at line 301 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_Ti50 = 0.12

Definition at line 302 of file constants.h.

Referenced by Galaxy::Galaxy().

const double ab_V51 = 0.0

Definition at line 303 of file constants.h.

Referenced by Galaxy::Galaxy().

const bool ADI = false

Controls whether to refine the solution with the ADI scheme.

Definition at line 161 of file constants.h.

Referenced by DRAGON::DRAGON(), DRAGON::particle_modulated_spectrum(), and DRAGON::PrintRatios().

const double AFR02MC2C = ALPHAf*Rele*Rele*C*(Mele*1.e6)

[eV*cm^3/s] = e^6/mc/hc

Definition at line 196 of file constants.h.

Referenced by TBremsstrahlungLoss::TBremsstrahlungLoss().

const double alpha = 0.09

Courant parameter for the ADI scheme

Definition at line 164 of file constants.h.

Referenced by TCREvolutorADI::TCREvolutorADI().

const double ALPHAf = 1.0/137.035989561

Fine structure constant.

Definition at line 188 of file constants.h.

const int antiproton_cs = 2

antiproton cross section option.

Definition at line 240 of file constants.h.

Referenced by TSpallationNetwork::TSpallationNetwork().

const double bet_e = sqrt(2.*BK*Te/Mele)

Some parameter

Definition at line 212 of file constants.h.

const double Bh = 5.

[muG] intensity of MF at Sun's position.

Definition at line 169 of file constants.h.

Referenced by DRAGON::Print(), and TBField::TBField().

const double BK = 1.38066e-23/1.60218e-13

[MeV/K] Bolzmann constant

Definition at line 210 of file constants.h.

const std::string BNLdata = "data/bnl_data.dat"

Nuclear properties database.

Definition at line 131 of file constants.h.

Referenced by TNucleiList::TNucleiList().

const double C = 2.99792458e10

[cm s^-1] Speed of light in vacuum.

Definition at line 134 of file constants.h.

Referenced by TParticle::FindNormalization(), TParticle::Print(), TParticle::PrintSpectrum(), DRAGON::Run(), TBremsstrahlungLoss::TBremsstrahlungLoss(), and TISRF::TISRF().

const double Clight = C * Myr

[cm Myr^-1]

Definition at line 142 of file constants.h.

Referenced by TGalpropXSec::GetXSec(), TGalpropXSec::Print(), TInelasticCrossSection::TInelasticCrossSection(), and TSpallationNetwork::TSpallationNetwork().

const int cross_section_option = 12

Option related to Galprop cross section model.

Definition at line 152 of file constants.h.

Referenced by TGalpropXSec::GetXSec().

const double cutoff_rig = 5.e3

[GV] Rigidity of cut off of primary electron injection spectra.

Definition at line 233 of file constants.h.

Referenced by TSpectrum::TSpectrum().

const double D_ref_rig = 3.

[GV] reference rigidity for D.

Definition at line 170 of file constants.h.

Referenced by DRAGON::Print(), and TDiffusionCoefficient::TDiffusionCoefficient().

const int diff_reacc = 1

Use diffusion reacceleration model: 1 = Seo & Ptuskin 94; 2 = Ptuskin 2003.

Definition at line 174 of file constants.h.

Referenced by TReaccelerationCoefficient::TReaccelerationCoefficient().

const DMreaction DMr = Annihilation

Set DM reaction channel.

Definition at line 113 of file constants.h.

Referenced by TDMSource::TDMSource().

const DMspectrum DMs = DarkSUSY

Set DM reaction model.

Definition at line 114 of file constants.h.

Referenced by TSpectrum::TSpectrum().

const double dtfactor = 0.25

Factor of reduction of timestep after

See also:

Nrept iterations.

Definition at line 157 of file constants.h.

Referenced by TCREvolutor::Run().

const double dtmax = 1e3

[Myr] Maximum (starting) value of timestep.

Definition at line 159 of file constants.h.

Referenced by TCREvolutor::Run().

const double dtmin = 1e-6

[Myr] Minimum value of timestep.

Definition at line 158 of file constants.h.

Referenced by TCREvolutor::Run().

const double dzzGal = 0.01

[kpc] value to average over the gas distribution, used in

See also:

GasType::Galprop.

Definition at line 173 of file constants.h.

Referenced by TDMSource::DM_profile_av(), TGas::TGas(), TH2Gas::TH2Gas(), THIGas::THIGas(), and THIIGas::THIIGas().

const double EH = 19.e-6

[MeV] H eff. ioniz. potential (in nucl. loss)

Definition at line 206 of file constants.h.

Referenced by TIonizationLoss::TIonizationLoss().

const double EHe = 44.e-6

[MeV] He eff. ioniz. potential (in nucl. loss)

Definition at line 207 of file constants.h.

Referenced by TIonizationLoss::TIonizationLoss().

const double EiH = 13.6e-6

[MeV] Hydrogen (in electron EL)

Definition at line 198 of file constants.h.

Referenced by TIonizationLoss::TIonizationLoss().

const double EiHe = 24.59e-6

[MeV] Helium (in electron EL)

Definition at line 199 of file constants.h.

Referenced by TIonizationLoss::TIonizationLoss().

const double EkDelta = 100.0

Ek in case of Delta injection for DM

Definition at line 116 of file constants.h.

Referenced by TSpectrum::TSpectrum().

const double Ekmax = 3.e3

[Gev/n] Maximal energy.

Definition at line 230 of file constants.h.

Referenced by DRAGON::Print(), and TGrid::TGrid().

const double Ekmin = 0.01

[GeV/n] Minimal energy.

Definition at line 231 of file constants.h.

Referenced by DRAGON::Print(), and TGrid::TGrid().

const std::string ElTablefile = "data/Electron_production.dat"

Tabulated electron yields.

Definition at line 127 of file constants.h.

Referenced by TSpallationNetwork::InitDataTables().

const double erg_to_eV = 1./eV_to_erg

Definition at line 147 of file constants.h.

Referenced by TBField::GetEnDensity(), and TISRF::TISRF().

const double eV_to_erg = 1.60217733e-12

Definition at line 146 of file constants.h.

Referenced by TICSLoss::TICSLoss().

const bool fullstore = true

Set whether to store full spatial density of propagated CRs.

Definition at line 122 of file constants.h.

Referenced by DRAGON::DRAGON().

const double gam1 = 100.

Boosts

Definition at line 218 of file constants.h.

Referenced by TBremsstrahlungLoss::TBremsstrahlungLoss().

const double gam2 = 800.

Boosts

Definition at line 219 of file constants.h.

Referenced by TBremsstrahlungLoss::TBremsstrahlungLoss().

const GasType gas_model = Galprop

Set Gas model.

Definition at line 181 of file constants.h.

Referenced by DRAGON::Print(), TGas::TGas(), TH2Gas::TH2Gas(), THIGas::THIGas(), and THIIGas::THIIGas().

const double H2PiC = 0.19732705359e-10

[MeV*cm] =hc/(2Pi) conversion const.

Definition at line 192 of file constants.h.

Referenced by TCoulombLoss::TCoulombLoss().

const double h_planck = 6.6260755e-27

[erg sec] Planck's constant.

Definition at line 135 of file constants.h.

Referenced by TICSLoss::TICSLoss(), and TISRF::TISRF().

const double He_abundance = 0.11

Helium abundance in the ISM.

Definition at line 185 of file constants.h.

Referenced by TXSecBase::GetHefactor(), TGalpropXSec::GetXSec(), TWebber03::Print(), DRAGON::Print(), TBremsstrahlungLoss::TBremsstrahlungLoss(), TInelasticCrossSection::TInelasticCrossSection(), TIonizationLoss::TIonizationLoss(), and TSpallationNetwork::TSpallationNetwork().

const double ind_low_nucl = 2.75

Source injection psectral index for primary nuclei below

See also:

sp_ref_rig.

Definition at line 237 of file constants.h.

Referenced by Galaxy::Galaxy().

const std::string ISRFfile = "data/MilkyWay_DR0.5_DZ0.1_DPHI10_RMAX20_ZMAX5_galprop_format.fits.gz"

ISRF model.

Definition at line 126 of file constants.h.

Referenced by Galaxy::Galaxy().

const double km = 1e5

Definition at line 140 of file constants.h.

Referenced by DRAGON::Print().

const double kpc = 3.08568e21

Definition at line 139 of file constants.h.

Referenced by TParticle::FindNormalization(), Utility::nH2_Gal(), TParticle::Print(), DRAGON::Print(), TParticle::PrintSpectrum(), DRAGON::Run(), and TDMSource::TDMSource().

const LeptonYields ly = GalpropTable

Set lepton yields.

Definition at line 120 of file constants.h.

Referenced by TSpallationNetwork::TSpallationNetwork().

const double Mele = 0.5109990615

[MeV/c^2] electron rest mass.

Definition at line 189 of file constants.h.

Referenced by TBremsstrahlungLoss::TBremsstrahlungLoss(), TCoulombLoss::TCoulombLoss(), TIonizationLoss::TIonizationLoss(), and TISRF::TISRF().

const double MeleGeV = 1e-3*Mele

[GeV/c^2] electron rest mass.

Definition at line 190 of file constants.h.

Referenced by TGrid::TGrid(), and TSpectrum::TSpectrum().

const double MH = 1.67e-24

[grams] the mass of H atom.

Definition at line 216 of file constants.h.

Referenced by TBremsstrahlungLoss::TBremsstrahlungLoss().

const double MHe = 1.66e-24 *4.

[grams] the mass of He atom.

Definition at line 217 of file constants.h.

Referenced by TBremsstrahlungLoss::TBremsstrahlungLoss().

const double mp = .938

Definition at line 144 of file constants.h.

Referenced by TGalpropXSec::antiproton_cc1(), TGalpropXSec::nucleon_cs(), DRAGON::particle_modulated_spectrum(), tan_ng(), TCoulombLoss::TCoulombLoss(), TGrid::TGrid(), TIonizationLoss::TIonizationLoss(), and TSpectrum::TSpectrum().

const double Myr = 1e6*year

Definition at line 138 of file constants.h.

Referenced by DRAGON::Print(), TBremsstrahlungLoss::TBremsstrahlungLoss(), TCoulombLoss::TCoulombLoss(), TDMSource::TDMSource(), TICSLoss::TICSLoss(), TIonizationLoss::TIonizationLoss(), TNucleiList::TNucleiList(), TParticle::TParticle(), and TSynchrotronLoss::TSynchrotronLoss().

const int Nrept = 20

Number of iterations for a single timestep value

See also:

dtfactor

dtmin

dtmat. Only used in class

TCREvolutor .

Definition at line 156 of file constants.h.

Referenced by TCREvolutor::Run().

const bool OpSplit = true

Controls whether to use the OS (Galprop) scheme.

Definition at line 162 of file constants.h.

Referenced by DRAGON::DRAGON().

const double p = 3.0

Number of differential operators to be updated.

Definition at line 151 of file constants.h.

Referenced by DRAGON::Print(), TCREvolutorADI::Run(), and TCREvolutor::Run().

const bool partialstore = true

Set whether to store only CR spectra at Solar system position.

Definition at line 123 of file constants.h.

Referenced by DRAGON::DRAGON().

const double Pi = 3.14159265358979312

Definition at line 145 of file constants.h.

Referenced by TGalpropXSec::antiproton_cc1(), TBField::GetEnDensity(), TGalpropXSec::nucleon_cs(), DRAGON::Run(), TCoulombLoss::TCoulombLoss(), and TISRF::TISRF().

const double PIR02C = Pi*Rele*Rele*C

[cm^3/s]

Definition at line 195 of file constants.h.

Referenced by TSynchrotronLoss::TSynchrotronLoss().

const double PIR02MC2C = Pi*Rele*Rele*C*(Mele*1.e6)

[eV*cm^3/s] = Pi*e^4/mc

Definition at line 194 of file constants.h.

Referenced by TCoulombLoss::TCoulombLoss(), and TIonizationLoss::TIonizationLoss().

const double PIR0H2C2 = Pi*Rele*H2PiC*H2PiC

[MeV^2*cm^3]

Definition at line 193 of file constants.h.

Referenced by TCoulombLoss::TCoulombLoss().

const std::string PosTablefile = "data/Positron_production.dat"

Tabulated positron yields.

Definition at line 128 of file constants.h.

Referenced by TSpallationNetwork::InitDataTables().

const bool prop_ap = false

Whether to propagate antiprotons.

Definition at line 246 of file constants.h.

Referenced by TNucleiList::TNucleiList(), and TSpallationNetwork::TSpallationNetwork().

const bool prop_el = false

Whether to propagate electrons and positrons.

Definition at line 247 of file constants.h.

Referenced by Galaxy::Galaxy(), TNucleiList::TNucleiList(), and TSpallationNetwork::TSpallationNetwork().

const bool prop_extracomp = false

Whether to propagate the extra component (electron/positron).

Warning:

Be careful! if == true then ONLY the extra component is propagated, and the other settings are overridden.

Definition at line 248 of file constants.h.

Referenced by Galaxy::Galaxy(), DRAGON::Run(), and TNucleiList::TNucleiList().

const double r0 = 2.

[kpc] radial position of maximum of D.

Definition at line 172 of file constants.h.

Referenced by TDiffusionCoefficient::GetProfile().

const double rB = 10.

[kpc] MF radial scale length.

Definition at line 167 of file constants.h.

Referenced by DRAGON::Print(), and TBField::TBField().

const double rd = 2.

[kpc] radial length scale of D.

Definition at line 171 of file constants.h.

Referenced by TDiffusionCoefficient::GetProfile().

const double Rele = 2.8179409238e-13

[cm] =e^2/mc^2 classical electron radius.

Definition at line 191 of file constants.h.

Referenced by TCoulombLoss::TCoulombLoss(), and TISRF::TISRF().

const double rhos = 0.36

DM solar system density GeV/cm^3

Definition at line 115 of file constants.h.

Referenced by TDMSource::TDMSource().

const double ringmax = 4.

[kpc] Outer radius of the ring if option

See also:

SNRType::Rings is selected.

Definition at line 184 of file constants.h.

Referenced by TAstrophysicalSource::TAstrophysicalSource().

const double ringmin = 3.

[kpc] Inner radius of the ring if option

See also:

SNRType::Rings is selected.

Definition at line 183 of file constants.h.

Referenced by TAstrophysicalSource::TAstrophysicalSource().

const double Rmax = 20.

[kpc] Radial end of the Galaxy.

Definition at line 180 of file constants.h.

Referenced by TParticle::FindNormalization(), DRAGON::Print(), TParticle::PrintSpectrum(), and TGrid::TGrid().

const double Rmin = 0.

[kpc]

Definition at line 179 of file constants.h.

Referenced by TParticle::FindNormalization(), DRAGON::Print(), TParticle::PrintSpectrum(), and TGrid::TGrid().

const double robs = 8.5

[kpc] value of r at Sun's position.

Definition at line 177 of file constants.h.

Referenced by TDMSource::DM_profile(), TParticle::FindNormalization(), TDiffusionCoefficient::GetProfile(), DRAGON::particle_modulated_spectrum(), DRAGON::Print(), TParticle::PrintSpectrum(), TAstrophysicalSource::TAstrophysicalSource(), and TBField::TBField().

const bool scaling = false

use the propagated He spectrum to compute He contribution to antiproton spectrum, or use a scaling relation from the proton spectrum.

Definition at line 241 of file constants.h.

Referenced by TSpallationNetwork::TSpallationNetwork().

const double sigmaT = 6.6524e-25/(kpc*kpc)

[kpc^2] Thomson cross section

Definition at line 221 of file constants.h.

const double small = 1e-45

Definition at line 117 of file constants.h.

const SNRType SNR_model = Ferriere

Set Source distribution model.

Definition at line 182 of file constants.h.

Referenced by DRAGON::Print(), and TAstrophysicalSource::TAstrophysicalSource().

const double sp_ref_rig = 9.

[GV] energy for abundance normalization or breaks in the sources.

Definition at line 232 of file constants.h.

Referenced by DRAGON::Print(), and TSpectrum::TSpectrum().

const double sp_ref_rig_break_el = 4.

[GV] energy for break in electron injection spectrum.

Definition at line 236 of file constants.h.

Referenced by TSpectrum::TSpectrum().

const double sp_ref_rig_el = 33.20

[GeV] energy for normalization of the final electron flux.

Definition at line 225 of file constants.h.

Referenced by TParticle::FindNormalization().

const double sp_ref_rig_norm = 1.e2

[GeV] energy for normalization of the final hadron flux.

Definition at line 224 of file constants.h.

Referenced by TParticle::FindNormalization(), and DRAGON::Print().

const SpallXSec spallationxsec = GalpropXSec

Set spallation cross section model.

Definition at line 119 of file constants.h.

Referenced by DRAGON::Run(), TInelasticCrossSection::TInelasticCrossSection(), and TSpallationNetwork::TSpallationNetwork().

const double spect_ind_Al = 2.75

Definition at line 332 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Ar = 2.75

Definition at line 337 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_B = 2.75

Definition at line 324 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Be = 2.75

Definition at line 323 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_C = 2.75

Definition at line 325 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Ca = 2.75

Definition at line 339 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Cl = 2.75

Definition at line 336 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Co = 2.75

Definition at line 346 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Cr = 2.75

Definition at line 343 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_el = 2.25

Source injection spectral index for primary electrons above

See also:

sp_ref_rig_break_el.

Definition at line 235 of file constants.h.

Referenced by TSpectrum::TSpectrum().

const double spect_ind_el_low = 2.25

Source injection spectral index for primary electrons below

See also:

sp_ref_rig_break_el.

Definition at line 234 of file constants.h.

Referenced by TSpectrum::TSpectrum().

const double spect_ind_F = 2.75

Definition at line 328 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Fe = 2.75

Definition at line 345 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_He = 2.75

Definition at line 322 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_K = 2.75

Definition at line 338 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Mg = 2.75

Definition at line 331 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Mn = 2.75

Definition at line 344 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_N = 2.75

Definition at line 326 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Na = 2.75

Definition at line 330 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Ne = 2.75

Definition at line 329 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Ni = 2.75

Definition at line 347 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_O = 2.75

Definition at line 327 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_P = 2.75

Definition at line 334 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_pr = 2.75

< Propagated spectral indexes above

See also:

sp_ref_rig.

Definition at line 321 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_S = 2.75

Definition at line 335 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Sc = 2.75

Definition at line 340 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Si = 2.75

Definition at line 333 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_Ti = 2.75

Definition at line 341 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_ind_V = 2.75

Definition at line 342 of file constants.h.

Referenced by Galaxy::Galaxy().

const double spect_norm = 5.e-2

[m^-2 sr^-1 s^-1 GeV^-1] normalization at

See also:

sp_ref_rig_norm from Biermann A&A 330, 389 (1998).

Definition at line 226 of file constants.h.

Referenced by TParticle::FindNormalization(), and DRAGON::Print().

const double spect_norm_el = 4.2e-3

[m^-2 sr^-1 s^-1 GeV^-1] normalization at

See also:

sp_ref_rig_el from Biermann A&A 330, 389 (1998).

Definition at line 227 of file constants.h.

Referenced by TParticle::FindNormalization().

const double t_half_limit = 1.e4

[yr] Lifetime below which nuclei are considered completely unstable.

Definition at line 153 of file constants.h.

Referenced by TGalpropXSec::GetXSec().

const double Te = 1.0e4

[K] temperature of warm ionized medium

Definition at line 211 of file constants.h.

const double TH = 62.8

[g/cm^2] Hydrogen

Definition at line 200 of file constants.h.

Referenced by TBremsstrahlungLoss::TBremsstrahlungLoss().

const double THe = 93.1

[g/cm^2] Helium

Definition at line 201 of file constants.h.

Referenced by TBremsstrahlungLoss::TBremsstrahlungLoss().

const double tol = 3e-4

Tolerance for convergence of the solution. Used only when ADI == true

Definition at line 163 of file constants.h.

Referenced by TCREvolutorADI::TCREvolutorADI().

const double u = 1./200.

[kpc] just a little value to regularize diffusion behavior near GC.

Definition at line 150 of file constants.h.

Referenced by DRAGON::Print(), TDiffusionCoefficient::TDiffusionCoefficient(), and TSpallationNetwork::TSpallationNetwork().

const std::string Webber03Data = "data/webber_xsec.dat"

Tabulated Webber 2003 spallation cross sections.

Definition at line 129 of file constants.h.

Referenced by TWebber03::TWebber03().

const std::string Webber03DataTotal = "data/webber_xsec_total.dat"

Tabulated Webber 2003 inelastic cross sections.

Definition at line 130 of file constants.h.

Referenced by TWebber03::TWebber03().

const double xm = pow(pow(3.*sqrt(Pi)/4.,1./3.) *bet_e,3)

Some other parameter

Definition at line 213 of file constants.h.

Referenced by TCoulombLoss::TCoulombLoss().

const double year = 365.25*24.*3600.

Definition at line 137 of file constants.h.

Referenced by TGalpropXSec::GetXSec(), TGalpropXSec::nucdata(), and TNucleiList::TNucleiList().

const double ZHe = 2.

He charge.

Definition at line 203 of file constants.h.

Referenced by TBremsstrahlungLoss::TBremsstrahlungLoss(), and TIonizationLoss::TIonizationLoss().

const int Zmax = 14

Maximal propagated charge.

Definition at line 244 of file constants.h.

Referenced by TNucleiList::TNucleiList().

const int Zmin = 1

Minimal propagated charge.

Definition at line 245 of file constants.h.

Referenced by TNucleiList::TNucleiList().

const double zobs = 0.

[kpc] value of z at Sun's position.

Definition at line 178 of file constants.h.

Referenced by TParticle::FindNormalization(), TDiffusionCoefficient::GetProfile(), DRAGON::particle_modulated_spectrum(), DRAGON::Print(), and TParticle::PrintSpectrum().

const double zr = 2.0

[kpc] regular MF vertical scale.

Definition at line 168 of file constants.h.

Referenced by DRAGON::Print(), and TBField::TBField().