lanlsourcerun_2.mac

/run/verbose 1
/run/printProgress 1000
#/event/verbose 1
#/tracking/verbose 1

/g4simple/setRandomSeed true

/g4simple/setReferencePhysList Shielding_LIV 

##Geometry File to run simulation with
/g4simple/setDetectorGDML gscanner.gdml

##Specify output file
/g4simple/setOutputFormat hdf5
/analysis/setFileName Sim_2

##Set Volume ID's for Geometry Volumes (volume.in.geometry.file_PV volID)
/g4simple/setVolID PPC_Germanium_Detector_PV 1
#/g4simple/setVolID Gold_Foil_PV 2
#/g4simple/setVolID LANL_Gadget_PV 4
#/g4simple/setVolID vacuum_volume_PV 3

##to track all steps, including steps for which Edep=0
#/g4simple/recordAllSteps

##initialize the run
/process/em/fluoBearden true
/process/deex/correlatedGamma true
/run/initialize
/process/em/auger true
/process/em/fluo true
/process/em/pixe true
/process/em/deexcitationIgnoreCut true 
/run/setCutForAGivenParticle e- 100 nm
/run/setCutForAGivenParticle gamma 100 nm
/run/setCutForAGivenParticle e+ 100 nm
/process/em/printParameters

/g4simple/listPhysVols

##Set Source Specifications; particle type, ion specifications (proton number, nucleon number, ionic charge, excitation energy in keV),  geometrical specifications (3 lines, comment these out to make point source), position of source, what volume in geometry file to confine the source to (will take overlap of source geometry and volume in geometry file to make new source geometry, comment this out to not confine source to any volume in geometry file).
#/gps/particle ion
#/gps/ion 95 241 0 0
#/grdm/nucleusLimits 241 241 95 95
#/gps/energy 0 MeV
#/gps/pos/type Volume
#/gps/pos/shape Sphere
#/gps/pos/radius 5 mm
#/gps/pos/centre 0 0 59.139 mm
#/gps/pos/confine LANL_Source_Disk_PV

#/gps/particle gamma
#/gps/energy 60 keV
#/gps/direction 0 0 -1
#/gps/pos/centre 0 0 41 mm

/gps/particle ion
/gps/ion 43 99 0 0
/gps/pos/type Volume
/gps/pos/shape Cylinder
/gps/pos/centre 14.5 0 50.501 mm
/gps/pos/radius 3.125 mm
/gps/pos/halfz 0.0005 mm

##set number of runs
/run/beamOn 1000000