[Gate-users] Optical Photon Tracking

Vesna Cuplov vesna.cuplov at gmail.com
Wed Nov 16 14:41:41 CET 2011


Dear Gate users:

I'm a Gate beginner and I'm trying to run the V6.1 version of the code
using optical photons
(I enabled this option in the env_gate.sh).

After a successful compilation, my macro crashes before the end of the
first event processing.

Start Run processing.
/grdm/analogueMC 1
/grdm/verbose 0
/grdm/allVolumes
/gate/decay/setDecayTimeLimit     1000 s
/gate/decay/setPrimaryDecayTimeGeneration 0
=====================================
  G4EventManager::ProcessOneEvent()
=====================================
1 vertices passed from G4Event.
G4PrimaryTransformer::PrimaryVertex (0(mm),0(mm),0(mm),1.75993e+07(nsec))
Primary particle (gamma) --- Transfered with momentum (0,3e-06,0)
A new track 0x10e8b00 (trackID 1, parentID 0) is passed to G4StackManager.
1 primaries are passed from G4EventTransformer.
!!!!!!! Now start processing an event !!!!!!!
### pop requested out of 1 stacked tracks.
Track 0x10e8b00 (trackID 1, parentID 0) is passed to G4TrackingManager.
Track (trackID 1, parentID 0) is processed with stopping code 2
### pop requested out of 0 stacked tracks.
### 0 waiting tracks are re-classified to
     0 urgent tracks and 0 waiting tracks.
NULL returned from G4StackManager.
Terminate current event processing.

 *** Break *** segmentation violation

Something is going bad with the tracking (
http://www-geant4.kek.jp/lxr/source//event/src/G4EventManager.cc#L257)
Below, I copy/paste the relevant pieces of my macro. I use a scanner system
geometry. My module is just an array
of silicon pixels. My phantom is just a water drop.

/gate/physics/addProcess Scintillation
/gate/physics/addProcess OpticalRayleigh
/gate/physics/addProcess OpticalAbsorption
/gate/physics/addProcess OpticalBoundary
/gate/physics/processList Enabled
/gate/physics/processList Initialized

/gate/digitizer/Singles/insert  opticaladder
/gate/digitizer/Singles/insert  readout
/gate/digitizer/Singles/readout/setDepth 1

/gate/output/analysis/disable
/gate/output/fastanalysis/enable

If you use/used optical photons or try/tried to use them, maybe you
encountered the same
problem. I understood from an old post (august 2011) on this mailing list,
the command

/gate/physics/optical/enable 1

doesn't exist anymore. Or the Gate wiki page explains that "To use the
optical photon tracking capabilities of Gate, this option has to be
switched on first." As my issue is related to theTracking, is there
something else I am supposed to do in order to switch on the optical
photon tracking?

Any help is welcome,
vesna


More information about the Gate-users mailing list