Message: Total energy deposited by each gamma ray Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Total energy deposited by each gamma ray 

Forum: Event and Track Management
Date: 23 Oct, 2004
From: Kiran Bharwani <Kiran Bharwani>

Hello,
  I am trying to plot a histogram (using AIDA) of energy deposited by 
each incident gamma ray in a detector.
So I want something like: 100 gamma rays deposited 500 keV of energy, 
200 gamma rays deposited 300 keV of energy etc.
  I have a physics list like below:

	//	Gamma
	G4ParticleDefinition* gammaParticle = G4Gamma::GammaDefinition();
	G4ProcessManager* pManager = gammaParticle->GetProcessManager();

	// Register processes to gamma’s process manager
	pManager->AddDiscreteProcess(&thePhotoElectricEffect);
	pManager->AddDiscreteProcess(&theComptonScattering);
	pManager->AddDiscreteProcess(&theGammaConversion);

	pManager->AddDiscreteProcess(&theLowEnPhoto);
	pManager->AddDiscreteProcess(&theLowEnCompton);
	pManager->AddDiscreteProcess(&theLowEnPair);
	pManager->AddDiscreteProcess(&theLowEnRayleigh);

	//	Electron
	G4ParticleDefinition* electronParticle = G4Electron::ElectronDefinition();
	pManager = electronParticle->GetProcessManager();

	// Register processes to electron’s process manager
	pManager->AddDiscreteProcess(&theElectronBremsStrahlung);  
	pManager->AddProcess(&theElectronIonisation, ordInActive, 2, 2);

	pManager->AddDiscreteProcess(&theLowEnBremss);
	pManager->AddProcess(&theLowEnIon, ordInActive, 2, 2);

	pManager->AddProcess(&theElectronMultipleScattering);
	pManager->SetProcessOrdering(&theElectronMultipleScattering, idxAlongStep,  1);
	pManager->SetProcessOrdering(&theElectronMultipleScattering, idxPostStep,  1);

How do I get the required information? Which function should I use ? How should I add up energies deposited in the detector. I check in my stepping action if the track volume is my detector but I am still confused. Sometimes the energy adds up to more than the energy of the incident gamma ray (how is that possible). Eg. if I use a 500 keV gamma, and add up all the deposited energy [ step->GetTotalEnergyDeposit() ] it comes out to be more than 500 keV. So, how do I add up the energy deposited by various processes (compton, photoelectric etc.) I just want the total energy deposited by each gamma ray .. does not matter by compton or photoelectric etc. If the gamma delects out of the detector, I want only the energy deposited before it left the detector.

Thank you for your help, KB

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

1 None: Re: Total energy deposited by each gamma ray   (michel maire - 25 Oct, 2004)
 Add Message Add Message
to: "Total energy deposited by each gamma ray"

 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 ]