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
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