The dnaphysics example shows how to simulate track structures in liquid water using the Geant4-DNA physics models and processes.
The Geant4-DNA processes and models are further described at: http://geant4-dna.org
Any report or published results obtained using the Geant4-DNA software shall cite the following Geant4-DNA collaboration publications: Med. Phys. 45, (2018) e722-e739 Phys. Med. 31 (2015) 861-874 Med. Phys. 37 (2010) 4692-4708 Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178
The geometry is a 100-micron side cube (World) made of liquid water (G4_WATER material). Particles are shot from the center of the volume.
The variable density feature of materials is illustrated in DetectorConstruction. The material can be changed directly in dnaphysics.in macro file.
Make sure $G4LEDATA points to the low energy electromagnetic data files.
In interactive mode, run:
./dnaphysics
In batch, the macro dnaphysics.in can be used. It shows how to shoot different particle types.
The PhysicsList uses Geant4-DNA Physics constructors or activator. Activator can be used to apply Geant4-DNA to a specific Region only.
1) Geant4-DNA Physics constructors can be selected using the command:
/dna/test/addPhysics DNA_OptX
where X is 0 to 8 (2, 4 or 6 are recommended).
2) Alternatively, Geant4-DNA activator for World Region can be selected using the command:
/process/em/AddDNARegion World DNA_OptX
where X = 0, 2, 4, or 6.
3) In addition to 1) or 2), to enable radioactive decay, one can use:
/dna/test/addPhysics raddecay
4) Warning regarding ions: when the incident particle type is ion (/gun/particle ion), specified with Z and A numbers (/gun/ion A Z), the Rudd ionisation extended model is used. The particles are tracked by default down to 0.5 MeV/u. This tracking cut can be bypassed using :
/dna/test/addIonsTrackingCut false
The output results consists in a dna.root file, containing for each simulation step:
This information is extracted from the SteppingAction class.
The ROOT file can be easily analyzed using for example the provided ROOT macro file plot.C; to do so : be sure to have ROOT installed on your machine be sure to be in the directory containing the ROOT files created by dnaphysics copy plot.C into this directory from there, launch ROOT by typing root under your ROOT session, type in : .X plot.C to execute the macro file alternatively you can type directly under your session : root plot.C
The naming scheme on the displayed ROOT plots is as follows (see SteppingAction.cc):
e-_G4DNAElectronSolvation: 10
e-_G4DNAElastic: 11
e-_G4DNAExcitation: 12
e-_G4DNAIonisation: 13
e-_G4DNAAttachment: 14
e-_G4DNAVibExcitation: 15
msc: 110
CoulombScat: 120
eIoni: 130
proton_G4DNAElastic: 21
proton_G4DNAExcitation: 22
proton_G4DNAIonisation: 23
proton_G4DNAChargeDecrease: 24
msc: 210
CoulombScat: 220
hIoni: 230
nuclearStopping: 240
hydrogen_G4DNAElastic: 31
hydrogen_G4DNAExcitation: 32
hydrogen_G4DNAIonisation: 33
hydrogen_G4DNAChargeIncrease: 35
alpha_G4DNAElastic: 41
alpha_G4DNAExcitation: 42
alpha_G4DNAIonisation: 43
alpha_G4DNAChargeDecrease: 44
msc: 410
CoulombScat: 420
ionIoni: 430
nuclearStopping: 440
alpha+_G4DNAElastic: 51
alpha+_G4DNAExcitation: 52
alpha+_G4DNAIonisation: 53
alpha+_G4DNAChargeDecrease: 54
alpha+_G4DNAChargeIncrease: 55
msc: 510
CoulombScat: 520
hIoni: 530
nuclearStopping: 540
helium_G4DNAElastic: 61
helium_G4DNAExcitation: 62
helium_G4DNAIonisation: 63
helium_G4DNAChargeIncrease: 65
GenericIon_G4DNAIonisation: 73
msc: 710
CoulombScat: 720
ionIoni: 730
nuclearStopping: 740
phot: 81
compt: 82
conv: 83
Rayl: 84
Should you have any enquiry, please do not hesitate to contact: incer.nosp@m.ti@c.nosp@m.enbg..nosp@m.in2p.nosp@m.3.fr or tran@.nosp@m.lp2i.nosp@m.b.in2.nosp@m.p3.f.nosp@m.r