Message: How to model carbon ions using DNA physics? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question How to model carbon ions using DNA physics? 

Forum: DNA/Very Low Energy
Date: 24 Jul, 2018
From: Ivan Konobeev <Ivan Konobeev>

Hello

Could you please help me. This question is probably very basic however I can't find the answer.

I'm using Geant4 v.10.04.p01 with Visual Studio 2017. I'm trying to model carbon ions (A = 12, Z = 6) in dnaphysics example (which attaches default DNA constructor to World region). To do that I made minor changes to dnaphysics.in macro file: I just replaced default particle which was 100 keV proton with 2.4 MeV carbon ion. These are the changes that I made to dnaphysics.in:

From:

  /tracking/verbose 0
  /gun/particle proton
  /gun/energy 100 keV
  /run/beamOn 2

Was changed to:

  /tracking/verbose 1
  /gun/particle ion
  /gun/ion 6 12
  /gun/energy 2400 keV
  /run/beamOn 10000000

Then I executed the program with:

RelWithDebInfo\dnaphysics dnaphysics.in

Output of the program:

*********************************************************************************************************
* G4Track Information:   Particle = C12,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0       2.4        0        0         0       World initStep
    1 1.91e-08        0        0         0      2.4 1.91e-08  1.91e-08       World GenericIon_G4DNAIonisation

*********************************************************************************************************
* G4Track Information:   Particle = C12,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0       2.4        0        0         0       World initStep
    1  9.1e-08        0        0         0      2.4  9.1e-08   9.1e-08       World GenericIon_G4DNAIonisation

So carbon ions with energy 2.4 MeV are immediately consumed by water medium having track length below 1 nanometer. I believe these results are incorrect. In G4EmDNAPhysics.cc file I saw commented instructions:

  genericIonsManager->GetIon("carbon");

  if ( particleName == "carbon" ) {
        ph->RegisterProcess(new G4DNAIonisation("carbon_G4DNAIonisation"), particle);
  }

Uncommenting those did't change anything. Also it is unclear from where the particle "carbon" can be "included" since there is no "carbon" particle in "geant4_10_04_p01\source\particles" directory.

So is it possible to use Geant4-DNA constructors to model carbon ions? How can I do that?

Thank you.

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: How to model carbon ions using DNA physics?   (Pedro Gomes - Apr 16, 04:52)
 Add Message Add Message
to: "How to model carbon ions using DNA physics?"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews