Dear Geant4 experts,
In the extended TestEm7 example in Geant4 I get Bragg curves for protons, but not with carbon ions. I will be grateful if you could help me.
Best regards
ionC12.mac
Macro file for “TestEm7.cc”
(can be run in batch, without graphic)
/control/verbose 2
/run/verbose 2
/tracking/verbose 0
/testem/det/setMat G4_WATER
/testem/det/setSizeX 20 cm
/testem/det/setSizeYZ 20 cm
/testem/det/tallyNumber 4
/testem/det/tallyDefinition 0 2. 2. 2. mm
/testem/det/tallyDefinition 1 2. 2. 2. mm
/testem/det/tallyDefinition 2 2. 2. 2. mm
/testem/det/tallyDefinition 3 2. 2. 2. mm
/testem/det/tallyPosition 0 -9.9 0. 0. cm
/testem/det/tallyPosition 1 -5. 0. 0. cm
/testem/det/tallyPosition 2 0. 0. 0. cm
/testem/det/tallyPosition 3 5. 0. 0. cm
/testem/phys/addPhysics emstandard_opt3
/testem/phys/addPhysics binary_ion
/testem/phys/addPhysics elastic
#/testem/phys/addPhysics ionGasModels
/run/initialize
/gun/particle ion
/gun/ion 6 12 6
/gun/energy 100 MeV
step limit
/run/setCut 0.5 mm
/run/printProgress 100
/analysis/setFileName c12
/analysis/h1/set 1 100 0 100 mm
/analysis/h1/set 2 100 0 100 mm
/analysis/h1/set 3 100 0 100 mm
/run/beamOn 1000
proton.mac
Macro file for “TestEm7.cc”
(can be run in batch, without graphic)
/control/verbose 2
/run/verbose 2
/testem/det/setMat Water
/testem/det/setSizeX 10 cm
/testem/det/setSizeYZ 10 cm
#/testem/phys/addPhysics local
/testem/phys/addPhysics emstandard_opt3
/testem/phys/addPhysics binary
/testem/phys/addPhysics elastic
/run/initialize
/gun/particle proton
/gun/energy 100 MeV
step limit
#/run/setCut 0.5 mm
/run/setCut 1 mm
#/testem/stepMax 0.5 mm
/run/printProgress 1000
/analysis/setFileName testem7
/analysis/h1/set 2 100 0 100 mm
/analysis/h1/set 3 100 0 100 mm
/run/beamOn 10000