Skip to content

nCompton counter not working correctly for low SetCutInRegion values #711

Description

@ghanima86

Describe the bug

When the following setting is used:
/gate/physics/Gamma/SetCutInRegion world 0.1mm
The nCrystalCompton counter doesn't work correctly, neither in Hits or in Singles (rootfiles)

Example:
******************* Hits *********************************************************************************

  • Row * eventID * edep * parentID * trackID * primaryID * PDGEncodi * nCrystalC * processNa *

  •    3 *       296 * 0.4689259 *         0 *         1 *         1 *        22 *         0 *     compt *
    
  •    4 *       296 * 0.1627675 *         0 *         1 *         1 *        22 *         0 *      phot *
    
  •    5 *       296 *                0 *          1 *         2 *         1 *        22 *         0 *      Rayl *
    
  •    6 *       296 * 0.0283064 *         1 *         2 *         1 *        22 *         0 *      phot *
    

When the cut is set to higher values, e.g. 1m, the problem doesn't occur and counting seems correct

Desktop (please complete the following information):

  • OS: EL9
  • Gate version or commit hash: tested in v9.2 and v9.4
  • Geant4 version: 11
  • Root version: views LCG_106:x86_64-el9-gcc13-opt (from cvmfs)
  • compiler version: cmake 3.27.5

Minimal example
attached
bug_minex9p4.txt

Expected behavior
nCrystalCompton should increase when a compton interaction occurs in the crystal

Additional context
Possible correlated with issue#539

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions