<div dir="ltr">Hello Gate users,<div>I am simulating a cylindricalPET system and storing data in ROOT file format. I have found lately Schmidtlein's code to read coincidences from the ROOT TTree and generate corresponding Michelogram. I am encountering very strange behaviour and can't seem to understand whether it is due to the simulation, the code, or whatever else. Trying to keep it simple, I am now simulating a point source in the central axial plane . Still I get strange output from the code. I will now include the parameters that I have changed in the code, according to my PET scanner specifications:</div><div><br></div><div><div>// DLS Parameters</div><div>#define  N_RINGS               45 // DLS: Number of rings</div><div>#define  N_SEG                 89 // DLS: Number of segments (2 N_RINGS - 1)</div><div>#define  N_DET                448 // DLS: Detectors per ring</div><div>#define  S_WIDTH              325 // DLS: Number of radial sinogram bins</div><div>//#define  N_RSEC                28 // DLS: Number of resector</div><div>//#define  N_MODULE               20 // DLS: Number of modules (2x3)</div><div>#define  N_MOD_xy               4 // DLS: Number of tangential modules</div><div>#define  N_MOD_z                5 // DLS: Number of axial modules</div><div>//#define  N_SUBMOD               1 // DLS: Number of submodules (1x1)</div><div>#define  N_SMOD_xy              1 // DLS: Number of tangential submodules</div><div>#define  N_SMOD_z               1 // DLS: Number of axial submodules</div><div>//#define  N_CRYSTAL             36 // DLS: Number of crystals (6x6)</div><div>#define  N_CRY_xy               4 // DLS: Number of tangential crystals</div><div>#define  N_CRY_z                9 // DLS: Number of axial crystals</div><div>#define  MAX_D_RING            1 // DLS: Maximum ring difference</div><div>//#define  N_PLANES             265 // DLS: Total number of sinograms</div><div>//#define  FOV                  550 // DLS: Width of the FOV (mm)</div><div>//#define  SLICES_PER_FOV        35 // DLS: Number of slices per FOV</div><div>#define  USE_OFFSET             0 // DLS: On/Off use of offset</div><div>#define  OFFSET               499 // DLS: Sets initial sinogram angle</div></div><div><br></div><div>I will also attach an image representing the supposed sinogram at the center of Z axis (Z=0) . If anybody has ever encountered similar behaviour and has the knowledge on how to overcome this situation, it would be greatly appreciated. </div><div>Kind regards,</div><div>Vincenzo</div></div>