Skip to content

nCompton counter not working correctly for low SetCutInRegion values #711

Open
@ghanima86

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

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions