MTDThresholdClusterizer Class Reference

An explicit threshold-based clustering algorithm. More...

#include <MTDThresholdClusterizer.h>

Inheritance diagram for MTDThresholdClusterizer:

Public Member Functions
void	clusterize (const FTLRecHitCollection &input, const MTDGeometry *geom, const MTDTopology *topo, FTLClusterCollection &output) override
	Cluster hits. This method operates on a matrix of hits and finds the largest contiguous cluster around each seed hit. Input and output data stored in DetSet. More...
	MTDThresholdClusterizer (edm::ParameterSet const &conf)
	Constructor: More...
	~MTDThresholdClusterizer () override
Public Member Functions inherited from MTDClusterizerBase
virtual	~MTDClusterizerBase ()

Static Public Member Functions
static void	fillDescriptions (edm::ParameterSetDescription &desc)

Private Member Functions
void	clear_buffer (RecHitIterator itr)
	Clear the internal buffer array. More...
void	copy_to_buffer (RecHitIterator itr)
	Copy FTLRecHit into the buffer, identify seeds. More...
FTLCluster	make_cluster (const FTLCluster::FTLHitPos &hit)
	The actual clustering algorithm: group the neighboring hits around the seed. More...
bool	setup (const MTDGeometry *geometry, const MTDTopology *topo, const DetId &id)

Private Attributes
bool	bufferAlreadySet
MTDArrayBuffer	theBuffer
	Data storage. More...
std::vector< FTLCluster >	theClusters
float	theClusterThreshold
DetId	theCurrentId
float	theHitThreshold
	Clustering-related quantities: More...
int	theNumOfCols
int	theNumOfRows
	Geometry-related information. More...
std::vector< FTLCluster::FTLHitPos >	theSeeds
float	theSeedThreshold
float	theTimeThreshold

Additional Inherited Members
Public Types inherited from MTDClusterizerBase
typedef FTLClusterCollection::const_iterator	ClusterIterator
typedef FTLRecHitCollection::const_iterator	RecHitIterator

Detailed Description

An explicit threshold-based clustering algorithm.

A threshold-based clustering algorithm which clusters FTLRecHits into FTLClusters for each DetUnit. The algorithm is straightforward and purely topological: the clustering process starts with seed hits and continues by adding adjacent hits above the hit threshold. Once the cluster is made, it has to be above the cluster threshold as well.

The clusterization is performed on a matrix with size equal to the size of the MTD detector, each cell containing the cahrge and time of the corresponding hit The matrix is reset after each clusterization.

The search starts from seed hits, i.e. hits with sufficiently large amplitudes

FTLCluster contains a barycenter, but it should be noted that that information is largely useless. One must use a PositionEstimator class to compute the RecHit position and its error for every given cluster.

Sets the MTDArrayBuffer dimensions and pixel thresholds. Makes clusters and stores them in theCache if the option useCache has been set.

Definition at line 51 of file MTDThresholdClusterizer.h.

Constructor & Destructor Documentation

MTDThresholdClusterizer()

MTDThresholdClusterizer::MTDThresholdClusterizer

(

edm::ParameterSet const &

conf

)

Constructor:

Definition at line 32 of file MTDThresholdClusterizer.cc.

    :  // Get energy thresholds
      theHitThreshold(conf.getParameter<double>("HitThreshold")),
      theSeedThreshold(conf.getParameter<double>("SeedThreshold")),
      theClusterThreshold(conf.getParameter<double>("ClusterThreshold")),
      theTimeThreshold(conf.getParameter<double>("TimeThreshold")),
      theNumOfRows(0),
      theNumOfCols(0),
      theCurrentId(0),
      theBuffer(theNumOfRows, theNumOfCols),
      bufferAlreadySet(false) {}

~MTDThresholdClusterizer()

MTDThresholdClusterizer::~MTDThresholdClusterizer ( )

override

Definition at line 44 of file MTDThresholdClusterizer.cc.

{}

Member Function Documentation

clear_buffer()

void MTDThresholdClusterizer::clear_buffer ( RecHitIterator  itr )

private

Clear the internal buffer array.

MTDs which are not part of recognized clusters are NOT ERASED during the cluster finding. Erase them now.

Definition at line 190 of file MTDThresholdClusterizer.cc.

{ theBuffer.clear(itr->row(), itr->column()); }

References MTDArrayBuffer::clear(), and theBuffer.

Referenced by clusterize().

clusterize()

void MTDThresholdClusterizer::clusterize ( const FTLRecHitCollection &  input, const MTDGeometry *  geom, const MTDTopology *  topo, FTLClusterCollection &  output  )

overridevirtual

Cluster hits. This method operates on a matrix of hits and finds the largest contiguous cluster around each seed hit. Input and output data stored in DetSet.

Implements MTDClusterizerBase.

Definition at line 94 of file MTDThresholdClusterizer.cc.

                                                                       {
  FTLRecHitCollection::const_iterator begin = input.begin();
  FTLRecHitCollection::const_iterator end = input.end();
 
  // Do not bother for empty detectors
  if (begin == end) {
    edm::LogInfo("MTDThresholdClusterizer") << "No hits to clusterize";
    return;
  }
 
  DEBUG("Input collection " << input.size());
  assert(output.empty());
 
  std::set<unsigned> geoIds;
  std::multimap<unsigned, unsigned> geoIdToIdx;
 
  unsigned index = 0;
  for (const auto& hit : input) {
    MTDDetId mtdId = MTDDetId(hit.detid());
    if (mtdId.subDetector() != MTDDetId::FastTime) {
      throw cms::Exception("MTDThresholdClusterizer")
          << "MTDDetId: " << std::hex << mtdId.rawId() << " is invalid!" << std::dec << std::endl;
    }
 
    if (mtdId.mtdSubDetector() == MTDDetId::BTL) {
      BTLDetId hitId(hit.detid());
      //for BTL topology gives different layout id
      DetId geoId = hitId.geographicalId(MTDTopologyMode::crysLayoutFromTopoMode(topo->getMTDTopologyMode()));
      geoIdToIdx.emplace(geoId, index);
      geoIds.emplace(geoId);
      ++index;
    } else if (mtdId.mtdSubDetector() == MTDDetId::ETL) {
      ETLDetId hitId(hit.detid());
      DetId geoId = hitId.geographicalId();
      geoIdToIdx.emplace(geoId, index);
      geoIds.emplace(geoId);
      ++index;
    } else {
      throw cms::Exception("MTDThresholdClusterizer")
          << "MTDDetId: " << std::hex << mtdId.rawId() << " is invalid!" << std::dec << std::endl;
    }
  }
 
  //cluster hits within geoIds (modules)
  for (unsigned id : geoIds) {
    //  Set up the clusterization on this DetId.
    if (!setup(geom, topo, DetId(id)))
      return;
 
    auto range = geoIdToIdx.equal_range(id);
    DEBUG("Matching Ids for " << std::hex << id << std::dec << " [" << range.first->second << ","
                              << range.second->second << "]");
 
    //  Copy MTDRecHits to the buffer array; select the seed hits
    //  on the way, and store them in theSeeds.
    for (auto itr = range.first; itr != range.second; ++itr) {
      const unsigned hitidx = itr->second;
      copy_to_buffer(begin + hitidx);
    }
 
    FTLClusterCollection::FastFiller clustersOnDet(output, id);
 
    for (unsigned int i = 0; i < theSeeds.size(); i++) {
      if (theBuffer.energy(theSeeds[i]) > theSeedThreshold) {  // Is this seed still valid?
        //  Make a cluster around this seed
        const FTLCluster& cluster = make_cluster(theSeeds[i]);
 
        //  Check if the cluster is above threshold
        if (cluster.energy() > theClusterThreshold) {
          DEBUG("putting in this cluster " << i << " #hits:" << cluster.size() << " E:" << cluster.energy()
                                           << " T:" << cluster.time() << " X:" << cluster.x() << " Y:" << cluster.y());
          clustersOnDet.push_back(cluster);
        }
      }
    }
 
    // Erase the seeds.
    theSeeds.clear();
    //  Need to clean unused hits from the buffer array.
    for (auto itr = range.first; itr != range.second; ++itr) {
      const unsigned hitidx = itr->second;
      clear_buffer(begin + hitidx);
    }
  }
}

References cms::cuda::assert(), begin, MTDDetId::BTL, clear_buffer(), copy_to_buffer(), MTDTopologyMode::crysLayoutFromTopoMode(), DEBUG, TauDecayModes::dec, end, MTDArrayBuffer::energy(), FTLCluster::energy(), MTDDetId::ETL, Exception, MTDDetId::FastTime, relativeConstraints::geom, MTDTopology::getMTDTopologyMode(), mps_fire::i, input, make_cluster(), MTDDetId::mtdSubDetector(), convertSQLitetoXML_cfg::output, edmNew::DetSetVector< T >::FastFiller::push_back(), FastTimerService_cff::range, DetId::rawId(), setup(), FTLCluster::size(), MTDDetId::subDetector(), theBuffer, theClusterThreshold, theSeeds, theSeedThreshold, FTLCluster::time(), FTLCluster::x(), and FTLCluster::y().

copy_to_buffer()

void MTDThresholdClusterizer::copy_to_buffer ( RecHitIterator  itr )

private

Copy FTLRecHit into the buffer, identify seeds.

Definition at line 195 of file MTDThresholdClusterizer.cc.

                                                               {
  int row = itr->row();
  int col = itr->column();
  float energy = itr->energy();
  float time = itr->time();
  float timeError = itr->timeError();
 
  DEBUG("ROW " << row << " COL " << col << " ENERGY " << energy << " TIME " << time);
  if (energy > theHitThreshold) {
    theBuffer.set(row, col, energy, time, timeError);
    if (energy > theSeedThreshold)
      theSeeds.push_back(FTLCluster::FTLHitPos(row, col));
    //sort seeds?
  }
}

References cuy::col, DEBUG, HCALHighEnergyHPDFilter_cfi::energy, MTDArrayBuffer::set(), theBuffer, theHitThreshold, theSeeds, theSeedThreshold, and ntuplemaker::time.

Referenced by clusterize().

fillDescriptions()

void MTDThresholdClusterizer::fillDescriptions ( edm::ParameterSetDescription &  desc )

static

Definition at line 47 of file MTDThresholdClusterizer.cc.

                                                                               {
  desc.add<double>("HitThreshold", 0.);
  desc.add<double>("SeedThreshold", 0.);
  desc.add<double>("ClusterThreshold", 0.);
  desc.add<double>("TimeThreshold", 10.);
}

References edm::ParameterSetDescription::add().

Referenced by MTDClusterProducer::fillDescriptions().

make_cluster()

FTLCluster MTDThresholdClusterizer::make_cluster ( const FTLCluster::FTLHitPos &  hit )

private

The actual clustering algorithm: group the neighboring hits around the seed.

Definition at line 214 of file MTDThresholdClusterizer.cc.

                                                                               {
  //First we acquire the seeds for the clusters
  float seed_energy = theBuffer.energy(hit.row(), hit.col());
  float seed_time = theBuffer.time(hit.row(), hit.col());
  float seed_time_error = theBuffer.time_error(hit.row(), hit.col());
  theBuffer.clear(hit);
 
  AccretionCluster acluster;
  acluster.add(hit, seed_energy, seed_time, seed_time_error);
 
  bool stopClus = false;
  //Here we search all hits adjacent to all hits in the cluster.
  while (!acluster.empty() && !stopClus) {
    //This is the standard algorithm to find and add a hit
    auto curInd = acluster.top();
    acluster.pop();
    for (auto c = std::max(0, int(acluster.y[curInd] - 1));
         c < std::min(int(acluster.y[curInd] + 2), int(theBuffer.columns())) && !stopClus;
         ++c) {
      for (auto r = std::max(0, int(acluster.x[curInd] - 1));
           r < std::min(int(acluster.x[curInd] + 2), int(theBuffer.rows())) && !stopClus;
           ++r) {
        if (theBuffer.energy(r, c) > theHitThreshold) {
          if (std::abs(theBuffer.time(r, c) - seed_time) >
              theTimeThreshold *
                  sqrt(theBuffer.time_error(r, c) * theBuffer.time_error(r, c) + seed_time_error * seed_time_error))
            continue;
          FTLCluster::FTLHitPos newhit(r, c);
          if (!acluster.add(newhit, theBuffer.energy(r, c), theBuffer.time(r, c), theBuffer.time_error(r, c))) {
            stopClus = true;
            break;
          }
          theBuffer.clear(newhit);
        }
      }
    }
  }  // while accretion
 
  FTLCluster cluster(theCurrentId,
                     acluster.isize,
                     acluster.energy.data(),
                     acluster.time.data(),
                     acluster.timeError.data(),
                     acluster.x.data(),
                     acluster.y.data(),
                     acluster.xmin,
                     acluster.ymin);
  return cluster;
}

References funct::abs(), MTDClusterizerBase::AccretionCluster::add(), HltBtagPostValidation_cff::c, MTDArrayBuffer::clear(), MTDArrayBuffer::columns(), MTDClusterizerBase::AccretionCluster::empty(), MTDArrayBuffer::energy(), MTDClusterizerBase::AccretionCluster::energy, MTDClusterizerBase::AccretionCluster::isize, SiStripPI::max, min(), MTDClusterizerBase::AccretionCluster::pop(), alignCSCRings::r, MTDArrayBuffer::rows(), mathSSE::sqrt(), theBuffer, theCurrentId, theHitThreshold, theTimeThreshold, MTDArrayBuffer::time(), MTDClusterizerBase::AccretionCluster::time, MTDArrayBuffer::time_error(), MTDClusterizerBase::AccretionCluster::timeError, MTDClusterizerBase::AccretionCluster::top(), MTDClusterizerBase::AccretionCluster::x, MTDClusterizerBase::AccretionCluster::xmin, MTDClusterizerBase::AccretionCluster::y, and MTDClusterizerBase::AccretionCluster::ymin.

Referenced by clusterize().

setup()

bool MTDThresholdClusterizer::setup ( const MTDGeometry *  geom, const MTDTopology *  topo, const DetId &  id  )

private

Prepare the Clusterizer to work on a particular DetUnit. Re-init the size of the panel/plaquette (so update nrows and ncols),

Definition at line 58 of file MTDThresholdClusterizer.cc.

                                                                                                     {
  theCurrentId = id;
  //using geopraphicalId here
  const auto& thedet = geom->idToDet(id);
  if (thedet == nullptr) {
    throw cms::Exception("MTDThresholdClusterizer")
        << "GeographicalID: " << std::hex << id.rawId() << " is invalid!" << std::dec << std::endl;
    return false;
  }
  const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
  const RectangularMTDTopology& topol = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());
 
  // Get the new sizes.
  unsigned int nrows = topol.nrows();     // rows in x
  unsigned int ncols = topol.ncolumns();  // cols in y
 
  theNumOfRows = nrows;  // Set new sizes
  theNumOfCols = ncols;
 
  DEBUG("Buffer size [" << theNumOfRows + 1 << "," << theNumOfCols + 1 << "]");
 
  if (nrows > theBuffer.rows() || ncols > theBuffer.columns()) {  // change only when a larger is needed
    // Resize the buffer
    theBuffer.setSize(nrows + 1, ncols + 1);  // +1 needed for MTD
    bufferAlreadySet = true;
  }
 
  return true;
}

References bufferAlreadySet, MTDArrayBuffer::columns(), DEBUG, TauDecayModes::dec, Exception, relativeConstraints::geom, triggerObjects_cff::id, hgcalPlots::ncols, RectangularMTDTopology::nrows(), MTDArrayBuffer::rows(), MTDArrayBuffer::setSize(), ProxyMTDTopology::specificTopology(), theBuffer, theCurrentId, theNumOfCols, and theNumOfRows.

Referenced by clusterize().

Member Data Documentation

bufferAlreadySet

bool MTDThresholdClusterizer::bufferAlreadySet

private

Definition at line 82 of file MTDThresholdClusterizer.h.

Referenced by setup().

theBuffer

MTDArrayBuffer MTDThresholdClusterizer::theBuffer

private

Data storage.

Definition at line 81 of file MTDThresholdClusterizer.h.

Referenced by clear_buffer(), clusterize(), copy_to_buffer(), make_cluster(), and setup().

theClusters

std::vector<FTLCluster> MTDThresholdClusterizer::theClusters

private

Definition at line 66 of file MTDThresholdClusterizer.h.

theClusterThreshold

float MTDThresholdClusterizer::theClusterThreshold

private

Definition at line 71 of file MTDThresholdClusterizer.h.

Referenced by clusterize().

theCurrentId

DetId MTDThresholdClusterizer::theCurrentId

private

Definition at line 78 of file MTDThresholdClusterizer.h.

Referenced by make_cluster(), and setup().

theHitThreshold

float MTDThresholdClusterizer::theHitThreshold

private

Clustering-related quantities:

Definition at line 69 of file MTDThresholdClusterizer.h.

Referenced by copy_to_buffer(), and make_cluster().

theNumOfCols

int MTDThresholdClusterizer::theNumOfCols

private

Definition at line 76 of file MTDThresholdClusterizer.h.

Referenced by setup().

theNumOfRows

int MTDThresholdClusterizer::theNumOfRows

private

Geometry-related information.

Definition at line 75 of file MTDThresholdClusterizer.h.

Referenced by setup().

theSeeds

std::vector<FTLCluster::FTLHitPos> MTDThresholdClusterizer::theSeeds

private

Definition at line 65 of file MTDThresholdClusterizer.h.

Referenced by clusterize(), and copy_to_buffer().

theSeedThreshold

float MTDThresholdClusterizer::theSeedThreshold

private

Definition at line 70 of file MTDThresholdClusterizer.h.

Referenced by clusterize(), and copy_to_buffer().

theTimeThreshold

float MTDThresholdClusterizer::theTimeThreshold

private

Definition at line 72 of file MTDThresholdClusterizer.h.

Referenced by make_cluster().