# Demonstrates Environment Source
# The environment is a radiation field that might be experienced, for
# example, by a spacecraft in a radiation belt, or a robot (or, indeed
# a human) in a damaged nuclear reactor.

# A notional radiation "cavity" is created enclosing all volumes in a
# component. The cavity is a sphere of radius R.  Primary particles are
# generated on the surface of the sphere, directed inwards, following a
# cosine angular distribution (Lambert's cosine law) relative to the inward
# direction. This produces an isotropic, homogeneous, "omnidirectional flux".

# Even if the radiation has some directional dependency it is often the case
# that the instrument (your detector) is rotating or moving about so the
# flux will average to isotropic over time.

# The basic definition of flux, f, which in principle can vary with
# direction and position, is defined by dN/dt = f*da*dOmega, the
# rate of flow of particles out of an element of area da perpendicular to
# the direction into an element of solid angle dOmega. If the flux is
# homogenous and isotropic, we can define the "omnidirectional flux"
# F = 4*pi*f per unit area. It is quoted, for example, as number per cm2
# per second.

# Fluence is simply F*T, the flux F over a time period T, so quoted, for
# example, as number per cm2.

# One can derive equivalent definitions of fluence:
# - the number of particles that enter a sphere of unit cross-sectional
#   area;
# - the track length per unit volume;
# - etc. For more details see TOPAS documentation.

# For N particles (histories), the fluence will be N/(pi*R^2). This is
# printed at the end of run. It is up to you to decide if this is enough
# for your application. Thus:
# - to simulate flux F for time T you need pi*R^2*F*T histories;
# - or, given N histories, you will have simulated a time period
#   T = N/(F*pi*R^2).

# A test sphere of radius r will attract N*r^2/R^2 particles.

# A thin test disc of radius r will attract (N/2)*r^2/R^2 particles.

# This parameter file implements:
# - a test sphere
# - a thin test disc
# - a test cube
# In all cases:
# - the fluence is "measured" by a "Fluence" scorer (scores the sum of all
#   track lengths divides by the volume);
# - the number of tracks crossing a surface is scored, and can be related
#   to the fluence by the above formulae.

# All materials are set to vacuum so that flux is not attenuated
# within the cavity.

# Note: the world must be bigger than the radiation cavity, which may be
# bigger than a box enclosing your geometry. TOPAS will tell you if you
# need to increase the world size.
s:Ge/World/Material = "Vacuum"
s:Ge/World/Type     = "TsBox"
d:Ge/World/HLX      = 5 m
d:Ge/World/HLY      = 5 m
d:Ge/World/HLZ      = 5 m

s:Ge/TestVolumes/Type    = "Group"
s:Ge/TestVolumes/Parent  = "World"
dc:Ge/TestVolumes/RotZ   = 0. deg
dc:Ge/TestVolumes/TransZ = 1 m

s:Ge/TestSphere/Type      = "TsSphere"
s:Ge/TestSphere/Parent    = "TestVolumes"
sc:Ge/TestSphere/Material = "Vacuum"
dc:Ge/TestSphere/RMax     = 50 cm
dc:Ge/TestSphere/TransY   = -2 m
dc:Ge/TestSphere/TransZ   = 0.5 m

s:Ge/TestBox/Type      = "TsBox"
s:Ge/TestBox/Parent    = "TestVolumes"
sc:Ge/TestBox/Material = "Vacuum"
dc:Ge/TestBox/HLX      = 0.5 m
dc:Ge/TestBox/HLY      = 0.5 m
dc:Ge/TestBox/HLZ      = 0.5 m
dc:Ge/TestBox/TransY   = -2 m
dc:Ge/TestBox/TransZ   = -1 m

s:Ge/TestDisc/Type     = "TsCylinder"
s:Ge/TestDisc/Parent   = "TestVolumes"
sc:Ge/TestDisc/Material = "Vacuum"
dc:Ge/TestDisc/RMax     = 50 cm
dc:Ge/TestDisc/HL       = 1 nm
dc:Ge/TestDisc/TransY   = 2 m
dc:Ge/TestDisc/TransZ   = 0.5 m

