<html><head><style type="text/css">body{font:12px Arial;margin:3px;overflow-y:auto;overflow-x:auto}p{margin:0px;}blockquote, ol, ul{margin-top:0px;margin-bottom:0px;}</style></head>
<body><div style="display: block; font-family: Arial; font-size: 12px;"><div style="display: block; font-family: Arial; font-size: 12px;">Greetings,<br>
<div style=""><div style="display: block; font-family: Arial; font-size: 12px;"><div style=""><div style="display: block; font-family: Arial; font-size: 12px;">
<br>
I'm trying to simulate a simple geometry to access only primary
interactions. It consists of a block of lead (9cm thick), and I shoot
gamma-ray (511keV, 10Bq for 1 sec) in one direction (Source-Surface
Distance=10cm). I have no phantom. I set electron cut to 1 km, X- and
Delta-Ray cuts to 1 GeV. I activate only photoelectric and Compton
effects (low energy for both).<br>
I finally obtain 14 events.<br>
<br>
I cannot understand why I have so many Compton effects (20 single
diffusions and 7 double) with nCrystalCompton (on the Hits branch)
(attached) whereas I have only 10 LowEnCompton on the processName sheet
(attached). Could someone explain it to me? The GATE display gives me
well 10 Compton.<br>
Moreover, on the Gate display, I can't understand this lines:<br>
<br>
*********************************************************************************************************<br>
* G4Track Information: Particle = gamma, Track ID = 1, Parent ID = 0<br>
*********************************************************************************************************<br>
<br>
Step# X(mm) Y(mm)
Z(mm) KinE(MeV) dE(MeV) StepLeng TrackLeng NextVolume
ProcName<br>
0
0
0 0
0.511
0
0
0 world_P initStep<br>
1 100 -1.22e-14
-6.12e-15
0.511
0 100
100 level1_P Transportation<br>
2 105 -1.29e-14
-6.46e-15 0.349
0.162 5.43
105 level1_P LowEnCompton<br>
3
108 -1.76
4.19
0.148
0 5.38
111 level1_P LowEnCompton<br>
:----- List of 2ndaries - #SpawnInStep=
1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 1 ---------------<br>
:
108 -1.76
4.19
0.202
e-<br>
:----------------------------------------------------------------- EndOf2ndaries Info ---------------<br>
4
108 -1.36
3.05 0
0.0158 1.39
112 level1_P LowEnPhotoElec<br>
:----- List of 2ndaries - #SpawnInStep=
1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 2 ---------------<br>
:
108 -1.36
3.05
0.132
e-<br>
:----------------------------------------------------------------- EndOf2ndaries Info ---------------<br>
<br>
*********************************************************************************************************<br>
* G4Track Information: Particle = e-, Track ID = 3, Parent ID = 1<br>
*********************************************************************************************************<br>
<br>
Step# X(mm) Y(mm)
Z(mm) KinE(MeV) dE(MeV) StepLeng TrackLeng NextVolume
ProcName<br>
0
108 -1.36
3.05
0.132
0
0 0
level1_P initStep<br>
1
108 -1.36
3.05
0 0.132 0.0428
0.0428 level1_P eIoni<br>
<br>
*********************************************************************************************************<br>
* G4Track Information: Particle = e-, Track ID = 2, Parent ID = 1<br>
*********************************************************************************************************<br>
<br>
Step# X(mm) Y(mm)
Z(mm) KinE(MeV) dE(MeV) StepLeng TrackLeng NextVolume
ProcName<br>
0
108 -1.76
4.19
0.202
0
0 0
level1_P initStep<br>
1
108 -1.76
4.19
0 0.202 0.0826
0.0826 level1_P eIoni<br>
>>> Event 3<br>
<br>
Why Gate tracks the electron #2 whereas it comes from a Compton effect? It is the only Compton electron tracked.<br>
<br>
Has anyone explanation for it? Please help me.<br>
<br>
Aurélien Vasseur<br>
<br>
PS: I am working with:<br>
Mac OSX 10.3.8<br>
CLHEP1.8.2.0<br>
GEANT4.6.2.p02<br>
GATE1.0.2-4.6.2<br>
Root3.10</div></div></div>
</div></div></br><p style="margin-top:11px;padding-top:3px;background-image: url(http://mail.lycos.co.uk/Images/Mail/_content/dot.gif);background-repeat: repeat-x;background-position: 0px 0px;"><img src="http://premiummail.caramail.lycos.fr/Images/Mail/_icons/premium.gif"><font color="#CC0000"><b> 300 Mo gratuits sur CaraMail : </b></font><a href="http://secure.caramail.lycos.fr/services/signin/mail.jsp?TARGETCODE=FR_footermail_link" target="_blank">Cliquez ici pour en profiter!</a></div></body></html>