CMS 3D CMS Logo
List of all members | Classes | Public Member Functions | Public Attributes | Private Member Functions | Private Attributes
cscdqm::StripClusterFinder Class Reference

Object used to find Strip Clusters. More...

#include <CSCDQM_StripClusterFinder.h>

Classes

struct  C1
 
struct  C2
 
class  Sort
 

Public Member Functions

void DoAction (int layerId, float *cathodes)
 
std::vector< StripCluster > getClusters ()
 
void printClusters (void)
 
 StripClusterFinder (int l, int s, int cf, int st, bool ME11=false)
 

Public Attributes

ClusterLocalMax localMaxTMP
 
std::vector< StripCluster > MEStripClusters
 
std::vector< StripClusterFitData > thePulseHeightMap
 

Private Member Functions

bool FindAndMatch (void)
 
void KillCluster (uint32_t ic1, uint32_t ic2, C1 const &, C2 const &)
 
void Match (void)
 
void RefindMax (void)
 
void SearchBorders (void)
 
void SearchMax (int32_t layerId)
 

Private Attributes

bool is7DCFEBs
 
bool isME11
 
int LayerNmb
 
int StripNmb
 
int TimeSliceNmb
 

Detailed Description

Object used to find Strip Clusters.

Definition at line 23 of file CSCDQM_StripClusterFinder.h.

Constructor & Destructor Documentation

◆ StripClusterFinder()

cscdqm::StripClusterFinder::StripClusterFinder ( int  l,
int  s,
int  cf,
int  st,
bool  ME11 = false 
)

Definition at line 7 of file CSCDQM_StripClusterFinder.cc.

References is7DCFEBs, isME11, MainPageGenerator::l, LayerNmb, Reference_intrackfit_cff::ME11, alignCSCRings::s, StripNmb, and TimeSliceNmb.

7  {
8  //
9  // Options
10  //
11  // fOpt = new CalibOptions();
12  LayerNmb = l;
13  TimeSliceNmb = s;
14  StripNmb = cf * 16;
15  isME11 = ME11;
16  if (cf == 7) {
17  is7DCFEBs = true;
18  isME11 = true;
19  } else {
20  is7DCFEBs = false;
21  }
22  }
Reference_intrackfit_cff.ME11
list ME11
Definition: Reference_intrackfit_cff.py:70

Member Function Documentation

◆ DoAction()

void cscdqm::StripClusterFinder::DoAction ( int  layerId,
float *  cathodes 
)

Definition at line 23 of file CSCDQM_StripClusterFinder.cc.

References cscdqm::StripClusterFitData::bx_, cscdqm::StripClusterFitData::channel_, cscdqm::StripClusterFitData::height_, mps_fire::i, dqmiolumiharvest::j, dqmdumpme::k, MainPageGenerator::l, Match(), MEStripClusters, RefindMax(), SearchBorders(), SearchMax(), jetUpdater_cfi::sort, StripNmb, thePulseHeightMap, and TimeSliceNmb.

Referenced by cscdqm::EventProcessor::processCSC().

