Message: Possible bug in extended example rdecay02 Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Possible bug in extended example rdecay02 

Forum: Documentation and Examples
Date: 06 Jan, 2014
From: Raymond Tsang <Raymond Tsang>

I'm using the example: extended/radioactivedecay/rdecay02/ in version 4.9.6.p02 as the starting point of my Germanium detector simulation. 

I noticed that line 215 in src/exrdmAnalysisManager.cc is a little peculiar:
214     for (size_t i = 1; i < fEdepo.size(); i++) {
215       if (std::fabs((fEdepo[i].GetTime()- Time)/second) <= fPulseWidth) {
216         if ( fEdepo[i].GetEnergy() > 0. ) {
217           TarE += fEdepo[i].GetEnergy();
218           TarW += fEdepo[i].GetEnergy()*fEdepo[i].GetWeight();

It seems to me, the "/second" bit is unnecessary. Because fPulseWidth has proper units:
 65 exrdmAnalysisManager::exrdmAnalysisManager()
 66 : fVerbose(0), fNEvt1(-1), fNEvt2(-2),
 67   fHistEMax (15.0*MeV), fHistEMin (0.),fHistNBin(100),
 68   fTargetThresE(10*keV), fDetectorThresE(10*keV),fPulseWidth(1.*microsecond)

Am I misunderstanding it?

Raymond Tsang

 Add Message Add Message
to: "Possible bug in extended example rdecay02"

 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 ]