def mt_with(self, *cands):
     myP4 = self.p4()
     lepP4 = TLorentzVector()
     for lep in cands: lepP4 += lep.p4()
     # squared transverse mass of system = 
     #     (met.Et + lep.Et)^2 - (met + lep).Pt)^2
     mt = math.sqrt(abs((lepP4.et() + myP4.et())**2 - ((lepP4 + myP4).pt())**2))
     return mt
Exemple #2
0
 def mt_with(self, *cands):
     myP4 = self.p4()
     lepP4 = TLorentzVector()
     for lep in cands:
         lepP4 += lep.p4()
     # squared transverse mass of system =
     #     (met.Et + lep.Et)^2 - (met + lep).Pt)^2
     mt = math.sqrt(
         abs((lepP4.et() + myP4.et())**2 - ((lepP4 + myP4).pt())**2))
     return mt