# Score surface track count of primary protons
# in energy bins in two phantoms

s:Ge/World/Material  = "Vacuum"
d:Ge/World/HLX       = 2.0 m
d:Ge/World/HLY       = 2.0 m
d:Ge/World/HLZ       = 2.0 m
b:Ge/World/Invisible = "TRUE"

s:Ge/WaterBox/Type     = "TsBox"
s:Ge/WaterBox/Parent   = "World"
s:Ge/WaterBox/Material = "G4_WATER"
d:Ge/WaterBox/HLX      = 30.0 cm
d:Ge/WaterBox/HLY      = 30.0 cm
d:Ge/WaterBox/HLZ      = 2.0 cm
d:Ge/WaterBox/TransZ   = -20. cm

s:Ge/LeadBox/Type     = "TsBox"
s:Ge/LeadBox/Parent   = "World"
s:Ge/LeadBox/Material = "Lead"
d:Ge/LeadBox/HLX      = 30.0 cm
d:Ge/LeadBox/HLY      = 30.0 cm
d:Ge/LeadBox/HLZ      = 10.0 cm
d:Ge/LeadBox/TransZ   = -40. cm
d:Ge/LeadBox/RotX     = 0. deg
d:Ge/LeadBox/RotY     = 0. deg
d:Ge/LeadBox/RotZ     = 0. deg

s:Sc/EnergyDepAtPhantom/Quantity                          = "SurfaceTrackCount"
s:Sc/EnergyDepAtPhantom/Surface                           = "WaterBox/ZPlusSurface"
b:Sc/EnergyDepAtPhantom/OutputToConsole                   = "True"
s:Sc/EnergyDepAtPhantom/OutputType                        = "csv"
s:Sc/EnergyDepAtPhantom/IfOutputFileAlreadyExists         = "Overwrite"
i:Sc/EnergyDepAtPhantom/EBins                             = 10
sv:Sc/EnergyDepAtPhantom/OnlyIncludeParticlesNamed        = 1 "proton"
s:Sc/EnergyDepAtPhantom/OnlyIncludeParticlesOfGeneration  = "primary"
d:Sc/EnergyDepAtPhantom/EBinMin                           = 0. MeV
d:Sc/EnergyDepAtPhantom/EBinMax                           = 100. MeV

s:Sc/EnergyDepAtPhantom2/Quantity                         = "SurfaceTrackCount"
s:Sc/EnergyDepAtPhantom2/Surface                          = "LeadBox/ZPlusSurface"
b:Sc/EnergyDepAtPhantom2/OutputToConsole                  = "True"
s:Sc/EnergyDepAtPhantom2/OutputType                       = "csv"
s:Sc/EnergyDepAtPhantom2/IfOutputFileAlreadyExists        = "Overwrite"
i:Sc/EnergyDepAtPhantom2/EBins                            = 10
sv:Sc/EnergyDepAtPhantom2/OnlyIncludeParticlesNamed       = 1 "proton"
s:Sc/EnergyDepAtPhantom2/OnlyIncludeParticlesOfGeneration = "primary"
d:Sc/EnergyDepAtPhantom2/EBinMin                          = 0. MeV
d:Sc/EnergyDepAtPhantom2/EBinMax                          = 100. MeV

s:Gr/ViewA/Type                             = "OpenGL"
i:Gr/ViewA/WindowSizeX                      = 900
i:Gr/ViewA/WindowSizeY                      = 900
d:Gr/ViewA/Theta                            = 55 deg
d:Gr/ViewA/Phi                              = 20 deg
s:Gr/ViewA/Projection                       = "Perspective"
d:Gr/ViewA/PerspectiveAngle                 = 30 deg
u:Gr/ViewA/Zoom                             = 1.3
b:Gr/ViewA/HiddenLineRemovalForTrajectories = "True"

s:So/Example/Type                     = "Beam"
s:So/Example/Component                = "BeamPosition"
s:So/Example/BeamParticle             = "proton"
d:So/Example/BeamEnergy               = 169.23 MeV
u:So/Example/BeamEnergySpread         = 0.757504
s:So/Example/BeamPositionDistribution = "Gaussian"
s:So/Example/BeamPositionCutoffShape  = "Ellipse"
d:So/Example/BeamPositionCutoffX      = 10. cm
d:So/Example/BeamPositionCutoffY      = 10. cm
d:So/Example/BeamPositionSpreadX      = 0.65 cm
d:So/Example/BeamPositionSpreadY      = 0.65 cm
s:So/Example/BeamAngularDistribution  = "Gaussian"
d:So/Example/BeamAngularCutoffX       = 90. deg
d:So/Example/BeamAngularCutoffY       = 90. deg
d:So/Example/BeamAngularSpreadX       = 0.0032 rad
d:So/Example/BeamAngularSpreadY       = 0.0032 rad
i:So/Example/NumberOfHistoriesInRun   = 1000

i:Ts/ShowHistoryCountAtInterval = 10
b:Ts/PauseBeforeQuit = "True"
