TTCluster.cc

Go to the documentation of this file.

#include "DataFormats/L1TrackTrigger/interface/TTCluster.h"

template< >
unsigned int TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findWidth() const
{
  int rowMin = 99999999;
  int rowMax = 0;
  for ( unsigned int i = 0; i < theHits.size(); i++ )
  {
    int row = 0;
    if ( this->getRows().empty() )
    {
      row = theHits[i]->row();
    }
    else
    {
      row = this->getRows()[i];
    }
    if ( row < rowMin )
      rowMin = row;
    if ( row > rowMax )
      rowMax = row;
  }
  return abs( rowMax - rowMin + 1 ); 
}

template< >
MeasurementPoint TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findHitLocalCoordinates( unsigned int hitIdx ) const
{
  if ( this->getRows().empty() || this->getCols().empty() )
  {
    MeasurementPoint mp( theHits[hitIdx]->row(), theHits[hitIdx]->column() );
    return mp;
  }
  else
  {
    int row = this->getRows()[hitIdx];
    int col = this->getCols()[hitIdx];
    MeasurementPoint mp( row, col );
    return mp;
  }
}

template< >
MeasurementPoint TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findAverageLocalCoordinates() const
{
  double averageCol = 0.0;
  double averageRow = 0.0;

  if ( !theHits.empty() )
  {
    if ( this->getRows().empty() || this->getCols().empty() )
    {
      typename std::vector< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::const_iterator hitIter;
      for ( hitIter = theHits.begin();
            hitIter != theHits.end();
            hitIter++ )
      {
        averageCol += (*hitIter)->column();
        averageRow += (*hitIter)->row();
      }
      averageCol /= theHits.size();
      averageRow /= theHits.size();
    }
    else
    {
      for ( unsigned int j = 0; j < theHits.size(); j++ )
      {
        averageCol += theCols[j];
        averageRow += theRows[j];
      }
      averageCol /= theHits.size();
      averageRow /= theHits.size();
    }
  }
  return MeasurementPoint( averageRow, averageCol );
}

template< >
MeasurementPoint TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findAverageLocalCoordinatesCentered() const
{
  double averageCol = 0.0;
  double averageRow = 0.0;

  if ( !theHits.empty() )
  {
    if ( this->getRows().empty() || this->getCols().empty() )
    {
      typename std::vector< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::const_iterator hitIter;
      for ( hitIter = theHits.begin();
            hitIter != theHits.end();
            hitIter++ )
      {
        averageCol += (*hitIter)->column()+0.5;
        averageRow += (*hitIter)->row()+0.5;
      }
      averageCol /= theHits.size();
      averageRow /= theHits.size();
    }
    else
    {
      for ( unsigned int j = 0; j < theHits.size(); j++ )
      {
        averageCol += theCols[j]+0.5;
        averageRow += theRows[j]+0.5;
      }
      averageCol /= theHits.size();
      averageRow /= theHits.size();
    }
  }
  return MeasurementPoint( averageRow, averageCol );
}

template< >
std::vector< int > TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findRows() const
{
  std::vector< int > temp;
  temp.reserve(theHits.size());
  for ( unsigned int i = 0; i < theHits.size(); i++ )
  {
    temp.push_back( theHits[i]->row() );
  }
  return temp;
}

template< >
std::vector< int > TTCluster< edm::Ref< edm::DetSetVector< Phase2TrackerDigi >, Phase2TrackerDigi > >::findCols() const
{
  std::vector< int > temp;
  temp.reserve(theHits.size());
  for ( unsigned int i = 0; i < theHits.size(); i++ )
  {
    temp.push_back( theHits[i]->column() );
  }
  return temp;
}