CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
HcalHaloAlgo Class Reference

#include <HcalHaloAlgo.h>

Public Member Functions

reco::HcalHaloData Calculate (const CaloGeometry &TheCaloGeometry, edm::Handle< HBHERecHitCollection > &TheHBHERecHits, edm::Handle< CaloTowerCollection > &TheCaloTowers, edm::Handle< EBRecHitCollection > &TheEBRecHits, edm::Handle< EERecHitCollection > &TheEERecHits, const edm::EventSetup &TheSetup)
 
reco::HcalHaloData Calculate (const CaloGeometry &TheCaloGeometry, edm::Handle< HBHERecHitCollection > &TheHBHERecHits, edm::Handle< EBRecHitCollection > &TheEBRecHits, edm::Handle< EERecHitCollection > &TheEERecHits, const edm::EventSetup &TheSetup)
 
std::vector< reco::HaloClusterCandidateHCAL > GetHaloClusterCandidateHB (edm::Handle< EcalRecHitCollection > &ebrechitcoll, edm::Handle< HBHERecHitCollection > &hbherechitcoll, float et_thresh_seedrh)
 
std::vector< reco::HaloClusterCandidateHCAL > GetHaloClusterCandidateHE (edm::Handle< EcalRecHitCollection > &eerechitcoll, edm::Handle< HBHERecHitCollection > &hbherechitcoll, float et_thresh_seedrh)
 
float GetHBRecHitEnergyThreshold ()
 
float GetHERecHitEnergyThreshold ()
 
float GetPhiWedgeEnergyThreshold ()
 
int GetPhiWedgeNHitsThreshold ()
 
bool HBClusterShapeandTimeStudy (reco::HaloClusterCandidateHCAL hcand, bool ishlt)
 
 HcalHaloAlgo ()
 
bool HEClusterShapeandTimeStudy (reco::HaloClusterCandidateHCAL hcand, bool ishlt)
 
void SetPhiWedgeEnergyThreshold (float SumE)
 
void SetPhiWedgeNHitsThreshold (int nhits)
 
void SetPhiWedgeThresholds (float SumE, int nhits)
 
void SetRecHitEnergyThresholds (float HB, float HE)
 
 ~HcalHaloAlgo ()
 

Private Member Functions

math::XYZPoint getPosition (const DetId &id, reco::Vertex::Point vtx)
 

Private Attributes

const CaloGeometry * geo_
 
float HBRecHitEnergyThreshold
 
float HERecHitEnergyThreshold
 
const HcalGeometry * hgeo_
 
int NHitsThreshold
 
float SumEnergyThreshold
 

Detailed Description

Definition at line 39 of file HcalHaloAlgo.h.

Constructor & Destructor Documentation

◆ HcalHaloAlgo()

HcalHaloAlgo::HcalHaloAlgo ( )

Definition at line 23 of file HcalHaloAlgo.cc.

23  : geo_(nullptr), hgeo_(nullptr) {
24  HBRecHitEnergyThreshold = 0.;
25  HERecHitEnergyThreshold = 0.;
26  SumEnergyThreshold = 0.;
27  NHitsThreshold = 0;
28 }

References HBRecHitEnergyThreshold, HERecHitEnergyThreshold, NHitsThreshold, and SumEnergyThreshold.

◆ ~HcalHaloAlgo()

HcalHaloAlgo::~HcalHaloAlgo ( )
inline

Definition at line 44 of file HcalHaloAlgo.h.

44 {}

Member Function Documentation

◆ Calculate() [1/2]

HcalHaloData HcalHaloAlgo::Calculate ( const CaloGeometry &  TheCaloGeometry,
edm::Handle< HBHERecHitCollection > &  TheHBHERecHits,
edm::Handle< CaloTowerCollection > &  TheCaloTowers,
edm::Handle< EBRecHitCollection > &  TheEBRecHits,
edm::Handle< EERecHitCollection > &  TheEERecHits,
const edm::EventSetup &  TheSetup 
)

Definition at line 39 of file HcalHaloAlgo.cc.

44  {
45  HcalHaloData TheHcalHaloData;
46  // ieta overlap geometrically w/ HB
47  const int iEtaOverlap = 22;
48  const int nPhiMax = 73;
49  // Store Energy sum of rechits as a function of iPhi (iPhi goes from 1 to 72)
50  float SumE[nPhiMax];
51  // Store Number of rechits as a function of iPhi
52  int NumHits[nPhiMax];
53  // Store minimum time of rechit as a function of iPhi
54  float MinTimeHits[nPhiMax];
55  // Store maximum time of rechit as a function of iPhi
56  float MaxTimeHits[nPhiMax];
57  for (unsigned int i = 0; i < nPhiMax; i++) {
58  SumE[i] = 0;
59  NumHits[i] = 0;
60  MinTimeHits[i] = 0.;
61  MaxTimeHits[i] = 0.;
62  }
63 
64  for (const auto& hit : (*TheHBHERecHits)) {
65  HcalDetId id = HcalDetId(hit.id());
66  switch (id.subdet()) {
67  case HcalBarrel:
68  if (hit.energy() < HBRecHitEnergyThreshold)
69  continue;
70  break;
71  case HcalEndcap:
72  if (hit.energy() < HERecHitEnergyThreshold)
73  continue;
74  break;
75  default:
76  continue;
77  }
78 
79  int iEta = id.ieta();
80  int iPhi = id.iphi();
81  if (iPhi < nPhiMax && std::abs(iEta) <= iEtaOverlap) {
82  SumE[iPhi] += hit.energy();
83  NumHits[iPhi]++;
84 
85  float time = hit.time();
86  MinTimeHits[iPhi] = time < MinTimeHits[iPhi] ? time : MinTimeHits[iPhi];
87  MaxTimeHits[iPhi] = time > MaxTimeHits[iPhi] ? time : MaxTimeHits[iPhi];
88  }
89  }
90 
91  for (int iPhi = 1; iPhi < nPhiMax; iPhi++) {
92  if (SumE[iPhi] >= SumEnergyThreshold && NumHits[iPhi] > NHitsThreshold) {
93  // Build PhiWedge and store to HcalHaloData if energy or #hits pass thresholds
94  PhiWedge wedge(SumE[iPhi], iPhi, NumHits[iPhi], MinTimeHits[iPhi], MaxTimeHits[iPhi]);
95 
96  // Loop over rechits again to calculate direction based on timing info
97  std::vector<const HBHERecHit*> Hits;
98  for (const auto& hit : (*TheHBHERecHits)) {
99  HcalDetId id = HcalDetId(hit.id());
100  if (id.iphi() != iPhi)
101  continue;
102  if (std::abs(id.ieta()) > iEtaOverlap)
103  continue; // has to overlap geometrically w/ HB
104  switch (id.subdet()) {
105  case HcalBarrel:
106  if (hit.energy() < HBRecHitEnergyThreshold)
107  continue;
108  break;
109  case HcalEndcap:
110  if (hit.energy() < HERecHitEnergyThreshold)
111  continue;
112  break;
113  default:
114  continue;
115  }
116  Hits.push_back(&(hit));
117  }
118 
119  std::sort(Hits.begin(), Hits.end(), CompareTime);
120  float MinusToPlus = 0.;
121  float PlusToMinus = 0.;
122  for (unsigned int i = 0; i < Hits.size(); i++) {
123  HcalDetId id_i = HcalDetId(Hits[i]->id());
124  int ieta_i = id_i.ieta();
125  for (unsigned int j = (i + 1); j < Hits.size(); j++) {
126  HcalDetId id_j = HcalDetId(Hits[j]->id());
127  int ieta_j = id_j.ieta();
128  if (ieta_i > ieta_j)
129  PlusToMinus += std::abs(ieta_i - ieta_j);
130  else
131  MinusToPlus += std::abs(ieta_i - ieta_j);
132  }
133  }
134  float PlusZOriginConfidence = (PlusToMinus + MinusToPlus) ? PlusToMinus / (PlusToMinus + MinusToPlus) : -1.;
135  wedge.SetPlusZOriginConfidence(PlusZOriginConfidence);
136  TheHcalHaloData.GetPhiWedges().push_back(wedge);
137  }
138  }
139 
140  // Don't use HF.
141  int maxAbsIEta = 29;
142 
143  std::map<int, float> iPhiHadEtMap;
144  std::vector<const CaloTower*> sortedCaloTowers;
145  for (const auto& tower : (*TheCaloTowers)) {
146  if (std::abs(tower.ieta()) > maxAbsIEta)
147  continue;
148 
149  int iPhi = tower.iphi();
150  if (!iPhiHadEtMap.count(iPhi))
151  iPhiHadEtMap[iPhi] = 0.0;
152  iPhiHadEtMap[iPhi] += tower.hadEt();
153 
154  if (tower.numProblematicHcalCells() > 0)
155  sortedCaloTowers.push_back(&(tower));
156  }
157 
158  // Sort towers such that lowest iphi and ieta are first, highest last, and towers
159  // with same iphi value are consecutive. Then we can do everything else in one loop.
160  std::sort(sortedCaloTowers.begin(), sortedCaloTowers.end(), CompareTowers);
161 
162  HaloTowerStrip strip;
163 
164  int prevIEta = -99, prevIPhi = -99;
165  float prevHadEt = 0.;
166  float prevEmEt = 0.;
167  std::pair<uint8_t, CaloTowerDetId> prevPair, towerPair;
168  bool wasContiguous = true;
169 
170  // Loop through and store a vector of pairs (problematicCells, DetId) for each contiguous strip we find
171  for (unsigned int i = 0; i < sortedCaloTowers.size(); i++) {
172  const CaloTower* tower = sortedCaloTowers[i];
173 
174  towerPair = std::make_pair((uint8_t)tower->numProblematicHcalCells(), tower->id());
175 
176  bool newIPhi = tower->iphi() != prevIPhi;
177  bool isContiguous = tower->ieta() == 1 ? tower->ieta() - 2 == prevIEta : tower->ieta() - 1 == prevIEta;
178 
179  isContiguous = isContiguous || (tower->ieta() == -maxAbsIEta);
180  if (newIPhi)
181  isContiguous = false;
182 
183  if (!wasContiguous && isContiguous) {
184  strip.cellTowerIds.push_back(prevPair);
185  strip.cellTowerIds.push_back(towerPair);
186  strip.hadEt += prevHadEt + tower->hadEt();
187  strip.emEt += prevEmEt + tower->emEt();
188  }
189 
190  if (wasContiguous && isContiguous) {
191  strip.cellTowerIds.push_back(towerPair);
192  strip.hadEt += tower->hadEt();
193  strip.emEt += tower->emEt();
194  }
195 
196  if ((wasContiguous && !isContiguous) || i == sortedCaloTowers.size() - 1) { //ended the strip, so flush it
197 
198  if (strip.cellTowerIds.size() > 3) {
199  int iPhi = strip.cellTowerIds.at(0).second.iphi();
200  int iPhiLower = (iPhi == 1) ? 72 : iPhi - 1;
201  int iPhiUpper = (iPhi == 72) ? 1 : iPhi + 1;
202 
203  float energyRatio = 0.0;
204  if (iPhiHadEtMap.count(iPhiLower))
205  energyRatio += iPhiHadEtMap[iPhiLower];
206  if (iPhiHadEtMap.count(iPhiUpper))
207  energyRatio += iPhiHadEtMap[iPhiUpper];
208  iPhiHadEtMap[iPhi] = std::max(iPhiHadEtMap[iPhi], 0.001F);
209 
210  energyRatio /= iPhiHadEtMap[iPhi];
211  strip.energyRatio = energyRatio;
212 
213  TheHcalHaloData.getProblematicStrips().push_back(strip);
214  }
215  strip = HaloTowerStrip();
216  }
217 
218  wasContiguous = isContiguous;
219  prevPair = towerPair;
220  prevEmEt = tower->emEt();
221  prevIPhi = tower->iphi();
222  prevIEta = tower->ieta();
223  prevHadEt = tower->hadEt();
224  }
225 
226  edm::ESHandle<CaloGeometry> pGeo;
227  TheSetup.get<CaloGeometryRecord>().get(pGeo);
228  geo_ = pGeo.product();
229  hgeo_ = dynamic_cast<const HcalGeometry*>(geo_->getSubdetectorGeometry(DetId::Hcal, 1));
230 
231  //Halo cluster building:
232  //Various clusters are built, depending on the subdetector.
233  //In barrel, one looks for deposits narrow in phi.
234  //In endcaps, one looks for localized deposits (dr condition in EE where r =sqrt(dphi*dphi+deta*deta)
235  //E/H condition is also applied.
236  //The halo cluster building step targets a large efficiency (ideally >99%) for beam halo deposits.
237  //These clusters are used as input for the halo pattern finding methods in HcalHaloAlgo and for the CSC-calo matching methods in GlobalHaloAlgo.
238 
239  //Et threshold hardcoded for now. Might one to get it from config
240 
241  std::vector<HaloClusterCandidateHCAL> haloclustercands_HB;
242  haloclustercands_HB = GetHaloClusterCandidateHB(TheEBRecHits, TheHBHERecHits, 5);
243 
244  std::vector<HaloClusterCandidateHCAL> haloclustercands_HE;
245  haloclustercands_HE = GetHaloClusterCandidateHE(TheEERecHits, TheHBHERecHits, 10);
246 
247  TheHcalHaloData.setHaloClusterCandidatesHB(haloclustercands_HB);
248  TheHcalHaloData.setHaloClusterCandidatesHE(haloclustercands_HE);
249 
250  return TheHcalHaloData;
251 }

References funct::abs(), CompareTime(), CompareTowers(), F(), geo_, edm::EventSetup::get(), get, GetHaloClusterCandidateHB(), GetHaloClusterCandidateHE(), reco::HcalHaloData::GetPhiWedges(), reco::HcalHaloData::getProblematicStrips(), CaloGeometry::getSubdetectorGeometry(), HBRecHitEnergyThreshold, DetId::Hcal, HcalBarrel, HcalEndcap, HERecHitEnergyThreshold, hgeo_, mps_fire::i, hit::id, L1TowerCalibrationProducer_cfi::iEta, LEDCalibrationChannels::ieta, HcalDetId::ieta(), LEDCalibrationChannels::iphi, dqmiolumiharvest::j, SiStripPI::max, NHitsThreshold, edm::ESHandle< T >::product(), reco::HcalHaloData::setHaloClusterCandidatesHB(), reco::HcalHaloData::setHaloClusterCandidatesHE(), reco::PhiWedge::SetPlusZOriginConfidence(), digitizers_cfi::strip, SumEnergyThreshold, ntuplemaker::time, and hgcalTowerProducer_cfi::tower.

Referenced by Calculate(), and reco::HcalHaloDataProducer::produce().

◆ Calculate() [2/2]

HcalHaloData HcalHaloAlgo::Calculate ( const CaloGeometry &  TheCaloGeometry,
edm::Handle< HBHERecHitCollection > &  TheHBHERecHits,
edm::Handle< EBRecHitCollection > &  TheEBRecHits,
edm::Handle< EERecHitCollection > &  TheEERecHits,
const edm::EventSetup &  TheSetup 
)

Definition at line 30 of file HcalHaloAlgo.cc.

34  {
35  edm::Handle<CaloTowerCollection> TheCaloTowers;
36  return Calculate(TheCaloGeometry, TheHBHERecHits, TheCaloTowers, TheEBRecHits, TheEERecHits, TheSetup);
37 }

References Calculate().

◆ GetHaloClusterCandidateHB()

std::vector< HaloClusterCandidateHCAL > HcalHaloAlgo::GetHaloClusterCandidateHB ( edm::Handle< EcalRecHitCollection > &  ebrechitcoll,
edm::Handle< HBHERecHitCollection > &  hbherechitcoll,
float  et_thresh_seedrh 
)

Definition at line 253 of file HcalHaloAlgo.cc.

256  {
257  std::vector<HaloClusterCandidateHCAL> TheHaloClusterCandsHB;
258 
259  reco::Vertex::Point vtx(0, 0, 0);
260 
261  for (size_t ihit = 0; ihit < hbherechitcoll->size(); ++ihit) {
262  HaloClusterCandidateHCAL clustercand;
263 
264  const HBHERecHit& rechit = (*hbherechitcoll)[ihit];
265  math::XYZPoint rhpos = getPosition(rechit.id(), vtx);
266  //Et condition
267  double rhet = rechit.energy() * sqrt(rhpos.perp2() / rhpos.mag2());
268  if (rhet < et_thresh_seedrh)
269  continue;
270  if (std::abs(rhpos.z()) > zseparation_HBHE)
271  continue;
272  double eta = rhpos.eta();
273  double phi = rhpos.phi();
274 
275  bool isiso = true;
276  double etcluster(0);
277  int nbtowerssameeta(0);
278  double timediscriminatorITBH(0), timediscriminatorOTBH(0);
279  double etstrip_phiseedplus1(0), etstrip_phiseedminus1(0);
280 
281  //Building the cluster
282  edm::RefVector<HBHERecHitCollection> bhrhcandidates;
283  for (size_t jhit = 0; jhit < hbherechitcoll->size(); ++jhit) {
284  const HBHERecHit& rechitj = (*hbherechitcoll)[jhit];
285  HBHERecHitRef rhRef(hbherechitcoll, jhit);
286  math::XYZPoint rhposj = getPosition(rechitj.id(), vtx);
287  double rhetj = rechitj.energy() * sqrt(rhposj.perp2() / rhposj.mag2());
288  if (rhetj < 2)
289  continue;
290  if (std::abs(rhposj.z()) > zseparation_HBHE)
291  continue;
292  double etaj = rhposj.eta();
293  double phij = rhposj.phi();
294  double deta = eta - etaj;
295  double dphi = deltaPhi(phi, phij);
296  if (std::abs(deta) > 0.4)
297  continue; //This means +/-4 towers in eta
298  if (std::abs(dphi) > 0.2)
299  continue; //This means +/-2 towers in phi
300  if (std::abs(dphi) > 0.1 && std::abs(deta) < 0.2) {
301  isiso = false;
302  break;
303  } //The strip should be isolated
304  if (std::abs(dphi) > 0.1)
305  continue;
306  if (std::abs(dphi) < 0.05)
307  nbtowerssameeta++;
308  if (dphi > 0.05)
309  etstrip_phiseedplus1 += rhetj;
310  if (dphi < -0.05)
311  etstrip_phiseedminus1 += rhetj;
312 
313  etcluster += rhetj;
314  //Timing discriminator
315  //We assign a weight to the rechit defined as:
316  //Log10(Et)*f(T,R,Z)
317  //where f(T,R,Z) is the separation curve between halo-like and IP-like times.
318  //The time difference between a deposit from a outgoing IT halo and a deposit coming from a particle emitted at the IP is given by:
319  //dt= ( - sqrt(R^2+z^2) + |z| )/c
320  // For OT beam halo, the time difference is:
321  //dt= ( 25 + sqrt(R^2+z^2) + |z| )/c
322  //only consider the central part of HB as things get hard at large z.
323  //The best fitted value for R leads to 240 cm (IT) and 330 cm (OT)
324  double rhtj = rechitj.time();
325  timediscriminatorITBH +=
326  std::log10(rhetj) *
327  (rhtj + 0.5 * (sqrt(240. * 240. + rhposj.z() * rhposj.z()) - std::abs(rhposj.z())) / c_cm_per_ns);
328  if (std::abs(rhposj.z()) < 300)
329  timediscriminatorOTBH +=
330  std::log10(rhetj) *
331  (rhtj - 0.5 * (25 - (sqrt(330. * 330. + rhposj.z() * rhposj.z()) + std::abs(rhposj.z())) / c_cm_per_ns));
332  bhrhcandidates.push_back(rhRef);
333  }
334  //Isolation conditions
335  if (!isiso)
336  continue;
337  if (etstrip_phiseedplus1 / etcluster > 0.2 && etstrip_phiseedminus1 / etcluster > 0.2)
338  continue;
339 
340  //Calculate E/H
341  double eoh(0);
342  for (size_t jhit = 0; jhit < ecalrechitcoll->size(); ++jhit) {
343  const EcalRecHit& rechitj = (*ecalrechitcoll)[jhit];
344  math::XYZPoint rhposj = getPosition(rechitj.id(), vtx);
345  double rhetj = rechitj.energy() * sqrt(rhposj.perp2() / rhposj.mag2());
346  if (rhetj < 2)
347  continue;
348  double etaj = rhposj.eta();
349  double phij = rhposj.phi();
350  if (std::abs(eta - etaj) > 0.2)
351  continue;
352  if (std::abs(deltaPhi(phi, phij)) > 0.2)
353  continue;
354  eoh += rhetj / etcluster;
355  }
356  //E/H condition
357  if (eoh > 0.1)
358  continue;
359 
360  clustercand.setClusterEt(etcluster);
361  clustercand.setSeedEt(rhet);
362  clustercand.setSeedEta(eta);
363  clustercand.setSeedPhi(phi);
364  clustercand.setSeedZ(rhpos.Z());
365  clustercand.setSeedR(sqrt(rhpos.perp2()));
366  clustercand.setSeedTime(rechit.time());
367  clustercand.setEoverH(eoh);
368  clustercand.setNbTowersInEta(nbtowerssameeta);
369  clustercand.setEtStripPhiSeedPlus1(etstrip_phiseedplus1);
370  clustercand.setEtStripPhiSeedMinus1(etstrip_phiseedminus1);
371  clustercand.setTimeDiscriminatorITBH(timediscriminatorITBH);
372  clustercand.setTimeDiscriminatorOTBH(timediscriminatorOTBH);
373  clustercand.setBeamHaloRecHitsCandidates(bhrhcandidates);
374 
375  bool isbeamhalofrompattern = HBClusterShapeandTimeStudy(clustercand, false);
376  clustercand.setIsHaloFromPattern(isbeamhalofrompattern);
377  bool isbeamhalofrompattern_hlt = HBClusterShapeandTimeStudy(clustercand, true);
378  clustercand.setIsHaloFromPattern_HLT(isbeamhalofrompattern_hlt);
379 
380  TheHaloClusterCandsHB.push_back(clustercand);
381  }
382 
383  return TheHaloClusterCandsHB;
384 }

References funct::abs(), SiPixelRawToDigiRegional_cfi::deltaPhi, CaloRecHit::energy(), EcalRecHit::energy(), PVValHelper::eta, getPosition(), HBClusterShapeandTimeStudy(), HBHERecHit::id(), EcalRecHit::id(), phi, edm::RefVector< C, T, F >::push_back(), reco::HaloClusterCandidateHCAL::setBeamHaloRecHitsCandidates(), reco::HaloClusterCandidateHCAL::setClusterEt(), reco::HaloClusterCandidateHCAL::setEoverH(), reco::HaloClusterCandidateHCAL::setEtStripPhiSeedMinus1(), reco::HaloClusterCandidateHCAL::setEtStripPhiSeedPlus1(), reco::HaloClusterCandidateHCAL::setIsHaloFromPattern(), reco::HaloClusterCandidateHCAL::setIsHaloFromPattern_HLT(), reco::HaloClusterCandidateHCAL::setNbTowersInEta(), reco::HaloClusterCandidateHCAL::setSeedEt(), reco::HaloClusterCandidateHCAL::setSeedEta(), reco::HaloClusterCandidateHCAL::setSeedPhi(), reco::HaloClusterCandidateHCAL::setSeedR(), reco::HaloClusterCandidateHCAL::setSeedTime(), reco::HaloClusterCandidateHCAL::setSeedZ(), reco::HaloClusterCandidateHCAL::setTimeDiscriminatorITBH(), reco::HaloClusterCandidateHCAL::setTimeDiscriminatorOTBH(), edm::SortedCollection< T, SORT >::size(), mathSSE::sqrt(), CaloRecHit::time(), and badGlobalMuonTaggersAOD_cff::vtx.

Referenced by Calculate().

◆ GetHaloClusterCandidateHE()

std::vector< HaloClusterCandidateHCAL > HcalHaloAlgo::GetHaloClusterCandidateHE ( edm::Handle< EcalRecHitCollection > &  eerechitcoll,
edm::Handle< HBHERecHitCollection > &  hbherechitcoll,
float  et_thresh_seedrh 
)

Definition at line 386 of file HcalHaloAlgo.cc.

389  {
390  std::vector<HaloClusterCandidateHCAL> TheHaloClusterCandsHE;
391 
392  reco::Vertex::Point vtx(0, 0, 0);
393 
394  for (size_t ihit = 0; ihit < hbherechitcoll->size(); ++ihit) {
395  HaloClusterCandidateHCAL clustercand;
396 
397  const HBHERecHit& rechit = (*hbherechitcoll)[ihit];
398  math::XYZPoint rhpos = getPosition(rechit.id(), vtx);
399  //Et condition
400  double rhet = rechit.energy() * sqrt(rhpos.perp2() / rhpos.mag2());
401  if (rhet < et_thresh_seedrh)
402  continue;
403  if (std::abs(rhpos.z()) < zseparation_HBHE)
404  continue;
405  double eta = rhpos.eta();
406  double phi = rhpos.phi();
407  double rhr = sqrt(rhpos.perp2());
408  bool isiso = true;
409  double etcluster(0), hdepth1(0);
410  int clustersize(0);
411  double etstrip_phiseedplus1(0), etstrip_phiseedminus1(0);
412 
413  //Building the cluster
414  edm::RefVector<HBHERecHitCollection> bhrhcandidates;
415  for (size_t jhit = 0; jhit < hbherechitcoll->size(); ++jhit) {
416  const HBHERecHit& rechitj = (*hbherechitcoll)[jhit];
417  HBHERecHitRef rhRef(hbherechitcoll, jhit);
418  math::XYZPoint rhposj = getPosition(rechitj.id(), vtx);
419  double rhetj = rechitj.energy() * sqrt(rhposj.perp2() / rhposj.mag2());
420  if (rhetj < 2)
421  continue;
422  if (std::abs(rhposj.z()) < zseparation_HBHE)
423  continue;
424  if (rhpos.z() * rhposj.z() < 0)
425  continue;
426  double phij = rhposj.phi();
427  double dphi = deltaPhi(phi, phij);
428  if (std::abs(dphi) > 0.4)
429  continue;
430  double rhrj = sqrt(rhposj.perp2());
431  if (std::abs(rhr - rhrj) > 50)
432  continue;
433  if (std::abs(dphi) > 0.2 || std::abs(rhr - rhrj) > 20) {
434  isiso = false;
435  break;
436  } //The deposit should be isolated
437  if (dphi > 0.05)
438  etstrip_phiseedplus1 += rhetj;
439  if (dphi < -0.05)
440  etstrip_phiseedminus1 += rhetj;
441  clustersize++;
442  etcluster += rhetj;
443  if (std::abs(rhposj.z()) < 405)
444  hdepth1 += rhetj;
445  //No timing condition for now in HE
446  bhrhcandidates.push_back(rhRef);
447  }
448  //Isolation conditions
449  if (!isiso)
450  continue;
451  if (etstrip_phiseedplus1 / etcluster > 0.1 && etstrip_phiseedminus1 / etcluster > 0.1)
452  continue;
453 
454  //Calculate E/H
455  double eoh(0);
456  for (size_t jhit = 0; jhit < ecalrechitcoll->size(); ++jhit) {
457  const EcalRecHit& rechitj = (*ecalrechitcoll)[jhit];
458  math::XYZPoint rhposj = getPosition(rechitj.id(), vtx);
459  double rhetj = rechitj.energy() * sqrt(rhposj.perp2() / rhposj.mag2());
460  if (rhetj < 2)
461  continue;
462  if (rhpos.z() * rhposj.z() < 0)
463  continue;
464  double etaj = rhposj.eta();
465  double phij = rhposj.phi();
466  double dr = sqrt((eta - etaj) * (eta - etaj) + deltaPhi(phi, phij) * deltaPhi(phi, phij));
467  if (dr > 0.3)
468  continue;
469 
470  eoh += rhetj / etcluster;
471  }
472  //E/H condition
473  if (eoh > 0.1)
474  continue;
475 
476  clustercand.setClusterEt(etcluster);
477  clustercand.setSeedEt(rhet);
478  clustercand.setSeedEta(eta);
479  clustercand.setSeedPhi(phi);
480  clustercand.setSeedZ(rhpos.Z());
481  clustercand.setSeedR(sqrt(rhpos.perp2()));
482  clustercand.setSeedTime(rechit.time());
483  clustercand.setEoverH(eoh);
484  clustercand.setH1overH123(hdepth1 / etcluster);
485  clustercand.setClusterSize(clustersize);
486  clustercand.setEtStripPhiSeedPlus1(etstrip_phiseedplus1);
487  clustercand.setEtStripPhiSeedMinus1(etstrip_phiseedminus1);
488  clustercand.setTimeDiscriminator(0);
489  clustercand.setBeamHaloRecHitsCandidates(bhrhcandidates);
490 
491  bool isbeamhalofrompattern = HEClusterShapeandTimeStudy(clustercand, false);
492  clustercand.setIsHaloFromPattern(isbeamhalofrompattern);
493  bool isbeamhalofrompattern_hlt = HEClusterShapeandTimeStudy(clustercand, true);
494  clustercand.setIsHaloFromPattern_HLT(isbeamhalofrompattern_hlt);
495 
496  TheHaloClusterCandsHE.push_back(clustercand);
497  }
498 
499  return TheHaloClusterCandsHE;
500 }

References funct::abs(), SiPixelRawToDigiRegional_cfi::deltaPhi, flavorHistoryFilter_cfi::dr, CaloRecHit::energy(), EcalRecHit::energy(), PVValHelper::eta, getPosition(), HEClusterShapeandTimeStudy(), HBHERecHit::id(), EcalRecHit::id(), phi, edm::RefVector< C, T, F >::push_back(), reco::HaloClusterCandidateHCAL::setBeamHaloRecHitsCandidates(), reco::HaloClusterCandidateHCAL::setClusterEt(), reco::HaloClusterCandidateHCAL::setClusterSize(), reco::HaloClusterCandidateHCAL::setEoverH(), reco::HaloClusterCandidateHCAL::setEtStripPhiSeedMinus1(), reco::HaloClusterCandidateHCAL::setEtStripPhiSeedPlus1(), reco::HaloClusterCandidateHCAL::setH1overH123(), reco::HaloClusterCandidateHCAL::setIsHaloFromPattern(), reco::HaloClusterCandidateHCAL::setIsHaloFromPattern_HLT(), reco::HaloClusterCandidateHCAL::setSeedEt(), reco::HaloClusterCandidateHCAL::setSeedEta(), reco::HaloClusterCandidateHCAL::setSeedPhi(), reco::HaloClusterCandidateHCAL::setSeedR(), reco::HaloClusterCandidateHCAL::setSeedTime(), reco::HaloClusterCandidateHCAL::setSeedZ(), reco::HaloClusterCandidateHCAL::setTimeDiscriminator(), edm::SortedCollection< T, SORT >::size(), mathSSE::sqrt(), CaloRecHit::time(), and badGlobalMuonTaggersAOD_cff::vtx.

Referenced by Calculate().

◆ GetHBRecHitEnergyThreshold()

float HcalHaloAlgo::GetHBRecHitEnergyThreshold ( )
inline

Definition at line 75 of file HcalHaloAlgo.h.

75 { return HBRecHitEnergyThreshold; }

References HBRecHitEnergyThreshold.

◆ GetHERecHitEnergyThreshold()

float HcalHaloAlgo::GetHERecHitEnergyThreshold ( )
inline

Definition at line 76 of file HcalHaloAlgo.h.

76 { return HERecHitEnergyThreshold; }

References HERecHitEnergyThreshold.

◆ GetPhiWedgeEnergyThreshold()

float HcalHaloAlgo::GetPhiWedgeEnergyThreshold ( )
inline

Definition at line 79 of file HcalHaloAlgo.h.

79 { return SumEnergyThreshold; }

References SumEnergyThreshold.

◆ GetPhiWedgeNHitsThreshold()

int HcalHaloAlgo::GetPhiWedgeNHitsThreshold ( )
inline

Definition at line 80 of file HcalHaloAlgo.h.

80 { return NHitsThreshold; }

References NHitsThreshold.

◆ getPosition()

math::XYZPoint HcalHaloAlgo::getPosition ( const DetId &  id,
reco::Vertex::Point  vtx 
)
private

Definition at line 558 of file HcalHaloAlgo.cc.

558  {
559  const GlobalPoint pos = ((id.det() == DetId::Hcal) ? hgeo_->getPosition(id) : GlobalPoint(geo_->getPosition(id)));
560  math::XYZPoint posV(pos.x() - vtx.x(), pos.y() - vtx.y(), pos.z() - vtx.z());
561  return posV;
562 }

References geo_, CaloGeometry::getPosition(), HcalGeometry::getPosition(), DetId::Hcal, hgeo_, and badGlobalMuonTaggersAOD_cff::vtx.

Referenced by GetHaloClusterCandidateHB(), and GetHaloClusterCandidateHE().

◆ HBClusterShapeandTimeStudy()

bool HcalHaloAlgo::HBClusterShapeandTimeStudy ( reco::HaloClusterCandidateHCAL  hcand,
bool  ishlt 
)

Definition at line 502 of file HcalHaloAlgo.cc.

502  {
503  //Conditions on the central strip size in eta.
504  //For low size, extra conditions on seed et, isolation and cluster timing
505  //Here we target both IT and OT beam halo. Two separate discriminators were built for the two cases.
506 
507  if (hcand.getSeedEt() < 10)
508  return false;
509 
510  if (hcand.getNbTowersInEta() < 3)
511  return false;
512  //Isolation criteria for very short eta strips
513  if (hcand.getNbTowersInEta() == 3 && (hcand.getEtStripPhiSeedPlus1() > 0.1 || hcand.getEtStripPhiSeedMinus1() > 0.1))
514  return false;
515  if (hcand.getNbTowersInEta() <= 5 && (hcand.getEtStripPhiSeedPlus1() > 0.1 && hcand.getEtStripPhiSeedMinus1() > 0.1))
516  return false;
517 
518  //Timing conditions for short eta strips
519  if (hcand.getNbTowersInEta() == 3 && hcand.getTimeDiscriminatorITBH() >= 0.)
520  return false;
521  if (hcand.getNbTowersInEta() <= 6 && hcand.getTimeDiscriminatorITBH() >= 5. && hcand.getTimeDiscriminatorOTBH() < 0.)
522  return false;
523 
524  //For HLT, only use conditions without timing
525  if (ishlt && hcand.getNbTowersInEta() < 7)
526  return false;
527 
528  hcand.setIsHaloFromPattern(true);
529 
530  return true;
531 }

References reco::HaloClusterCandidateHCAL::getEtStripPhiSeedMinus1(), reco::HaloClusterCandidateHCAL::getEtStripPhiSeedPlus1(), reco::HaloClusterCandidateHCAL::getNbTowersInEta(), reco::HaloClusterCandidateHCAL::getSeedEt(), reco::HaloClusterCandidateHCAL::getTimeDiscriminatorITBH(), reco::HaloClusterCandidateHCAL::getTimeDiscriminatorOTBH(), and reco::HaloClusterCandidateHCAL::setIsHaloFromPattern().

Referenced by GetHaloClusterCandidateHB().

◆ HEClusterShapeandTimeStudy()

bool HcalHaloAlgo::HEClusterShapeandTimeStudy ( reco::HaloClusterCandidateHCAL  hcand,
bool  ishlt 
)

Definition at line 533 of file HcalHaloAlgo.cc.

533  {
534  //Conditions on H1/H123 to spot halo interacting only in one HCAL layer.
535  //For R> about 170cm, HE has only one layer and this condition cannot be applied
536  //Note that for R>170 cm, the halo is in CSC acceptance and will most likely be spotted by the CSC-calo matching method
537  //A method to identify halos interacting in both H1 and H2/H3 at low R is still missing.
538 
539  if (hcand.getSeedEt() < 20)
540  return false;
541  if (hcand.getSeedR() > 170)
542  return false;
543 
544  if (hcand.getH1overH123() > 0.02 && hcand.getH1overH123() < 0.98)
545  return false;
546 
547  //This method is one of the ones with the highest fake rate: in JetHT dataset, it happens in around 0.1% of the cases that a low pt jet (pt= 20) leaves all of its energy in only one HCAL layer.
548  //At HLT, one only cares about large deposits from BH that would lead to a MET/SinglePhoton trigger to be fired.
549  //Rising the seed Et threshold at HLT has therefore little impact on the HLT performances but ensures that possible controversial events are still recorded.
550  if (ishlt && hcand.getSeedEt() < 50)
551  return false;
552 
553  hcand.setIsHaloFromPattern(true);
554 
555  return true;
556 }

References reco::HaloClusterCandidateHCAL::getH1overH123(), reco::HaloClusterCandidateHCAL::getSeedEt(), reco::HaloClusterCandidateHCAL::getSeedR(), and reco::HaloClusterCandidateHCAL::setIsHaloFromPattern().

Referenced by GetHaloClusterCandidateHE().

◆ SetPhiWedgeEnergyThreshold()

void HcalHaloAlgo::SetPhiWedgeEnergyThreshold ( float  SumE)
inline

Definition at line 67 of file HcalHaloAlgo.h.

67 { SumEnergyThreshold = SumE; }

References SumEnergyThreshold.

◆ SetPhiWedgeNHitsThreshold()

void HcalHaloAlgo::SetPhiWedgeNHitsThreshold ( int  nhits)
inline

Definition at line 68 of file HcalHaloAlgo.h.

68 { NHitsThreshold = nhits; }

References nhits, and NHitsThreshold.

◆ SetPhiWedgeThresholds()

void HcalHaloAlgo::SetPhiWedgeThresholds ( float  SumE,
int  nhits 
)
inline

Definition at line 69 of file HcalHaloAlgo.h.

69  {
70  SumEnergyThreshold = SumE;
71  NHitsThreshold = nhits;
72  }

References nhits, NHitsThreshold, and SumEnergyThreshold.

Referenced by reco::HcalHaloDataProducer::produce().

◆ SetRecHitEnergyThresholds()

void HcalHaloAlgo::SetRecHitEnergyThresholds ( float  HB,
float  HE 
)
inline

Definition at line 61 of file HcalHaloAlgo.h.

61  {
62  HBRecHitEnergyThreshold = HB;
63  HERecHitEnergyThreshold = HE;
64  }

References HB, HBRecHitEnergyThreshold, HE, and HERecHitEnergyThreshold.

Referenced by reco::HcalHaloDataProducer::produce().

Member Data Documentation

◆ geo_

const CaloGeometry* HcalHaloAlgo::geo_
private

Definition at line 102 of file HcalHaloAlgo.h.

Referenced by Calculate(), and getPosition().

◆ HBRecHitEnergyThreshold

float HcalHaloAlgo::HBRecHitEnergyThreshold
private

Definition at line 95 of file HcalHaloAlgo.h.

Referenced by Calculate(), GetHBRecHitEnergyThreshold(), HcalHaloAlgo(), and SetRecHitEnergyThresholds().

◆ HERecHitEnergyThreshold

float HcalHaloAlgo::HERecHitEnergyThreshold
private

Definition at line 96 of file HcalHaloAlgo.h.

Referenced by Calculate(), GetHERecHitEnergyThreshold(), HcalHaloAlgo(), and SetRecHitEnergyThresholds().

◆ hgeo_

const HcalGeometry* HcalHaloAlgo::hgeo_
private

Definition at line 103 of file HcalHaloAlgo.h.

Referenced by Calculate(), and getPosition().

◆ NHitsThreshold

int HcalHaloAlgo::NHitsThreshold
private

Definition at line 100 of file HcalHaloAlgo.h.

Referenced by Calculate(), GetPhiWedgeNHitsThreshold(), HcalHaloAlgo(), SetPhiWedgeNHitsThreshold(), and SetPhiWedgeThresholds().

◆ SumEnergyThreshold

float HcalHaloAlgo::SumEnergyThreshold
private

Definition at line 99 of file HcalHaloAlgo.h.

Referenced by Calculate(), GetPhiWedgeEnergyThreshold(), HcalHaloAlgo(), SetPhiWedgeEnergyThreshold(), and SetPhiWedgeThresholds().

edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
reco::HaloClusterCandidateHCAL::setSeedZ
void setSeedZ(double x)
Definition: HaloClusterCandidateHCAL.h:52
EcalRecHit
Definition: EcalRecHit.h:15
mps_fire.i
i
Definition: mps_fire.py:355
HcalHaloAlgo::HBRecHitEnergyThreshold
float HBRecHitEnergyThreshold
Definition: HcalHaloAlgo.h:95
hit::id
unsigned int id
Definition: SiStripHitEffFromCalibTree.cc:92
CaloRecHit::energy
constexpr float energy() const
Definition: CaloRecHit.h:29
EcalRecHit::id
DetId id() const
get the id
Definition: EcalRecHit.h:77
reco::HaloClusterCandidateHCAL::setTimeDiscriminator
void setTimeDiscriminator(double x)
Definition: HaloClusterCandidateHCAL.h:68
reco::HaloClusterCandidateHCAL::setBeamHaloRecHitsCandidates
void setBeamHaloRecHitsCandidates(edm::RefVector< HBHERecHitCollection > x)
Definition: HaloClusterCandidateHCAL.h:60
CaloGeometry::getPosition
GlobalPoint getPosition(const DetId &id) const
Get the position of a given detector id.
Definition: CaloGeometry.cc:50
reco::HaloClusterCandidateHCAL::setSeedEt
void setSeedEt(double x)
Definition: HaloClusterCandidateHCAL.h:49
CaloGeometryRecord
Definition: CaloGeometryRecord.h:30
reco::HcalHaloData::setHaloClusterCandidatesHB
void setHaloClusterCandidatesHB(const std::vector< HaloClusterCandidateHCAL > &x)
Definition: HcalHaloData.h:58
digitizers_cfi.strip
strip
Definition: digitizers_cfi.py:19
HBHERecHit
Definition: HBHERecHit.h:13
reco::HaloClusterCandidateHCAL::setClusterEt
void setClusterEt(double x)
Definition: HaloClusterCandidateHCAL.h:48
pos
Definition: PixelAliasList.h:18
reco::HaloClusterCandidateHCAL::setSeedR
void setSeedR(double x)
Definition: HaloClusterCandidateHCAL.h:53
DetId::Hcal
Definition: DetId.h:28
reco::HaloClusterCandidateHCAL::setNbTowersInEta
void setNbTowersInEta(double x)
Definition: HaloClusterCandidateHCAL.h:62
hgcalTowerProducer_cfi.tower
tower
Definition: hgcalTowerProducer_cfi.py:3
CaloGeometry::getSubdetectorGeometry
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:34
reco::HcalHaloData
Definition: HcalHaloData.h:37
reco::HaloClusterCandidateHCAL::getNbTowersInEta
int getNbTowersInEta() const
Definition: HaloClusterCandidateHCAL.h:40
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:215
reco::HaloClusterCandidateHCAL::getTimeDiscriminatorOTBH
double getTimeDiscriminatorOTBH() const
Definition: HaloClusterCandidateHCAL.h:42
reco::HaloClusterCandidateHCAL::setClusterSize
void setClusterSize(int x)
Definition: HaloClusterCandidateHCAL.h:67
reco::HaloClusterCandidateHCAL::setTimeDiscriminatorOTBH
void setTimeDiscriminatorOTBH(double x)
Definition: HaloClusterCandidateHCAL.h:64
reco::HaloClusterCandidateHCAL::setIsHaloFromPattern_HLT
void setIsHaloFromPattern_HLT(bool x)
Definition: HaloClusterCandidateHCAL.h:56
HcalHaloAlgo::GetHaloClusterCandidateHE
std::vector< reco::HaloClusterCandidateHCAL > GetHaloClusterCandidateHE(edm::Handle< EcalRecHitCollection > &eerechitcoll, edm::Handle< HBHERecHitCollection > &hbherechitcoll, float et_thresh_seedrh)
Definition: HcalHaloAlgo.cc:386
edm::RefVector
Definition: EDProductfwd.h:27
EcalRecHit::energy
float energy() const
Definition: EcalRecHit.h:68
HcalBarrel
Definition: HcalAssistant.h:33
HB
Definition: GlobalHaloAlgo.cc:15
edm::Handle
Definition: AssociativeIterator.h:50
edm::Ref
Definition: AssociativeIterator.h:58
F
static uInt32 F(BLOWFISH_CTX *ctx, uInt32 x)
Definition: blowfish.cc:163
LEDCalibrationChannels.iphi
iphi
Definition: LEDCalibrationChannels.py:64
HcalHaloAlgo::HBClusterShapeandTimeStudy
bool HBClusterShapeandTimeStudy(reco::HaloClusterCandidateHCAL hcand, bool ishlt)
Definition: HcalHaloAlgo.cc:502
HcalHaloAlgo::geo_
const CaloGeometry * geo_
Definition: HcalHaloAlgo.h:102
reco::PhiWedge
Definition: PhiWedge.h:13
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
SiPixelRawToDigiRegional_cfi.deltaPhi
deltaPhi
Definition: SiPixelRawToDigiRegional_cfi.py:9
PVValHelper::eta
Definition: PVValidationHelpers.h:69
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
edm::ESHandle< CaloGeometry >
reco::HaloClusterCandidateHCAL::getH1overH123
double getH1overH123() const
Definition: HaloClusterCandidateHCAL.h:44
HcalHaloAlgo::HERecHitEnergyThreshold
float HERecHitEnergyThreshold
Definition: HcalHaloAlgo.h:96
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
Point3DBase< float, GlobalTag >
nhits
Definition: HIMultiTrackSelector.h:42
badGlobalMuonTaggersAOD_cff.vtx
vtx
Definition: badGlobalMuonTaggersAOD_cff.py:5
CaloRecHit::time
constexpr float time() const
Definition: CaloRecHit.h:31
reco::HaloClusterCandidateHCAL::getSeedR
double getSeedR() const
Definition: HaloClusterCandidateHCAL.h:31
reco::HcalHaloData::getProblematicStrips
const std::vector< HaloTowerStrip > & getProblematicStrips() const
Definition: HcalHaloData.h:49
LEDCalibrationChannels.ieta
ieta
Definition: LEDCalibrationChannels.py:63
HcalDetId::ieta
constexpr int ieta() const
get the cell ieta
Definition: HcalDetId.h:155
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
HcalGeometry::getPosition
GlobalPoint getPosition(const DetId &id) const
Definition: HcalGeometry.cc:178
reco::HaloClusterCandidateHCAL::setH1overH123
void setH1overH123(double x)
Definition: HaloClusterCandidateHCAL.h:66
SiStripPI::max
Definition: SiStripPayloadInspectorHelper.h:169
reco::HcalHaloData::setHaloClusterCandidatesHE
void setHaloClusterCandidatesHE(const std::vector< HaloClusterCandidateHCAL > &x)
Definition: HcalHaloData.h:59
reco::HaloClusterCandidateHCAL::getSeedEt
double getSeedEt() const
Definition: HaloClusterCandidateHCAL.h:27
reco::HaloClusterCandidateHCAL::setEtStripPhiSeedMinus1
void setEtStripPhiSeedMinus1(double x)
Definition: HaloClusterCandidateHCAL.h:59
reco::HaloClusterCandidateHCAL::setIsHaloFromPattern
void setIsHaloFromPattern(bool x)
Definition: HaloClusterCandidateHCAL.h:55
HcalDetId
Definition: HcalDetId.h:12
CaloTower
Definition: CaloTower.h:26
HE
Definition: GlobalHaloAlgo.cc:15
reco::HaloClusterCandidateHCAL::setTimeDiscriminatorITBH
void setTimeDiscriminatorITBH(double x)
Definition: HaloClusterCandidateHCAL.h:63
HcalHaloAlgo::getPosition
math::XYZPoint getPosition(const DetId &id, reco::Vertex::Point vtx)
Definition: HcalHaloAlgo.cc:558
CompareTime
bool CompareTime(const HBHERecHit *x, const HBHERecHit *y)
Definition: HcalHaloAlgo.cc:18
reco::HcalHaloData::GetPhiWedges
const std::vector< PhiWedge > & GetPhiWedges() const
Definition: HcalHaloData.h:45
get
#define get
edm::RefVector::push_back
void push_back(value_type const &ref)
Add a Ref<C, T> to the RefVector.
Definition: RefVector.h:67
reco::HaloClusterCandidateHCAL::setSeedTime
void setSeedTime(double x)
Definition: HaloClusterCandidateHCAL.h:54
reco::HaloClusterCandidateHCAL::getTimeDiscriminatorITBH
double getTimeDiscriminatorITBH() const
Definition: HaloClusterCandidateHCAL.h:41
HcalHaloAlgo::hgeo_
const HcalGeometry * hgeo_
Definition: HcalHaloAlgo.h:103
DDAxes::phi
reco::Vertex::Point
math::XYZPoint Point
point in the space
Definition: Vertex.h:40
CompareTowers
bool CompareTowers(const CaloTower *x, const CaloTower *y)
Definition: HcalHaloAlgo.cc:19
reco::HaloClusterCandidateHCAL
Definition: HaloClusterCandidateHCAL.h:21
HcalEndcap
Definition: HcalAssistant.h:34
reco::HaloClusterCandidateHCAL::getEtStripPhiSeedMinus1
double getEtStripPhiSeedMinus1() const
Definition: HaloClusterCandidateHCAL.h:37
flavorHistoryFilter_cfi.dr
dr
Definition: flavorHistoryFilter_cfi.py:37
reco::HaloClusterCandidateHCAL::getEtStripPhiSeedPlus1
double getEtStripPhiSeedPlus1() const
Definition: HaloClusterCandidateHCAL.h:36
HaloTowerStrip
Definition: HcalHaloData.h:16
L1TowerCalibrationProducer_cfi.iEta
iEta
Definition: L1TowerCalibrationProducer_cfi.py:60
HcalHaloAlgo::NHitsThreshold
int NHitsThreshold
Definition: HcalHaloAlgo.h:100
HcalHaloAlgo::SumEnergyThreshold
float SumEnergyThreshold
Definition: HcalHaloAlgo.h:99
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
HBHERecHit::id
constexpr HcalDetId id() const
get the id
Definition: HBHERecHit.h:39
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
ntuplemaker.time
time
Definition: ntuplemaker.py:310
reco::HaloClusterCandidateHCAL::setEtStripPhiSeedPlus1
void setEtStripPhiSeedPlus1(double x)
Definition: HaloClusterCandidateHCAL.h:58
HcalHaloAlgo::HEClusterShapeandTimeStudy
bool HEClusterShapeandTimeStudy(reco::HaloClusterCandidateHCAL hcand, bool ishlt)
Definition: HcalHaloAlgo.cc:533
reco::HaloClusterCandidateHCAL::setEoverH
void setEoverH(double x)
Definition: HaloClusterCandidateHCAL.h:57
reco::HaloClusterCandidateHCAL::setSeedPhi
void setSeedPhi(double x)
Definition: HaloClusterCandidateHCAL.h:51
HcalHaloAlgo::Calculate
reco::HcalHaloData Calculate(const CaloGeometry &TheCaloGeometry, edm::Handle< HBHERecHitCollection > &TheHBHERecHits, edm::Handle< CaloTowerCollection > &TheCaloTowers, edm::Handle< EBRecHitCollection > &TheEBRecHits, edm::Handle< EERecHitCollection > &TheEERecHits, const edm::EventSetup &TheSetup)
Definition: HcalHaloAlgo.cc:39
reco::HaloClusterCandidateHCAL::setSeedEta
void setSeedEta(double x)
Definition: HaloClusterCandidateHCAL.h:50
hit
Definition: SiStripHitEffFromCalibTree.cc:88
HcalHaloAlgo::GetHaloClusterCandidateHB
std::vector< reco::HaloClusterCandidateHCAL > GetHaloClusterCandidateHB(edm::Handle< EcalRecHitCollection > &ebrechitcoll, edm::Handle< HBHERecHitCollection > &hbherechitcoll, float et_thresh_seedrh)
Definition: HcalHaloAlgo.cc:253
Generated for CMSSW Reference Manual by doxygen 1.8.16