Message: How to change average ionization potential I ? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question How to change average ionization potential I ? 

Forum: Medical Applications
Date: 22 Oct, 2013
From: Asya Demir <Asya Demir>

Dear experts,

I inject C12 ions on water and I am able to obtain the Bragg peak by making changes in the TestEm7 example. Below are parts of the changes I made in SteppingAction.cc:

 G4Track* fTrack = step->GetTrack();
  G4double edep = step->GetTotalEnergyDeposit();
  if (edep <= 0.) return;  

  G4double niel = step->GetNonIonizingEnergyDeposit();

  runAction->FillEdep(edep, niel);

  if (step->GetTrack()->GetTrackID() == 1) runAction->AddPrimaryStep();  

  //Bragg curve
  //	
  G4StepPoint* prePoint  = step->GetPreStepPoint();
  G4StepPoint* postPoint = step->GetPostStepPoint();
  G4int ih;G4double weight;
  G4double x1 = prePoint->GetPosition().x();
  G4double x2 = postPoint->GetPosition().x();  
  G4double x  = x1 + G4UniformRand()*(x2-x1) + 0.5*(detector->GetAbsorSizeX());
  // histoManager->FillHisto(1,x, edep);

  histoManager->FillHisto(2,x, edep);
  histoManager->FillHisto(1,x, edep);
  // if (fTrack->GetDefinition()->GetParticleName()=="proton")
  //{

Now, I would like to choose a different ionization potential value (I) instead of the default one to see how it alters the Bragg curve. Unfortunately I cannot find the location of the formula including "I". Hoping that my question is not silly, I would appreciate any help.

Asya

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

1 None: Re: How to change average ionization potential I ?   (michel maire - 26 Oct, 2013)
 Add Message Add Message
to: "How to change average ionization potential I ?"

 Subscribe Subscribe

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


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]