23  {
24  int TimeId, StripId;
25  MEStripClusters.clear();
26  StripClusterFitData PulseHeightMapTMP;
27 
28  thePulseHeightMap.clear();
29 
30  // fill
31  //===========================================================================
32 
33  for (StripId = 0; StripId < StripNmb; StripId++) {
34  for (TimeId = 0; TimeId < TimeSliceNmb; TimeId++) {
35  PulseHeightMapTMP.height_[TimeId] = *(Cathodes + StripNmb * (TimeSliceNmb * LayerId + TimeId) + StripId);
36  }
37  PulseHeightMapTMP.bx_ = 0;
38  PulseHeightMapTMP.channel_ = StripId; // was StripId
39  thePulseHeightMap.push_back(PulseHeightMapTMP);
40  }
41  sort(thePulseHeightMap.begin(), thePulseHeightMap.end(), Sort());
42  //===========================================================================
43 
44  if (thePulseHeightMap.empty())
45  return;
46  SearchMax(LayerId);
47  SearchBorders();
48  Match();
49  RefindMax();
50  /*
51  int val;
52  for(i=0;i<MEStripClusters.size();i++){
53  val=MEStripClusters[i].LFTBNDStrip;
54  MEStripClusters[i].LFTBNDStrip=thePulseHeightMap[val].channel_;
55  val=MEStripClusters[i].IRTBNDStrip;
56  MEStripClusters[i].IRTBNDStrip=thePulseHeightMap[val].channel_;
57  for(j=0;j<MEStripClusters[i].localMax.size();j++){
58  val=MEStripClusters[i].localMax[j].Strip;
59  MEStripClusters[i].localMax[j].Strip=thePulseHeightMap[val].channel_;
60  }
61  }
62  */
63 
64  float sumstrip;
65  float sumtime;
66  float sumheight;
67 
68  for (uint32_t i = 0; i < MEStripClusters.size(); i++) {
69  MEStripClusters[i].ClusterPulseMapHeight.clear();
70  for (uint32_t j = 0; j < thePulseHeightMap.size(); j++) {
71  if (thePulseHeightMap[j].channel_ >= MEStripClusters[i].LFTBNDStrip &&
72  thePulseHeightMap[j].channel_ <= MEStripClusters[i].IRTBNDStrip)
73  MEStripClusters[i].ClusterPulseMapHeight.push_back(thePulseHeightMap[j]);
74  }
75  sumstrip = 0;
76  sumtime = 0;
77  sumheight = 0;
78  for (uint32_t k = 0; k < MEStripClusters[i].ClusterPulseMapHeight.size(); k++) {
79  for (int l = 0; l < 16; l++) {
80  sumstrip += MEStripClusters[i].ClusterPulseMapHeight[k].height_[l] *
81  (MEStripClusters[i].ClusterPulseMapHeight[k].channel_ + 1);
82  sumtime += MEStripClusters[i].ClusterPulseMapHeight[k].height_[l] * (l + 1);
83  sumheight += MEStripClusters[i].ClusterPulseMapHeight[k].height_[l];
84  }
85  }
86  if (sumheight) {
87  MEStripClusters[i].Mean[0] = sumstrip / sumheight;
88  MEStripClusters[i].Mean[1] = sumtime / sumheight;
89  }
90  }
91  // printClusters();
92  return;
93  }
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:30
mps_fire.i
i
Definition: mps_fire.py:429
cscdqm::StripClusterFinder::thePulseHeightMap
std::vector< StripClusterFitData > thePulseHeightMap
Definition: CSCDQM_StripClusterFinder.h:28
cscdqm::StripClusterFinder::MEStripClusters
std::vector< StripCluster > MEStripClusters
Definition: CSCDQM_StripClusterFinder.h:35

◆ FindAndMatch()

bool cscdqm::StripClusterFinder::FindAndMatch ( void  )
private

Definition at line 192 of file CSCDQM_StripClusterFinder.cc.

References alignmentValidation::c1, cscdqm::StripClusterFinder::C2::IC2MAX, cscdqm::StripClusterFinder::C2::IC2MIN, cscdqm::StripClusterFinder::C2::JC2MAX, cscdqm::StripClusterFinder::C2::JC2MIN, KillCluster(), and MEStripClusters.

Referenced by Match().

