#set number of threads
/run/numberOfThreads 4

#set Verbosity
#/control/verbose 2
#/tracking/verbose 2

#Initialize the kernel
/run/initialize

#Primary particles
#source shape
#/gps/pos/type Volume		
#/gps/pos/shape Cylinder
#/gps/pos/radius 0.74 cm
#/gps/pos/halfz 1.6 cm
#position
/gps/pos/centre 0. 0. -26.6 cm	    
#volume confinement
/gps/pos/confine VialSolventPhysical #confinement must be after the source definition
#direction
/gps/ang/type iso      				#isotropic emission
#particle
/gps/particle ion					#Am241->(95, 241), Co57->(27, 57), Cs137->(55, 137), Na22->(11, 22)	
/gps/ion 27 57 0 0             		#Z, A, Q, E [checked->it works! (expect Am241->more lines than expected)]
#kinetic energy
/gps/ene/type Mono
/gps/ene/mono 0. keV            	#particle at rest

#/gps/particle gamma
#/gps/ene/mono 59.54 keV

#/control/execute mymacro.mac

#set output name (without extention)
/run/setfilenamesave output

# Start the simulation
/run/printProgress 10000
/run/beamOn 100000