[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