192  {
193  // Find clusters to match
194  for (uint32_t ic1 = 0; ic1 < MEStripClusters.size(); ic1++) {
195  C1 c1;
196  c1.IC1MIN = MEStripClusters[ic1].LFTBNDStrip;
197  c1.IC1MAX = MEStripClusters[ic1].IRTBNDStrip;
198  c1.JC1MIN = MEStripClusters[ic1].LFTBNDTime;
199  c1.JC1MAX = MEStripClusters[ic1].IRTBNDTime;
200  for (uint32_t ic2 = ic1 + 1; ic2 < MEStripClusters.size(); ic2++) {
201  C2 c2;
202  c2.IC2MIN = MEStripClusters[ic2].LFTBNDStrip;
203  c2.IC2MAX = MEStripClusters[ic2].IRTBNDStrip;
204  c2.JC2MIN = MEStripClusters[ic2].LFTBNDTime;
205  c2.JC2MAX = MEStripClusters[ic2].IRTBNDTime;
206  if ((c2.IC2MIN >= c1.IC1MIN && c2.IC2MIN <= c1.IC1MAX && c2.JC2MIN >= c1.JC1MIN && c2.JC2MIN <= c1.JC1MAX) ||
207  (c2.IC2MIN >= c1.IC1MIN && c2.IC2MIN <= c1.IC1MAX && c2.JC2MAX >= c1.JC1MIN && c2.JC2MAX <= c1.JC1MAX) ||
208  (c2.IC2MAX >= c1.IC1MIN && c2.IC2MAX <= c1.IC1MAX && c2.JC2MIN >= c1.JC1MIN && c2.JC2MIN <= c1.JC1MAX) ||
209  (c2.IC2MAX >= c1.IC1MIN && c2.IC2MAX <= c1.IC1MAX && c2.JC2MAX >= c1.JC1MIN && c2.JC2MAX <= c1.JC1MAX)) {
210  KillCluster(ic1, ic2, c1, c2);
211  return true;
212  } else {
213  if ((c1.IC1MIN >= c2.IC2MIN && c1.IC1MIN <= c2.IC2MAX && c1.JC1MIN >= c2.JC2MIN && c1.JC1MIN <= c2.JC2MAX) ||
214  (c1.IC1MIN >= c2.IC2MIN && c1.IC1MIN <= c2.IC2MAX && c1.JC1MAX >= c2.JC2MIN && c1.JC1MAX <= c2.JC2MAX) ||
215  (c1.IC1MAX >= c2.IC2MIN && c1.IC1MAX <= c2.IC2MAX && c1.JC1MIN >= c2.JC2MIN && c1.JC1MIN <= c2.JC2MAX) ||
216  (c1.IC1MAX >= c2.IC2MIN && c1.IC1MAX <= c2.IC2MAX && c1.JC1MAX >= c2.JC2MIN && c1.JC1MAX <= c2.JC2MAX)) {
217  KillCluster(ic1, ic2, c1, c2);
218  return true;
219  }
220  }
221  }
222  }
223  return false;
224  }
cscdqm::StripClusterFinder::KillCluster
void KillCluster(uint32_t ic1, uint32_t ic2, C1 const &, C2 const &)
Definition: CSCDQM_StripClusterFinder.cc:225
cscdqm::StripClusterFinder::MEStripClusters
std::vector< StripCluster > MEStripClusters
Definition: CSCDQM_StripClusterFinder.h:35

◆ getClusters()

std::vector<StripCluster> cscdqm::StripClusterFinder::getClusters ( )
inline

Definition at line 37 of file CSCDQM_StripClusterFinder.h.

References MEStripClusters.

Referenced by cscdqm::EventProcessor::processCSC().

37 { return MEStripClusters; }
cscdqm::StripClusterFinder::MEStripClusters
std::vector< StripCluster > MEStripClusters
Definition: CSCDQM_StripClusterFinder.h:35

◆ KillCluster()

void cscdqm::StripClusterFinder::KillCluster ( uint32_t  ic1,
uint32_t  ic2,
C1 const &  c1,
C2 const &  c2 
)
private

Definition at line 225 of file CSCDQM_StripClusterFinder.cc.

References alignmentValidation::c1, cscdqm::StripClusterFinder::C2::IC2MAX, cscdqm::StripClusterFinder::C2::IC2MIN, cscdqm::StripClusterFinder::C2::JC2MAX, cscdqm::StripClusterFinder::C2::JC2MIN, and MEStripClusters.

Referenced by FindAndMatch().

225  {
226  // Match Clusters and kill one of clusters.
227  if (c1.IC1MIN < c2.IC2MIN)
228  MEStripClusters[ic1].LFTBNDStrip = c1.IC1MIN;
229  else
230  MEStripClusters[ic1].LFTBNDStrip = c2.IC2MIN;
231  if (c1.JC1MIN < c2.JC2MIN)
232  MEStripClusters[ic1].LFTBNDTime = c1.JC1MIN;
233  else
234  MEStripClusters[ic1].LFTBNDTime = c2.JC2MIN;
235  if (c1.IC1MAX > c2.IC2MAX)
236  MEStripClusters[ic1].IRTBNDStrip = c1.IC1MAX;
237  else
238  MEStripClusters[ic1].IRTBNDStrip = c2.IC2MAX;
239  if (c1.JC1MAX > c2.JC2MAX)
240  MEStripClusters[ic1].IRTBNDTime = c1.JC1MAX;
241  else
242  MEStripClusters[ic1].IRTBNDTime = c2.JC2MAX;
243 
244  MEStripClusters.erase(MEStripClusters.begin() + ic2);
245  return;
246  }
cscdqm::StripClusterFinder::MEStripClusters
std::vector< StripCluster > MEStripClusters
Definition: CSCDQM_StripClusterFinder.h:35

