constants.h

Go to the documentation of this file.

#ifndef _CONSTANTS_H
#define _CONSTANTS_H

#define REAC

#define CONVECTION


#include "../config.h"

#include <string>
#include <cmath>

// DM

typedef enum {
  Decay, 
  Annihilation 
} DMreaction;

typedef enum {
  Delta, 
  DarkSUSY 
} DMspectrum;

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

typedef enum {
  BronfFerr, 
  NS, 
  Galprop 
} GasType; 

typedef enum {
  Galprop_, 
  Ferriere, 
  PointSource, 
  Ring 
} SNRType; 

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

typedef enum {
  GalpropXSec, 
  Webber03 
} SpallXSec;

typedef enum {
  GalpropTable, 
  Kamae 
} LeptonYields;

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; 
// Data files
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"; 
// Fundamental constants
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 ; //272029          ; /**< [GeV] Proton mass. */
const double Pi        = 3.14159265358979312 ;
const double eV_to_erg = 1.60217733e-12 ;
const double erg_to_eV = 1./eV_to_erg ;

// Utils
const double u                 = 1./200. ; 
const double p                 = 3.0 ;     
const int cross_section_option = 12 ; 
const double t_half_limit      = 1.e4 ; 
// Algorithm
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; 
// Galactic Magnetic Field
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 ; 
// Galaxy properties
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 ;
// Energy losses
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.; 
// Ionization
const double EH = 19.e-6; 
const double EHe= 44.e-6; 
// Coulomb
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); 
// Bremsstrahlung
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) ; 
// CR parameters
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 ; 
// Parameter space
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; 
// Antiprotons
const int antiproton_cs = 2 ; 
const bool scaling = false ; 
// Limits for the nuclear chain
const int Zmax     = 14 ; 
const int Zmin     = 1 ; 
const bool prop_ap = false ; 
const bool prop_el = false ; 
const bool prop_extracomp = false;  
// Source abundances (GALPROP-like)
const double ab_pr   = 1.06e6;  // 1e5;
const double ab_He3   = 9.033;   // 8511.;
const double ab_He4   = 7.199e4;   // 8511.;
const double ab_Li   = 0.;       // 1.12e-6;
const double ab_Be9   = 0.;       // 2.4e-6;
const double ab_Be10   = 0.;       // 2.4e-6;
const double ab_B10    = 0.;       // 5.01e-5;
const double ab_B11    = 0.;       // 5.01e-5;
//const double ab_C    = 2819.;
const double ab_C13  = 5.268e-07;
const double ab_C14  = 0.0;
//const double ab_N    = 182.8;    
const double ab_N15  = 5.961e-5;      // 0.;
const double ab_O    = 3822.;    // 3450.; //45.7*0.54;
const double ab_O17  = 6.713e-7;      // 0.;
const double ab_O18  = 1.286;      // 0.;
const double ab_F19  = 2.664e-8;      // 0.0036;      
const double ab_Ne20 = 312.5;    // 6.;
const double ab_Ne21 = 0.003556;      // 0.;
const double ab_Ne22 = 100.1;     // 0.91;
const double ab_Na22 = 0.0;
const double ab_Na23 = 22.84;     // 0.15;     
const double ab_Mg24 = 658.1;    // 3.39;      
const double ab_Mg25 = 82.5;     // 0.;
const double ab_Mg26 = 104.7;     // 0.;
const double ab_Al26 = 0.0;
const double ab_Al27 = 76.42;     // 0.23;     
const double ab_Si28 = 725.7 ;   // 3.23;      
const double ab_Si29 = 35.02;     // 0.;
const double ab_Si30 = 24.68;     // 0.;
const double ab_Si32 = 24.68;     // 0.;
const double ab_P31  = 7.28;     // 0.02;
const double ab_S32  = 94.40;    // 1.38;     
const double ab_S33  = 0.8;      // 0.;
const double ab_S34  = 4.45;     // 0.;
const double ab_S36  = 0.02;     // 0.;
const double ab_Cl35 = 2.11;     // 0.03;      
const double ab_Cl37 = 0.72;     // 0.;
const double ab_Ar36 = 11.17;    // 0.15;      
const double ab_Ar38 = 2.22;     // 0.;
const double ab_K39  = 4.47;     // 0.01;   
const double ab_Ca40 = 39.87;    // 0.2;
const double ab_Ca42 = 0.28;     // 0.;
const double ab_Ca44 = 0.93;     // 0.;
const double ab_Ca48 = 0.09;     // 0.;
const double ab_Sc45 = 0.121;    // .0001;     
const double ab_Ti46 = 0.17;     // 0.;
const double ab_Ti47 = 0.16;     // 0.;
const double ab_Ti48 = 1.64;     // 0.008;
const double ab_Ti49 = 0.124;    // 0.;
const double ab_Ti50 = 0.12;     // 0.;
const double ab_V51  = 0.0;      // 0.001;   
const double ab_Cr50 = 0.72 ;    // 0.;     
const double ab_Cr52 = 14.51;    // 0.04;
const double ab_Cr53 = 1.69;     // 0.;
const double ab_Cr54 = 0.43;     // 0.;
const double ab_Mn55 = 16.15;    // 0.02;     
const double ab_Fe54 = 38.00;    // 0.;     
const double ab_Fe56 = 620.2;    // 2.82;
const double ab_Fe57 = 15.07;    // 0.;
const double ab_Fe58 = 2.31;     // 0.;
const double ab_Co59 = 1.24;     // 0.008;    
const double ab_Ni58 = 26.21;    // 0.12;    
const double ab_Ni60 = 10.44;    // 0.05;
const double ab_Ni61 = 0.48;     // 0.;
const double ab_Ni62 = 1.51;     // 0.;
const double ab_Ni64 = 0.46;     // 0.;

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;

/*
// Observed spectra at Earth (Wiebel-Sooth, Biermann...)
const double spect_ind_pr = 2.77;
const double spect_ind_He = 2.64;
const double spect_ind_C = 2.66;
const double spect_ind_N = 2.72;
const double spect_ind_O = 2.68;
const double spect_ind_F = 2.69;
const double spect_ind_Ne = 2.64;
const double spect_ind_Na = 2.66;
const double spect_ind_Mg = 2.64;
const double spect_ind_Al = 2.66;
const double spect_ind_Si = 2.75;
const double spect_ind_P = 2.69;
const double spect_ind_S = 2.55;
const double spect_ind_Cl = 2.68;
const double spect_ind_Ar = 2.64;
const double spect_ind_K = 2.65;
const double spect_ind_Ca = 2.70;
const double spect_ind_Sc = 2.64;
const double spect_ind_Ti = 2.61;
const double spect_ind_V = 2.63;
const double spect_ind_Cr = 2.67;
const double spect_ind_Mn = 2.46;
const double spect_ind_Fe = 2.60;
const double spect_ind_Co = 2.72;
const double spect_ind_Ni = 2.751;
*/

#endif