# ---------------------------------------------------------
# THIS MACRO PERMIT TO RECONSTRUCT THE SPREAD OUT BRAGG PEAK
# MODULATING THE PROTON BEAM ENERGY BY MEANS OF A ROTATING
# MODULATOR WHEEL.
# HERE THE WHEEL IS ROTATED OF 1 DEGREE AT EACH RUN AND 
# ONE THOUSAND PROTON HYSTORIES ARE SIMULATED IN EACH 
# RUN
# IN ORDER TO INCREASE THE RESOLUTION MORE EVENT CAN ADDED
# OR A SMALLER ANGLE CAN BE CHOSEN.
# 
# NOTE THAT THE MODULATOR MATERIAL IS POLTMETHYLMETHACRILATE 
# (PMMA) FOR DEFAULT. IF ONE WANT CARRY OUT A SIMULATION WITHOUT 
# THE MODULATOR HE/SHE MUST SET "Air" the <<ModMater>> in the
# <<GetMater>> function of the HadrontherapyModulator.cc class
#
# USERS SHOULD GIVE A LOOK TO THE HELP OF THE IDLE TO KNOW
# THE ACTIVATED MESSSENGERS FOR THE GEOMETRY
#
# ADDITIONAL INFORMATIONS ON THE MESSENGER AVAILABLE CAN BE FOUND
# INSIDE THE HADRONTHERAPY DOCUMENTATION (http://www.ge.infn.it/geant4/examples/).
# 
# ANYWAY SEND ME AN E-MAIL FOR ANY QUESTION: cirrone@lns.infn.it.   
# --------------------------------------------------------------------------------  

/control/verbose 1
/tracking/verbose 0
/run/verbose 0
/event/verbose 0


# SETTING FOR THE PHYSICS MODELS
/physics/addPhysics Decay
/physics/addPhysics EM-Photon-EPDL
/physics/addPhysics EM-Electron-EEDL
/physics/addPhysics EM-Positron-Standard
/physics/addPhysics EM-HadronIon-LowE
/physics/addPhysics EM-Muon-Standard
/physics/addPhysics HadronicEl-HadronIon-LElastic
/physics/addPhysics HadronicInel-ProtonNeutron-LEP
/physics/addPhysics HadronicInel-Ion-LEP
/physics/addPhysics HadronicInel-Pion-LEP
/physics/addPhysics HadronicAtRest-MuonMinus-Capture

/run/initialize	

# SET OF THE VISUALISATION CHARACTERISTICS
/vis/scene/create
/vis/open OGLIX
/vis/viewer/zoom 3.5
/vis/viewer/set/viewpointThetaPhi  50 130 deg
/vis/viewer/pan -60 40 cm
/tracking/storeTrajectory 1
/vis/scene/endOfEventAction accumulate

# IF YOU WANT A RANGE SHIFTER PLEASE CHANGE
# THE FOLLOWIG PARAMETERS
# SETTING FOR THE RANGE SHIFTER

/beamLine/RangeShifter/RSMat Air
/beamLine/RangeShifter/thickness 1 cm


# START OF THE MODULATOR ROTATION
# IT CONSISTS IN A COMPLETE ROTATION 
# OF THE WHHEEL (360 DEG)

/run/beamOn 10
/vis/viewer/flush

/modulator/angle 1 deg
/run/beamOn 10
/vis/viewer/flush

/modulator/angle 1 deg
/run/beamOn 10
/vis/viewer/flush