◆ Match()

void cscdqm::StripClusterFinder::Match ( void  )
private

Definition at line 182 of file CSCDQM_StripClusterFinder.cc.

References FindAndMatch().

Referenced by DoAction().

182  {
183  // MATCHING THE OVERLAPING CLASTERS
184  bool find2match = true;
185  do {
186  find2match = FindAndMatch();
187  } while (find2match);
188 
189  return;
190  }

◆ printClusters()

void cscdqm::StripClusterFinder::printClusters ( void  )

Definition at line 329 of file CSCDQM_StripClusterFinder.cc.

References gather_cfg::cout, mps_fire::i, dqmiolumiharvest::j, dqmdumpme::k, MainPageGenerator::l, and MEStripClusters.

329  {
330  int iS, jT;
331  std::cout << "====================================================================" << std::endl;
332  std::cout << "debug information from StripClusterFinder" << std::endl;
333  for (uint32_t i = 0; i < MEStripClusters.size(); i++) {
334  if (MEStripClusters[i].localMax.empty())
335  continue;
336  std::cout << " Cluster: " << i + 1 << " Number of local Maximums " << MEStripClusters[i].localMax.size()
337  << std::endl;
338  for (uint32_t j = 0; j < MEStripClusters[i].localMax.size(); j++) {
339  iS = MEStripClusters[i].localMax[j].Strip;
340  jT = MEStripClusters[i].localMax[j].Time;
341 
342  // std::cout << "Local Max: " << j << " Strip: " << iS << " Time: " << jT << std::endl;
343  for (uint32_t k = 0; k < MEStripClusters[i].ClusterPulseMapHeight.size(); k++) {
344  if (MEStripClusters[i].ClusterPulseMapHeight[k].channel_ == iS)
345  std::cout << "Local Max: " << j + 1 << " Strip: " << iS + 1 << " Time: " << jT + 1
346  << " Height: " << MEStripClusters[i].ClusterPulseMapHeight[k].height_[jT] << std::endl;
347  }
348  }
349  for (uint32_t k = 0; k < MEStripClusters[i].ClusterPulseMapHeight.size(); k++) {
350  std::cout << "Strip: " << MEStripClusters[i].ClusterPulseMapHeight[k].channel_ + 1;
351  for (int l = 0; l < 16; l++)
352  std::cout << " " << MEStripClusters[i].ClusterPulseMapHeight[k].height_[l];
353  std::cout << std::endl;
354  }
355 
356  std::cout << " Left Top corner strip: " << MEStripClusters[i].LFTBNDStrip + 1 << " "
357  << " time: " << MEStripClusters[i].LFTBNDTime + 1 << std::endl;
358  std::cout << " Right Bottom corner strip: " << MEStripClusters[i].IRTBNDStrip + 1 << " "
359  << " time: " << MEStripClusters[i].IRTBNDTime + 1 << std::endl;
360  }
361  std::cout << "======================================================================" << std::endl;
362  return;
363  }
mps_fire.i
i
Definition: mps_fire.py:429
cscdqm::StripClusterFinder::MEStripClusters
std::vector< StripCluster > MEStripClusters
Definition: CSCDQM_StripClusterFinder.h:35
gather_cfg.cout
cout
Definition: gather_cfg.py:144

◆ RefindMax()

void cscdqm::StripClusterFinder::RefindMax ( void  )
private

Definition at line 247 of file CSCDQM_StripClusterFinder.cc.

References dqmMemoryStats::float, mps_fire::i, is7DCFEBs, isME11, dqmiolumiharvest::j, localMaxTMP, MEStripClusters, cscdqm::ClusterLocalMax::Strip, thePulseHeightMap, and cscdqm::ClusterLocalMax::Time.

Referenced by DoAction().

