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.

00054              {
00055   ISO=0, 
00056   NFW=1, 
00057   Kra=2, 
00058   Moore=3, 
00059   Einasto=4 
00060 } 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.

00036              {
00037   Decay, 
00038   Annihilation 
00039 } 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.

00045              {
00046   Delta, 
00047   DarkSUSY 
00048 } 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.

00087              {
00088   Constant, 
00089   Exp, 
00090   Expr, 
00091   Radial, 
00092   ExpRadial 
00093 } 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.

00066              {
00067   BronfFerr, 
00068   NS, 
00069   Galprop 
00070 } GasType; 

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.

00108              {
00109   GalpropTable, 
00110   Kamae 
00111 } 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.

00076              {
00077   Galprop_, 
00078   Ferriere, 
00079   PointSource, 
00080   Ring 
00081 } 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.

00099              {
00100   GalpropXSec, 
00101   Webber03 
00102 } 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
const double Clight = C * Myr
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
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
const LeptonYields ly = GalpropTable

Set lepton yields.

Definition at line 120 of file constants.h.

Referenced by TSpallationNetwork::TSpallationNetwork().

const double Mele = 0.5109990615
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
const double Myr = 1e6*year
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
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.
const double robs = 8.5
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.
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.
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().

 All Classes Namespaces Files Functions Variables Enumerations Enumerator
Generated on Mon Sep 27 12:59:56 2010 for DRAGON by  doxygen 1.6.3