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: Particles
Date: 21 Jul, 2006
From: Jianglai Liu <Jianglai Liu>

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 G4Scintillation.cc, 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 G4Scintillation.cc. 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!

Regards,
Jianglai

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

 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 ]