247  {
248  // SEARCHING EXTREMUMS IN THE CLUSTERS
249 
250  for (uint32_t i = 0; i < MEStripClusters.size(); i++) {
251  MEStripClusters[i].localMax.clear();
252  int ilocal = 0;
253  int iLS = MEStripClusters[i].LFTBNDStrip;
254  int iRS = MEStripClusters[i].IRTBNDStrip;
255  int iLT = MEStripClusters[i].LFTBNDTime;
256  int iRT = MEStripClusters[i].IRTBNDTime;
257 
258  for (int iS = iLS; iS <= iRS; iS++) {
259  if (isME11 && (thePulseHeightMap[iS].channel_ == 63 || thePulseHeightMap[iS].channel_ == 64))
260  continue;
261  for (int jT = iLT; jT <= iRT; jT++) {
262  if (iS == 0 || jT == 0 || (!is7DCFEBs && (iS == 79)) || (is7DCFEBs && (iS == 111)) || jT == 7)
263  continue;
264  if (thePulseHeightMap[iS].height_[jT] > thePulseHeightMap[iS - 1].height_[jT] &&
265  thePulseHeightMap[iS].height_[jT] > thePulseHeightMap[iS + 1].height_[jT] &&
266  thePulseHeightMap[iS].height_[jT] > thePulseHeightMap[iS].height_[jT - 1] &&
267  thePulseHeightMap[iS].height_[jT] > thePulseHeightMap[iS].height_[jT + 1] &&
268  thePulseHeightMap[iS].height_[jT] > thePulseHeightMap[iS - 1].height_[jT - 1] &&
269  thePulseHeightMap[iS].height_[jT] > thePulseHeightMap[iS - 1].height_[jT + 1] &&
270  thePulseHeightMap[iS].height_[jT] > thePulseHeightMap[iS + 1].height_[jT - 1] &&
271  thePulseHeightMap[iS].height_[jT] > thePulseHeightMap[iS + 1].height_[jT + 1]) {
272  localMaxTMP.Strip = iS;
273  localMaxTMP.Time = jT;
274  MEStripClusters[i].localMax.push_back(localMaxTMP);
275  ilocal++;
276  }
277  }
278  }
279  // kill local maximums rellated to noise, maximums with pulse height less then 10% of Global max of clust.
280  //fing Global Max
281  float GlobalMax = 0;
282  if (!MEStripClusters[i].localMax.empty()) {
283  //std::cout << "Cluster: " << i << " Number of local maximums before erase: "
284  // << MEStripClusters[i].localMax.size() << std::endl;
285  for (uint32_t j = 0; j < MEStripClusters[i].localMax.size(); j++) {
286  int iS = MEStripClusters[i].localMax[j].Strip;
287  int jT = MEStripClusters[i].localMax[j].Time;
288  /*
289  std::cout << "Current Max:"
290  << " " << iS
291  << " " << jT
292  << " " << thePulseHeightMap[iS].height_[jT] << std::endl;
293  */
294  if (thePulseHeightMap[iS].height_[jT] > GlobalMax)
295  GlobalMax = thePulseHeightMap[iS].height_[jT];
296  }
297  GlobalMax = (float)(GlobalMax / 10.);
298  //erase noise localMaximums
299  bool Erased;
300  do {
301  Erased = false;
302  for (uint32_t j = 0; j < MEStripClusters[i].localMax.size(); j++) {
303  int iS = MEStripClusters[i].localMax[j].Strip;
304  int jT = MEStripClusters[i].localMax[j].Time;
305  if (thePulseHeightMap[iS].height_[jT] < GlobalMax) {
306  MEStripClusters[i].localMax.erase(MEStripClusters[i].localMax.begin() + j);
307  Erased = true;
308  break;
309  }
310  }
311  } while (Erased);
312 
313  //debug outs
314  //std::cout << "Cluster: " << i << " Number of local maximums: "
315  // << MEStripClusters[i].localMax.size() << std::endl;
316  /*
317  for(j=0;j<MEStripClusters[i].localMax.size();j++){
318  iS=MEStripClusters[i].localMax[j].Strip;
319  jT=MEStripClusters[i].localMax[j].Time;
320  std::cout << "Local Max: " << j << " Strip: " << iS << " Time: " << jT
321  << " Height: " << thePulseHeightMap[iS].height_[jT]
322  << " Cut Value: " << GlobalMax << std::endl;
323  }
324  */
325  }
326  }
327  return;
328  }
mps_fire.i
i
Definition: mps_fire.py:429
cscdqm::StripClusterFinder::thePulseHeightMap
std::vector< StripClusterFitData > thePulseHeightMap
Definition: CSCDQM_StripClusterFinder.h:28
cscdqm::StripClusterFinder::MEStripClusters
std::vector< StripCluster > MEStripClusters
Definition: CSCDQM_StripClusterFinder.h:35

