Phase2TrackerClusterizerAlgorithm.cc

Go to the documentation of this file.

#include "Phase2TrackerClusterizerAlgorithm.h"

#include "Geometry/CommonTopologies/interface/PixelTopology.h"

/*
 * Initialise the clusterizer algorithm
 */

Phase2TrackerClusterizerAlgorithm::Phase2TrackerClusterizerAlgorithm(unsigned int maxClusterSize,
                                                                     unsigned int maxNumberClusters)
    : maxClusterSize_(maxClusterSize), maxNumberClusters_(maxNumberClusters), nrows_(0), ncols_(0) {}

/* 
 * Change the size of the 2D matrix for this module (varies from pixel to strip modules)
 */

void Phase2TrackerClusterizerAlgorithm::setup(const PixelGeomDetUnit* pixDet) {
  const PixelTopology& topol(pixDet->specificTopology());
  nrows_ = topol.nrows();
  ncols_ = topol.ncolumns();
  matrix_ = decltype(matrix_)(nrows_, ncols_);
}

/* 
 * Go over the Digis and create clusters
 */

void Phase2TrackerClusterizerAlgorithm::clusterizeDetUnit(const edm::DetSet<Phase2TrackerDigi>& digis,
                                                          Phase2TrackerCluster1DCollectionNew::FastFiller& clusters) {
  // Fill the 2D matrix with the hit information : (hit or not)
  fillMatrix(digis.begin(), digis.end());

  // Number of clusters
  unsigned int numberClusters(0);
  Phase2TrackerDigi firstDigi;
  unsigned int sizeCluster(0);
  bool closeCluster(false);
  bool HIPbit(false);

  // Loop over the Digis
  // for the S modules, 1 column = 1 strip, so adjacent digis are along the rows
  // same for P modules
  for (unsigned int col(0); col < ncols_; ++col) {
    for (unsigned int row(0); row < nrows_; ++row) {
      // If the Digi is hit
      if (matrix_(row, col)) {
        // No cluster is open, create a new one
        if (sizeCluster == 0) {
          // Define first digi
          firstDigi = Phase2TrackerDigi(row, col);
          sizeCluster = 1;
        }
        // A cluster is open, increase its size
        else
          ++sizeCluster;
        // Check if we reached the maximum size of the cluster and need to close it
        closeCluster = ((maxClusterSize_ != 0 and sizeCluster >= maxClusterSize_) ? true : false);
      }
      // Otherwise check if we need to close a cluster (end of cluster)
      else
        closeCluster = ((sizeCluster != 0) ? true : false);

      // update the HIP bit
      HIPbit |= (matrix_(row, col) == 2);

      // Always close a cluster if we reach the end of the loop
      if (sizeCluster != 0 and row == (nrows_ - 1))
        closeCluster = true;

      // If we have to close a cluster, do it
      if (closeCluster) {
        // Add the cluster to the list
        clusters.push_back(Phase2TrackerCluster1D(firstDigi, sizeCluster, HIPbit));
        // Reset the variables
        sizeCluster = 0;
        // Increase the number of clusters
        ++numberClusters;
        HIPbit = false;
      }

      // Check if we hit the maximum number of clusters per module
      if (maxNumberClusters_ != 0 and numberClusters > maxNumberClusters_)
        return;
    }
  }

  // Reset the matrix
  clearMatrix(digis.begin(), digis.end());
}

/* 
 * Copy the value of the Digis to the 2D matrix (hit or not). 
 */

void Phase2TrackerClusterizerAlgorithm::fillMatrix(edm::DetSet<Phase2TrackerDigi>::const_iterator begin,
                                                   edm::DetSet<Phase2TrackerDigi>::const_iterator end) {
  for (edm::DetSet<Phase2TrackerDigi>::const_iterator di(begin); di != end; ++di)
    matrix_.set(di->row(), di->column(), true, di->overThreshold());
}

/*
 * Clear the array of hits
 */

void Phase2TrackerClusterizerAlgorithm::clearMatrix(edm::DetSet<Phase2TrackerDigi>::const_iterator begin,
                                                    edm::DetSet<Phase2TrackerDigi>::const_iterator end) {
  for (edm::DetSet<Phase2TrackerDigi>::const_iterator di(begin); di != end; ++di)
    matrix_.set(di->row(), di->column(), false, false);
}