def pixelize(jet_csts, edges, cutoff=0.1):
"""Return eta-phi histogram of transverse energy deposits.
Optionally set all instensities below cutoff to zero.
"""
image, _, _ = np.histogram2d(
jet_csts['eta'], jet_csts['phi'],
bins=(edges[0], edges[1]),
weights=jet_csts['ET'] * (jet_csts['ET'] > cutoff))
return image
评论列表
文章目录