Message: optical photons in plastic scintillator Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question optical photons in plastic scintillator  

Forum: Particles
Date: 08 Nov, 2013
From: Balint Radics <Balint Radics>

Hello,

I am using geant4-09-05-patch-02. I have a plastic scintillator implementation which has scintillation yield 8400/MeV. I get sometimes several tens of thousands of optical photons in an event in my scintillator and in such a case subsequently the Geant4 runs are extremely slow...(~several seconds per event, while previously I had around 100 times faster event rate). I tried to find the reason for it, I first tried my main guesses: I removed any step limiters I have, I turned off the I/O writing to my output root file - none of them improved the speed.

Following example, in my Physics List I specify the process as:

 G4OpticalPhysics* opticalPhysics = new G4OpticalPhysics();
 fPList->RegisterPhysics( opticalPhysics );

 opticalPhysics->SetWLSTimeProfile("delta");

 opticalPhysics->SetScintillationYieldFactor(1.0);
 opticalPhysics->SetScintillationExcitationRatio(0.0);

 opticalPhysics->SetMaxNumPhotonsPerStep(100);
 opticalPhysics->SetMaxBetaChangePerStep(10.0);

 opticalPhysics->SetTrackSecondariesFirst(kScintillation,true);
 opticalPhysics->SetTrackSecondariesFirst(kCerenkov,true);

Is there any other parameter that might speed up the code? Perhaps some additional/coarser step limits can be configured for these guys? Or other suggestions how to speed up the optical photon processsing in Geant?

Thanks, Balint

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

1 Feedback: Re: optical photons in plastic scintillator   (Gumplinger Peter - 08 Nov, 2013)
(_ Agree: Re: optical photons in plastic scintillator   (Balint Radics - 08 Nov, 2013)
 Add Message Add Message
to: "optical photons in plastic scintillator "

 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 ]