Acceptance, Friday March 1st

[A-A-Abdelhamid]

The acceptance of the 400 GeV 1 ns smuon events is 0.682 in the paper, but I am getting 0.360 (after I redefined the d0 function)

Here is how I calculate d0:

def CalcD0(particle):

  ver= particle.production_vertex.position
  vx= ver.x
  vy= ver.y
  vz= ver.z

  momentum =particle.momentum
  pt = momentum.pt()
  px=momentum.px
  py=momentum.py
  pz=momentum.pz
  v0 = ROOT.TVector3(vx, vy, vz)  # Position vector from origin to production vertex
  p = ROOT.TVector3(px, py, pz)   # Momentum vector
  cross_product = v0.Cross(p)

  d0 = cross_product.Mag() / p.Mag()

  return d0

[A-A-Abdelhamid]

How I graphed the proper lifetime tau

 pos=vertex.position
 Tpos= pos.t
 lifetime= 1E9*(Tpos/(3E11))                                
 proper=lifetime/gamma                            
 hist.Fill(proper)                            

[A-A-Abdelhamid]