s:Sc/TestSphereFluenceScorer/Quantity = "Fluence"
s:Sc/TestSphereFluenceScorer/Component = "TestSphere"
b:Sc/TestSphereFluenceScorer/PropagateToChildren = "True"
b:Sc/TestSphereFluenceScorer/OutputToConsole = "True"
b:Sc/TestSphereFluenceScorer/OutputAfterRun = "True"
s:Sc/TestSphereFluenceScorer/OutputType = "csv"
s:Sc/TestSphereFluenceScorer/IfOutputFileAlreadyExists = "Overwrite"

s:Sc/TestSphereTrackScorer/Quantity = "SurfaceTrackCount"
s:Sc/TestSphereTrackScorer/Component = "TestSphere"
s:Sc/TestSphereTrackScorer/Surface = "TestSphere/OuterCurvedSurface"
b:Sc/TestSphereTrackScorer/PropagateToChildren = "True"
b:Sc/TestSphereTrackScorer/OutputToConsole = "True"
b:Sc/TestSphereTrackScorer/OutputAfterRun = "True"
s:Sc/TestSphereTrackScorer/OutputType = "csv"
s:Sc/TestSphereTrackScorer/IfOutputFileAlreadyExists = "Overwrite"

s:Sc/TestBoxFluenceScorer/Quantity = "Fluence"
s:Sc/TestBoxFluenceScorer/Component = "TestBox"
b:Sc/TestBoxFluenceScorer/PropagateToChildren = "True"
b:Sc/TestBoxFluenceScorer/OutputToConsole = "True"
b:Sc/TestBoxFluenceScorer/OutputAfterRun = "True"
s:Sc/TestBoxFluenceScorer/OutputType = "csv"
s:Sc/TestBoxFluenceScorer/IfOutputFileAlreadyExists = "Overwrite"

s:Sc/TestBoxTrackScorer/Quantity = "SurfaceTrackCount"
s:Sc/TestBoxTrackScorer/Component = "TestBox"
s:Sc/TestBoxTrackScorer/Surface = "TestBox/XPlusSurface"
b:Sc/TestBoxTrackScorer/PropagateToChildren = "True"
b:Sc/TestBoxTrackScorer/OutputToConsole = "True"
b:Sc/TestBoxTrackScorer/OutputAfterRun = "True"
s:Sc/TestBoxTrackScorer/OutputType = "csv"
s:Sc/TestBoxTrackScorer/IfOutputFileAlreadyExists = "Overwrite"

s:Sc/TestDiscFluenceScorer/Quantity = "Fluence"
s:Sc/TestDiscFluenceScorer/Component = "TestDisc"
b:Sc/TestDiscFluenceScorer/PropagateToChildren = "True"
b:Sc/TestDiscFluenceScorer/OutputToConsole = "True"
b:Sc/TestDiscFluenceScorer/OutputAfterRun = "True"
s:Sc/TestDiscFluenceScorer/OutputType = "csv"
s:Sc/TestDiscFluenceScorer/IfOutputFileAlreadyExists = "Overwrite"

s:Sc/TestDiscTrackScorer/Quantity = "SurfaceTrackCount"
s:Sc/TestDiscTrackScorer/Surface = "TestDisc/ZMinusSurface"
s:Sc/TestDiscTrackScorer/Component = "TestDisc"
b:Sc/TestDiscTrackScorer/PropagateToChildren = "True"
b:Sc/TestDiscTrackScorer/OutputToConsole = "True"
b:Sc/TestDiscTrackScorer/OutputAfterRun = "True"
s:Sc/TestDiscTrackScorer/OutputType = "csv"
s:Sc/TestDiscTrackScorer/IfOutputFileAlreadyExists = "Overwrite"

s:So/Example/Type                     = "Environment"
s:So/Example/Component                = "TestVolumes"
s:So/Example/BeamParticle             = "proton"
dc:So/Example/BeamEnergy              = 200 MeV
ic:So/Example/NumberOfHistoriesInRun  = 1000000

#i:Ts/Seed = 1
i:Ts/NumberOfThreads = 0  # All available threads
i:Ts/ShowHistoryCountAtInterval = 100000

#s:Gr/ViewA/Type              = "OpenGL"
#b:Gr/ViewA/IncludeStepPoints = "True"
#b:Ts/PauseBeforeQuit = "True"
#b:Ts/UseQt = "True"