def alphaT(pt0, phi0, pt1, phi1):
return pt1/sqrt(2 * pt0 * pt1 * (1.0 - cos(deltaPhi(phi0, phi1))) )
chain.GetEntry(ievt)
mindphi_c50 = 3.141593
mindphi_c30 = 3.141593
mindphi_c25 = 3.141593
mindphi_c20 = 3.141593
mindphi_cf50 = 3.141593
mindphi_cf30 = 3.141593
mindphi_cf25 = 3.141593
mindphi_cf20 = 3.141593
mindphi_c20f30 = 3.141593
metdphisig_sumet = METtype1corr.et + 0.1
for ihj in xrange(chain.nhJets):
dphi = abs(deltaPhi(METtype1corr.phi, chain.hJet_phi[ihj]))
if mindphi_c50 > dphi and chain.hJet_ptReg[ihj] > 50. and abs(chain.hJet_eta[ihj]) < 2.5:
mindphi_c50 = dphi
if mindphi_c30 > dphi and chain.hJet_ptReg[ihj] > 30. and abs(chain.hJet_eta[ihj]) < 2.5:
mindphi_c30 = dphi
if mindphi_c25 > dphi and chain.hJet_ptReg[ihj] > 25. and abs(chain.hJet_eta[ihj]) < 2.5:
mindphi_c25 = dphi
if mindphi_c20 > dphi and chain.hJet_ptReg[ihj] > 20. and abs(chain.hJet_eta[ihj]) < 2.5:
mindphi_c20 = dphi
if mindphi_cf50 > dphi and chain.hJet_ptReg[ihj] > 50. and abs(chain.hJet_eta[ihj]) < 5:
mindphi_cf50 = dphi
if mindphi_cf30 > dphi and chain.hJet_ptReg[ihj] > 30. and abs(chain.hJet_eta[ihj]) < 5:
mindphi_cf30 = dphi
if mindphi_cf25 > dphi and chain.hJet_ptReg[ihj] > 25. and abs(chain.hJet_eta[ihj]) < 5:
mindphi_cf25 = dphi
if mindphi_cf20 > dphi and chain.hJet_ptReg[ihj] > 20. and abs(chain.hJet_eta[ihj]) < 5: