[Gate-users] Scintillation does not work

Fri Jun 24 11:28:42 CEST 2022

Dear Users,

currently I am facing a lot of problems with the scintillation process. I use Gate9.1, Geant4.10.7 and installed everything with the precompiled Version for the VM. I recompiled Gate with Optical = On and LUT_Davis = On.

When I run my simulation I see in root Optical Data and NumScintillation and I should have produced some Scintillationphotons, but neither do I see them in the Hits file nor during the Simulation.

[cid:4cbf07eb-4d42-4bf4-83b7-eea5801404e0][cid:7069f9d8-397b-433c-a703-341078c136cf]

       0      10   1   0     0     0     0     0     3    14   1.19057513718741776178808e-04  2.910e-04  2.849e+00 -2.135e+00 -2.675e-01 -1.347e+00      22     1     0   1    0    0 Compton NULL NULL
       0      10   1   0     0     0     0     0     3    14   1.19057514463877581450263e-04  3.153e-02  1.678e-01 -2.087e+00 -3.181e-01 -1.500e+00      11     2     1   1    0    0 Transportation NULL NULL
       0      24   1   0     0     0     0     0     6    17   2.96347702829952576019518e-04  2.910e-04  2.963e+00  4.084e-01  3.499e-01 -1.463e+00      22     1     0   1    0    0 Compton NULL NULL
       0      24   1   0     0     0     0     0     6    17   2.96347702830890302478500e-04  4.815e-05  1.439e-04  4.085e-01  3.500e-01 -1.463e+00      11     2     1   1    0    0 Transportation NULL NULL
       0      24   1   0     0     0     0     0     0     7   2.96347703032492383861535e-04  8.383e-02  3.094e-02  4.271e-01  3.681e-01 -1.480e+00      11     2     1   1    0    0 ElectronIonisation NULL NULL
       0      24   1   0     0     0     0     0     0     7   2.96347703032492383861535e-04  0.000e+00  0.000e+00  4.271e-01  3.681e-01 -1.480e+00      11     2     1   1    0    0 Scintillation NULL NULL

Thats a part of my Hits ascii file. In the last row is a scintillationprocess, but somehow the particle is an electron(electron = 11, photon = 22). How does that make sense, I have no scintillation photons, only secondary produced photons by compton interactions.

I've seen the physicslist "optical" a lot but somehow it does not exist anymore.

I use a block consisting of LYSO and a plastic. I've defined their scintillationyield and so on in the materials.xml. When I run the macro, you can see that those parameters are read-in. But somehow theres no scintillation. And I should note that the Singles output is allways empty.

Of course every layer is defined as crystalSD. And I use a gamma source.

Solutions I allready tried and checked:

https://github.com/OpenGATE/Gate/commit/e274e87dafe32868d705131313b92c97f1e7dac4

https://github.com/OpenGATE/Gate/issues/531

https://github.com/OpenGATE/Gate/commit/42e158e7e6de774b655af663fc3ba949bdde3e1e

<https://github.com/OpenGATE/Gate/commit/42e158e7e6de774b655af663fc3ba949bdde3e1e>http://lists.opengatecollaboration.org/pipermail/gate-users/2009-May/003338.html

Thanks in advance and kind regards,

Paul

I run these physics:

/gate/physics/addProcess  PhotoElectric
/gate/physics/processes/PhotoElectric/setModel  PenelopeModel
/gate/physics/addProcess  Compton
/gate/physics/processes/Compton/setModel  PenelopeModel
/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  OpticalAbsorption
/gate/physics/addProcess   PositronAnnihilation
/gate/physics/addProcess  Scintillation
/gate/physics/addProcess  OpticalBoundary
/gate/physics/addProcess  OpticalRayleigh
/gate/physics/addProcess  RadioactiveDecay
/gate/physics/processList Enabled
/gate/physics/processList  Initialized

Surfaces:

/gate/LYSO/surfaces/name             surface1
/gate/LYSO/surfaces/insert             EJ232
/gate/LYSO/surfaces/surface1/setSurface     RoughTeflon_LUT

/gate/EJ232/surfaces/name             surface2
/gate/EJ232/surfaces/insert             LYSO
/gate/EJ232/surfaces/surface2/setSurface     RoughTeflon_LUT

/gate/PMT/surfaces/name             surface3
/gate/PMT/surfaces/insert             LYSO
/gate/PMT/surfaces/surface3/setSurface     Detector_LUT

/gate/PMT/surfaces/name             surface4
/gate/PMT/surfaces/insert             EJ232
/gate/PMT/surfaces/surface4/setSurface     Detector_LUT

Digitizer:

/gate/digitizer/Singles/insert adder
/gate/digitizer/name OpticalSingles
/gate/digitizer/insert singleChain
/gate/digitizer/OpticalSingles/insert opticaladder
/gate/digitizer/Singles/readout/setDepth    5

Data output:

/gate/output/root/enable
/gate/output/root/setFileName             LYSOEJ232PMTTAOpticalAdder
#/gate/output/root/setRootSinglesAdderFlag 1
#/gate/output/root/setRootSinglesReadoutFlag 1
/gate/output/root/setRootHitFlag 1
/gate/output/root/setRootSinglesFlag 1
/gate/output/root/setRootOpticalFlag  1

#/gate/output/tree/addCollection Singles

/gate/output/ascii/enable
/gate/output/ascii/setFileName LYSOEJ232PMTTAOpticalAdder
/gate/output/ascii/setOutFileSinglesFlag 1
/gate/output/ascii/setOutFileOpticalSinglesFlag 1

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20220624/ab9712bd/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pastedImage.png
Type: image/png
Size: 22007 bytes
Desc: pastedImage.png
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20220624/ab9712bd/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pastedImage.png
Type: image/png
Size: 23522 bytes
Desc: pastedImage.png
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20220624/ab9712bd/attachment-0003.png>