#
#    **********************************
#    *                                *
#    *   Macro 1.0MeV   *
#    *                                *
#    **********************************
#


#---------------------------------------------------------------
# Choosing the physics model 
#---------------------------------------------------------------

/UseEM/phys/addPhysics emlivermore


#---------------------------------------------------------------
# Confirm cutoff
#---------------------------------------------------------------

/cuts/setLowEdge 1.0* keV

#---------------------------------------------------------------
# Initialize kernel
#---------------------------------------------------------------

/run/initialize

#---------------------------------------------------------------
# Control the amount of information printed.
#---------------------------------------------------------------

/control/verbose  2
/tracking/verbose 2
/run/verbose      2
/event/verbose    2
/gps/verbose      2

#-------------------------------------------------------------------------------------------
# Define the particle source
#-------------------------------------------------------------------------------------------
# 1. Define the particle type
# 2. Define the source distribution type 
# 3. Define the source shape type
# 4. Define the center coordinates of the source 
# 5. Define the half length in x
# 6. Define the half length in y (infinite)
# 7. Define the half length in z (infinite)
# 8./9./10. Define Parallelepiped-forming angles
# 11. Set the angular distribution type (here isotropic) 
# 12. Set the source energy spectrum (here monoenergetic)
# 13. Set the gamma energy for mono-energetic sources
# 14. Set the number of particles emitted per event (events are selected by /run/beamOn)
#-------------------------------------------------------------------------------------------
/gps/particle   e-
/gps/pos/type   Plane  #Beam Surface Plane Point or Volume 
/gps/pos/shape  Rectangle
/gps/pos/centre -.1400 0.0 0.0 cm
/gps/pos/halfx  0.0 cm
/gps/pos/halfy  .1400 cm
/gps/pos/halfz  .1400 cm
/gps/pos/paralp 0.0   rad
/gps/pos/parthe 0.0   rad
/gps/pos/parphi 0.0   rad
/gps/ene/type  Mono
/gps/direction  1 0 0
/gps/ene/mono   1.0 MeV
/gps/number     1  
/gps/source/show

#-------------------------------------------------------------------------------
# Define scoring mesh
#-------------------------------------------------------------------------------
/score/create/boxMesh BoxMeshGlobal
/score/mesh/boxSize .1400 2500.00 2500.00 cm 
/score/mesh/nBin    200    1      1     
/score/quantity/energyDeposit eDep J
/score/quantity/doseDeposit   dDep Gy
/score/quantity/cellFlux      flux 
/score/quantity/cellCharge    ChargeDep
/score/list
/score/close


#-------------------------------------------------------------------------------
# Control on secondary phenomenons
#-------------------------------------------------------------------------------
/process/em/deexcitation region true true true 
/process/em/fluo            true true true true true 
/process/em/auger           true true true true true
#/process/em/pixe            true    //Particle-induced X-ray emission
#/process/em/pixeXSmodel     Empirical
#/process/em/pixeElecXSmodel Empirical

#-------------------------------------------
# Auger Cascade and Cutoff
#-------------------------------------------
/process/em/augerCascade true
/process/em/deexcitationIgnoreCut false

#---------------------------------------------------------------------------------
# Running Events
#---------------------------------------------------------------------------------
/run/beamOn 15000000

#---------------------------------------------------------------------------------
# Scoring Saving
#---------------------------------------------------------------------------------
/score/dumpQuantityToFile BoxMeshGlobal eDep       E2_inf_mobetron_uni_1.0MeV_200mesh_emlivermore_energy_15M.txt
/score/dumpQuantityToFile BoxMeshGlobal dDep       E2_inf_mobetron_uni_1.0MeV_200mesh_emlivermore_dose_15M.txt
/score/dumpQuantityToFile BoxMeshGlobal flux       E2_inf_mobetron_uni_1.0MeV_200mesh_emlivermore_flux_15M.txt
/score/dumpQuantityToFile BoxMeshGlobal ChargeDep  E2_inf_mobetron_uni_1.0MeV_200mesh_emlivermore_ChargeDep_15M.txt


#-----------------------------
# Visualization
#-----------------------------
#/vis/ogl/printEPS
#/vis/ogl/set/printMode pixmap
#/vis/ogl/set/printSize 3000 3000
#/vis/ogl/set/printFilename benchmark_120_deg
#/vis/ogl/export