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

[David Sarrut]

Hello,

unsure, but please use build-in physic list instead of
"/gate/physics/addProcess" stuff (which is deprecated):
https://opengate.readthedocs.io/en/latest/setting_up_the_physics.html#new-physics-list-mechanism

also try to remove "/gate/Phantom0/attachPhantomSD" to check if those
events come from the phantom (it should be anyway indicated in the root
file and remove from the analysis)

David


On Tue, Aug 24, 2021 at 4:33 AM 升级 <hsjbit at 163.com> wrote:

> 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
>
>
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users



-- 
David Sarrut, Phd
Directeur de recherche CNRS
CREATIS, UMR CNRS 5220, Inserm U1294
Centre de lutte contre le cancer Léon Bérard
28 rue Laënnec, 69373 Lyon cedex 08
Tel : 04 78 78 51 51 / 06 74 72 05 42
http://dsarrut.github.io <http://www.creatis.insa-lyon.fr/~dsarrut>
_________________________________
 "2 + 2 = 5,  for extremely large values of 2"
_________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210825/b50a81d6/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/20210825/b50a81d6/attachment-0001.png>