Dear experts,
I am working on the interaction between 0.511 photon and LYSO crystals. When using mgdraw.f
to dump the information of interactions and energy deposition, the results shows like:
Pri. 10 Gen. 2 Eng.Dep. 211 LocalDe. 1.60 keV in Reg. CST5 0.000000 0.000000 0.106425
Pri. 10 Gen. 1 Interac. 7 Comp. 511.00 keV in Reg. CST5 0.000000 0.000000 0.106425
secondary 7 7.65 keV
secondary 3 58.26 keV
secondary 7 443.50 keV
Pri. 10 Gen. 3 Eng.Dep. 211 LocalDe. 2.26 keV in Reg. CST5 0.001348 0.002773 0.102942
Pri. 10 Gen. 2 Interac. 7 Ph.e. 7.65 keV in Reg. CST5 0.001348 0.002773 0.102942
secondary 3 5.39 keV
Pri. 10 Gen. 3 Eng.Dep. 3 Bel.Th.2 5.39 keV in Reg. CST5 0.001348 0.002773 0.102942
Pri. 10 Gen. 2 Eng.Dep. 3 Bel.Th.2 9.93 keV in Reg. CST5 -0.000013 0.000181 0.107021
Pri. 10 Gen. 2 Eng.Dep. 7 Escape 443.50 keV in Reg. BLKBODY ******************************
Pri. 12 Gen. 1 Interac. 7 Rayl. 511.00 keV in Reg. CST5 0.000000 0.000000 0.129740
secondary 7 511.00 keV
Pri. 12 Gen. 2 Eng.Dep. 7 Escape 511.00 keV in Reg. BLKBODY ******************************
I really don’t understand the interaction and energy deposition caused by primary 10 (generation 1). My question is:
- what causes the local energy deposition (ICODE=20 and value=1.6 keV) in line1, as the first interaction of the photon is Compton scattering?
- why the Compton scattering reaction generated 3 secondaries (2 photon and 1 electron)?
- the local energy deposition also occurred in the photoelectric process, is it because of the characteristic X ray or Auger electron emitted by the ionized atom, and their energy is lower than the production or transportation threshold?
Besides, the threshold for e+/e- and photon are 1e-5 and 1e-6, respectively:
And there are the input file, mgdraw.f and part of the dumped data.
inters.inp (3.9 KB)
mgdraw.f (9.5 KB)
inters001_Dump.dat (9.3 KB)
Best regards,
Shan