◆ SearchBorders()

void cscdqm::StripClusterFinder::SearchBorders ( void  )
private

Definition at line 124 of file CSCDQM_StripClusterFinder.cc.

References mps_fire::i, isME11, MEStripClusters, and thePulseHeightMap.

Referenced by DoAction().

124  {
125  uint32_t iS, iT, iL, jL, iR, jR;
126 
127  // SEARCHING PARAMETERS OF THE CLASTERS (LEFT DOWN & RIGHT UP)
128 
129  for (uint32_t i = 0; i < MEStripClusters.size(); i++) {
130  if (MEStripClusters[i].localMax.empty()) {
131  edm::LogWarning("NoLocalMax") << "Cluster " << i << " has no local Maxima";
132  continue;
133  }
134  iS = MEStripClusters[i].localMax[0].Strip;
135  iT = MEStripClusters[i].localMax[0].Time;
136  // LEFT DOWN
137  // strip
138  MEStripClusters[i].LFTBNDStrip = 0;
139  for (iL = iS - 1; iL > 0; iL--) {
140  if (isME11 && (thePulseHeightMap[iL].channel_ == 64)) {
141  MEStripClusters[i].LFTBNDStrip = iL;
142  break;
143  }
144  if (thePulseHeightMap[iL].height_[iT] == 0.) {
145  MEStripClusters[i].LFTBNDStrip = iL + 1;
146  break;
147  }
148  }
149  //time
150  MEStripClusters[i].LFTBNDTime = 0;
151  for (jL = iT - 1; jL > 0; jL--) {
152  if (thePulseHeightMap[iS].height_[jL] == 0.) {
153  MEStripClusters[i].LFTBNDTime = jL + 1;
154  break;
155  }
156  }
157  // RIGHT UP
158  //strip
159  MEStripClusters[i].IRTBNDStrip = thePulseHeightMap.size() - 1;
160  for (iR = iS + 1; iR < thePulseHeightMap.size(); iR++) {
161  if (isME11 && (thePulseHeightMap[iR].channel_ == 63)) {
162  MEStripClusters[i].IRTBNDStrip = iR;
163  break;
164  }
165  if (thePulseHeightMap[iR].height_[iT] == 0.) {
166  MEStripClusters[i].IRTBNDStrip = iR - 1;
167  break;
168  }
169  }
170  //time
171  MEStripClusters[i].IRTBNDTime = 15;
172  for (jR = iT + 1; jR < 16; jR++) {
173  if (thePulseHeightMap[iS].height_[jR] == 0.) {
174  MEStripClusters[i].IRTBNDTime = jR - 1;
175  break;
176  }
177  }
178  }
179  return;
180  }
mps_fire.i
i
Definition: mps_fire.py:429
cscdqm::StripClusterFinder::thePulseHeightMap
std::vector< StripClusterFitData > thePulseHeightMap
Definition: CSCDQM_StripClusterFinder.h:28
cscdqm::StripClusterFinder::MEStripClusters
std::vector< StripCluster > MEStripClusters
Definition: CSCDQM_StripClusterFinder.h:35
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122

◆ SearchMax()

void cscdqm::StripClusterFinder::SearchMax ( int32_t  layerId)
private

Definition at line 95 of file CSCDQM_StripClusterFinder.cc.

References mps_fire::i, cscdqm::StripCluster::IRTBNDStrip, cscdqm::StripCluster::IRTBNDTime, isME11, dqmiolumiharvest::j, cscdqm::StripCluster::LayerId, cscdqm::StripCluster::LFTBNDStrip, cscdqm::StripCluster::LFTBNDTime, cscdqm::StripCluster::localMax, localMaxTMP, MEStripClusters, cscdqm::ClusterLocalMax::Strip, thePulseHeightMap, and cscdqm::ClusterLocalMax::Time.

Referenced by DoAction().

