@@ -133,22 +133,24 @@ void Clusterizer::findClusters(const gsl::span<const Digit>& digits, std::vector
133133void Clusterizer::addDigitToCluster (Cluster& cluster, int row, int col, const gsl::span<const Digit>& digits)
134134{
135135 auto & geo = Geometry::instance ();
136- if (row < 0 || row >= geo.getNrows () || col < 0 || col >= geo.getNcols ())
136+ if (row < 0 || row >= geo.getNrows () || col < 0 || col >= geo.getNcols ()) {
137137 return ;
138+ }
138139 int digitIndex = mDigitIndices [row][col];
139140 LOGP (debug, " checking row={} and col={} digitIndex={}" , row, col, digitIndex);
140- if (digitIndex < 0 )
141+ if (digitIndex < 0 ) {
141142 return ;
142-
143+ }
143144 const Digit& digit = digits[digitIndex];
144145 if (cluster.getMultiplicity () > 0 ) {
145146 // check if new digit is in the same chamber and sector
146147 const Digit& digit2 = digits[cluster.getDigitIndex (0 )];
147148 auto [sector1, ch1] = geo.getSectorChamber (digit.getTower ());
148149 auto [sector2, ch2] = geo.getSectorChamber (digit2.getTower ());
149150 LOGP (debug, " checking if sector and chamber are the same: ({},{}) ({},{})" , sector1, ch1, sector2, ch2);
150- if (sector1 != sector2 || ch1 != ch2)
151+ if (sector1 != sector2 || ch1 != ch2) {
151152 return ;
153+ }
152154 }
153155
154156 mDigitIndices [row][col] = -1 ;
@@ -179,11 +181,13 @@ void Clusterizer::makeClusters(const gsl::span<const Digit>& digits, std::vector
179181 auto [row, col] = geo.globalRowColFromIndex (digit.getTower ());
180182 bool isCrystal = geo.isCrystal (digit.getTower ());
181183 if (isCrystal) {
182- if (digit.getEnergy () < mCrystalDigitThreshold )
184+ if (digit.getEnergy () < mCrystalDigitThreshold ) {
183185 continue ;
186+ }
184187 } else {
185- if (digit.getEnergy () < mSamplingDigitThreshold )
188+ if (digit.getEnergy () < mSamplingDigitThreshold ) {
186189 continue ;
190+ }
187191 }
188192 mDigitIndices [row][col] = i;
189193 }
@@ -192,10 +196,12 @@ void Clusterizer::makeClusters(const gsl::span<const Digit>& digits, std::vector
192196 for (int i = 0 ; i < nDigits; i++) {
193197 const Digit& digitSeed = digits[i];
194198 auto [row, col] = geo.globalRowColFromIndex (digitSeed.getTower ());
195- if (mDigitIndices [row][col] < 0 )
199+ if (mDigitIndices [row][col] < 0 ) {
196200 continue ; // digit was already added in one of the clusters
197- if (digitSeed.getEnergy () < mClusteringThreshold )
201+ }
202+ if (digitSeed.getEnergy () < mClusteringThreshold ) {
198203 continue ;
204+ }
199205 LOGP (debug, " starting new cluster at row={} and col={}" , row, col);
200206 auto & cluster = clusters.emplace_back ();
201207 addDigitToCluster (cluster, row, col, digits);
@@ -382,8 +388,9 @@ void Clusterizer::evalClusters(std::vector<Cluster>& clusters)
382388 double xi, yi, zi;
383389 geo.detIdToGlobalPosition (towerId, xi, yi, zi);
384390 double r = std::sqrt ((x - xi) * (x - xi) + (y - yi) * (y - yi) + (z - zi) * (z - zi));
385- if (r > 2.2 )
391+ if (r > 2.2 ) {
386392 continue ;
393+ }
387394 double frac = fCrystalShowerShape ->Eval (r);
388395 double rms = fCrystalRMS ->Eval (r);
389396 chi2 += std::pow ((energy / ee - frac) / rms, 2 .);
@@ -409,10 +416,10 @@ void Clusterizer::evalClusters(std::vector<Cluster>& clusters)
409416 bool isEdge = 0 ;
410417 for (size_t i = 0 ; i < cluster.getMultiplicity (); i++) {
411418 int towerId = cluster.getDigitTowerId (i);
412- if (! geo.isAtTheEdge (towerId))
413- continue ;
414- isEdge = 1 ;
415- break ;
419+ if (geo.isAtTheEdge (towerId)) {
420+ isEdge = 1 ;
421+ break ;
422+ }
416423 }
417424 cluster.setEdgeFlag (isEdge);
418425
@@ -434,8 +441,9 @@ int Clusterizer::getNumberOfLocalMax(Cluster& clu, int* maxAt, float* maxAtEnerg
434441 for (int i = 0 ; i < n; i++) {
435442 isLocalMax[i] = false ;
436443 float en1 = clu.getDigitEnergy (i);
437- if (en1 > mClusteringThreshold )
444+ if (en1 > mClusteringThreshold ) {
438445 isLocalMax[i] = true ;
446+ }
439447 }
440448
441449 for (int i = 0 ; i < n; i++) {
0 commit comments