Message: How to input GCR spectra? (multi source GPS) Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question How to input GCR spectra? (multi source GPS) 

Forum: Particles
Date: 11 May, 2017
From: <Norikonbu>

Hi, now I am trying to input multi-source GCR spectra (souce:CREME 96) by GPS command like this:

======================================================================

# =============================
# CREME96 GCR source
# Particle type: proton
# =============================
/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/pos/centre 0. 0. 0. cm
/gps/pos/radius 300 cm
#/gps/source/multiplevertex true 9.664E+11
/gps/particle proton
/gps/ene/type Arb
/gps/ene/diffspec true
/gps/hist/type arb
/gps/ene/min 3.0000 MeV
/gps/ene/max 85675.0 MeV
/gps/source/intensity 3725.431908
gps/hist/point 3.0237 9.29999e+08
gps/hist/point 7.55925 5.70937e+06

(omitting)

gps/hist/point 25197.5 9590.15
gps/hist/point 35276.5 2780.31
gps/hist/point 45355.5 1070.65
gps/hist/point 55434.5 508.812
gps/hist/point 65513.5 295.687
gps/hist/point 75592.5 149.791
gps/hist/point 85671.5 107.317


/gps/hist/inter Lin

/gps/ang/type cos
/gps/ang/mintheta 0.0E+00 deg
/gps/ang/maxtheta 9.0E+01 deg

# =============================
# CREME96 GCR source
# Particle type: He
# =============================
/gps/source/add 0.094
/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/pos/centre 0. 0. 0. cm # Value changed from the ADM example
/gps/pos/radius 300 cm
/gps/particle ion
/gps/ion 2 4
#/gps/source/multiplevertex true
/gps/ene/type Arb
/gps/ene/diffspec true
/gps/hist/type arb
/gps/ene/min 12.0000 MeV
/gps/ene/max 340225.00 MeV
#/gps/source/intensity 3.289E+13
gps/hist/point 12.01 6.10314e+07
gps/hist/point 30.02 2.32241e+06
gps/hist/point 60.04 2.02840e+06
gps/hist/point 120.08 2.02216e+06
gps/hist/point 200.13 9.63785e+05
gps/hist/point 280.18 7.02448e+05

(omitting)

gps/hist/point 220143.00 2.70177e+0
gps/hist/point 260169.00 1.55823e+01
gps/hist/point 300195.00 7.91681e+00
gps/hist/point 340221.00 5.65487e+00




/gps/hist/inter Lin

/gps/ang/type cos
/gps/ang/mintheta 0.0E+00 deg
/gps/ang/maxtheta 9.0E+01 deg

/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/pos/centre 0. 0. 0. cm # Value changed from the ADM example
/gps/pos/radius 300 cm

# =============================
# CREME96 GCR source
# Particle type: Li
# =============================
/gps/source/add 0.00036
/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/pos/centre 0. 0. 0. cm # Value changed from the ADM example
/gps/pos/radius 300 cm
/gps/particle ion
/gps/ion 3 7
/gps/ene/type Arb
/gps/ene/diffspec true
/gps/hist/type arb
/gps/ene/min 36.0000 MeV
/gps/ene/max 300270 MeV
#/gps/source/intensity 1.828E+10
gps/hist/point 20.82 6.92024e+04
gps/hist/point 52.06 2.93978e+03
gps/hist/point 104.11 1.88659e+03
gps/hist/point 208.23 2.74274e+03
gps/hist/point 347.05 2.38359e+03
gps/hist/point 485.87 2.14319e+03

(omitting)

gps/hist/point 58998.50 1.00531e+01
gps/hist/point 65939.50 6.78584e+00
gps/hist/point 86762.50 1.38230e+01
gps/hist/point 121467.50 4.02124e+00
gps/hist/point 173525.00 2.63894e+00




/gps/hist/inter Lin

/gps/ang/type cos
/gps/ang/mintheta 0.0E+00 deg
/gps/ang/maxtheta 9.0E+01 deg

==========================================================

I've confirmed that these commands work properly inserting another spectra(histogram) and different values of "/gps/source/add ".
However,using above values of spectra and "/gps/source/add ", Geant4 always aborts.

Please help me with the problem.
Is the problem merely due to values (eg numerical format etc)? Or is it due to a more fundamental cause?

Thanks, Noriko

 Add Message Add Message
to: "How to input GCR spectra? (multi source GPS)"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]