s.mLL = sqrt((p1 + p2)*(p1 + p2) - p*p)
s.jzb = getVarValue(c, "mht") - s.ptZ
s.pdg1 = l1["pdg"]
s.pdg2 = l2["pdg"]
s.phi1 = l1["phi"]
s.phi2 = l2["phi"]
s.eta1 = l1["eta"]
s.eta2 = l2["eta"]
s.pt1 = l1["pt"]
s.pt2 = l2["pt"]
s.relIso1 = l1["relIso"]
s.relIso2 = l2["relIso"]
metPhi = getVarValue(c, "rawMetphi")
s.deltaPhiLL = deltaPhi(l1["phi"], l2["phi"])
s.deltaEtaLL = l1["eta"] - l2["eta"]
s.deltaRLL = sqrt(s.deltaPhiLL**2 + s.deltaEtaLL**2)
s.minLepMetIso = minAbsPiMinusDeltaPhi(metPhi,[l1["phi"], l2["phi"]])
jetPhis=[]
bjetPhis=[]
bjetEtas=[]
bjets3Vec = []
bjets = []
for jet in goodJets:
jetPhis.append(jet["phi"])
if jet["btag"]>=0.679:
bjetPhis.append(jet["phi"])
bjetEtas.append(jet["eta"])
px = jet["pt"]*cos(jet["phi"]);
py = jet["pt"]*sin(jet["phi"]);
pz = jet["pt"]*sinh(jet["eta"]);
