#	C O N F I G U R A T I O N

/control/execute mac/configuration.mac

/control/execute mac/verbosity.mac
#===============
#vis
#===============
/vis/disable
#/control/execute mac/visu.mac
#/vis/ogl/export Simulation.jpg

#===============
# GEOMETRY
#===============
/gate/geometry/setMaterialDatabase mac/GateMaterials.db

#===============
# GEOMETRY
#===============
# world
/gate/world/geometry/setXLength 2 m
/gate/world/geometry/setYLength 2 m
/gate/world/geometry/setZLength 2 m
/gate/world/setMaterial Vacuum

#=======================
# PHANTOM——ladder model
#=======================
/gate/world/daughters/name phantom
/gate/world/daughters/insert box
/gate/phantom/setMaterial Li
/gate/phantom/geometry/setXLength   100 mm
/gate/phantom/geometry/setYLength   100 mm
/gate/phantom/geometry/setZLength   5 mm
/gate/phantom/placement/setTranslation 0.0 0.0 2.5 mm
/gate/phantom/vis/forceWireframe
/gate/phantom/vis/setColor blue

# ============================
# ==  PHYSICS             ==
# ===========================

#/gate/physics/addPhysicsList emstandard_opt4
/gate/physics/addPhysicsList QGSP_BERT_HP

/gate/physics/Gamma/SetCutInRegion world 1 mm
/gate/physics/Electron/SetCutInRegion world 1 mm
/gate/physics/Positron/SetCutInRegion world 1 mm

/gate/physics/Gamma/SetCutInRegion phantom 1 um
/gate/physics/Electron/SetCutInRegion phantom 1 um
/gate/physics/Positron/SetCutInRegion phantom 1 um

/gate/physics/SetMinKineticEnergyInRegion phantom 0.1 eV

/gate/physics/SetMaxStepSizeInRegion phantom 1 um
/gate/physics/ActivateStepLimiter proton
/gate/physics/ActivateStepLimiter neutron
/gate/physics/ActivateStepLimiter deuteron
/gate/physics/ActivateStepLimiter triton
/gate/physics/ActivateStepLimiter alpha
/gate/physics/ActivateStepLimiter GenericIon

/gate/physics/setEMin 0.1 eV

#/gate/physics/setDEDXBinning 100
#/gate/physics/setLambdaBinning 100

/gate/physics/processList Enabled
/gate/physics/processList Initialized
/gate/physics/displayCuts

#=============================
#== ACTOR
#=============================

/gate/actor/addActor PhaseSpaceActor          gemactor
/gate/actor/gemactor/save                     output/gemactor.root
/gate/actor/gemactor/attachTo                 phantom
/gate/actor/gemactor/useVolumeFrame           false
/gate/actor/gemactor/enableEkine              true
/gate/actor/gemactor/enableXPosition          true
/gate/actor/gemactor/enableYPosition          true
/gate/actor/gemactor/enableZPosition          true
/gate/actor/gemactor/enableXDirection         true
/gate/actor/gemactor/enableYDirection         true
/gate/actor/gemactor/enableZDirection         true
/gate/actor/gemactor/enableParticleName       true
/gate/actor/gemactor/enableWeight             true
/gate/actor/gemactor/enableTime               true
#/gate/actor/gemactor/storeOutgoingParticles   true
/gate/actor/gemactor/storeSecondaries         true
#/gate/actor/gemactor/addFilter                       particleFilter
#/gate/actor/gemactor/particleFilter/addParticle      neutron

/gate/actor/addActor  DoseActor             DoseDose
/gate/actor/DoseDose/save                   output/Depth_Dose.txt
/gate/actor/DoseDose/attachTo               phantom
/gate/actor/DoseDose/stepHitType            random
/gate/actor/DoseDose/setVoxelSize           100 100 0.001 mm
/gate/actor/DoseDose/enableEdep             true

/gate/actor/addActor     SimulationStatisticActor stat 
/gate/actor/stat/save              output/stats_adder.txt 
/gate/actor/stat/saveEveryNSeconds 60

# ============================
# == INITIALIZE             ==
# ===========================

/gate/run/initialize

# ============================
# == SOURCE                 ==
# ============================

/gate/source/addSource PBS PencilBeam
/gate/source/PBS/setParticleType proton
/gate/source/PBS/setEnergy 2.5 MeV
/gate/source/PBS/setSigmaEnergy 0.0 MeV
/gate/source/PBS/setPosition 0 0 -500 mm
/gate/source/PBS/setSigmaX 2 mm
/gate/source/PBS/setSigmaY 2 mm
/gate/source/PBS/setSigmaTheta 2 mrad
/gate/source/PBS/setSigmaPhi 2 mrad
/gate/source/PBS/setEllipseXThetaEmittance 3.8 mm*mrad
/gate/source/PBS/setEllipseXThetaRotationNorm negative
/gate/source/PBS/setEllipseYPhiEmittance 3.8 mm*mrad
/gate/source/PBS/setEllipseYPhiRotationNorm negative

# ============================
# == OUTPUT                ==
# ============================

/gate/application/noGlobalOutput

#====================================================
#  R A N D O M   E N G I N E  A N D  S E E D
#====================================================
#/control/execute                   mac/SetSeed_MersenneTwister.mac
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto

#=====================================================
#   M E A S U R E M E N T   S E T T I N G S   
#=====================================================

/gate/application/setTotalNumberOfPrimaries 1000000

/gate/application/startDAQ
/gate/application/describe
exit