C12 ion Bragg peak

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