# run1.mac
# Run in batch: ./exampleB1 run1.mac
# Or interactively: /control/execute run1.mac

# Print commands
/control/verbose 2

#############################
##Target 1###
/score/create/cylinderMesh boxMesh_1
/score/mesh/cylinderSize 0.75 0.15 cm
#/score/mesh/translate/xyz 0.1 0.1 0.0 cm
/score/mesh/nBin 30 30 30
/score/quantity/energyDeposit eDep
/score/quantity/nOfStep nOfStepGamma
/score/filter/particle gammaFilter gamma
/score/quantity/nOfStep nOfStepEMinus
/score/filter/particle eMinusFilter e-
/score/filter/particle muMinusFilter mu-
/score/quantity/nOfStep nOfStepEPlus
/score/filter/particle ePlusFilter e+
/score/filter/particle muPlusFilter mu+
/score/close
/score/list

##Target 2###
/score/create/cylinderMesh boxMesh_2
/score/mesh/cylinderSize 0.75 0.15 cm
/score/mesh/translate/xyz 0. 0. 2.3 cm
/score/mesh/nBin 30 30 30
/score/quantity/energyDeposit eDep
/score/quantity/nOfStep nOfStepGamma
/score/filter/particle gammaFilter gamma
/score/quantity/nOfStep nOfStepEMinus
/score/filter/particle eMinusFilter e-
/score/filter/particle muMinusFilter mu-
/score/quantity/nOfStep nOfStepEPlus
/score/filter/particle ePlusFilter e+
/score/filter/particle muPlusFilter mu+
/score/close
/score/list

# positron 45 GeV to the direction (0.,0.,1.)
#/gun/particle e+
#/gun/energy 45 GeV
/gps/ene/type Mono
/gps/particle e+
/gps/energy 45 GeV
/gps/pos/radius 0 mm
/gps/direction 0 0 1
/gps/position 0 0 -100 cm 

#/run/beamOn 10000000
#/run/beamOn 100000
/run/beamOn 1000

/score/dump
/score/dumpAllQuantitiesToFile boxMesh_1 boxMesh_1.dump
/score/dumpAllQuantitiesToFile boxMesh_2 boxMesh_2.dump

# 
# proton 210 MeV to the direction (0.,0.,1.)
#
#/gun/particle proton
#/gun/energy 210 MeV
#/tracking/verbose 2
#
#/run/beamOn 1