#include <math.h>#include <CLHEP/Random/RandomEngine.h>#include <CLHEP/Random/RandGeneral.h>#include <CLHEP/Random/JamesRandom.h>#include "CrElectronPrimary.h"Go to the source code of this file.
Typedefs | |
| typedef double | G4double |
| CrElectronPrimary: The primary particle generator for primary cosmic-ray electrons. More... | |
Functions | |
| G4double | beta (G4double E ) |
| G4double | rigidity (G4double E ) |
| G4double | energy (G4double rigidity ) |
| G4double | geomag_cut (G4double E , G4double cor ) |
| G4double | org_spec (G4double E ) |
| G4double | mod_spec (G4double E , G4double phi ) |
| G4double | primaryCRspec (G4double E , G4double cor , G4double phi ) |
| G4double | primaryCRenvelope1 (G4double E , G4double cor , G4double phi ) |
| G4double | primaryCRenvelope1_integral (G4double E , G4double cor , G4double phi ) |
| G4double | primaryCRenvelope1_integral_rev (G4double value, G4double cor , G4double phi ) |
| G4double | primaryCRenvelope2 (G4double E , G4double cor , G4double phi ) |
| G4double | primaryCRenvelope2_integral (G4double E , G4double cor , G4double phi ) |
| G4double | primaryCRenvelope2_integral_rev (G4double value, G4double cor , G4double phi ) |
| G4double | primaryCRenergy (HepRandomEngine* engine, G4double cor, G4double solarPotential) |
Variables | |
| const G4double | pi = 3.14159265358979323846264339 |
| const G4double | restE = 5.11e-4 |
| const G4double | lowE_primary = 1.0 |
| const G4double | highE_primary = 100.0 |
| const G4double | INTEGRAL_primary = 4.3643 |
| const G4double | A_primary = 7.23e-1 |
| Generate a random distribution of primary cosmic ray protons j(E) = mod_spec(E, phi) * beta(E) * geomag_cut(E, CutOff) mod_spec(E, phi) = org_spec(E+phi*1e-3) * ((E+restE)**2 - restE**2)/((E+restE+phi*1e-3)**2-restE**2) org_spec(E) = A * rigidity(E)**-a a = 2.79 rigidity(E) = sqrt((E+restE)**2 - restE**2) beta(E) = sqrt(1 - (E/restE+1)**-2) geomag_cut(E, CutOff) = 1/(1 + (rigidity(E)/CutOff)**-12.0) CutOff = 4.44 for Theta_M = 0.735 and altitude = 50km (balloon experiment) phi = 540, 1100 [MV] for Solar minimum, maximum E: [GeV] j: [m**-2 s**-1 sr**-1 MeV**-1 ]. More... | |
| const G4double | a_primary = 3.33 |
|
|
CrElectronPrimary: The primary particle generator for primary cosmic-ray electrons. Ver 1.0 on 2001-04-18 by Masanobu Ozaki <ozaki@astro.isas.ac.jp> separate the primary component from the all-in-one code. 2001-04-26 comments are added by T. Mizuno
Definition at line 25 of file CrElectronPrimary.cxx. |
|
|
Definition at line 45 of file CrElectronPrimary.cxx. |
|
|
Definition at line 67 of file CrElectronPrimary.cxx. |
|
|
Definition at line 112 of file CrElectronPrimary.cxx. |
|
|
Definition at line 126 of file CrElectronPrimary.cxx. |
|
|
Definition at line 119 of file CrElectronPrimary.cxx. |
|
|
Definition at line 208 of file CrElectronPrimary.cxx. |
|
|
Definition at line 151 of file CrElectronPrimary.cxx. |
|
|
Definition at line 163 of file CrElectronPrimary.cxx. |
|
|
Definition at line 175 of file CrElectronPrimary.cxx. Referenced by primaryCRenergy(). |
|
|
Definition at line 187 of file CrElectronPrimary.cxx. |
|
|
Definition at line 194 of file CrElectronPrimary.cxx. |
|
|
Definition at line 201 of file CrElectronPrimary.cxx. Referenced by primaryCRenergy(). |
|
|
Definition at line 134 of file CrElectronPrimary.cxx. |
|
|
Definition at line 56 of file CrElectronPrimary.cxx. |
|
|
Generate a random distribution of primary cosmic ray protons j(E) = mod_spec(E, phi) * beta(E) * geomag_cut(E, CutOff) mod_spec(E, phi) = org_spec(E+phi*1e-3) * ((E+restE)**2 - restE**2)/((E+restE+phi*1e-3)**2-restE**2) org_spec(E) = A * rigidity(E)**-a a = 2.79 rigidity(E) = sqrt((E+restE)**2 - restE**2) beta(E) = sqrt(1 - (E/restE+1)**-2) geomag_cut(E, CutOff) = 1/(1 + (rigidity(E)/CutOff)**-12.0) CutOff = 4.44 for Theta_M = 0.735 and altitude = 50km (balloon experiment) phi = 540, 1100 [MV] for Solar minimum, maximum E: [GeV] j: [m**-2 s**-1 sr**-1 MeV**-1 ]. References: org_spec: Komori, Y. et al. 1999, Proc. of Dai-Kikyu (large balloon) Sympo. Heisei-11yr, 33-36 (Fig. 2) mod_spec: Gleeson, L. J. and Axford, W. I. 1968, ApJ, 154, 1011-1026 (Eq. 11) geomag_cut formula: an eyeball fitting function of the ratio of BESS98 and AMS98 proton spectra. CutOff: calculated as (Rc/GV) = 14.9 * (1+h/R)^-2 * (cos(theta_M))^4, where h gives altitude from earth surface, R means an radius of earth, and theta_M is geomagnetic lattitude. See references below. "Handbook of space astronomy and astrophysics" 2nd edition, p225 (Zombeck, 1990, Cambridge University Press) "High Energy Astrophysics" 2nd edition, p325-330 (M. S. Longair, 1992, Cambridge University Press) Definition at line 90 of file CrProtonPrimary.cxx. |
|
|
Definition at line 42 of file CrProtonPrimary.cxx. |
|
|
Definition at line 91 of file CrProtonPrimary.cxx. |
|
|
Definition at line 33 of file CrProtonPrimary.cxx. |
|
|
Definition at line 32 of file CrProtonPrimary.cxx. |
|
|
Definition at line 29 of file CrProtonSplash.cxx. |
|
|
Definition at line 31 of file CrProtonSplash.cxx. |
1.2.3 written by Dimitri van Heesch,
© 1997-2000