201 def ProcessHits(self, step, rohist):
202 preStepPoint = step.GetPreStepPoint()
203 track = step.GetTrack()
204 part = track.GetDynamicParticle()
205 pid = part.GetPDGcode()
206 if abs(pid) in leptons :
207 mom = track.GetMomentum()
208 pos = track.GetPosition()
209
210
211 w = track.GetWeight()
212 parentid = int(w)/100000-10000
213 pythiaid = int(w)%100000-10000
214 h['ntuple'].Fill(float(pid), float(mom.x/GeV),float(mom.y/GeV),float(mom.z/GeV),
215 float(pos.x/m),float(pos.y/m),float(pos.z/m),pythiaid,parentid)
216
217
218