Message: How to change average ionization potential I ? Not Logged In (login)

## 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: Outline Depth: Add message:

1 Re: How to change average ionization potential I ?   (michel maire - 26 Oct, 2013)
 to: "How to change average ionization potential I ?"
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 ]