Message: Re: simplest example of reverse monte carlo for GRAS Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Idea Re: simplest example of reverse monte carlo for GRAS 

Forum: Space Applications
Re: None simplest example of reverse monte carlo for GRAS (Svetlana Shasharina)
Re: None Re: simplest example of reverse monte carlo for GRAS (Tevfik Kaplanoglu)
Re: None Re: simplest example of reverse monte carlo for GRAS (Svetlana Shasharina)
Date: May 11, 10:22
From: Leonardo Ghizoni <Leonardo Ghizoni>

Hi,

I'm not sure if you already figured this out, but here's an example of one of my past simulations using GRAS RMC. It's a trapped protons file of a SSO 800km orbit (the levels were taken from SPENVIS):

/gps/ene/type Arb
/gps/hist/type arb
/gps/hist/point   1.000000E-01   4.949600E+04
/gps/hist/point   1.500000E-01   3.489700E+04
/gps/hist/point   2.000000E-01   2.289800E+04
/gps/hist/point   3.000000E-01   1.044100E+04
/gps/hist/point   4.000000E-01   5.505700E+03
/gps/hist/point   5.000000E-01   2.978600E+03
/gps/hist/point   6.000000E-01   1.823800E+03
/gps/hist/point   7.000000E-01   1.159600E+03
/gps/hist/point   1.000000E+00   5.049900E+02
/gps/hist/point   1.500000E+00   1.599600E+02
/gps/hist/point   2.000000E+00   8.276700E+01
/gps/hist/point   3.000000E+00   3.253700E+01
/gps/hist/point   4.000000E+00   1.883800E+01
/gps/hist/point   5.000000E+00   1.211600E+01
/gps/hist/point   6.000000E+00   9.357900E+00
/gps/hist/point   7.000000E+00   7.460900E+00
/gps/hist/point   1.000000E+01   4.596300E+00
/gps/hist/point   1.500000E+01   2.429900E+00
/gps/hist/point   2.000000E+01   1.706700E+00
/gps/hist/point   3.000000E+01   1.148200E+00
/gps/hist/point   4.000000E+01   9.444400E-01
/gps/hist/point   5.000000E+01   8.249200E-01
/gps/hist/point   6.000000E+01   7.412400E-01
/gps/hist/point   7.000000E+01   6.907000E-01
/gps/hist/point   1.000000E+02   5.362600E-01
/gps/hist/point   1.500000E+02   3.412200E-01
/gps/hist/point   2.000000E+02   2.116100E-01
/gps/hist/point   3.000000E+02   8.763900E-02
/gps/hist/point   4.000000E+02   7.265500E-03

/gras/analysis/adjoint/addUserDefinedPrimSpectrum proton 1.091733E+12 1/cm2 Arb
#
#Definition of the adjoint source
##################################
/adjoint/DefineAdjSourceOnExtSurfaceOfAVolume {TargetName}
/adjoint/SetAdjSourceEmin 1.0000e-03 MeV
/adjoint/SetAdjSourceEmax 4.0000e+02 MeV
#
#Definition of the external source
###################################
/adjoint/DefineSphericalExtSource 0 0 0 20 cm            #<X Y Z R length_unit>
/adjoint/SetExtSourceEmax  4.0000e+02 MeV

Observe that I set the SetAdjSourceEmin and SetAdjSourceEmax as the minimum and maximum energies provided by the SPENVIS histogram.

Hope it helps.

Regards,
Leonardo.

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: simplest example of reverse monte carlo for GRAS   (Svetlana G. Shasharina - May 11, 10:28)
 Add Message Add Message
to: "Re: simplest example of reverse monte carlo for GRAS"

 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 ]