[Gate-users] Optical Photon Tracking

[Vesna Cuplov]

Hello

I was able to make Gate run on my macro enabling the Optical Photons. The
wiki page that I follow (
http://wiki.opengatecollaboration.org/index.php/Users_Guide_V6.1:Generating_and_tracking_optical_photons)
needs to be updated as it has some typos or obsolete commands.

/gate/output/fastanalysis enable  should be used replacing
/gate/output/fastanalysis/enable 1 (which makes the run crash)

vesna







On Wed, Nov 16, 2011 at 2:41 PM, Vesna Cuplov <vesna.cuplov at gmail.com>wrote:

> 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
>



-- 
Dr. Vesna Cuplov
Rice University, Houston, USA