Double electrons from Compton process

Hello,

Can someone, please, help me understand this following simulated event?

I have a photon generating two electrons (at the same initial position and time) via a Compton process.
One electron doesn’t move in space, lives for a short time and has low initial kinetic energy, while the second electron moves longer in space and time and has larger energy.

--- Printing MCTrack information for track (0x7f1da05eecb0) ----------------------
Particle type (PDG ID):      22
Production process:      annihil (G4 process type: 2)
Production in G4Volume:  water_block_phys
Initial position:          0.420922 mm | -0.529835 mm | -0.052366 mm
Final position:            -15.9854 mm |   6.48652 mm |   56.0835 mm
Initial time:            0.00412973 ns
Final time:                0.278189 ns
Initial kinetic energy:    0.510999 MeV | Final kinetic energy:      0.150391 MeV   
Initial total energy:      0.510999 MeV | Final total energy:        0.150391 MeV   
Linked parent: 0x7f1da05eeae0
----------------------------------------------------------------------------------


--- Printing MCTrack information for track (0x7f1da05eed98) ----------------------
Particle type (PDG ID):      11
Production process:       compt (G4 process type: 2)
Production in G4Volume:  water_block_phys
Initial position:          -7.74191 mm |  -9.23105 mm |  -7.92486 mm
Final position:            -7.74191 mm |  -9.23105 mm |  -7.92486 mm
Initial time:             0.0566245 ns
Final time:               0.0566252 ns
Initial kinetic energy:  0.00050894 MeV | Final kinetic energy:             0 MeV   
Initial total energy:      0.511508 MeV | Final total energy:        0.510999 MeV   
Linked parent: 0x7f1da05eecb0
----------------------------------------------------------------------------------


--- Printing MCTrack information for track (0x7f1da05eee80) ----------------------
Particle type (PDG ID):      11
Production process:       compt (G4 process type: 2)
Production in G4Volume:  water_block_phys
Initial position:          -7.74191 mm |  -9.23105 mm |  -7.92486 mm
Final position:            -7.72954 mm |  -9.24719 mm |  -8.14298 mm
Initial time:             0.0566245 ns
Final time:               0.0584198 ns
Initial kinetic energy:     0.18219** MeV | Final kinetic energy:             0 MeV   
Initial total energy:      0.693189 MeV | Final total energy:        0.510999 MeV   
Linked parent: 0x7f1da05eecb0
----------------------------------------------------------------------------------

Both electrons are from the same photon.
Is this an event that electrons of different orbitals in the same atom are being excited?

Thank you for the help!
Best,
Mateus

Ciao,

my impression is that one electron is the Compton electron, directly kicked out by the incoming gamma, while the other is an Auger electron coming from the filling of the atomic vacancy… Which model are you using for the Compton scattering?

Ciao,
Luciano

Ciao Luciano,
Thank you for the quick reply :smile:

I’m using FTFP_BERT_LIV…
Would you indicate another physics list?
(I am simulating a PET scanner using silicon pixel sensors for the detection layers)

Thanks!
Mateus

Hello,

please, use UI command “/tracking/verbose 2” - history of simulation will allow to understand what happens. Log file may be long if you will require many events but you may limit the output by this particular event.

VI

Thank you for your tip!
Just to let you know, I will try to implement this as soon as possible (I am using Allpix-squared for the Geant4 simulation, and I will have to hard-code this UI command in the framework code).
I will post the result here as soon as I have it.
Cheers,
M

Hello,

This is what I got:

