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;
}