[Gate-users] Question about coincidence events in root file analysis

[升级]

Dear GATE expert,



              First, I got "True" coincidence event from a root file by filtering：
index = np.where((comptonPhantom1 == 0) & (comptonPhantom2 == 0)
& (eventID1 == eventID2)
& ((RayleighPhantom1 == 0) & (RayleighPhantom2 == 0))
                )


              Then, I see True LORs in attachment(LORs.png),  (as below,  red "*" is source_xyz, green line and dots were LOR and positions)
      




          The left subplot is my question!  I can’t understand why this is included in True coincidence events


        Could you do me a favor, please?  I'm very confused.  Is it possible that it is a comptonCrystal or RayleighCrystal scattering?


        Any advice would be greatly appreciated. Thanks in advance.




The following is the content of my macro file： （phantom， physics，source）


#=====================================================
# Phantom
#=====================================================
/gate/world/daughters/name Phantom0 
/gate/world/daughters/insert cylinder
/gate/Phantom0/setMaterial Water
/gate/Phantom0/vis/forceSolid
/gate/Phantom0/vis/setColor blue
/gate/Phantom0/geometry/setRmax   1.25 mm
/gate/Phantom0/geometry/setHeight 70. cm
/gate/Phantom0/placement/setTranslation  0. 0. 0. cm 
/gate/Phantom0/attachPhantomSD


#=====================================================
# PHYSICS
#=====================================================
/gate/physics/addProcess PhotoElectric
/gate/physics/processes/PhotoElectric/setModel StandardModel


/gate/physics/addProcess Compton
/gate/physics/processes/Compton/setModel StandardModel


/gate/physics/addProcess RayleighScattering
/gate/physics/processes/RayleighScattering/setModel PenelopeModel


/gate/physics/addProcess ElectronIonisation
/gate/physics/processes/ElectronIonisation/setModel StandardModel e-
/gate/physics/processes/ElectronIonisation/setModel StandardModel e+


/gate/physics/addProcess Bremsstrahlung
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+


/gate/physics/addProcess PositronAnnihilation


/gate/physics/addProcess eMultipleScattering
/gate/physics/processes/eMultipleScattering/setModel UrbanModel e-
/gate/physics/processes/eMultipleScattering/setModel UrbanModel e+


/gate/physics/processList Enabled
/gate/physics/processList Initialized


/gate/physics/Gamma/SetCutInRegion      crystal 1.0 cm
/gate/physics/Electron/SetCutInRegion   crystal 1.0 cm
/gate/physics/Positron/SetCutInRegion   crystal 1.0 cm



#=====================================================
# Source
#=====================================================
/gate/source/addSource                  TwoGamma0
/gate/source/TwoGamma0/setActivity      0.15e-3 Ci


/gate/source/TwoGamma0/gps/particle     gamma
/gate/source/TwoGamma0/setType          backtoback
/gate/source/TwoGamma0/gps/energytype   Mono
/gate/source/TwoGamma0/gps/monoenergy   511 keV


/gate/source/TwoGamma0/setForcedUnstableFlag true
/gate/source/TwoGamma0/setForcedHalfLife 6588. s


/gate/source/TwoGamma0/gps/type         Volume
/gate/source/TwoGamma0/gps/shape        Cylinder
/gate/source/TwoGamma0/gps/radius       1 mm
/gate/source/TwoGamma0/gps/halfz        35. cm
/gate/source/TwoGamma0/gps/angtype      iso
/gate/source/TwoGamma0/gps/centre       0. 0. 0. cm
/gate/source/TwoGamma0/attachTo         Phantom0
/gate/source/list










My best regards,
ConnorHe

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210824/97b539c0/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: LORs.png
Type: image/png
Size: 79691 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210824/97b539c0/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: LORs.png
Type: image/png
Size: 79691 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210824/97b539c0/attachment-0003.png>