CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Private Member Functions | Private Attributes
CosmicClusterAlgo Class Reference

#include <CosmicClusterAlgo.h>

Public Types

enum  EcalPart { barrel = 0, endcap = 1 }
 
typedef math::XYZPoint Point
 point in the space More...
 
enum  VerbosityLevel { pDEBUG = 0, pWARNING = 1, pINFO = 2, pERROR = 3 }
 

Public Member Functions

 CosmicClusterAlgo ()
 
 CosmicClusterAlgo (double ebst, double ebSt, double ebDt, double ebSp, double ecst, double ecSt, double ecDt, double ecSp, const PositionCalc &posCalc, VerbosityLevel the_verbosity=pERROR)
 
std::vector< reco::BasicClustermakeClusters (const EcalRecHitCollection *hits, const EcalUncalibratedRecHitCollection *uncalibhits, const CaloSubdetectorGeometry *geometry, const CaloSubdetectorTopology *topology_p, const CaloSubdetectorGeometry *geometryES_p, EcalPart ecalPart, bool regional=false, const std::vector< RectangularEtaPhiRegion > &regions=std::vector< RectangularEtaPhiRegion >())
 
void setVerbosity (VerbosityLevel the_verbosity)
 
virtual ~CosmicClusterAlgo ()
 

Private Member Functions

void addCrystal (const DetId &det, const bool in9)
 
bool checkMaxima (CaloNavigator< DetId > &navigator)
 
void mainSearch (const CaloSubdetectorGeometry *geometry_p, const CaloSubdetectorTopology *topology_p, const CaloSubdetectorGeometry *geometryES_p, EcalPart ecalPart)
 
void makeCluster (const CaloSubdetectorGeometry *geometry_p, const CaloSubdetectorGeometry *geometryES_p, DetId seedId)
 
void prepareCluster (CaloNavigator< DetId > &navigator, const CaloSubdetectorGeometry *geometry)
 

Private Attributes

std::set< DetIdcanSeed_s
 
std::vector< reco::BasicClusterclusters_v
 
std::vector< DetIdcurrent_v25
 
std::vector< std::pair< DetId, float > > current_v25Sup
 
std::vector< DetIdcurrent_v9
 
double ecalBarrelSecondThreshold
 
double ecalBarrelSeedThreshold
 
double ecalBarrelSingleThreshold
 
double ecalBarrelSupThreshold
 
double ecalEndcapSecondThreshold
 
double ecalEndcapSeedThreshold
 
double ecalEndcapSingleThreshold
 
double ecalEndcapSupThreshold
 
bool inEB
 
PositionCalc posCalculator_
 
const EcalRecHitCollectionrecHits_
 
std::vector< EcalRecHitseeds
 
const EcalUncalibratedRecHitCollectionuncalibRecHits_
 
std::set< DetIdused_s
 
VerbosityLevel verbosity
 

Detailed Description

Definition at line 27 of file CosmicClusterAlgo.h.

Member Typedef Documentation

◆ Point

point in the space

Definition at line 72 of file CosmicClusterAlgo.h.

Member Enumeration Documentation

◆ EcalPart

Enumerator
barrel 
endcap 

Definition at line 29 of file CosmicClusterAlgo.h.

29 { barrel = 0, endcap = 1 };

◆ VerbosityLevel

Enumerator
pDEBUG 
pWARNING 
pINFO 
pERROR 

Definition at line 30 of file CosmicClusterAlgo.h.

30 { pDEBUG = 0, pWARNING = 1, pINFO = 2, pERROR = 3 };

Constructor & Destructor Documentation

◆ CosmicClusterAlgo() [1/2]

CosmicClusterAlgo::CosmicClusterAlgo ( )
inline

Definition at line 32 of file CosmicClusterAlgo.h.

32 {}

◆ CosmicClusterAlgo() [2/2]

CosmicClusterAlgo::CosmicClusterAlgo ( double  ebst,
double  ebSt,
double  ebDt,
double  ebSp,
double  ecst,
double  ecSt,
double  ecDt,
double  ecSp,
const PositionCalc posCalc,
VerbosityLevel  the_verbosity = pERROR 
)
inline

Definition at line 34 of file CosmicClusterAlgo.h.

52  verbosity(the_verbosity) {
53  posCalculator_ = posCalc;
54  }

