#include <SiStripCluster.h>

Public Types
typedef std::vector < SiStripDigi > ::const_iterator	SiStripDigiIter
typedef std::pair < SiStripDigiIter, SiStripDigiIter >	SiStripDigiRange
Public Member Functions
const std::vector< uint8_t > &	amplitudes () const
float	barycenter () const
uint16_t	firstStrip () const
uint32_t	geographicalId () const
	SiStripCluster (uint32_t detid, const SiStripDigiRange &range)
	SiStripCluster (const uint32_t &detid, const uint16_t &firstStrip, std::vector< uint16_t >::const_iterator begin, std::vector< uint16_t >::const_iterator end)
	SiStripCluster ()
Private Attributes
std::vector< uint8_t >	amplitudes_
uint32_t	detId_
uint16_t	firstStrip_

Detailed Description

Definition at line 8 of file SiStripCluster.h.

Member Typedef Documentation

typedef std::vector<SiStripDigi>::const_iterator SiStripCluster::SiStripDigiIter

Definition at line 11 of file SiStripCluster.h.

typedef std::pair<SiStripDigiIter,SiStripDigiIter> SiStripCluster::SiStripDigiRange

Definition at line 12 of file SiStripCluster.h.

Constructor & Destructor Documentation

SiStripCluster::SiStripCluster ( ) [inline]

Construct from a range of digis that form a cluster and from a DetID. The range is assumed to be non-empty.

Definition at line 18 of file SiStripCluster.h.

: detId_(0) {}

SiStripCluster::SiStripCluster	(	uint32_t	detid,
		const SiStripDigiRange &	range
	)

check if digis consecutive

Definition at line 4 of file SiStripCluster.cc.

References amplitudes_, i, and j.

                                                                             :
  detId_(detid), firstStrip_(range.first->strip())
{

  amplitudes_.reserve( range.second - range.first);
  
  uint16_t lastStrip=0;
  bool firstInloop = true;
  for (SiStripDigiIter i=range.first; i!=range.second; i++) {
    
    if (!firstInloop && i->strip() != lastStrip + 1) {
      for (int j=0; j < i->strip()-(lastStrip+1); j++) {
        amplitudes_.push_back( 0);
      }
    }
    lastStrip = i->strip();
    firstInloop = false;
    
    amplitudes_.push_back(i->adc()); 
  }
}

SiStripCluster::SiStripCluster	(	const uint32_t &	detid,
		const uint16_t &	firstStrip,
		std::vector< uint16_t >::const_iterator	begin,
		std::vector< uint16_t >::const_iterator	end
	)

Definition at line 27 of file SiStripCluster.cc.

                                                                      :

  detId_(detid),
  firstStrip_(firstStrip),
  amplitudes_(begin,end)
{}

Member Function Documentation

const std::vector<uint8_t>& SiStripCluster::amplitudes ( ) const [inline]

The amplitudes of the strips forming the cluster. The amplitudes are on consecutive strips; if a strip is missing the amplitude is set to zero. A strip may be missing in the middle of a cluster because of a clusterizer that accepts holes. A strip may also be missing anywhere in the cluster, including the edge, to record a dead/noisy channel.

Definition at line 45 of file SiStripCluster.h.

References amplitudes_.

Referenced by TrackerHitAssociator::associateSimpleRecHitCluster(), SimpleCosmicBONSeeder::checkCharge(), DeDxDiscriminatorProducer::ClusterSaturatingStrip(), SiStripGainFromData::ComputeChargeOverPath(), DeDxEstimatorProducer::getCharge(), DeDxDiscriminatorProducer::GetProbability(), SiStripRecHitConverterAlgorithm::isMasked(), TkStripMeasurementDet::isMasked(), CalibrationTrackSelector::isOkCharge(), DeDxDiscriminatorLearner::Learn(), StripCPEgeometric::localParameters(), StripCPEfromTrackAngle::localParameters(), ShallowGainCalibration::produce(), ShallowClustersProducer::produce(), jptJetAnalysis::StripSignalOverNoiseCalculator::signalOverNoise(), SiStripClusterInfo::stripCharges(), and SiStripClusterInfo::width().

{return amplitudes_;}

float SiStripCluster::barycenter ( ) const

The barycenter of the cluster, not corrected for Lorentz shift; should not be used as position estimate for tracking.

Definition at line 37 of file SiStripCluster.cc.

References amplitudes_, firstStrip_, and i.

Referenced by SiStripClusterInfo::baryStrip(), SiStripFineDelayHit::closestCluster(), SiStripRecHitConverterAlgorithm::isMasked(), TkStripMeasurementDet::isMasked(), StripCPE::localParameters(), StripCPEfromTrackAngle::localParameters(), StripClusterAboveU::operator()(), and ShallowClustersProducer::produce().

                                      {
  int sumx = 0;
  int suma = 0;
  size_t asize = amplitudes_.size();
  for (size_t i=0;i<asize;++i) {
    sumx += (firstStrip_+i)*(amplitudes_[i]);
    suma += amplitudes_[i];
  }
  
  // strip centers are offcet by half pitch w.r.t. strip numbers,
  // so one has to add 0.5 to get the correct barycenter position
  return sumx / static_cast<float>(suma) + 0.5;
}