# Write phase space in TOPAS ASCII form
b:Ge/World/Invisible = "TRUE"
s:Ge/VacFilm/Type = "TsBox"
s:Ge/VacFilm/Parent = "World"
s:Ge/VacFilm/Material = "G4_WATER"
d:Ge/VacFilm/HLX = 50.0 cm
d:Ge/VacFilm/HLY = 50.0 cm
d:Ge/VacFilm/HLZ = 1.0 cm
d:Ge/VacFilm/TransX = 0. cm
d:Ge/VacFilm/TransY = 0. cm
d:Ge/VacFilm/TransZ = 0. cm
d:Ge/VacFilm/RotX = 0. deg
d:Ge/VacFilm/RotY = 0. deg
d:Ge/VacFilm/RotZ = 0. deg
s:Ge/VacFilm/Color = "skyblue"
s:Sc/PhaseSpaceAtVacFilm/Quantity = "PhaseSpace"
b:Sc/PhaseSpaceAtVacFilm/OutputToConsole = "True"
s:Sc/PhaseSpaceAtVacFilm/Surface = "VacFilm/ZMinusSurface"
s:Sc/PhaseSpaceAtVacFilm/OutputType = "ROOT" # ASCII, Binary, Limited or ROOT
s:Sc/PhaseSpaceAtVacFilm/OutputFile = "ROOTOutput"
i:Sc/PhaseSpaceAtVacFilm/OutputBufferSize = 1000
#s:Sc/PhaseSpaceAtVacFilm/OnlyIncludeParticlesGoing = "In"
b:Sc/PhaseSpaceAtVacFilm/IncludeTOPASTime = "True"
b:Sc/PhaseSpaceAtVacFilm/IncludeTimeOfFlight = "True"
b:Sc/PhaseSpaceAtVacFilm/IncludeRunID = "True"
b:Sc/PhaseSpaceAtVacFilm/IncludeEventID = "True"
b:Sc/PhaseSpaceAtVacFilm/IncludeTrackID = "True"
b:Sc/PhaseSpaceAtVacFilm/IncludeParentID = "True"
b:Sc/PhaseSpaceAtVacFilm/IncludeVertexInfo = "True"
b:Sc/PhaseSpaceAtVacFilm/IncludeSeed = "True"
#sv:Sc/PhaseSpaceAtVacFilm/OnlyIncludeParticlesNamed = 1 "Proton"
s:Sc/PhaseSpaceAtVacFilm/IfOutputFileAlreadyExists = "Overwrite"
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 = 100
i:Ts/ShowHistoryCountAtInterval = 10
b:Ts/PauseBeforeQuit = "True"