<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
<div>Dear Gate users </div><div><br></div>I try to follow all the steps to run CT scanner, I have an image of 512 x 512 pixels. I'm voxelez the phantom and the source (position = -256, -256, 0) but I have not seen any interaction and I have this message.<div><br></div><div><br></div><div>Can you help me</div><div><br></div><div>Thank you <br><br><div><br></div><div><br></div><div>/tracking/storeTrajectory 1</div><div>Attributes available for modeling and filtering with</div><div>  "/vis/modeling/trajectories/create/drawByAttribute" and</div><div>  "/vis/filtering/trajectories/create/attributeFilter" commands:</div><div>G4Trajectory:</div><div>  Charge (Ch): unit: e+ (G4double)</div><div>  Track ID (ID): G4int</div><div>  Initial kinetic energy (IKE): G4BestUnit (G4double)</div><div>  Initial momentum magnitude (IMag): G4BestUnit (G4double)</div><div>  Initial momentum (IMom): G4BestUnit (G4ThreeVector)</div><div>  No. of points (NTP): G4int</div><div><br></div> <br>   PDG Encoding (PDG): G4int<div>  Parent ID (PID): G4int</div><div>  Particle Name (PN): G4String</div><div>G4TrajectoryPoint:</div><div>  Position (Pos): G4BestUnit (G4ThreeVector)</div><div><font color="#ff0000">WARNING</font>: Trajectory storing has been requested.  This action may be</div><div>  reversed with "/tracking/storeTrajectory 0".</div><div><font color="#ff0000">WARNING:</font> The vis manager will keep up to 100 events.</div><div>  This may use a lot of memory.</div><div>  It may be changed with, e.g., "/vis/scene/endOfEventAction accumulate 10".</div><div> CREATED THE COMMAND /gate/systems/cylindricalPET/addAnewRsector</div><div> Call VoxelParameterized::InsertReader</div><div> Range Translator</div><div>        Range   0 : imageValue in [    0.000e+00 , 1.500e+03 ]  ---> material Air, visibility 1, coulour (0.000e+00,1.000e+00,0.000e+00,2.000e-01)</div><div>        Range   1 : imageValue in [    1.501e+03 , 2.250e+03 ]  ---> material Water, visibility 1, coulour (1.000e+00,0.000e+00,0.000e+00,2.000e-01)</div><div>        Range   2 : imageValue in [    2.251e+03 , 2.750e+03 ]  ---> material Air, visibility 1, coulour (0.000e+00,1.000e+00,0.000e+00,2.000e-01)</div><div>        Range   3 : imageValue in [    2.751e+03 , 3.848e+03 ]  ---> material RibBone, visibility 1, coulour (0.000e+00,0.000e+00,1.000e+00,2.000e-01)</div><div>', assuming default BIGENDIAN'LITTLEENDIAN</div><div><br></div><div>'nrecognised type name 'unsigned integer</div><div> Header read from       'CT-IM00120_uint16.h33'</div><div>'Data file name         'm001-IM00120.i33</div><div> Nb of planes:           1</div><div> Nb of pixels per plane: 512 512</div><div> Pixel size:             4.590e-01 4.590e-01</div><div> Slice thickness:        2.179e+00</div><div> Matrix size:            2.350e+02 2.350e+02</div><div> Data type:              </div><div> Data byte order:        BIGENDIAN</div><div><br></div><div><br></div><div>'<font color="#ff0000">!ror:</font> Could not open header file 'm001-IM00120.i33</div><div> <!> *** WARNING *** <!>  Unknown process: PositronAnnihilationStd</div><div> <!> *** WARNING *** <!>  Unknown process: MultipleScattering</div><div> <!> *** WARNING *** <!>  Unknown process: MultipleScattering</div><div><br></div> <br> <br> List of Enabled processes:<div><br></div><div><br></div><div>   ===  Process: PhotoElectric  ===</div><div>Particle: gamma</div><div>    * Model(s):</div><div>       - StandardModel</div><div><br></div><div>   ===  Process: Compton  ===</div><div>Particle: gamma</div><div>    * Model(s):</div><div>       - StandardModel</div><div><br></div><div>   ===  Process: Bremsstrahlung  ===</div><div>Particle: e+</div><div>    * Model(s):</div><div>       - StandardModel</div><div>Particle: e-</div><div>    * Model(s):</div><div>       - StandardModel</div><div><br></div><div>   ===  Process: ElectronIonisation  ===</div><div>Particle: e+</div><div>    * Model(s):</div><div>       - StandardModel</div><div>Particle: e-</div><div>    * Model(s):</div><div>       - StandardModel</div><div><br></div><div>List of initialized processes:</div><div><br></div><div>[Core-0] Initialization of geometry</div><div>G4PhysicalVolumeModel::Validate(): A volume of the same name and</div><div>  copy number ("world_phys", copy 0) still exists and is being used.</div><div>  But it is not the same volume you originally specified</div><div>  in /vis/scene/add/.</div><div>[Core-0] Initialization of physics</div><div>[Core-0] Initialization of actors</div><div><br></div><div>GATE object:        'systems/cylindricalPET'</div> <br> Components:    <div><br></div><div>GATE object:        'systems/cylindricalPET/base'</div><div>Attached to volume: cylindricalPET</div><div>Nb of children:       1</div><div><br></div><div>        GATE object:        'systems/cylindricalPET/rsector'</div><div>        Attached to volume: head</div><div>        Nb of children:       1</div><div><br></div><div>                GATE object:        'systems/cylindricalPET/module'</div><div>                Attached to volume: ---</div><div>                Nb of children:       1</div><div><br></div><div>                        GATE object:        'systems/cylindricalPET/submodule'</div><div>                        Attached to volume: block</div><div>                        Nb of children:       1</div><div><br></div><div>                                GATE object:        'systems/cylindricalPET/crystal'</div><div>                                Attached to volume: crystal</div><div>                                Nb of children:       4</div><div><br></div><div>                                        GATE object:        'systems/cylindricalPET/layer0'</div><div>                                        Attached to volume: LYSO</div><div>                                        Nb of children:       0</div><div><br></div><div>                                        GATE object:        'systems/cylindricalPET/layer1'</div><div>                                        Attached to volume: LSO</div><div>                                        Nb of children:       0</div><div><br></div><div>                                        GATE object:        'systems/cylindricalPET/layer2'</div><div>                                        Attached to volume: ---</div><div>                                        Nb of children:       0</div><div><br></div><div>                                        GATE object:        'systems/cylindricalPET/layer3'</div><div>                                        Attached to volume: ---</div><div>                                        Nb of children:       0</div><div>geometrical design type: 1</div><div>ring diameter: -1.750e+01 cm </div><div>rsector axial pitch: 0.000e+00 fm </div><div>rsector azimuthal pitch: -4.500e+01 degree</div><div>rsector tangential size: 2.000e+00 cm </div><div>rsector axial size: 2.530e+01 cm </div><div>module axial size: 2.530e+01 cm </div><div>module tangential size: 2.000e+00 cm </div> module axial pitch: 0.000e+00 fm <div>module tangential pitch: 0.000e+00 fm </div><div>submodule axial size: 1.900e+00 cm </div><div>submodule tangential size: 2.000e+00 cm </div><div>submodule axial pitch: 1.950e+00 cm </div><div>submodule tangential pitch: 0.000e+00 fm </div><div>crystal radial size: 2.200e+00 mm </div><div>crystal axial size: 2.200e+00 mm </div><div>crystal tangential size: 2.000e+00 cm </div><div>crystal axial pitch: 2.400e+00 mm </div><div>crystal tangential pitch: 0.000e+00 fm </div><div>layer0 radial size: 2.200e+00 mm </div><div>in layer0 interaction length: 1.100e+00 mm </div><div>layer1 radial size: 2.200e+00 mm </div><div>in layer1 interaction length: 1.100e+00 mm </div><div>Axial nb of rsectors: 1</div><div>Azimuthal nb of rsectors: 8</div><div>Axial nb of modules: 1</div><div>Tangential nb of modules: 1</div><div>Axial nb of submodules: 13</div><div>Tangential nb of submodules: 1</div><div>Axial nb of crystals: 8</div><div>Tangential nb of crystals: 1</div><div>Radial nb of layers: 2</div><div> Range Translator</div><div>        Range   0 : imageValue in [    0.000e+00 , 1.500e+03 ]  ---> activity (Bq) 0.000e+00</div><div>        Range   1 : imageValue in [    1.501e+03 , 2.250e+03 ]  ---> activity (Bq) 1.000e+01</div><div>        Range   2 : imageValue in [    2.251e+03 , 2.750e+03 ]  ---> activity (Bq) 2.000e+01</div><div>        Range   3 : imageValue in [    2.751e+03 , 3.848e+03 ]  ---> activity (Bq) 3.000e+01</div><div>GateSourceVoxelImageReader::ReadFile : fileName: CT-IM00120_uint16.h33</div><div>', assuming default BIGENDIAN'LITTLEENDIAN</div><div><br></div><div>'nrecognised type name 'unsigned integer</div><div> Header read from       'CT-IM00120_uint16.h33'</div><div>'Data file name         'm001-IM00120.i33</div><div> Nb of planes:           1</div><div> Nb of pixels per plane: 512 512</div><div> Pixel size:             4.590e-01 4.590e-01</div><div> Slice thickness:        2.179e+00</div><div> Matrix size:            2.350e+02 2.350e+02</div><div> Data type:              </div><div> Data byte order:        BIGENDIAN</div><div><br></div><div><br></div><div>'!ror: Could not open header file 'm001-IM00120.i33</div><div> **** Error: Volume does not exist **** </div><div> Ignoring confine condition</div><div>Source ---------------> voxel_brain</div><div>  ID                  : 0</div><div>  type                : backtoback</div><div>  startTime (s)       : 0.000e+00</div><div>  time (s)            : 0.000e+00</div><div>  forcedUnstable      : 0</div><div>  verboseLevel        : 0</div><div>----------------------- </div><div>  Voxel reader ---------->   number of voxels       : 0</div><div>  total activity (Bq)    : 0.000e+00</div><div>  position  (mm)         : 0.000e+00 0.000e+00 0.000e+00</div><div>  voxel size  (mm)    : 1.000e+00 1.000e+00 1.000e+00</div><div>Source --------------> voxel_brain</div><div>  ID                 : 0</div><div>  type               : backtoback</div><div>  activity (Bq)      : 0.000e+00</div><div>  startTime (s)      : 0.000e+00</div><div>  time (s)           : 0.000e+00</div><div>  forcedUnstable     : 0</div><div>  forcedHalfLife (s) : -6.931e-01</div><div>  verboseLevel       : 0</div><div>  relative to vol    : World</div><div>---------------------- </div><div><br></div><div>    GPS info ----------------> </div><div>    particle                 : gamma</div><div>    SourcePosType            : Point</div><div>    Shape                    : NULL</div><div>    halfx,halfy,halfz  (mm)  : 0.000e+00 0.000e+00 0.000e+00</div><div>    Radius             (mm)  : 0.000e+00</div><div>    CentreCoords       (mm)  : 0.000e+00 0.000e+00 0.000e+00</div><div>    EnergyDisType            : Mono</div><div>    AngleDisType            : iso</div><div>    MinTheta, MaxTheta (deg) : 0.000e+00 1.800e+02</div><div>    MinPhi, MaxPhi     (deg) : 0.000e+00 3.600e+02</div><div>    -------------------------- </div><div><br></div><div>GateSourceMgr::ListSources: List of the sources in the source manager</div><div>Source ---------------> voxel_brain</div><div>  ID                  : 0</div><div>  type                : backtoback</div><div>  startTime (s)       : 0.000e+00</div><div>  time (s)            : 0.000e+00</div><div>  forcedUnstable      : 0</div><div><br></div><div><div> verboseLevel        : 0</div><div>----------------------- </div><div>***** COMMAND NOT FOUND </gate/output/root/setRootdelayFlag 1> *****</div><div><br></div><div>***** Batch is interrupted!! *****</div><div>[Acquisition-0]                                      </div><div>[Acquisition-0] ============= Source initialization =============</div><div>[Acquisition-0]                                      </div><div>[Acquisition-0] ============= Acquisition starts! =============</div><div>[Acquisition-0] Simulation start time = 0 sec</div><div>[Acquisition-0] Simulation end time   = 2.5 sec</div><div>[Acquisition-0] Simulation will have  = 1 run(s)</div><div>[Acquisition-0] </div><div>[Acquisition-0] Slice 0 from 0 to 2.5 s [slice=2.5 s], final stop at 2.5 s.</div><div><br></div><div>PhotoElectric:   for  gamma    SubType= 12</div><div>      ===== EM models for the G4Region  DefaultRegionForTheWorld ======</div><div>       PhotoElectric :  Emin= 0.000e+00 eV    Emax= 1.000e+01 TeV</div><div><br></div><div>Compton:   for  gamma    SubType= 13</div><div>      Lambda tables from 1.000e+02 eV  to 1.000e+01 TeV in 77 bins, spline: 1</div><div>      ===== EM models for the G4Region  DefaultRegionForTheWorld ======</div><div>       Klein-Nishina :  Emin= 0.000e+00 eV    Emax= 1.000e+01 TeV</div><div><br></div><div>Bremsstrahlung:   for  e-    SubType= 3</div><div>      dE/dx and range tables from 1.000e+02 eV  to 1.000e+01 TeV in 77 bins</div><div>      Lambda tables from threshold to 1.000e+01 TeV in 77 bins, spline: 1</div><div>      LPM flag: 1 for E > 1.000e+00 GeV</div><div>      ===== EM models for the G4Region  DefaultRegionForTheWorld ======</div><div>             eBremSB :  Emin= 0.000e+00 eV    Emax= 1.000e+00 GeV   AngularGenUrban</div><div>            eBremLPM :  Emin= 1.000e+00 GeV   Emax= 1.000e+01 TeV   AngularGenUrban</div><div><br></div><div>ElectronIonisation:   for  e-    SubType= 2</div><div>      dE/dx and range tables from 1.000e+02 eV  to 1.000e+01 TeV in 77 bins</div><div>      Lambda tables from threshold to 1.000e+01 TeV in 77 bins, spline: 1</div><div>      finalRange(mm)= 1.000e+00, dRoverRange= 2.000e-01, integral: 1, fluct: 1, linLossLimit= 1.000e-02</div><div>      ===== EM models for the G4Region  DefaultRegionForTheWorld ======</div><div>        MollerBhabha :  Emin= 0.000e+00 eV    Emax= 1.000e+01 TeV</div><div><br></div><div>Bremsstrahlung:   for  e+    SubType= 3</div><div>      dE/dx and range tables from 1.000e+02 eV  to 1.000e+01 TeV in 77 bins</div><div>      Lambda tables from threshold to 1.000e+01 TeV in 77 bins, spline: 1</div><div>      LPM flag: 1 for E > 1.000e+00 GeV</div><div>      ===== EM models for the G4Region  DefaultRegionForTheWorld ======</div><div>             eBremSB :  Emin= 0.000e+00 eV    Emax= 1.000e+00 GeV   AngularGenUrban</div><div>            eBremLPM :  Emin= 1.000e+00 GeV   Emax= 1.000e+01 TeV   AngularGenUrban</div><div><br></div><div>ElectronIonisation:   for  e+    SubType= 2</div><div>      dE/dx and range tables from 1.000e+02 eV  to 1.000e+01 TeV in 77 bins</div><div>      Lambda tables from threshold to 1.000e+01 TeV in 77 bins, spline: 1</div><div>      finalRange(mm)= 1.000e+00, dRoverRange= 2.000e-01, integral: 1, fluct: 1, linLossLimit= 1.000e-02</div><div>      ===== EM models for the G4Region  DefaultRegionForTheWorld ======</div><div>        MollerBhabha :  Emin= 0.000e+00 eV    Emax= 1.000e+01 TeV</div><div>1 event has been kept for refreshing and/or reviewing.</div><div>     User simulation time (sec)   := 4.900e-01</div><div>     Real simulation time (sec)   := 2.860e+00</div><div>     System simulation time (sec) := 6.000e-02</div><div>[Core-0] End of macro hrrt.mac</div><div>GateOutputMgr deleting...</div><div>Graphics systems deleted.</div><div>Visualization Manager deleting...</div><div>WARNING - Attempt to delete the physical volume store while geometry closed !</div><div>WARNING - Attempt to delete the logical volume store while geometry closed !</div><div>WARNING - Attempt to delete the solid store while geometry closed !</div><div>WARNING - Attempt to delete the region store while geometry closed !</div></div><div><br></div><br> <br><br></div>                                           </div></body>
</html>