etaMax = abs(ak7Def[1].Rapidity())
mjetReco = (ak7Def[0].M() + ak7Def[1].M()) * 0.5
ptAvg = (ak7Def[0].Perp() + ak7Def[1].Perp()) * 0.5
mjetFinal.append(mjetReco)
ptAsymmetry = (ak7Def[0].Perp() - ak7Def[1].Perp()) / (ak7Def[0].Perp() + ak7Def[1].Perp())
ptThreshold = max( ptAvg * 0.1, vetoPtCut )
njetsAboveThreshold = 0
for ijet in ak7Def :
if ijet.Perp() > ptThreshold :
njetsAboveThreshold += 1
[passTrig,iTrigHist] = trigHelper.passEventData(event, ptAvg )
passEvent = passTrig and njetsAboveThreshold == 2 and ptAsymmetry < 0.3
if not passEvent :
continue
# Here we have two reconstructed jets. Get the
# average jet mass, dijet mass, etaMax, average pt,
# and trigger bin
if options.verbose :
print 'mjjReco = ' + str(mjjReco) + ', etaMax = ' + str(etaMax)