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:
s.mbl1MinDeltaHT = invMassOfLightObjects(bjets3Vec[1], [px1, py1, pz1])
s.mbl2MinDeltaHT = invMassOfLightObjects(bjets3Vec[0], [px2, py2, pz2])
Ellbb = sqrt(bjets3Vec[0][0]**2 + bjets3Vec[0][1]**2 + bjets3Vec[0][2]**2) + sqrt(bjets3Vec[1][0]**2 + bjets3Vec[1][1]**2 + bjets3Vec[1][2]**2) + sqrt(px1**2 + py1**2 + pz1**2) + sqrt(px2**2 + py2**2 + pz2**2)
pxllbb = bjets3Vec[0][0] + bjets3Vec[1][0] + px1 + px2
pyllbb = bjets3Vec[0][1] + bjets3Vec[1][1] + py1 + py2
pzllbb = bjets3Vec[0][2] + bjets3Vec[1][2] + pz1 + pz2
s.mllbb = sqrt(Ellbb**2 - pxllbb**2 - pyllbb**2 - pzllbb**2)
s.minJetsMetIso = minAbsPiMinusDeltaPhi(metPhi, jetPhis)
s.minbJetsMetIso = minAbsDeltaPhi(metPhi, bjetPhis)
# print "mbb", s.mbb, "phibb", s.phibb, "minJetsMetIso", s.minJetsMetIso , "minbJetsMetIso", s.minbJetsMetIso, "mbl", s.mbl, "minLepbJetIso", s.minLepbJetIso, "minLepJetIso", s.minLepJetIso
# print s.deltaPhiLL, s.deltaEtaLL , s.deltaRLL, s.deltaPhibb, s.deltaEtabb, s.deltaRbb
# print "allgoodmuons", allGoodLeptons
# print "s.ptZ", s.ptZ, "s.mLL", s.mLL, "s.jzb",s.jzb
t.Fill()
del elist
else:
print "Zero entries in", bin, sample["name"]
del c
if not small:
