bjetEtas.append(jet["eta"]) px = jet["pt"]*cos(jet["phi"]); py = jet["pt"]*sin(jet["phi"]); pz = jet["pt"]*sinh(jet["eta"]); bjets3Vec.append([px, py, pz]) bjets.append(jet) if len(bjets)>0: s.btag0pt = bjets[0]["pt"] if len(bjets)>1: s.btagpt = bjets[1]["pt"] if len(bjets3Vec)>0: s.mbl = invMassOfLightObjects(bjets3Vec[0], [px1, py1, pz1]) #Mass of the hardest b-jet and the hardes lepton s.minLepbJetIso = min(minAbsDeltaPhi(l1["phi"], bjetPhis) , minAbsDeltaPhi(l2["phi"], bjetPhis)) s.minLepJetIso = min(minAbsDeltaPhi(l1["phi"], jetPhis) , minAbsDeltaPhi(l2["phi"], jetPhis)) if len(bjets3Vec)>1: s.mbb = invMassOfLightObjects(bjets3Vec[0], bjets3Vec[1]) s.deltaPhibb = deltaPhi(bjetPhis[0], bjetPhis[1]) s.deltaEtabb = bjetEtas[0] - bjetEtas[1] s.deltaRbb = sqrt(s.deltaPhibb**2 + s.deltaEtabb**2) pairL1B0 = True deltaHT1 = abs((s.pt1 + sqrt(bjets3Vec[0][0]**2 + bjets3Vec[0][1]**2)) - (s.pt2 + sqrt(bjets3Vec[1][0]**2 + bjets3Vec[1][1]**2))) deltaHT2 = abs((s.pt2 + sqrt(bjets3Vec[0][0]**2 + bjets3Vec[0][1]**2)) - (s.pt1 + sqrt(bjets3Vec[1][0]**2 + bjets3Vec[1][1]**2))) s.minDeltaHT = deltaHT1 if deltaHT1>deltaHT2: pairL1B0 = False s.minDeltaHT = deltaHT2 if pairL1B0: s.mbl1MinDeltaHT = invMassOfLightObjects(bjets3Vec[0], [px1, py1, pz1]) s.mbl2MinDeltaHT = invMassOfLightObjects(bjets3Vec[1], [px2, py2, pz2]) else: