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: