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

None 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)
Re: Idea Re: simplest example of reverse monte carlo for GRAS (Leonardo Ghizoni)
Date: May 11, 10:28
From: Svetlana Shasharina <Svetlana Shasharina>

Thank you, Leonardo!  Yes, I figured this out.  I appreciate your help.

Sveta


On 5/11/18 11:22 AM, Leonardo Ghizoni wrote:
> *** Discussion title: Space Applications
>
> 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.
>
> -------------------------------------------------------------
> Visit this GEANT4 at hypernews.slac.stanford.edu message (to reply or unsubscribe) at:
> http://hypernews.slac.stanford.edu/HyperNews/geant4/get/space_app/130/1/1/1.html
>

-- 
Svetlana Shasharina                            | sveta@txcorp.com
Vice-President of Computational Infrastructure | Tech-X Corporation

 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 ]