Message: How to flag a scattered track Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question How to flag a scattered track 

Keywords: scattered photons, flag track
Forum: Event and Track Management
Date: 05 Jun, 2008
From: Francisco <Francisco>

Hi
I am trying to do this by myself, but I give up.
I have my detector and a water phantom underneath. I want to get both primary 
and backscattered spectra for the energy deposition. The scattered spectrum 
should be the one due to photons that have scattered in the water phantom, so
I need to have two outputs, one for energy deposition in my detector 
due to primary or scattered photons in any volume apart from the phantom, and
other for scattered photons coming from the phantom. So when I do the following 

for stepping action:

if (volume==epitaxial){

  G4double EdepStep = aStep->GetTotalEnergyDeposit();
  if (EdepStep > 0.) {  
     runAction->AddEdep(EdepStep);
     eventAction->AddEdep(EdepStep);
  } 

 ..............
 ..............
}


for event action (I retrieve the deposited energy):

totE += (*CHC)[i]->GetEdep();

 .............
 .............

  ofstream outfile3("EnergyDeposited.dat", ios::app);
  outfile3.precision(3);
  outfile3.setf(ios::scientific);
  outfile3 << TotalEnergyDeposit/keV << G4endl;


Geant4 obviously does not differentiate between photons coming from my phantom
and those from anywhere else. So my question is:

How can I make Geant4 to put a flag in the tracks scattering in my water phantom
in such a way that when they hit my detector (in the volume called epitaxial)
the energy deposited by them can be retrieved apart from any other interaction in 
any other volume.  

Thanks

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

1 Feedback: Re: How to flag a scattered track   (Peter Gumplinger - 05 Jun, 2008)
 Add Message Add Message
to: "How to flag a scattered track"

 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 ]