# Systematic over-correction with scatter correction in TOF list-mode reconstruction **URL:** \ **Category:** General Discussions\ **Created:** [September 18, 2026, 8:31am UTC](https://castor-project.discourse.group/t/systematic-over-correction-with-scatter-correction-in-tof-list-mode-reconstruction/1501 "2026-09-18T08:31:13Z")\ **Posts on this page:** 1\ **Page:** 1
**Author:** ![D.Liang](https://avatars.discourse-cdn.com/v4/letter/d/58f4c7/32.png) [@D.Liang](https://castor-project.discourse.group/u/D.Liang)\ **Post date:** [September 18, 2026, 8:31am UTC](https://castor-project.discourse.group/t/systematic-over-correction-with-scatter-correction-in-tof-list-mode-reconstruction/1501/1 "2026-09-18T08:31:13Z")
Dear CASToR developers and users, I am experiencing a systematic apparent over-correction when applying scatter correction in TOF list-mode PET reconstruction with CASToR. I have already checked several common causes discussed in previous CASToR forum posts, but the problem remains. I am using simulated PET data with only the true and scatter events. Normalization, attenuation, and scatter corrections are applied. For non-TOF reconstruction, the correction works as expected: the contrast recovery coefficients (CRCs) of 37mm ROIs obtained from true+scatter events with scatter correction (85%) are very close to those obtained from true-only reconstruction (90%). However, when TOF information is enabled, the scatter-corrected reconstruction shows substantially higher CRCs, typically around 115–134%. For comparison, the corresponding TOF true-only reconstruction gives a CRC of approximately 90%. Therefore, the CRC values above 100% appear to be specifically associated with the TOF-dependent scatter correction rather than with TOF reconstruction itself. For clarity, I am using TOF list-mode data. The discrete TOF bins described below are used only during scatter estimation and are not reconstruction TOF bins. ### Scatter estimation procedure I first estimate the scatter event rate between detector-module pairs rather than directly estimating the scatter rate for every individual crystal-to-crystal LOR. The estimated scatter rate for each module pair is then distributed uniformly among all crystal-pair LORs connecting the corresponding pair of modules. Therefore, the crystal-level scatter estimate is based on a module-level estimate and is spatially smoothed by construction, rather than directly using noisy Monte Carlo scatter counts for each individual crystal-pair LOR. For the TOF dependence, I divide the TOF range into a number of bins and sort the simulated scatter events according to their TOF values. For each LOR and TOF bin, I calculate the scatter rate density as scatter density = N\_scatter / (acquisition duration × spatial TOF-bin width), where spatial TOF-bin width = c × Δt / 2. The resulting quantity therefore has units of counts/s/mm. This scatter correction is then written into the CASToR list-mode CDF for reconstruction. ### Checks already performed I have already investigated the following possible explanations: 1. Data mode 2. Number of TOF bins used for scatter estimation 3. Normalization and attenuation 4. Acquisition duration 5. TOF sign convention and detector ordering 6. True-only TOF reconstruction ### Main questions At this stage, I would particularly appreciate clarification regarding the expected normalization and units of the scatter term for continuous TOF list-mode reconstruction in CASToR. 1. For TOF list-mode reconstruction, is the scatter value stored for each event expected to represent a continuous scatter-rate density in counts/s/mm? 2. If a scatter rate is estimated within a finite TOF interval, is dividing the scatter rate by the spatial TOF-bin width 3. Is there any dependence on the TOF resolution or TOF projector normalization that needs to be considered when defining the scatter correction term? 4. Could the stronger over-correction observed when increasing the number of scatter-estimation TOF bins from 4 to 10 indicate that an additional TOF-bin-width-dependent scaling is being applied somewhere? 5. Is there a recommended numerical consistency test for verifying the normalization of the TOF scatter term? The most puzzling observation is that the same general normalization, attenuation, and scatter correction framework works well for non-TOF reconstruction, while the large quantitative bias appears only after making the scatter correction TOF-dependent. In addition, the fact that the apparent over-correction becomes stronger when increasing the number of scatter-estimation TOF bins from 4 to 10 makes me particularly suspicious of a TOF-width-dependent normalization issue. If useful, I can provide the CDF header, reconstruction command, TOF resolution, TOF range, scatter-estimation bin widths, scatter fraction, representative scatter correction values. Any clarification regarding the expected physical units and normalization of the TOF list-mode scatter correction would be greatly appreciated. Kinds regards, Da