References posCalculator_.

◆ ~CosmicClusterAlgo()

virtual CosmicClusterAlgo::~CosmicClusterAlgo ( )
inlinevirtual

Definition at line 56 of file CosmicClusterAlgo.h.

56 {}

Member Function Documentation

◆ addCrystal()

void CosmicClusterAlgo::addCrystal ( const DetId det,
const bool  in9 
)
private

Definition at line 359 of file CosmicClusterAlgo.cc.

359  {
361  if ((thisIt != recHits_->end()) && (thisIt->id() != DetId(0))) {
362  if ((used_s.find(thisIt->id()) == used_s.end())) {
364  used_s.insert(det);
365  if ((thisIt->energy() >= -1.) && !(thisItu->chi2() < -1.)) {
366  if (in9)
367  current_v9.push_back(det);
368  current_v25.push_back(det);
369  }
370  }
371  }
372 }

References current_v25, current_v9, edm::SortedCollection< T, SORT >::end(), edm::SortedCollection< T, SORT >::find(), recHits_, uncalibRecHits_, and used_s.

Referenced by prepareCluster().

◆ checkMaxima()

bool CosmicClusterAlgo::checkMaxima ( CaloNavigator< DetId > &  navigator)
private

Definition at line 275 of file CosmicClusterAlgo.cc.

275  {
276  bool maxima = true;
279  double thisEnergy = 0.;
280  double seedEnergy = seedHit->energy();
281 
282  std::vector<DetId> swissCrossVec;
283  swissCrossVec.clear();
284 
285  swissCrossVec.push_back(navigator.west());
286  navigator.home();
287  swissCrossVec.push_back(navigator.east());
288  navigator.home();
289  swissCrossVec.push_back(navigator.north());
290  navigator.home();
291  swissCrossVec.push_back(navigator.south());
292  navigator.home();
293 
294  for (unsigned int i = 0; i < swissCrossVec.size(); ++i) {
295  thisHit = recHits_->find(swissCrossVec[i]);
296  if ((swissCrossVec[i] == DetId(0)) || thisHit == recHits_->end())
297  thisEnergy = 0.0;
298  else
299  thisEnergy = thisHit->energy();
300  if (thisEnergy > seedEnergy) {
301  maxima = false;
302  break;
303  }
304  }
305 
306  return maxima;
307 }

References edm::SortedCollection< T, SORT >::end(), edm::SortedCollection< T, SORT >::find(), mps_fire::i, HLT_FULL_cff::navigator, and recHits_.

Referenced by mainSearch().

◆ mainSearch()

void CosmicClusterAlgo::mainSearch ( const CaloSubdetectorGeometry geometry_p,
const CaloSubdetectorTopology topology_p,
const CaloSubdetectorGeometry geometryES_p,
EcalPart  ecalPart 
)
private

Definition at line 134 of file CosmicClusterAlgo.cc.

137  {
138  if (verbosity < pINFO) {
139  std::cout << "Building clusters............" << std::endl;
140  }
141 
142  // Loop over seeds:
143  std::vector<EcalRecHit>::iterator it;
144  for (it = seeds.begin(); it != seeds.end(); it++) {
145  // check if this crystal is able to seed
146  // (event though it is already used)
147  bool usedButCanSeed = false;
148  if (canSeed_s.find(it->id()) != canSeed_s.end())
149  usedButCanSeed = true;
150 
151  // make sure the current seed does not belong to a cluster already.
152  if ((used_s.find(it->id()) != used_s.end()) && (usedButCanSeed == false)) {
153  if (it == seeds.begin()) {
154  if (verbosity < pINFO) {
155  std::cout << "##############################################################" << std::endl;
156  std::cout << "DEBUG ALERT: Highest energy seed already belongs to a cluster!" << std::endl;
157  std::cout << "##############################################################" << std::endl;
158  }
159  }
160 
161  // seed crystal is used or is used and cannot seed a cluster
162  // so continue to the next seed crystal...
163  continue;
164  }
165 
166  // clear the vector of hits in current clusterconst EcalRecHitCollection* hits,
167  current_v9.clear();
168  current_v25.clear();
169  current_v25Sup.clear();
170 
171  // Create a navigator at the seed
172  CaloNavigator<DetId> navigator(it->id(), topology_p);
173  DetId seedId = navigator.pos();
174  navigator.setHome(seedId);
175 
176  // Is the seed a local maximum?
177  bool localMaxima = checkMaxima(navigator);
178 
179  if (localMaxima) {
180  // build the 5x5 taking care over which crystals //JHaupt 4-27-08 3x3 is a good option...
181  // can seed new clusters and which can't
182  prepareCluster(navigator, geometry_p);
183  }
184 
185  // If some crystals in the current vector then
186  // make them into a cluster
187  if (!current_v25.empty()) {
188  makeCluster(geometry_p, geometryES_p, seedId);
189  }
190 
191  } // End loop on seed crystals
192 }

