[Gate-users] benchSPECT run errors

Thibault Frisson frisson at creatis.insa-lyon.fr
Tue Apr 12 11:41:13 CEST 2011


Hi,

To remove the warning message, you should add a model to the e+
Bremsstrahlung:
/gate/physics/addProcess Bremsstrahlung
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+


Thibault


> hi, Didier
>
>    Thank you for your help, I think you conclusion is right, I check
> diable the visualization, it can run properly and the results is
> reasonable.
> But  there is a warning of e+:
>
> Particle: e+
>     * Model(s):
>        - <!> *** Warning *** <!> No model selected!
>
> I don't know if  it will effect the programs.
>
> I checked the run informations with visualization, there is no obvious
> mistake(also have the ablove warning information)
>
> [G4]
> [G4] *************************************************************
> [G4]  Geant4 version Name: geant4-09-03-patch-02    (17-September-2010)
> [G4]                       Copyright : Geant4 Collaboration
> [G4]                       Reference : NIM A 506 (2003), 250-303
> [G4]                             WWW : http://cern.ch/geant4
> [G4] *************************************************************
> [G4]
> [Core-0] Initialization of geometry
> [Core-0] Initialization of physics
> [Core-0] Initialization of actors
> [Core-0]
> [Core-0]
> **********************************************************************
> [Core-0]  GATE version name: gate_v6.1
> [Core-0]                     Copyright : OpenGATE Collaboration
> [Core-0]                     Reference : Phys. Med. Biol. 49 (2004)
> 4543-4561
> [Core-0]                     Reference : Phys. Med. Biol. 56 (2011)
> 881-901
> [Core-0]                     WWW : http://www.opengatecollaboration.org/
> [Core-0]
> **********************************************************************
> [Core-0]
> [Core-0] Starting macro benchSPECT.mac
> [G4] Got standard cmap
> [G4] Window name: viewer-0 (OpenGLStoredX)
> [G4] /tracking/storeTrajectory 1
> [G4]  GateDetectorConstruction::insertARFSD  entered
> [G4]  created a ARF Sensivitive Detector
> [G4]  created command /gate/systems/SPECThead/arf/setARFStage
> [G4]  GetDirectoryName()+ setHeight = /gate/hole/geometry/setHeight
> [G4] Translation speed = 0 0 0.04 cm/s
> [G4] Translation speed = 0 0 0.04 cm/s
> [G4]
> [G4] GATE object:        'systems/SPECThead'
> [G4] Components:
> [G4]
> [G4] GATE object:        'systems/SPECThead/base'
> [G4] Attached to volume: SPECThead
> [G4] Nb of children:       1
> [G4]    GATE object:        'systems/SPECThead/crystal'
> [G4]    Attached to volume: crystal
> [G4]    Nb of children:       1
> [G4]            GATE object:        'systems/SPECThead/pixel'
> [G4]            Attached to volume: ---
> List of Enabled processes:
>    ===  Process: RayleighScattering  ===
> Particle: gamma
>     * Model(s):
>        - PenelopeModel
>    ===  Process: PhotoElectric  ===
> Particle: gamma
>     * Model(s):
>        - StandardModel
>    ===  Process: MultipleScattering  ===
> Particle: e-
>    ===  Process: ElectronIonisation  ===
> Particle: e+
>     * Model(s):
>        - <!> *** Warning *** <!> No model selected!
> Particle: e-
>     * Model(s):
>        - StandardModel
>    ===  Process: Compton  ===
> Particle: gamma
>     * Model(s):
>        - PenelopeModel
>    ===  Process: Bremsstrahlung  ===
> Particle: e+
>     * Model(s):
>        - <!> *** Warning *** <!> No model selected!
> Particle: e-
>     * Model(s):
>        - StandardModel
> List of initialized processes:
> [Core-0] Initialization of geometry
> [Core-0] Initialization of physics
> [Core-0] Initialization of actors
> [G4]            Nb of children:       0
> Confirming confinement to volume 'movsource_phys'...
> [Acquisition-0]
> [Acquisition-0] ============= Source initialization =============
> [Acquisition-0]
> [Acquisition-0] ============= Acquisition starts! =============
> [Acquisition-0] Simulation start time = 0 sec
> [Acquisition-0] Simulation end time   = 600 sec
> [Acquisition-0] Simulation will have  = 16 run(s)
> --------- Ranlux engine status ---------
>  Initial seed = 123456789
>  randoms[] = 0.366829
> 0.959976
> 0.382801
> 0.71216
> 0.91131
> 0.29693
> 0.806733
> 0.788568
> 0.684343
> 0.334206
> 0.0928262
> 0.66146
>  carry = 0, index = 11
>  luxury = 1 pDiscard = 87
> ----------------------------------------
> [Acquisition-0]
> [Acquisition-0] Slice 0 from 0 to 37.5 s [slice=37.5 s], final stop at 600
> s.
> [G4]
> [G4] RayleighScattering:   for  gamma    SubType= 11
> [G4]       Lambda tables from 100 eV  to 10 TeV in 200 bins, spline: 1
> [G4]       ===== EM models for the G4Region  DefaultRegionForTheWorld
> ======
> [G4]          PenRayleigh :     Emin=          0 eV         Emax=   100
> GeV
> [G4]
> [G4] PhotoElectric:   for  gamma    SubType= 12
> [G4]       ===== EM models for the G4Region  DefaultRegionForTheWorld
> ======
> [G4]        PhotoElectric :     Emin=          0 eV         Emax=   10
> TeV
> [G4]
> [G4] Compton:   for  gamma    SubType= 13
> [G4]       Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
> [G4]       ===== EM models for the G4Region  DefaultRegionForTheWorld
> ======
> [G4]           PenCompton :     Emin=          0 eV         Emax=   10
> TeV
> [G4]
> [G4] MultipleScattering:   for e-    SubType= 10
> [G4]       Lambda tables from 100 eV  to 10 TeV in 77 bins, spline: 1
> [G4]       RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1, skin=
> 3, geomFactor= 2.5
> [G4]       ===== EM models for the G4Region  DefaultRegionForTheWorld
> ======
> [G4]           UrbanMsc90 :     Emin=          0 eV         Emax=   10
> TeV
> [G4]
> [G4] ElectronIonisation:   for  e-    SubType= 2
> [G4]       dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
> [G4]       Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
> [G4]       finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1,
> linLossLimit= 0.01
> [G4]       ===== EM models for the G4Region  DefaultRegionForTheWorld
> ======
> [G4]         MollerBhabha :     Emin=          0 eV         Emax=   10
> TeV
> [G4]
> [G4] Bremsstrahlung:   for  e-    SubType= 3
> [G4]       dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
> [G4]       Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
> [G4]       LPM flag: 1 for E > 1 GeV
> [G4]       ===== EM models for the G4Region  DefaultRegionForTheWorld
> ======
> [G4]                eBrem :     Emin=          0 eV         Emax=   1 GeV
> [G4]             eBremRel :     Emin=          1 GeV        Emax=   10
> TeV
> [G4]
> [G4] ElectronIonisation:   for  e+    SubType= 2
> [G4]       dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
> [G4]       Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
> [G4]       finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1,
> linLossLimit= 0.01
> [G4]       ===== EM models for the G4Region  DefaultRegionForTheWorld
> ======
> [G4]         MollerBhabha :     Emin=          0 eV         Emax=   10
> TeV
> [G4]
> [G4] Bremsstrahlung:   for  e+    SubType= 3
> [G4]       dE/dx and range tables from 100 eV  to 10 TeV in 77 bins
> [G4]       Lambda tables from threshold to 10 TeV in 77 bins, spline: 1
> [G4]       LPM flag: 1 for E > 1 GeV
> [G4]       ===== EM models for the G4Region  DefaultRegionForTheWorld
> ======
> [G4]                eBrem :     Emin=          0 eV         Emax=   1 GeV
> [G4]             eBremRel :     Emin=          1 GeV        Emax=   10
> TeV
> [G4]  clearing up the projection data matrices ... done
>
>
> But I found when I run Geant4 examples, there are output information about
> visualization of these:
>
> You have successfully registered the following graphics systems.
> Current available graphics systems are:
>   ASCIITree (ATree)
>   DAWNFILE (DAWNFILE)
>   G4HepRep (HepRepXML)
>   G4HepRepFile (HepRepFile)
>   RayTracer (RayTracer)
>   VRML1FILE (VRML1FILE)
>   VRML2FILE (VRML2FILE)
>   gMocrenFile (gMocrenFile)
>   FukuiRenderer (DAWN)
>   OpenGLImmediateX (OGLIX)
>   OpenGLStoredX (OGLSX)
>   RayTracerX (RayTracerX)
> Registering model factories...
> You have successfully registered the following model factories.
> Registered model factories:
>   generic
>   drawByCharge
>   drawByParticleID
>   drawByOriginVolume
>   drawByAttribute
> Registered filter factories:
>   chargeFilter
>   particleFilter
>   originVolumeFilter
>   attributeFilter
>
> But Gate programs do not have, does it matter?
>
>
> 2011-04-12
>
>
>
> xcao
>
>
>
> ·¢¼þÈË£º Didier Benoit
> ·¢ËÍʱ¼ä£º 2011-04-11  18:13:20
> ÊÕ¼þÈË£º xcao
> ³­ËÍ£º gate-users
> Ö÷Ì⣺ Re: [Gate-users] benchSPECT run errors
>
> Hi Cao,
>
> Your detector seems good on your screenshot, and it's normal for your
> gammas. You have to disable the visualization. The visualization is just
> to check the position of your detector and source.
> Do you have some errors when you launch GATE?
>
> Kind Regards,
> Didier
>
> On 04/06/2011 05:19 AM, xcao wrote:
> hi, gate users!
>
>      I have installed Gate6.1 on Scientific Linux CERN 5.6 (SLC5.6) with
> geant4.9.3.p02 and root-5.26¡¢clhep-2.0.4.5.
> and everything goes well, but when I run benchSPECT for test, it returns
> unproper results:
>     1. The detector is not constructed completely
>     2. There are so many gammas, and it will not stop.
> just like this :
>
>
> I run Geant4 examples as test, and it runs ok, no errors return, and the
> results is correct.
> My PATH and LD_LIBRARY_PATH is :
> $echo $PATH
> /opt/software_simulation/ROOT/root/bin:/opt/software_simulation/Gate/gate_v6.1/bin/Linuxg++:
> /opt/software_simulation/Dawn/dawn_3_90b:/opt/software_simulation/ROOT/root/bin:/usr/sue/sbin:
> /usr/sue/bin:/usr/kerberos/sbin:/usr/kerberos/bin:/usr/local/sbin:/usr/local/bin:/sbin:/bin:/usr/sbin:/usr/bin:
> /usr/X11R6/bin:/root/bin:/opt/software_simulation/Geant4/g4work/bin/Linux-g++
>
> $echo $LD_LIBRARY_PATH
> /opt/software_simulation/Gate/gate_v6.1/tmp/Linux-g++/Gate:/opt/software_simulation/ROOT/root/lib:
> /opt/software_simulation/CLHEP/lib:/opt/software_simulation/CLHEP/lib/:
> /opt/software_simulation/ROOT/root/lib::/opt/software_simulation/Geant4/geant4.9.3.p02/lib/Linux-g++:
> /opt/software_simulation/CLHEP/lib:/opt/software_simulation/lmf/lib
>
> Can anyone give me some helps or any suggestion.
> Thank you!
>
>
> 2011-04-06
>
>
>
> xcao
>
>
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>
>
>
> --
> ----------------------------------------------------------------
> Didier BENOIT, GATE Engineer
>
> Imagerie et Mod¨¦lisation en Neurobiologie et Canc¨¦rologie (IMNC)
> IMNC ¨C UMR 8165 CNRS ¨C Universit¨¦ Paris 7 ¨C Universit¨¦ Paris 11
> Campus d¡¯Orsay
> Batiment 440
> 91405 ORSAY Cedex
>
> T¨¦l : (33) 1 69 15 68 76
> E-mail : benoit at imnc.in2p3.fr
> http://www.imnc.in2p3.fr
> ----------------------------------------------------------------
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users


-- 
Thibault Frisson, Phd
Post-doctorant

Centre de lutte contre le cancer Léon Bérard
Service de radiothérapie
28 rue Laënnec
69373 Lyon cedex 08
Tel : 04 78 78 51 51
      06 70 37 05 21

Laboratoire CREATIS (UMR CNRS 5220)
INSA - Bâtiment Blaise Pascal
7, avenue Jean Capelle
69621 Villeurbanne cedex

http://www.creatis.insa-lyon.fr/rio


More information about the Gate-users mailing list