95  {
96  StripCluster tmpCluster;
97  for (uint32_t i = 1; i < (thePulseHeightMap.size() - 1); i++) {
98  if (isME11 && (thePulseHeightMap[i].channel_ == 63 || thePulseHeightMap[i].channel_ == 64))
99  continue;
100  for (uint32_t j = 1; j < 15; j++) {
101  if (thePulseHeightMap[i].height_[j] > thePulseHeightMap[i - 1].height_[j] &&
102  thePulseHeightMap[i].height_[j] > thePulseHeightMap[i + 1].height_[j] &&
103  thePulseHeightMap[i].height_[j] > thePulseHeightMap[i].height_[j - 1] &&
104  thePulseHeightMap[i].height_[j] > thePulseHeightMap[i].height_[j + 1] &&
105  thePulseHeightMap[i].height_[j] > thePulseHeightMap[i - 1].height_[j - 1] &&
106  thePulseHeightMap[i].height_[j] > thePulseHeightMap[i - 1].height_[j + 1] &&
107  thePulseHeightMap[i].height_[j] > thePulseHeightMap[i + 1].height_[j - 1] &&
108  thePulseHeightMap[i].height_[j] > thePulseHeightMap[i + 1].height_[j + 1]) {
109  tmpCluster.localMax.clear();
110  localMaxTMP.Strip = i;
111  localMaxTMP.Time = j;
112  tmpCluster.localMax.push_back(localMaxTMP);
113  tmpCluster.LayerId = layerId;
114  tmpCluster.LFTBNDTime = -100;
115  tmpCluster.LFTBNDStrip = -100;
116  tmpCluster.IRTBNDTime = -100;
117  tmpCluster.IRTBNDStrip = -100;
118  MEStripClusters.push_back(tmpCluster);
119  }
120  }
121  }
122  return;
123  }
mps_fire.i
i
Definition: mps_fire.py:429
cscdqm::StripClusterFinder::thePulseHeightMap
std::vector< StripClusterFitData > thePulseHeightMap
Definition: CSCDQM_StripClusterFinder.h:28
cscdqm::StripClusterFinder::MEStripClusters
std::vector< StripCluster > MEStripClusters
Definition: CSCDQM_StripClusterFinder.h:35

Member Data Documentation

◆ is7DCFEBs

bool cscdqm::StripClusterFinder::is7DCFEBs
private

Definition at line 53 of file CSCDQM_StripClusterFinder.h.

Referenced by RefindMax(), and StripClusterFinder().

◆ isME11

bool cscdqm::StripClusterFinder::isME11
private

Definition at line 54 of file CSCDQM_StripClusterFinder.h.

Referenced by RefindMax(), SearchBorders(), SearchMax(), and StripClusterFinder().

◆ LayerNmb

int cscdqm::StripClusterFinder::LayerNmb
private

Definition at line 46 of file CSCDQM_StripClusterFinder.h.

Referenced by StripClusterFinder().

◆ localMaxTMP

ClusterLocalMax cscdqm::StripClusterFinder::localMaxTMP

Definition at line 36 of file CSCDQM_StripClusterFinder.h.

Referenced by RefindMax(), and SearchMax().

◆ MEStripClusters

std::vector<StripCluster> cscdqm::StripClusterFinder::MEStripClusters

Definition at line 35 of file CSCDQM_StripClusterFinder.h.

Referenced by DoAction(), FindAndMatch(), getClusters(), KillCluster(), printClusters(), RefindMax(), SearchBorders(), and SearchMax().

◆ StripNmb

int cscdqm::StripClusterFinder::StripNmb
private

Definition at line 46 of file CSCDQM_StripClusterFinder.h.

Referenced by DoAction(), and StripClusterFinder().

◆ thePulseHeightMap

std::vector<StripClusterFitData> cscdqm::StripClusterFinder::thePulseHeightMap

Definition at line 28 of file CSCDQM_StripClusterFinder.h.

Referenced by DoAction(), RefindMax(), SearchBorders(), and SearchMax().

◆ TimeSliceNmb

int cscdqm::StripClusterFinder::TimeSliceNmb
private

Definition at line 46 of file CSCDQM_StripClusterFinder.h.

Referenced by DoAction(), and StripClusterFinder().

Generated for CMSSW Reference Manual by doxygen 1.8.14