References canSeed_s, checkMaxima(), gather_cfg::cout, current_v25, current_v25Sup, current_v9, makeCluster(), HLT_FULL_cff::navigator, pINFO, prepareCluster(), seeds, used_s, and verbosity.

Referenced by makeClusters().

◆ makeCluster()

void CosmicClusterAlgo::makeCluster ( const CaloSubdetectorGeometry geometry_p,
const CaloSubdetectorGeometry geometryES_p,
DetId  seedId 
)
private

Definition at line 194 of file CosmicClusterAlgo.cc.

196  {
197  double energy = 0;
198  double energySecond = 0.; //JHaupt 4-27-08 Added for the second crystal stream
199  double energyMax = 0.; //JHaupt 4-27-08 Added for the max crystal stream
200  DetId detFir;
201  DetId detSec;
202  //bool goodCluster = false; //JHaupt 4-27-08 Added so that some can be earased.. used another day Might not be needed as seeds are energy ordered...
203 
204  std::vector<DetId>::iterator it;
205  for (it = current_v9.begin(); it != current_v9.end(); it++) {
206  // Martina - check recoFlag for crystals sorrounding the good one
208  // double-check that iterator can be dereferenced
209  if (itt == recHits_->end())
210  continue;
211  // if rechit affected by features other than these, do not allow 2 crystal seeding either
212  // features are hard coded, for now; add severity?
213  uint32_t rhFlag = (*itt).recoFlag();
214  if (!(rhFlag == EcalRecHit::kGood || rhFlag == EcalRecHit::kOutOfTime || rhFlag == EcalRecHit::kPoorCalib))
215  continue;
216 
218 
220  EcalUncalibratedRecHit uhit_p = *ittu;
221 
222  if (uhit_p.amplitude() > energySecond) {
223  energySecond = uhit_p.amplitude();
224  detSec = uhit_p.id();
225  }
226  if (energySecond > energyMax) {
227  std::swap(energySecond, energyMax);
228  std::swap(detFir, detSec);
229  }
230  }
231 
232  //if ((detFir.det() == DetId::Ecal) && (detFir.subdetId() == EcalEndcap)) inEB = false;
233  //We ignore the possiblity that the cluster may span into the EE
234 
235  if ((energyMax < (inEB ? ecalBarrelSingleThreshold : ecalEndcapSingleThreshold)) &&
237  return;
238 
239  for (it = current_v25.begin(); it != current_v25.end(); it++) {
241  EcalRecHit hit_p = *itt;
242 
244  EcalUncalibratedRecHit uhit_p = *ittu;
245 
247  // energy fraction = 1
248  current_v25Sup.push_back(std::pair<DetId, float>(hit_p.id(), 1.));
249  energy += hit_p.energy(); //Keep the fully corrected energy
250  }
251  }
252 
253  Point position;
255 
256  //don't write empty clusters
257  if (energy == 0 && position == Point(0, 0, 0))
258  return;
259 
260  if (verbosity < pINFO) {
261  std::cout << "JH******** NEW CLUSTER ********" << std::endl;
262  std::cout << "JHNo. of crystals = " << current_v25Sup.size() << std::endl;
263  std::cout << "JH Energy = " << energy << std::endl;
264  std::cout << "JH Phi = " << position.phi() << std::endl;
265  std::cout << "JH Eta = " << position.eta() << std::endl;
266  std::cout << "JH*****************************" << std::endl;
267  // specialize this
268  // std::cout << "JH****Emax**** "<<energyMax << " ieta " <<detFir.ieta() <<" iphi "<<detFir.ieta() << std::endl;
269  // std::cout << "JH****Esec**** "<<energySecond << " ieta " <<detSec.ieta() <<" iphi "<<detSec.ieta() << std::endl;
270  }
271  clusters_v.push_back(
273 }

References EcalUncalibratedRecHit::amplitude(), PositionCalc::Calculate_Location(), clusters_v, gather_cfg::cout, current_v25, current_v25Sup, current_v9, ecalBarrelSecondThreshold, ecalBarrelSingleThreshold, ecalBarrelSupThreshold, ecalEndcapSecondThreshold, ecalEndcapSingleThreshold, ecalEndcapSupThreshold, edm::SortedCollection< T, SORT >::end(), HCALHighEnergyHPDFilter_cfi::energy, EcalRecHit::energy(), edm::SortedCollection< T, SORT >::find(), EcalUncalibratedRecHit::id(), EcalRecHit::id(), inEB, EcalRecHit::kGood, EcalRecHit::kOutOfTime, EcalRecHit::kPoorCalib, reco::CaloCluster::multi5x5, pINFO, posCalculator_, position, recHits_, std::swap(), uncalibRecHits_, and verbosity.

Referenced by mainSearch().

◆ makeClusters()

std::vector< reco::BasicCluster > CosmicClusterAlgo::makeClusters ( const EcalRecHitCollection hits,
const EcalUncalibratedRecHitCollection uncalibhits,
const CaloSubdetectorGeometry geometry,
const CaloSubdetectorTopology topology_p,
const CaloSubdetectorGeometry geometryES_p,
EcalPart  ecalPart,
bool  regional = false,
const std::vector< RectangularEtaPhiRegion > &  regions = std::vector<RectangularEtaPhiRegion>() 
)

Definition at line 16 of file CosmicClusterAlgo.cc.

23  {
24  seeds.clear();
25  used_s.clear();
26  canSeed_s.clear();
27  clusters_v.clear();
28 
29  recHits_ = hits;
30  uncalibRecHits_ = uncalibhits;
31 
32  inEB = true;
33  double threshold = 0;
34  std::string ecalPart_string;
35  if (ecalPart == endcap) {
37  ecalPart_string = "EndCap";
38  inEB = false;
39  }
40  if (ecalPart == barrel) {
42  ecalPart_string = "Barrel";
43  inEB = true;
44  }
45 
46  if (verbosity < pINFO) {
47  std::cout << "-------------------------------------------------------------" << std::endl;
48  std::cout << "Island algorithm invoked for ECAL" << ecalPart_string << std::endl;
49  std::cout << "Looking for seeds, threshold used = " << threshold << " ADC" << std::endl;
50  }
51 
52  int nregions = 0;
53  if (regional)
54  nregions = regions.size();
55 
56  if (!regional || nregions) {
58  for (it = recHits_->begin(); it != recHits_->end(); it++) {
59  if (!uncalibRecHits_) {
60  if (verbosity < pINFO) {
61  std::cout << "-------------------------------------------------------------" << std::endl;
62  std::cout << "No Uncalibrated RecHits no Uncalibrated rec hit collection available" << std::endl;
63  }
64  continue;
65  }
66 
67  // if rechit affected by features other than these, do not allow if seeding
68  // features are hard coded, for now; add severity?
69  uint32_t rhFlag = (*it).recoFlag();
70  if (!(rhFlag == EcalRecHit::kGood || rhFlag == EcalRecHit::kOutOfTime || rhFlag == EcalRecHit::kPoorCalib))
71  continue;
72 
74 
75  if (itt == uncalibRecHits_->end()) {
76  if (verbosity < pINFO) {
77  std::cout << "-------------------------------------------------------------" << std::endl;
78  std::cout << "No Uncalibrated RecHit associated with the RecHit Probably no Uncalibrated rec hit collection "
79  "available"
80  << std::endl;
81  }
82  continue;
83  }
84 
85  EcalUncalibratedRecHit uhit_p = *itt;
86 
87  // looking for cluster seeds
89  continue; //
90 
91  auto thisCell = geometry_p->getGeometry(it->id());
92 
93  // Require that RecHit is within clustering region in case
94  // of regional reconstruction
95  bool withinRegion = false;
96  if (regional) {
97  std::vector<RectangularEtaPhiRegion>::const_iterator region;
98  for (region = regions.begin(); region != regions.end(); region++) {
99  if (region->inRegion(thisCell->etaPos(), thisCell->phiPos())) {
100  withinRegion = true;
101  break;
102  }
103  }
104  }
105 
106  if (!regional || withinRegion) {
107  seeds.push_back(
108  *it); // JHaupt 4-27-2008 Et -> energy, most likely not needed as there is already a threshold requirement.
109  }
110  }
111  }
112 
113  sort(seeds.begin(), seeds.end(), [](auto const &x, auto const &y) { return x.energy() > y.energy(); });
114 
115  if (verbosity < pINFO) {
116  std::cout << "JH Total number of seeds found in event = " << seeds.size() << std::endl;
117  for (EcalRecHitCollection::const_iterator ji = seeds.begin(); ji != seeds.end(); ++ji) {
118  //std::cout << "JH Seed Energy " << ji->energy() << " hashed " << ((EBDetId)ji->id()).hashedIndex() << std::endl;
119  }
120  }
121 
122  mainSearch(geometry_p, topology_p, geometryES_p, ecalPart);
123  std::sort(clusters_v.rbegin(), clusters_v.rend(), isClusterEtLess);
124 
125  if (verbosity < pINFO) {
126  std::cout << "---------- end of main search. clusters have been sorted ----" << std::endl;
127  }
128 
129  return clusters_v;
130 }

References EcalUncalibratedRecHit::amplitude(), barrel, edm::SortedCollection< T, SORT >::begin(), canSeed_s, clusters_v, gather_cfg::cout, ecalBarrelSeedThreshold, ecalEndcapSeedThreshold, edm::SortedCollection< T, SORT >::end(), endcap, edm::SortedCollection< T, SORT >::find(), CaloSubdetectorGeometry::getGeometry(), hfClusterShapes_cfi::hits, inEB, isClusterEtLess(), EcalRecHit::kGood, EcalRecHit::kOutOfTime, EcalRecHit::kPoorCalib, mainSearch(), pINFO, recHits_, HLT_FULL_cff::region, seeds, jetUpdater_cfi::sort, AlCaHLTBitMon_QueryRunRegistry::string, remoteMonitoring_LED_IterMethod_cfg::threshold, uncalibRecHits_, used_s, verbosity, x, and y.

Referenced by CosmicClusterProducer::clusterizeECALPart().

◆ prepareCluster()

void CosmicClusterAlgo::prepareCluster ( CaloNavigator< DetId > &  navigator,
const CaloSubdetectorGeometry geometry 
)
private

Definition at line 309 of file CosmicClusterAlgo.cc.

309  {
310  DetId thisDet;
311  std::set<DetId>::iterator setItr;
312 
313  // now add the 5x5 taking care to mark the edges
314  // as able to seed and where overlapping in the central
315  // region with crystals that were previously able to seed
316  // change their status so they are not able to seed
317  //std::cout << std::endl;
318  for (int dx = -2; dx < 3; ++dx) //for (int dx = -2; dx < 3; ++dx)
319  {
320  for (int dy = -2; dy < 3; ++dy) //for (int dy = -2; dy < 3; ++ dy)
321  {
322  // navigate in free steps forming
323  // a full 5x5
324  thisDet = navigator.offsetBy(dx, dy);
325  navigator.home();
326 
327  // add the current crystal
328  //std::cout << "adding " << dx << ", " << dy << std::endl;
329 
330  // now consider if we are in an edge (outer 16)
331  // or central (inner 9) region
332  if ((abs(dx) > 1) || (abs(dy) > 1)) {
333  // this is an "edge" so should be allowed to seed
334  // provided it is not already used
335  //std::cout << " setting can seed" << std::endl;
336  addCrystal(thisDet, false); //These are in the V25
337  canSeed_s.insert(thisDet);
338  } // end if "edge"
339  else {
340  // or else we are in the central 3x3
341  // and must remove any of these crystals from the canSeed set
342  setItr = canSeed_s.find(thisDet);
343  addCrystal(thisDet, true); //These are in the V9
344  if (setItr != canSeed_s.end()) {
345  //std::cout << " unsetting can seed" << std::endl;
346  canSeed_s.erase(setItr);
347  }
348  } // end if "centre"
349 
350  } // end loop on dy
351 
352  } // end loop on dx
353 
354  //std::cout << "*** " << std::endl;
355  //std::cout << " current_v contains " << current_v.size() << std::endl;
356  //std::cout << "*** " << std::endl;
357 }

References funct::abs(), addCrystal(), canSeed_s, PVValHelper::dx, PVValHelper::dy, and HLT_FULL_cff::navigator.

Referenced by mainSearch().

◆ setVerbosity()

void CosmicClusterAlgo::setVerbosity ( VerbosityLevel  the_verbosity)
inline

Definition at line 58 of file CosmicClusterAlgo.h.

58 { verbosity = the_verbosity; }

References verbosity.

Member Data Documentation

◆ canSeed_s

std::set<DetId> CosmicClusterAlgo::canSeed_s
private

Definition at line 99 of file CosmicClusterAlgo.h.

Referenced by mainSearch(), makeClusters(), and prepareCluster().

◆ clusters_v

std::vector<reco::BasicCluster> CosmicClusterAlgo::clusters_v
private

Definition at line 111 of file CosmicClusterAlgo.h.

Referenced by makeCluster(), and makeClusters().

◆ current_v25

std::vector<DetId> CosmicClusterAlgo::current_v25
private

Definition at line 107 of file CosmicClusterAlgo.h.

Referenced by addCrystal(), mainSearch(), and makeCluster().

◆ current_v25Sup

std::vector<std::pair<DetId, float> > CosmicClusterAlgo::current_v25Sup
private

Definition at line 108 of file CosmicClusterAlgo.h.

Referenced by mainSearch(), and makeCluster().

◆ current_v9

std::vector<DetId> CosmicClusterAlgo::current_v9
private

Definition at line 106 of file CosmicClusterAlgo.h.

Referenced by addCrystal(), mainSearch(), and makeCluster().

◆ ecalBarrelSecondThreshold

double CosmicClusterAlgo::ecalBarrelSecondThreshold
private

Definition at line 81 of file CosmicClusterAlgo.h.

Referenced by makeCluster().

◆ ecalBarrelSeedThreshold

double CosmicClusterAlgo::ecalBarrelSeedThreshold
private

Definition at line 79 of file CosmicClusterAlgo.h.

Referenced by makeClusters().

◆ ecalBarrelSingleThreshold

double CosmicClusterAlgo::ecalBarrelSingleThreshold
private

Definition at line 80 of file CosmicClusterAlgo.h.

Referenced by makeCluster().

◆ ecalBarrelSupThreshold

double CosmicClusterAlgo::ecalBarrelSupThreshold
private

Definition at line 82 of file CosmicClusterAlgo.h.

Referenced by makeCluster().

◆ ecalEndcapSecondThreshold

double CosmicClusterAlgo::ecalEndcapSecondThreshold
private

Definition at line 86 of file CosmicClusterAlgo.h.

Referenced by makeCluster().

◆ ecalEndcapSeedThreshold

double CosmicClusterAlgo::ecalEndcapSeedThreshold
private

Definition at line 84 of file CosmicClusterAlgo.h.

Referenced by makeClusters().

◆ ecalEndcapSingleThreshold

double CosmicClusterAlgo::ecalEndcapSingleThreshold
private

Definition at line 85 of file CosmicClusterAlgo.h.

Referenced by makeCluster().

◆ ecalEndcapSupThreshold

double CosmicClusterAlgo::ecalEndcapSupThreshold
private

Definition at line 87 of file CosmicClusterAlgo.h.

Referenced by makeCluster().

◆ inEB

bool CosmicClusterAlgo::inEB
private

Definition at line 103 of file CosmicClusterAlgo.h.

Referenced by makeCluster(), and makeClusters().

◆ posCalculator_

PositionCalc CosmicClusterAlgo::posCalculator_
private

Definition at line 76 of file CosmicClusterAlgo.h.

Referenced by CosmicClusterAlgo(), and makeCluster().

◆ recHits_

const EcalRecHitCollection* CosmicClusterAlgo::recHits_
private

Definition at line 90 of file CosmicClusterAlgo.h.

Referenced by addCrystal(), checkMaxima(), makeCluster(), and makeClusters().

◆ seeds

std::vector<EcalRecHit> CosmicClusterAlgo::seeds
private

Definition at line 95 of file CosmicClusterAlgo.h.

Referenced by mainSearch(), and makeClusters().

◆ uncalibRecHits_

const EcalUncalibratedRecHitCollection* CosmicClusterAlgo::uncalibRecHits_
private

Definition at line 92 of file CosmicClusterAlgo.h.

Referenced by addCrystal(), makeCluster(), and makeClusters().

◆ used_s

std::set<DetId> CosmicClusterAlgo::used_s
private

Definition at line 98 of file CosmicClusterAlgo.h.

Referenced by addCrystal(), mainSearch(), and makeClusters().

◆ verbosity

VerbosityLevel CosmicClusterAlgo::verbosity
private

Definition at line 114 of file CosmicClusterAlgo.h.

Referenced by mainSearch(), makeCluster(), makeClusters(), and setVerbosity().

CosmicClusterAlgo::clusters_v
std::vector< reco::BasicCluster > clusters_v
Definition: CosmicClusterAlgo.h:111
DDAxes::y
EcalRecHit
Definition: EcalRecHit.h:15
CosmicClusterAlgo::Point
math::XYZPoint Point
point in the space
Definition: CosmicClusterAlgo.h:72
CosmicClusterAlgo::current_v9
std::vector< DetId > current_v9
Definition: CosmicClusterAlgo.h:106
mps_fire.i
i
Definition: mps_fire.py:428
edm::SortedCollection< EcalRecHit >::const_iterator
std::vector< EcalRecHit >::const_iterator const_iterator
Definition: SortedCollection.h:80
CosmicClusterAlgo::mainSearch
void mainSearch(const CaloSubdetectorGeometry *geometry_p, const CaloSubdetectorTopology *topology_p, const CaloSubdetectorGeometry *geometryES_p, EcalPart ecalPart)
Definition: CosmicClusterAlgo.cc:134
hfClusterShapes_cfi.hits
hits
Definition: hfClusterShapes_cfi.py:5
EcalRecHit::id
DetId id() const
get the id
Definition: EcalRecHit.h:77
HLT_FULL_cff.navigator
navigator
Definition: HLT_FULL_cff.py:13126
EcalUncalibratedRecHit::amplitude
float amplitude() const
Definition: EcalUncalibratedRecHit.h:28
geometry
Definition: geometry.py:1
gather_cfg.cout
cout
Definition: gather_cfg.py:144
CosmicClusterAlgo::ecalEndcapSeedThreshold
double ecalEndcapSeedThreshold
Definition: CosmicClusterAlgo.h:84
CosmicClusterAlgo::endcap
Definition: CosmicClusterAlgo.h:29
CosmicClusterAlgo::ecalEndcapSupThreshold
double ecalEndcapSupThreshold
Definition: CosmicClusterAlgo.h:87
DDAxes::x
CosmicClusterAlgo::prepareCluster
void prepareCluster(CaloNavigator< DetId > &navigator, const CaloSubdetectorGeometry *geometry)
Definition: CosmicClusterAlgo.cc:309
EcalRecHit::energy
float energy() const
Definition: EcalRecHit.h:68
CosmicClusterAlgo::ecalBarrelSingleThreshold
double ecalBarrelSingleThreshold
Definition: CosmicClusterAlgo.h:80
CosmicClusterAlgo::ecalBarrelSupThreshold
double ecalBarrelSupThreshold
Definition: CosmicClusterAlgo.h:82
CosmicClusterAlgo::recHits_
const EcalRecHitCollection * recHits_
Definition: CosmicClusterAlgo.h:90
CosmicClusterAlgo::ecalBarrelSecondThreshold
double ecalBarrelSecondThreshold
Definition: CosmicClusterAlgo.h:81
DetId
Definition: DetId.h:17
std::swap
void swap(edm::DataFrameContainer &lhs, edm::DataFrameContainer &rhs)
Definition: DataFrameContainer.h:209
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:262
PositionCalc::Calculate_Location
math::XYZPoint Calculate_Location(const HitsAndFractions &iDetIds, const edm::SortedCollection< HitType > *iRecHits, const CaloSubdetectorGeometry *iSubGeom, const CaloSubdetectorGeometry *iESGeom=nullptr)
Definition: PositionCalc.h:65
reco::CaloCluster::multi5x5
Definition: CaloCluster.h:38
reco::CaloCluster
Definition: CaloCluster.h:31
HCALHighEnergyHPDFilter_cfi.energy
energy
Definition: HCALHighEnergyHPDFilter_cfi.py:5
EcalRecHit::kGood
Definition: EcalRecHit.h:21
CosmicClusterAlgo::ecalEndcapSecondThreshold
double ecalEndcapSecondThreshold
Definition: CosmicClusterAlgo.h:86
CosmicClusterAlgo::seeds
std::vector< EcalRecHit > seeds
Definition: CosmicClusterAlgo.h:95
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
CosmicClusterAlgo::ecalEndcapSingleThreshold
double ecalEndcapSingleThreshold
Definition: CosmicClusterAlgo.h:85
Point
Structure Point Contains parameters of Gaussian fits to DMRs.
Definition: DMRtrends.cc:57
CosmicClusterAlgo::current_v25
std::vector< DetId > current_v25
Definition: CosmicClusterAlgo.h:107
EcalRecHit::kOutOfTime
Definition: EcalRecHit.h:23
CosmicClusterAlgo::verbosity
VerbosityLevel verbosity
Definition: CosmicClusterAlgo.h:114
EcalUncalibratedRecHit::id
DetId id() const
Definition: EcalUncalibratedRecHit.h:38
CosmicClusterAlgo::canSeed_s
std::set< DetId > canSeed_s
Definition: CosmicClusterAlgo.h:99
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:267
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:29
position
static int position[264][3]
Definition: ReadPGInfo.cc:289
HLT_FULL_cff.region
region
Definition: HLT_FULL_cff.py:88271
PVValHelper::dy
Definition: PVValidationHelpers.h:50
CosmicClusterAlgo::checkMaxima
bool checkMaxima(CaloNavigator< DetId > &navigator)
Definition: CosmicClusterAlgo.cc:275
reco::CaloID
Definition: CaloID.h:17
CosmicClusterAlgo::pERROR
Definition: CosmicClusterAlgo.h:30
isClusterEtLess
bool isClusterEtLess(const reco::CaloCluster &x, const reco::CaloCluster &y)
Definition: ClusterEtLess.h:7
EcalUncalibratedRecHit
Definition: EcalUncalibratedRecHit.h:8
CosmicClusterAlgo::used_s
std::set< DetId > used_s
Definition: CosmicClusterAlgo.h:98
edm::SortedCollection::find
iterator find(key_type k)
Definition: SortedCollection.h:240
CosmicClusterAlgo::ecalBarrelSeedThreshold
double ecalBarrelSeedThreshold
Definition: CosmicClusterAlgo.h:79
CaloNavigator
Definition: CaloNavigator.h:7
CosmicClusterAlgo::pWARNING
Definition: CosmicClusterAlgo.h:30
CosmicClusterAlgo::inEB
bool inEB
Definition: CosmicClusterAlgo.h:103
CosmicClusterAlgo::uncalibRecHits_
const EcalUncalibratedRecHitCollection * uncalibRecHits_
Definition: CosmicClusterAlgo.h:92
CosmicClusterAlgo::posCalculator_
PositionCalc posCalculator_
Definition: CosmicClusterAlgo.h:76
EcalRecHit::kPoorCalib
Definition: EcalRecHit.h:26
CosmicClusterAlgo::current_v25Sup
std::vector< std::pair< DetId, float > > current_v25Sup
Definition: CosmicClusterAlgo.h:108
AlignmentPI::regions
regions
Definition: AlignmentPayloadInspectorHelper.h:76
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
remoteMonitoring_LED_IterMethod_cfg.threshold
threshold
Definition: remoteMonitoring_LED_IterMethod_cfg.py:430
CosmicClusterAlgo::pDEBUG
Definition: CosmicClusterAlgo.h:30
CosmicClusterAlgo::makeCluster
void makeCluster(const CaloSubdetectorGeometry *geometry_p, const CaloSubdetectorGeometry *geometryES_p, DetId seedId)
Definition: CosmicClusterAlgo.cc:194
CosmicClusterAlgo::barrel
Definition: CosmicClusterAlgo.h:29
PVValHelper::dx
Definition: PVValidationHelpers.h:49
CosmicClusterAlgo::pINFO
Definition: CosmicClusterAlgo.h:30
CosmicClusterAlgo::addCrystal
void addCrystal(const DetId &det, const bool in9)
Definition: CosmicClusterAlgo.cc:359