|14:58:45.793|    (INFO) (Event 1) [Geant4] *********************************************************************************************************
|14:58:45.793|    (INFO) (Event 1) [Geant4] * G4Track Information:   Particle = gamma,   Track ID = 1,   Parent ID = 0
|14:58:45.793|    (INFO) (Event 1) [Geant4] *********************************************************************************************************
|14:58:45.793|    (INFO) (Event 1) [Geant4]
|14:58:45.793|    (INFO) (Event 1) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
|14:58:45.793|    (INFO) (Event 1) [Geant4]     0        0     12.8    -15.7     0.511        0        0         0       World initStep
|14:58:45.793|    (INFO) (Event 1) [Geant4]     1    0.348     11.6    -18.6     0.511        0     3.11      3.11 chip_Plane_1_1_a_phys Transportation
|14:58:45.793|    (INFO) (Event 1) [Geant4]     2     0.35     11.6    -18.6     0.511        0   0.0196      3.13 sensor_Plane_1_1_a_phys Transportation
|14:58:45.793|    (INFO) (Event 1) [Geant4]     3     0.38     11.5    -18.8     0.511        0    0.272      3.41 Support_1_a_phys Transportation
|14:58:45.793|    (INFO) (Event 1) [Geant4]     4    0.395     11.5    -18.9     0.511        0    0.131      3.54 chip_Plane_2_1_a_phys Transportation
|14:58:45.793|    (INFO) (Event 1) [Geant4]     5    0.397     11.4      -19     0.511        0   0.0196      3.56 sensor_Plane_2_1_a_phys Transportation
|14:58:45.794|    (INFO) (Event 1) [Geant4]     6    0.427     11.3    -19.2     0.511        0    0.272      3.83 Support_2_a_phys Transportation
|14:58:45.794|    (INFO) (Event 1) [Geant4]     7    0.442     11.3    -19.3     0.511        0    0.131      3.96 chip_Plane_3_1_a_phys Transportation
|14:58:45.794|    (INFO) (Event 1) [Geant4]     8    0.444     11.3    -19.4     0.511        0   0.0196      3.98 sensor_Plane_3_1_a_phys Transportation
|14:58:45.794|    (INFO) (Event 1) [Geant4]     9    0.475     11.2    -19.6     0.511        0    0.272      4.25 Support_3_a_phys Transportation
|14:58:45.794|    (INFO) (Event 1) [Geant4]    10    0.489     11.1    -19.7     0.511        0    0.131      4.38 chip_Plane_4_1_a_phys Transportation
|14:58:45.794|    (INFO) (Event 1) [Geant4]    11    0.491     11.1    -19.7     0.511        0   0.0196       4.4 sensor_Plane_4_1_a_phys Transportation
|14:58:45.794|    (INFO) (Event 1) [Geant4]    12    0.522       11      -20     0.511        0    0.272      4.67 Support_4_a_phys Transportation
|14:58:45.794|    (INFO) (Event 1) [Geant4]    13    0.536       11    -20.1     0.511        0    0.131       4.8 chip_Plane_5_1_a_phys Transportation
|14:58:45.794|    (INFO) (Event 1) [Geant4]    14    0.539       11    -20.1     0.511        0   0.0196      4.82 sensor_Plane_5_1_a_phys Transportation
      :              :        :       :         :      :          :       :        :          :      :          :         :
|14:58:45.796|    (INFO) (Event 1) [Geant4]    71     1.44      7.9    -27.5     0.511        0   0.0196      12.9 sensor_Plane_24_1_a_phys Transportation
|14:58:45.796|    (INFO) (Event 1) [Geant4]    72     1.45     7.84    -27.6     0.464 4.09e-05    0.161        13 sensor_Plane_24_1_a_phys compt
|14:58:45.796|    (INFO) (Event 1) [Geant4]     :----- List of 2ndaries - #SpawnInStep=  3(Rest= 0,Along= 0,Post= 3), #SpawnTotal=  3 ---------------
|14:58:45.796|    (INFO) (Event 1) [Geant4]     :      1.45      7.84     -27.6    0.0471                 e-
|14:58:45.796|    (INFO) (Event 1) [Geant4]     :      1.45      7.84     -27.6  2.99e-05                 e-
|14:58:45.796|    (INFO) (Event 1) [Geant4]     :      1.45      7.84     -27.6  8.07e-05                 e-
|14:58:45.796|    (INFO) (Event 1) [Geant4]     :----------------------------------------------------------------- EndOf2ndaries Info ---------------
|14:58:45.796|    (INFO) (Event 1) [Geant4]    73     1.42      7.8    -27.7     0.464        0    0.116      13.1 Support_24_a_phys Transportation
|14:58:45.796|    (INFO) (Event 1) [Geant4]    74     1.37     7.75    -27.9     0.464        0    0.135      13.3 chip_Plane_25_1_a_phys Transportation
|14:58:45.796|    (INFO) (Event 1) [Geant4]    75     1.37     7.74    -27.9     0.464        0   0.0203      13.3 sensor_Plane_25_1_a_phys Transportation
      :              :        :       :         :      :          :       :        :          :      :          :         :
