JuliaHEP / LorentzVectorHEP.jl

x, y, z, t and pt, eta, phi, M Lorentz Vectors as used in collider experiments at CERN and other places

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Incorrect transverse mass `mt2` definition

Moelf opened this issue · comments

mt2(lv::LorentzVector) = lv.t^2 - lv.z^2
mt(lv::LorentzVector) = mt2(lv)<0 ? -sqrt(-mt2(lv)) : sqrt(mt2(lv))

idk why we didn't do the naive px^2 + py^2, because that seems to be correct

here our mt^2 is effectively

E^2 - pz^2 = m^2 + px^2 + py^2 + pz^2 - pz^2

which only gives me correct plot if particle is massless

I think we're missing a 2-particle version of mt2 and mt because that's apparently different from mt(vec1 + vec2), makes sense