Message: RandomSeed saving in Sensitive Detector Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question RandomSeed saving in Sensitive Detector 

Forum: Event and Track Management
Date: 25 Apr, 2012
From: Nicole Ackerman <Nicole Ackerman>

Hello,

I seem to be quite stuck on retrieving specific events from the Ranecu seeds, even after consulting example N03 (and the other hypernews threads).

I have a Sensitive Detector (SD) implemented with relatively rare interactions (secondary particle passing through and energy deposits) in it - I would like to review the events that did result in an interaction. For the events triggering my SD, there are multiple secondary particles entering it (which I do expect) so my SD object is called multiple times (which I expect).

If I use CLHEP::HepRandom::showEngineStatus(); in ProcessHits() I see a different set of seeds for each secondary track. Is this correct? If so, is one of these the correct seed to use to recall the event?

I've tried CLHEP::HepRandom::saveEngineStatus(fileN); and G4RunManager::GetRunManager()->rndmSaveThisEvent();, both of which save the last seed that I see printed on the screen. I have G4RunManager::GetRunManager()->SetRandomNumberStore(true); in BeginOfRunAction. Using /random/resetEngineFrom those files does not reproduce the event I want. I have tried creating a file with the *first* seed I see printed, but that also does not reproduce the event I want.

Is there something wrong with trying to save the event from within a SD object? I can use the seeds to retrieve a specific event, but it hasn't been the event that I *want*. Can you recommend any further debugging I can do?

The output I am seeing is pasted below, which shows how I am seeing different seeds from within the same event (and my debug output showing event number, SD copy, track). I'm unsure what else would be useful to share.

Thank you,
Nicole Ackerman
------------------------------------------

Event number 70000 processed.
e 70726 cellDet55 track: 2

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 2058228980, 106044651
----------------------------------------
e 70726 cellDet55 track: 2

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 116055488, 1006125497
----------------------------------------
e 70726 cellDet55 track: 2

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 1692605599, 161040593
----------------------------------------
e 70726 cellDet55 track: 2

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 1696737136, 804001109
----------------------------------------
e 70726 cellDet55 track: 2

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 1118943336, 921282989
----------------------------------------
e 70726 cellDet55 track: 2

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 1161310476, 1545743581
----------------------------------------
e 70726 cellDet55 track: 2

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 1428050470, 1756524741
----------------------------------------
e 70726 cellDet55 track: 7

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 344358947, 1998522247
----------------------------------------
e 70726 cellDet55 track: 6

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 924365050, 818438193
----------------------------------------
e 70726 cellDet55 track: 5

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 402010132, 1942593024
----------------------------------------
e 70726 cellDet55 track: 4

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 1381534978, 1278898817
----------------------------------------
e 70726 cellDet55 track: 3

--------- Ranecu engine status ---------
 Initial seed (index) = 0
 Current couple of seeds = 1304554019, 112121197
----------------------------------------
Event number 80000 processed.

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

1 Feedback: Re: RandomSeed saving in Sensitive Detector   (Gumplinger Peter - 25 Apr, 2012)
(_ Ok: Re: RandomSeed saving in Sensitive Detector   (Nicole Ackerman - 01 May, 2012)
 Add Message Add Message
to: "RandomSeed saving in Sensitive Detector"

 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 ]