[Gate-users] Sinogram Output

arda arda-konik at uiowa.edu
Wed Jul 2 23:13:17 CEST 2008


 Hello,
 I am a beginner in GATE. 
(CLHEP 2.0.3.1 / Geant4 9.0 / root 5.14 / gate3.1.2)

 I am trying to obtain sinogram data from "PET_Ecat_System.mac" example
and visualize them on IDL.  I read the  "yourSinogram_1.ima" and
"yourSinogram_1_del.ima" (both 162 Mb) data on IDL and the entire data
are zero.
 Attached the log file and the macro I used

1) Could you tell me what I am doing wrong ?
2) Also, what are the differences between "yourSinogram_1.ima",
"yourSinogram_1_del.ima" and "yourSinogram_1.sct.ima" output data. Which
one I need to use on IDL to read the actuall sinogram?

Thanks,
Arda

 


>         
-------------- next part --------------

*************************************************************
 Geant4 version Name: geant4-09-00-patch-01    (28-August-2007)
                      Copyright : Geant4 Collaboration
                      Reference : NIM A 506 (2003), 250-303
                            WWW : http://cern.ch/geant4
*************************************************************

Visualization Manager instantiating...
Visualization Manager initialising...
Registering graphics systems...

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

Time set to (s) 0
Got standard cmap
Window name: viewer-0 (OpenGLStoredX)
WARNING: The vis manager will keep up to 100 events.
  This may use a lot of memory.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
