Message: cannot reproduce rare events with Random Status saved Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Warning cannot reproduce rare events with Random Status saved 

Forum: Event and Track Management
Date: 05 Feb, 2013
From: Ilya Kudinov <Ilya Kudinov>

Hello,

I have a problem when trying to replay selected rare events by Random Status saving/restoring. My physics list consists of default G4DecayPhysics, G4RadioactiveDecayPhysics and G4EmPenelopePhysics. Radioactive decay of a radioisotope is used as Primary Generator (and random numbers are directly invoked at this stage for the position selection). I'm interested in events when emitted electrons reach detector with a specified energy. Two schemes of Random Status storing were tried (both with no success):

1. Simple and usual one: add

runManager->SetRandomNumberStore(true);

in main.cc and then in Stepping Action call:

G4RunManager::GetRunManager()->rndmSaveThisEvent();

2. A little bit more complicated (as I wasn't sure that the previous one saves Random Status before Primory Generator is invoked): add

runManager->StoreRandomNumberStatusToG4Event(1);

in main.cc and then in Stepping Action execute:

G4String rnd = G4RunManager::GetRunManager()->GetCurrentEvent()->GetRandomNumberStatus();

std::ofstream of; of.open(filename); of << rnd.substr(rnd.find("\n")+1,rnd.npos); of.close();

Doesn't matter which of the above methods was utilized, it didn't provide the selected event reconstruction (after UI command /random/restoreEngineFrom filename). That was clear since no electron penetrated the detector.

What is the reason for such behavior?

Thanks in advance, Ilya Kudinov

P.S. Geant4 version is 9.6, the latest release.

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

1 Feedback: Re: cannot reproduce rare events with Random Status saved   (Gumplinger Peter - 05 Feb, 2013)
(_ Ok: Re: cannot reproduce rare events with Random Status saved   (Ilya Kudinov - 05 Feb, 2013)
(_ Feedback: Re: cannot reproduce rare events with Random Status saved   (Gumplinger Peter - 06 Feb, 2013)
(_ Feedback: Re: cannot reproduce rare events with Random Status saved   (Gumplinger Peter - 06 Feb, 2013)
 Add Message Add Message
to: "cannot reproduce rare events with Random Status saved"

 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 ]