/tracking/verbose 0
/run/verbose 0

#/Physics/addPhysics QGSP_BIC_EMY
#/Physics/setCuts 0.0001 mm

#
# Tell the messenger that we want to construct the iThemba model
#
/Beamline/iThemba/Construct  
#
# remove old beamline elements
#
/Beamline/iThemba/RemoveMultiLayerFaradayCup                       Yes          # Yes (default)/ No 
/Beamline/iThemba/RemoveRangeMonitor                               Yes          # Yes (default)/ No 
/Beamline/iThemba/RemoveOccludingRingAndCentralStopper             Yes          # Yes (default)/ No 
#
# new beam line elements
#
/Beamline/iThemba/RemoveVacuumTube                                 No           # Yes / No (default)
/Beamline/iThemba/RemoveReferenceIonizationChamber                 No           # Yes / No (default)
/Beamline/iThemba/RemoveYSteeringMagnets                           No           # Yes / No (default)
/Beamline/iThemba/RemoveMultiWireIonizationChamber                 No           # Yes / No (default)
/Beamline/iThemba/RemoveXSteeringMagnets                           No           # Yes / No (default)
/Beamline/iThemba/RemoveRangeTrimmerPlate                          No           # Yes / No (default)
/Beamline/iThemba/RemoveShieldingCollimators                       No           # Yes / No (default)
/Beamline/iThemba/RemoveQuadrantAndMonitorIonizationChamber        No           # Yes / No (default)
/Beamline/iThemba/RemoveBeamFinalCollimator                        No           # Yes / No (default)
/Beamline/iThemba/RemoveNewRangeMonitoriThemba                     No           # Yes / No (default)
/Beamline/iThemba/RemoveLeadPlate                                  No           # Yes / No (default)
/Beamline/iThemba/RemoveWedge                                      Yes          # Yes (default)/ No 
/Beamline/iThemba/RemoveRangeModulator                             Yes          # Yes(default) / No 
#
# Set the range of the protons and the collimator diameter   
/Beamline/iThemba/SetRange 60.0 mm                     
/Beamline/iThemba/TreatmentCollimator/Diameter 10 cm
#
# set the strip tracker parameters
#
/StripTracker/Construct              false

#
# Set parameters for the phantom
#
/Phantom/Compensator/Construct       false
/Phantom/Construct                   false


#
# Paramters for single layer RT
#
/RangeTelescope/setNoModuleRT        1
/RangeTelescope/setEpiThickness      14 um
/RangeTelescope/setFrontFacePosition 0 0 20 cm
/RangeTelescope/Construct            true

# these do not work on multithreaded for some reason
# need to fix before we can use these
#/Pixel/SetSensorLength               100 mm
#/Pixel/SetPixelPitch                 50 um
#/Pixel/SetThreshold                  950              # number of electrons #

/DoubleDynamite/Construct            false
#
# Initialize Geometry
#
/Beamline/Update


/run/initialize

/gps/particle proton
/gps/pos/type Beam

/gps/pos/shape Ellipse
/gps/pos/halfx 0.721249 mm
/gps/pos/halfy 0.5863502 mm

/gps/pos/centre 0. 0. -7278.0 mm

/gps/ene/type Gauss
/gps/energy 200.64 MeV
/gps/ene/sigma 2.0 MeV        
/gps/direction 0 0 1          # the beam is travelling along Z axis

/histo/setFileName singleDynamite_iThemba2013_30mmRange

#/vis/open OGLI
##/vis/open DAWNFILE
##/vis/viewer/reset
#/vis/viewer/set/viewpointThetaPhi 80 20 deg
#/vis/viewer/set/style w
#/vis/scene/create
##/vis/scene/add/axes
##/vis/scene/add/volume
##/vis/scene/add/trajectories smooth
#/vis/sceneHandler/attach
#/vis/scene/endOfEventAction accumulate 10
#/vis/viewer/set/autoRefresh true

/run/initialize
/run/beamOn 100000
#