[GateSystemListManager::CheckScannerAutoCreation:
	Creating new system based on volume inserter 'ecat'
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
   GateMaterialDatabase: did not load properties table for material 'BGO'.
     This is only a problem when OPTICAL PHOTONS are transported in this material.'
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
   GateMaterialDatabase: did not load properties table for material 'Tungsten'.
     This is only a problem when OPTICAL PHOTONS are transported in this material.'
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
   GateMaterialDatabase: did not load properties table for material 'Water'.
     This is only a problem when OPTICAL PHOTONS are transported in this material.'
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
Nof Higher level 2

GATE object:        'digitizer/Coincidences'
Is enabled?         Yes
Coincidence window: 10 ns 
Coincidence window jitter: 0 ps 
Coincidence offset: 0 ps 
Coincidence offset jitter: 0 ps 
Min sector diff.:   2
Input:              'Singles'
Output:             'Coincidences'

GATE object:        'digitizer/delay'
Is enabled?         Yes
Coincidence window: 10 ns 
Coincidence window jitter: 0 ps 
Coincidence offset: 500 ns 
Coincidence offset jitter: 0 ps 
Min sector diff.:   2
Input:              'Singles'
Output:             'delay'

GATE object:        'digitizer/finalCoinc'
Is enabled?         Yes
Nb of elements:     0
Attached to:        'digitizer'
Output:             'finalCoinc'

Output module: 'sinogram'
 >> Job:                                 build a set of 2D sinograms from a PET simulation
 >> Is enabled ?                         Yes
 >> Number of crystals per crystal ring: 0
 >> Number of crystal rings:             0
 >> Number of radial sinogram bins:      0
 >> Filled ?                             No
 >> Attached to system:                  systems/ecat
 >> Input data:                          finalCoincGateSourceMgr::ListSources: List of the sources in the source managaer
Source --------------> F18CylinderSource
  ID                 : 0
  type               : 
  activity (Bq)      : 10000
  startTime (s)      : 0
  time (s)           : 0
  forcedUnstable     : 1
  forcedHalfLife (s) : 6586
  verboseLevel       : 0
---------------------- 

Time set to (s) 0
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
 >> entering [GateToSinogram::RecordBeginOfAcquisition]
    Acquisition start time: 5.93605e+228
    Number of frames: 0
    Number of crystals per crystal rings: 576
    Number of crystal rings:              32
    Number of radial sinogram bins:    288
    Number of azimuthal sinogram bins: 288
    Crystal location blurring in tangential direction: 1.8 mm
    Crystal location blurring in axial direction: 1.8 mm
    Prompt coincidences are recorded in data = 0
    Delayed coincidences (if simulated) are recorded in data = 1
    True scattered coincidences (eventID1 = eventID2) are recorded in data = 4
 >> leaving [GateToSinogram::RecordBeginOfAcquisition]
G4VRangeToEnergyConverter::ConvertCutToKineticEnergy   for e-
The cut in range [1000 (mm)]   is too big   for material  idx=2
The cut in energy is set1.79769e+305GeV 

phot:  Total cross sections from Sandia parametrisation. 
      Sampling according PhotoElectric model

compt:  Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV
      Sampling according Klein-Nishina model
      tables are built for  gamma
      Lambda tables from 100 eV  to 100 GeV in 90 bins.

eBrem:   tables are built for  e-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Total cross sections and sampling from StandBrem model (based on the EEDL data library) 
      Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1

eIoni:   tables are built for  e-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Delta cross sections and sampling from MollerBhabha model
      Good description from 1 KeV to 100 GeV.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1

msc:  Model variant of multiple scattering for e-
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      LateralDisplacementFlag=  1   Skin= 0
      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 1

eBrem:   tables are built for  e+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Total cross sections and sampling from StandBrem model (based on the EEDL data library) 
      Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1

annihil:       Sampling according eplus2gg model
      tables are built for  e+
      Lambda tables from 100 eV  to 100 TeV in 120 bins.

annihil:       Sampling according eplus2gg model
      tables are built for  e+
      Lambda tables from 100 eV  to 100 TeV in 120 bins.

eIoni:   tables are built for  e+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Delta cross sections and sampling from MollerBhabha model
      Good description from 1 KeV to 100 GeV.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1

hIoni:   tables are built for  proton
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Scaling relation is used from proton dE/dx and range.
      Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
      Parametrisation from Bragg for protons below.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1

msc:  Model variant of multiple scattering for proton
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      LateralDisplacementFlag=  1   Skin= 0
      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 1

msc:  Model variant of multiple scattering for GenericIon
      LateralDisplacementFlag=  0   Skin= 0
      Boundary/stepping algorithm is active with RangeFactor= 0.2  Step limit type 1

hIoni:   tables are built for  anti_proton
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Scaling relation is used from proton dE/dx and range.
      Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV
      Parametrisation from Bragg for protons below.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1

hIoni:   tables are built for  pi+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Scaling relation is used from proton dE/dx and range.
      Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
      Parametrisation from Bragg for protons below.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1

hIoni:   tables are built for  pi-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Scaling relation is used from proton dE/dx and range.
      Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV
      Parametrisation from Bragg for protons below.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1

msc:  Model variant of multiple scattering for pi-
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      LateralDisplacementFlag=  1   Skin= 0
      Boundary/stepping algorithm is active with RangeFactor= 0.02  Step limit type 1

========= Table of registered couples ==============================

Index : 0     used in the geometry : Yes     recalculation needed : No 
 Material : Air
 Range cuts        :  gamma 1 mm     e- 1 m      e+ 1 mm 
 Energy thresholds :  gamma 990 eV     e- 378.379 keV    e+ 990 eV 
 Region(s) which use this couple : 
    DefaultRegionForTheWorld

Index : 1     used in the geometry : Yes     recalculation needed : No 
 Material : BGO
 Range cuts        :  gamma 1 mm     e- 1 m      e+ 1 mm 
 Energy thresholds :  gamma 78.3228 keV    e- 5.58977 TeV    e+ 976.329 keV
 Region(s) which use this couple : 
    DefaultRegionForTheWorld

Index : 2     used in the geometry : Yes     recalculation needed : No 
 Material : Tungsten
 Range cuts        :  gamma 1 mm     e- 1 m      e+ 1 mm 
 Energy thresholds :  gamma 105.401 keV    e- 2.88022e+295 J      e+ 2.14666 MeV
 Region(s) which use this couple : 
    DefaultRegionForTheWorld

Index : 3     used in the geometry : Yes     recalculation needed : No 
 Material : Water
 Range cuts        :  gamma 1 mm     e- 1 m      e+ 1 mm 
 Energy thresholds :  gamma 2.90186 keV    e- 317.969 MeV    e+ 338.695 keV
 Region(s) which use this couple : 
    DefaultRegionForTheWorld

====================================================================


----- HepJamesRandom engine status -----
 Initial seed = 9876
 u[] = 0.174685 0.883674 0.958512 0.690013 0.0250537 0.142294 0.783649 0.356081 0.214774 0.69287 0.700302 0.465109 0.22215 0.768752 0.461089 0.397602 0.522215 0.121546 0.960227 0.510892 0.0116568 0.0286447 0.238777 0.585219 0.780612 0.709182 0.082106 0.828601 0.61446 0.0798209 0.688002 0.00978529 0.370831 0.376435 0.710184 0.864017 0.167338 0.21114 0.349235 0.522479 0.943713 0.0702708 0.778742 0.639458 0.506555 0.151751 0.990123 0.564688 0.430497 0.889832 0.995398 0.725368 0.710512 0.0550886 0.0540335 0.0159176 0.728379 0.947141 0.34922 0.20551 0.351029 0.672539 0.118976 0.655784 0.195375 0.0380308 0.293556 0.237497 0.79723 0.158375 0.029761 0.602215 0.616657 0.995638 0.734827 0.943425 0.130007 0.277064 0.306344 0.973337 0.183926 0.547536 0.914344 0.666572 0.210074 0.300351 0.83602 0.267525 0.248549 0.533332 0.1754 0.299529 0.538308 0.555472 0.283894 0.485254 0.222273 
 c = 0.0216029, cd = 0.456233, cm = 1
 i97 = 96, u[i97] = 0.222273
 j97 = 32, u[j97] = 0.370831
----------------------------------------
 >> entering [GateToSinogram::RecordBeginOfRun]
    Frame ID = 1
 >> leaving [GateToSinogram::RecordBeginOfRun]
WARNING: G4VisManager::EndOfEvent: Event keeping suspended.
  The number of events exceeds the maximum that may be kept, 100.

----- HepJamesRandom engine status -----
 Initial seed = 9876
 u[] = 0.572855 0.313373 0.0667574 0.668544 0.538844 0.908549 0.560988 0.721999 0.728548 0.506417 0.189487 0.585622 0.889681 0.000579596 0.448573 0.36662 0.36138 0.473628 0.942809 0.91608 0.915743 0.123247 0.805196 0.643545 0.335681 0.304751 0.0588319 0.76266 0.18071 0.65194 0.1687 0.234811 0.263546 0.0220037 0.655546 0.0708445 0.293801 0.609984 0.0930434 0.700128 0.546337 0.594176 0.0117456 0.750611 0.816423 0.97325 0.842786 0.44602 0.900362 0.17387 0.891011 0.23293 0.432541 0.210821 0.324536 0.553205 0.559053 0.708254 0.54576 0.572873 0.97042 0.729443 0.309828 0.105378 0.584256 0.864154 0.822712 0.2178 0.410484 0.119888 0.0876788 0.71059 0.829143 0.783047 0.95268 0.470871 0.771367 0.973968 0.979513 0.976906 0.807165 0.23837 0.390288 0.134977 0.420943 0.859059 0.609108 0.444528 0.324008 0.666146 0.754101 0.831177 0.478049 0.451316 0.606791 0.783734 0.674524 
 c = 0.0966497, cd = 0.456233, cm = 1
 i97 = 72, u[i97] = 0.829143
 j97 = 8, u[j97] = 0.728548
----------------------------------------
 >> entering [GateToSinogram::RecordEndOfRun]
    Number of prompt coincidences for all ring combinations        276
      Number of true coincidences for all ring combinations          276
        Number of scattered coincidences for all ring combinations     113
      Number of random coincidences for all ring combinations        0
    Number of delayed coincidences for all ring combinations       0
    sinograms 1,1,1,0,0 written to the raw file yourSinogram_1.ima
    sinograms 1,1,1,1,0 written to the raw file yourSinogram_1_del.ima
    sinograms 1,1,1,4,0 written to the raw file yourSinogram_1_sct.ima
 >> leaving [GateToSinogram::RecordEndOfRun]
WARNING: G4VisManager::EndOfRun: Event keeping was suspended.
  The number of events in the run exceeded the maximum to be kept, 100.
  The number of events to be kept can be changed with
  "/vis/scene/endOfEventAction accumulate <N>", where N is the
  maximum number you wish to allow.  N < 0 means "unlimited".
100 events have been kept for refreshing and/or reviewing.
Time set to (s) 1
G4PhysicalVolumeModel::Validate() called.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  Be warned that this does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
 Current Output Root File: 0x89a5f38 named Your_PET_Simu.root
GateToRoot: ROOT: files writing...
GateToRoot: ROOT: files closing...
Graphics systems deleted.
Visualization Manager deleting...
-------------- next part --------------
#  *********************************************************************
#  * DISCLAIMER                                                        *
#  *                                                                   *
#  * Neither the authors of this software system, nor their employing  *
#  * institutes, nor the agencies providing financial support for this *
#  * work make any representation or warranty, express or implied,     *
#  * regarding this software system or assume any liability for its    *
#  * use.                                                              *
#  *                                                                   *
#  * This code implementation is the intellectual property of the      *
#  * OpenGATE collaboration.                                           *
#  * By copying, distributing or modifying the Program (or any work    *
#  * based on the Program) you indicate your acceptance of this        *
#  * statement, and all its terms.                                     *
#  *********************************************************************
#
#  This macro has been prepared by:  
#     Sebastien JAN

#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
# P E T   S C A N N E R   G E O M E T R Y   W I T H   E C A T     #
# S Y S T E M   D E F I N I T I O N   A N D   E C A T 7           #
#                    F O R M A T   O U T P U T                    #
#                                                                 #
#            Y O U   N E E D   T O   I N S T A L L                #
#              T H E   E C A T 7   L I B R A R Y                  #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#

###############################################################
#                                                             #
#                                                             #
# D E F I N E   Y O U R   P E T   S C A N N E R   A N D       #
#                                                             #
# Y O U R   C O M P L E T E   G A T E   S I M U L A T I O N   #
#                                                             #
#                                                             #
###############################################################


#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
# R U N   T H E   V I S U A L I S A T I O N   A N D               #
# C O N S T R U C T   Y O U R   E X P E R I M E N T A L   A R E A #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#




#
#	V I S U A L I S A T I O N
#
/vis/open OGLSX
/vis/viewer/reset
/vis/drawView 45 45
/vis/viewer/zoom 4
/vis/viewer/set/style surface
/vis/drawVolume
/tracking/storeTrajectory 1
/vis/scene/endOfEventAction accumulate
/vis/viewer/update
/gate/geometry/enableAutoUpdate

# OR disable the visu
#/vis/disable


#
#     W O R L D
#
/gate/world/geometry/setXLength 400. cm
/gate/world/geometry/setYLength 400. cm
/gate/world/geometry/setZLength 400. cm


#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
#     D E F I N I T I O N   A N D   D E S C R I T I O N           #
#        O F   Y O U R   P E T   D E V I C E                      #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#


#     E C A T
/gate/world/daughters/name ecat
/gate/world/daughters/insert cylinder
/gate/ecat/setMaterial Air
/gate/ecat/geometry/setRmax 44.2 cm
/gate/ecat/geometry/setRmin 41.2 cm
/gate/ecat/geometry/setHeight 15.52 cm
/gate/ecat/vis/forceWireframe


#     B L O C K
/gate/ecat/daughters/name block
/gate/ecat/daughters/insert box
/gate/block/placement/setTranslation 427.0 0.0 0.0 mm
/gate/block/geometry/setXLength 30.0 mm
/gate/block/geometry/setYLength 35.8594 mm
/gate/block/geometry/setZLength 38.7 mm
/gate/block/setMaterial Air
/gate/block/vis/forceWireframe



#	C R Y S T A L
/gate/block/daughters/name crystal
/gate/block/daughters/insert box
/gate/crystal/geometry/setXLength 30.0 mm
/gate/crystal/geometry/setYLength 4.4 mm
/gate/crystal/geometry/setZLength 4.75 mm
/gate/crystal/setMaterial BGO
/gate/crystal/vis/setColor yellow




#	R E P E A T    C R Y S T A L
/gate/crystal/repeaters/insert cubicArray
/gate/crystal/cubicArray/setRepeatNumberX 1
/gate/crystal/cubicArray/setRepeatNumberY 8
/gate/crystal/cubicArray/setRepeatNumberZ 8
/gate/crystal/cubicArray/setRepeatVector 0. 4.4942 4.85 mm


#	R E P E A T    BLOCK
/gate/block/repeaters/insert linear
/gate/block/linear/setRepeatNumber 4
/gate/block/linear/setRepeatVector 0. 0. 38.8 mm
/gate/block/repeaters/insert ring
/gate/block/ring/setRepeatNumber 72


#      T U N G S T E N    S H I E L D 
/gate/world/daughters/name   carter
/gate/world/daughters/insert cylinder
/gate/carter/setMaterial Air
/gate/carter/geometry/setRmax   44.  cm
/gate/carter/geometry/setRmin   28.  cm
/gate/carter/geometry/setHeight 9. cm
/gate/carter/placement/setTranslation 0.0 0.0 12.5 cm
/gate/carter/vis/forceWireframe

/gate/carter/daughters/name carter1
/gate/carter/daughters/insert cylinder
/gate/carter1/setMaterial Tungsten
/gate/carter1/geometry/setRmax 30.0 cm
/gate/carter1/geometry/setRmin 29.0 cm
/gate/carter1/geometry/setHeight 3.0 cm
/gate/carter1/vis/setColor grey
/gate/carter1/placement/setTranslation 0.0 0.0 -3 cm

/gate/carter/daughters/name carter2
/gate/carter/daughters/insert cylinder
/gate/carter2/setMaterial Tungsten
/gate/carter2/geometry/setRmax 40.7 cm
/gate/carter2/geometry/setRmin 30.0 cm
/gate/carter2/geometry/setHeight 1.0 cm
/gate/carter2/vis/setColor grey
/gate/carter2/placement/setTranslation 0.0 0.0 -2 cm

/gate/carter/daughters/name carter3
/gate/carter/daughters/insert cylinder
/gate/carter3/setMaterial Tungsten
/gate/carter3/geometry/setRmax 40.7 cm
/gate/carter3/geometry/setRmin 39.7 cm
/gate/carter3/geometry/setHeight 3.0 cm
/gate/carter3/vis/setColor grey
/gate/carter3/placement/setTranslation 0.0 0.0 0 cm

/gate/carter/daughters/name carter4
/gate/carter/daughters/insert cylinder
/gate/carter4/setMaterial Tungsten
/gate/carter4/geometry/setRmax 43.7 cm
/gate/carter4/geometry/setRmin 40.7 cm
/gate/carter4/geometry/setHeight 1.0 cm
/gate/carter4/vis/setColor grey
/gate/carter4/placement/setTranslation 0.0 0.0 1 cm


#     A T T A C H    S Y S T E M 
/gate/systems/ecat/block/attach block
/gate/systems/ecat/crystal/attach crystal

#	A T T A C H    C R Y S T A L  SD
/gate/crystal/attachCrystalSD



#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
#     D E F I N I T I O N   A N D   D E S C R I T I O N           #
#             O F   Y O U R   P H A N T O M                       #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#


/gate/world/daughters/name phantom
/gate/world/daughters/insert cylinder
/gate/phantom/setMaterial Water
#/gate/phantom/vis/forceWireframe
/gate/phantom/vis/setColor green
/gate/phantom/geometry/setRmax   10.0 cm
/gate/phantom/geometry/setHeight 20. cm
/gate/phantom/attachPhantomSD



#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
#            S E C O N D   P A R T   O F   Y O U R                #
#                P E T   S I M U L A T I O N                      #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#



#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
#     D E F I N I T I O N   A N D   D E S C R I T I O N           #
#             O F   Y O U R   P H Y S I C S                       #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#

#       EM  P R O C E S S
/gate/physics/gamma/selectRayleigh lowenergy
/gate/physics/gamma/selectPhotoelectric standard
/gate/physics/gamma/selectCompton standard

#	I N A C T I V E    S E C O N D A R Y   E L E C T R O N S 
/gate/physics/setElectronCut 1. m

#	I N A C T I V E    X - R A Y S
/gate/physics/setXRayCut 1. GeV
/gate/physics/setDeltaRayCut 1. GeV



#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
#     I N I T I A L I Z A T I O N   O F   Y O U R                 #
#                 S I M U L A T I O N                             #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#


/gate/geometry/enableAutoUpdate
/run/initialize




#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
#  D E F I N I T I O N   O F   Y O U R   A C Q U I S I T I O N    #
#  D I G I T I Z E R   &   C O I N C I D E N C E   S H O R T E R  #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#

#	A D D E R
/gate/digitizer/Singles/insert adder


#       R E A D O U T 
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 1


#      E N E R G Y   B L U R R I N G
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.26
/gate/digitizer/Singles/blurring/setEnergyOfReference 511. keV


#       E N E R G Y   C U T
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 250. keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 750. keV


#       D E A D    T I M E 
/gate/digitizer/Singles/insert deadtime
/gate/digitizer/Singles/deadtime/setDeadTime 3000000. ps
/gate/digitizer/Singles/deadtime/setMode paralysable
/gate/digitizer/Singles/deadtime/chooseDTVolume block


#	C O I N C I    S O R T E R 
/gate/digitizer/Coincidences/setWindow 10. ns
/gate/digitizer/Coincidences/setOffset 0. ns
/gate/digitizer/Coincidences/describe

/gate/digitizer/name delay
/gate/digitizer/insert coincidenceSorter
/gate/digitizer/delay/setWindow 10. ns
/gate/digitizer/delay/setOffset 500. ns
/gate/digitizer/delay/describe


/gate/digitizer/name finalCoinc
/gate/digitizer/insert coincidenceChain
/gate/digitizer/finalCoinc/addInputName delay
/gate/digitizer/finalCoinc/addInputName Coincidences
/gate/digitizer/finalCoinc/usePriority true
/gate/digitizer/finalCoinc/describe



#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
#                   D E F I N I T I O N   O F                     #
#                Y O U R   O U T P U T   F I L E                  #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#

#       C H A N G E    T H E     S E E D (1)  O R   N O T (0)
/gate/output/root/setSaveRndmFlag 1

#       S E T U P  - R O O T  F I L E
/gate/output/root/setFileName Your_PET_Simu
/gate/output/root/setRootSinglesAdderFlag 1
/gate/output/root/setRootSinglesReadoutFlag 1
/gate/output/root/setRootHitFlag 1
/gate/output/root/setRootSinglesFlag 1
/gate/output/root/setRootNtupleFlag 1
/gate/output/root/setRootCoincidencesFlag 1

#       S E T U P  - A S C I I  F I L E
/gate/output/ascii/setOutFileSinglesAdderFlag 0
/gate/output/ascii/setOutFileSinglesReadoutFlag 0
/gate/output/ascii/setOutFileSinglesThresholderFlag 0
/gate/output/ascii/setOutFileSinglesUpholderFlag 0
/gate/output/ascii/setOutFileSinglesFlag 0
/gate/output/ascii/setOutFileHitsFlag 0
/gate/output/ascii/setOutFileCoincidencesFlag 0

###################################################
#                E C A T 7   O U T P U T
#
#  W I T H   N E W   F U N C T I O N A L I T I E S
#  *-*-* according to the coincidence sorter *-*-*
#
#             01  J U N E  2 0 0 6
###################################################
/gate/output/sinogram/enable
/gate/output/sinogram/setTangCrystalBlurring 1.8 mm
/gate/output/sinogram/setAxialCrystalBlurring 1.8 mm
/gate/output/sinogram/verbose 2
#/gate/output/sinogram/RawOutputEnable false
/gate/output/sinogram/RawOutputEnable
/gate/output/sinogram/setFileName yourSinogram
/gate/output/sinogram/StoreDelayeds
/gate/output/sinogram/StoreScatters
/gate/output/sinogram/setInputDataName finalCoinc
/gate/output/sinogram/describe


#/gate/output/ecat7/enable
#/gate/output/ecat7/verbose 2
#/gate/output/ecat7/setFileName yourSinogram
#/gate/output/ecat7/describe
#/gate/output/ecat7/mashing 2
#/gate/output/ecat7/span 9
#/gate/output/ecat7/maxringdiff 22
#/gate/output/ecat7/system 962
#/gate/output/ecat7/IsotopeCode F-18
#/gate/output/ecat7/IsotopeHalflife 6586.2 second
#/gate/output/ecat7/IsotopeBranchingFraction 1.0


#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
#                   D E F I N I T I O N   O F                     #
#            Y O U R   V E R B O S I T Y   L E V E L              #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#

/control/verbose 0
/grdm/verbose 0
/run/verbose 0
/event/verbose 0
/tracking/verbose 0
/gate/application/verbose 0
/gate/generator/verbose 0
/gate/stacking/verbose 0
/gate/event/verbose 0
/gate/source/verbose 0


#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
#                   D E F I N I T I O N   O F                     #
#                    Y O U R   S O U R C E S                      #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#


# F A S T   P R O D U C T I O N   M O D U L E   F O R   F - 1 8   I S O T O P E S
/gate/source/addSource F18CylinderSource
/gate/source/F18CylinderSource/setActivity 10000. becquerel
/gate/source/F18CylinderSource/gps/particle e+
/gate/source/F18CylinderSource/setForcedUnstableFlag true
/gate/source/F18CylinderSource/setForcedHalfLife 6586 s
/gate/source/F18CylinderSource/gps/energytype Fluor18
/gate/source/F18CylinderSource/gps/type Volume
/gate/source/F18CylinderSource/gps/shape Cylinder
/gate/source/F18CylinderSource/gps/radius 0.5 mm
/gate/source/F18CylinderSource/gps/halfz 0.5 mm
/gate/source/F18CylinderSource/gps/angtype iso
/gate/source/F18CylinderSource/gps/centre 0.0 0.0 0. cm


/gate/source/list


#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#                                                                 #
#               S T A R T   A C Q U I S I T I O N                 #
#                 T I M E   P A R A M E T E R S                   #
#                                                                 #
#-------------------oooooOOOOO00000OOOOOooooo---------------------#

#
# EXPERIMENT
#
/gate/application/setTimeSlice     1.  s
/gate/application/setTimeStart     0.   s
/gate/application/setTimeStop      1.  s

#
# LET'S RUN THE SIMULATION!
#
/gate/application/startDAQ







More information about the Gate-users mailing list