|14:58:45.800|    (INFO) (Event 1) [Geant4]   179    -3.46     2.56    -41.4     0.464        0    0.135      28.6 chip_Plane_60_1_a_phys Transportation
|14:58:45.800|    (INFO) (Event 1) [Geant4]   180    -3.46     2.55    -41.5     0.464        0   0.0203      28.6 sensor_Plane_60_1_a_phys Transportation
|14:58:45.800|    (INFO) (Event 1) [Geant4]   181    -3.55     2.46    -41.7     0.464        0    0.282      28.9 Support_60_a_phys Transportation
|14:58:45.800|    (INFO) (Event 1) [Geant4]   182     -3.6     2.41    -41.8     0.464        0    0.135        29       World Transportation
|14:58:45.800|    (INFO) (Event 1) [Geant4]   183    -9.95    -4.42    -59.7     0.464        0     20.2      49.2  OutOfWorld Transportation
|14:58:45.801|    (INFO) (Event 1) [Geant4]
|14:58:45.801|    (INFO) (Event 1) [Geant4] *********************************************************************************************************
|14:58:45.801|    (INFO) (Event 1) [Geant4] * G4Track Information:   Particle = e-,   Track ID = 4,   Parent ID = 1
|14:58:45.801|    (INFO) (Event 1) [Geant4] *********************************************************************************************************
|14:58:45.801|    (INFO) (Event 1) [Geant4]
|14:58:45.801|    (INFO) (Event 1) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
|14:58:45.801|    (INFO) (Event 1) [Geant4]     0     1.45     7.84    -27.6  8.07e-05        0        0         0 sensor_Plane_24_1_a_phys initStep
|14:58:45.801|    (INFO) (Event 1) [Geant4]     1     1.45     7.84    -27.6         0 8.07e-05 2.53e-07  2.53e-07 sensor_Plane_24_1_a_phys msc
|14:58:45.801|    (INFO) (Event 1) [Geant4]
|14:58:45.801|    (INFO) (Event 1) [Geant4] *********************************************************************************************************
|14:58:45.801|    (INFO) (Event 1) [Geant4] * G4Track Information:   Particle = e-,   Track ID = 3,   Parent ID = 1
|14:58:45.801|    (INFO) (Event 1) [Geant4] *********************************************************************************************************
|14:58:45.801|    (INFO) (Event 1) [Geant4]
|14:58:45.801|    (INFO) (Event 1) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
|14:58:45.801|    (INFO) (Event 1) [Geant4]     0     1.45     7.84    -27.6  2.99e-05        0        0         0 sensor_Plane_24_1_a_phys initStep
|14:58:45.801|    (INFO) (Event 1) [Geant4]     1     1.45     7.84    -27.6         0 2.99e-05 1.42e-07  1.42e-07 sensor_Plane_24_1_a_phys msc
|14:58:45.801|    (INFO) (Event 1) [Geant4]
|14:58:45.801|    (INFO) (Event 1) [Geant4] *********************************************************************************************************
|14:58:45.801|    (INFO) (Event 1) [Geant4] * G4Track Information:   Particle = e-,   Track ID = 2,   Parent ID = 1
|14:58:45.801|    (INFO) (Event 1) [Geant4] *********************************************************************************************************
|14:58:45.801|    (INFO) (Event 1) [Geant4]
|14:58:45.801|    (INFO) (Event 1) [Geant4] Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
|14:58:45.801|    (INFO) (Event 1) [Geant4]     0     1.45     7.84    -27.6    0.0471        0        0         0 sensor_Plane_24_1_a_phys initStep
|14:58:45.801|    (INFO) (Event 1) [Geant4]     1     1.46     7.84    -27.6    0.0425  0.00464  0.00171   0.00171 sensor_Plane_24_1_a_phys msc
|14:58:45.801|    (INFO) (Event 1) [Geant4]     2     1.46     7.84    -27.6    0.0216   0.0209   0.0151    0.0168 sensor_Plane_24_1_a_phys eIoni
|14:58:45.801|    (INFO) (Event 1) [Geant4]     3     1.46     7.84    -27.7         0   0.0216  0.00548    0.0223 sensor_Plane_24_1_a_phys eIoni

Does it say something?
Thanks for the help!

Best,
Mateus