Gabriele Cosmo
60 lines
1.5 KiB
Plaintext
60 lines
1.5 KiB
Plaintext
# S.E. Mazzaglia Nov 14, 2009
|
|
#
|
|
# Example macro file showing how to calculate stopping powers for particles and ions
|
|
#
|
|
|
|
#######################
|
|
# Set of the verboses #
|
|
#######################
|
|
#
|
|
/control/verbose 1
|
|
/tracking/verbose 0
|
|
/run/verbose 1
|
|
/event/verbose 0
|
|
|
|
##########################
|
|
# Set of the physic models
|
|
#
|
|
/Physics/addPhysics HADRONTHERAPY_1
|
|
|
|
# Initialize geometry and physic
|
|
/run/initialize
|
|
|
|
# Set a very high time threshold to allow all decays to happen
|
|
/process/had/rdm/thresholdForVeryLongDecayTime 1.0e+60 year
|
|
|
|
###############################
|
|
# Stopping Powers calculation #
|
|
###############################
|
|
|
|
###############################
|
|
# List of available materials #
|
|
###############################
|
|
|
|
# Show all available materials (list with additional informations)
|
|
/parameter/nist all
|
|
|
|
# Show all available elements (list with additional informations)
|
|
#/parameter/nist simple
|
|
|
|
# Show all available elements and materials (simple full list)
|
|
#/parameter/nist
|
|
|
|
###############
|
|
# "Particles" #
|
|
###############
|
|
|
|
# Calculate stopping power for proton into liquid water at 1 MeV
|
|
/parameter/getstopping G4_WATER
|
|
|
|
# Calculate stopping powers for protons into Hydrogen from 1 keV up to 200 MeV with 100 data points
|
|
# Print data to user terminal
|
|
/parameter/getstopping G4_H 0.001 200 100
|
|
|
|
# Calculate stopping powers for alphas into Copper from 1 keV up to 200 MeV with 100 data points
|
|
/parameter/getstopping G4_Cu 0.001 200 100 alpha stopAlpha_Cu_1keV_200MeV.out
|
|
|
|
|
|
|
|
|