From 502bf13915b34af12ce244c612914fe84a01cb27 Mon Sep 17 00:00:00 2001 From: Yi Chen Date: Mon, 13 Dec 2021 10:41:47 +0100 Subject: [PATCH] Update HiPuRhoProducer.cc One-line bug fix in the jetty tower exclusion part of the code. The bug uses ieta instead of N(ieta) to calculate the number of towers. --- RecoHI/HiJetAlgos/plugins/HiPuRhoProducer.cc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/RecoHI/HiJetAlgos/plugins/HiPuRhoProducer.cc b/RecoHI/HiJetAlgos/plugins/HiPuRhoProducer.cc index 190bc65570a94..c6dd86b735b38 100644 --- a/RecoHI/HiJetAlgos/plugins/HiPuRhoProducer.cc +++ b/RecoHI/HiJetAlgos/plugins/HiPuRhoProducer.cc @@ -464,7 +464,7 @@ void HiPuRhoProducer::calculateOrphanInput(std::vector& orph for (auto const& im : towermap_) { double dr2 = reco::deltaR2(im.eta, im.phi, jet_etaphi.first, jet_etaphi.second); if (dr2 < radiusPU_ * radiusPU_ && !excludedTowers[std::pair(im.ieta, im.iphi)] && - (im.ieta - ntowersWithJets_[im.ieta]) > minimumTowersFraction_ * im.ieta) { + (geomtowers_[im.ieta] - ntowersWithJets_[im.ieta]) > minimumTowersFraction_ * geomtowers_[im.ieta]) { ntowersWithJets_[im.ieta]++; excludedTowers[std::pair(im.ieta, im.iphi)] = 1; }