Message: How to speed up tracking for individual optical photons Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None How to speed up tracking for individual optical photons 

Forum: Event and Track Management
Date: 21 Jul, 2006
From: Jianglai Liu <Jianglai Liu>

Below is a question I asked in the "particle" forum. I am posting
it here as well to look for guidance.

Dear Geanters,

I am working on a MC for a detector with large volume of liquid 
scintillator and PMTs. Apparently, tracking individual optical 
photons is very time consuming. I wonder if there is a simple 
modification to, to only track say 1/N
of the photons, while keeping the mean photoelectrons
and the energy resolution of the detector unchanged. 

I thought about one dirty fix in One can do 
something like:

MeanTrack = MeanNumberOfPhotons/N;
G4double sigma = ResolutionScale * sqrt(MeanTrack)/sqrt(N);
NumTrack = G4int(G4RandGauss::shoot(MeanTrack,sigma)+0.5);

Then go on with the tracking ...

Then in the "digitization" routine, for each track that makes a 
photoelectron, I assign a random weight Poisson(N), and count the 
weight as the true number of PE. 

I am not sure this is an air-tight approach. Would you please share 
your smarter ideas, if any? Thanks a lot in advance!


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

1 Feedback: Re: How to speed up tracking for individual optical photons   (Peter Gumplinger - 02 Aug, 2006)
 Add Message Add Message
to: "How to speed up tracking for individual optical photons"

 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 ]