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: f = ROOT.TFile(ofile, "recreate") t.Write() f.Close()