Message: Tracking and Readout geometry cell number Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Tracking and Readout geometry cell number 

Forum: Event and Track Management
Date: 19 Nov, 2013
From: Anastasia Makarova <Anastasia Makarova>

Hello, dear colleagues

My geometry is pretty simple so far. There is a rectangular water phantom with thin cylindrical readout geometry sensitive volume in the center divided into even segments as the depth increases. A proton beam is of energy 160 MeV and incident perpendicular to the phantom surface, so that the cylindrical sensitive volume coincides with the beam in the water. The uppermost plane of the cylinder coincides with the phantom surface. The step limit is such, that one cell contains exactly one step. So PreStepPoint and PostStepPoint are on the boundaries of cells. In sensitive detector, in ProcessHits(...) I check the coordinates of PreStepPoint and PostStepPoint of G4Step and how they coincide with G4TouchableHistory.GetReplicaNumber().

So, here is the problem:

When the PreStepPoint is on the Phantom boundary and PostStepPoint is on the boundary of the next cell, G4TouchableHistory.GetReplicaNumber() returns 0 (so let's call this cell 0th). But for the next step(PreStepPoint in on the boundary between oth and 1st and PostStepPoint is on the boundary between 1st and 2nd) the G4TouchableHistory.GetReplicaNumber() also returns 0. For next step G4TouchableHistory.GetReplicaNumber() is 1 and for the one after that it's 2. But the 0 cell receives two steps. I use G4TouchableHistory.GetReplicaNumber() to accumulate EnergyDeposited by cell and this results in cell 0 getting twice the energy which is physically wrong.

Please, see code attached.

How can this problem be solved? I will greatly appreciate any help!

Sincerely yours, Anastasia

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2013/11/19/09.42-13148-BWDetectorConstruction.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2013/11/19/09.42-14829-PBWDetectorROGeometry.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2013/11/19/09.42-12882-PBWDetectorSD.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2013/11/19/09.42-56495-BWDetectorConstruction.hh
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2013/11/19/09.42-55082-PBWDetectorROGeometry.hh
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2013/11/19/09.42-16417-PBWDetectorSD.hh

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

1 Idea: Re: Tracking and Readout geometry cell number   (Gumplinger Peter - 19 Nov, 2013)
(_ Ok: Re: Tracking and Readout geometry cell number   (Anastasia Makarova - 21 Nov, 2013)
 Add Message Add Message
to: "Tracking and Readout geometry cell number"

 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 ]