#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

int	charge () const

uint16_t	firstStrip () const

float	getSplitClusterError () const

bool	isMerged () const

void	setMerged (bool mergedState)

void	setSplitClusterError (float errx)

	SiStripCluster ()

	SiStripCluster (const SiStripDigiRange &range)

template<typename Iter >
	SiStripCluster (const uint16_t &firstStrip, Iter begin, Iter end)

template<typename Iter >
	SiStripCluster (const uint16_t &firstStrip, Iter begin, Iter end, bool merged)

Static Public Attributes
static const uint16_t	mergedValueMask = 0x8000

static const uint16_t	stripIndexMask = 0x7FFF

Private Attributes
std::vector< uint8_t >	amplitudes_

float	error_x = -99999.9

uint16_t	firstStrip_ = 0

Detailed Description

Definition at line 9 of file SiStripCluster.h.

Member Typedef Documentation

◆ SiStripDigiIter

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

Definition at line 11 of file SiStripCluster.h.

◆ SiStripDigiRange

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

Definition at line 12 of file SiStripCluster.h.

Constructor & Destructor Documentation

◆ SiStripCluster() [1/4]

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 21 of file SiStripCluster.h.

21 {}

◆ SiStripCluster() [2/4]

SiStripCluster::SiStripCluster ( const SiStripDigiRange & range )

explicit

check if digis consecutive

Definition at line 4 of file SiStripCluster.cc.

                                                             : firstStrip_(range.first->strip()), error_x(-99999.9) {
   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());
   }
 }

References amplitudes_, mps_fire::i, dqmiolumiharvest::j, and FastTimerService_cff::range.

◆ SiStripCluster() [3/4]

template<typename Iter >

SiStripCluster::SiStripCluster	(	const uint16_t &	firstStrip,
		Iter	begin,
		Iter	end
	)

inline

Definition at line 26 of file SiStripCluster.h.

26 : amplitudes_(begin, end), firstStrip_(firstStrip) {}

◆ SiStripCluster() [4/4]

template<typename Iter >

SiStripCluster::SiStripCluster	(	const uint16_t &	firstStrip,
		Iter	begin,
		Iter	end,
		bool	merged
	)

inline

Definition at line 29 of file SiStripCluster.h.

       : amplitudes_(begin, end), firstStrip_(firstStrip) {
     if (merged)
       firstStrip_ |= mergedValueMask;  // if this is a candidate merged cluster
   }

References firstStrip_, and mergedValueMask.

Member Function Documentation

◆ amplitudes()

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.

You can find the special meanings of values { 0, 254, 255} in section 3.4.1 of http://www.te.rl.ac.uk/esdg/cms-fed/firmware/Documents/FE_FPGA_Technical_Description.pdf

Definition at line 51 of file SiStripCluster.h.

51 { return amplitudes_; }

References amplitudes_.

Referenced by TrackerHitAssociator::associateSimpleRecHitCluster(), charge(), SimpleCosmicBONSeeder::checkCharge(), ClusterVariables::clusterSize(), SiStripGainFromData::ComputeChargeOverPath(), DeDxTools::getCharge(), ClusterShapeHitFilter::getSizes(), SiStripRecHitConverterAlgorithm::isMasked(), StMeasurementConditionSet::isMasked(), CalibrationTrackSelector::isOkCharge(), StripCPEfromTrackAngle::localParameters(), PrintRecoObjects::print(), ClusterTPAssociationProducer::produce(), ShallowGainCalibration::produce(), DeDxTools::shapeSelection(), ClusterVariables::stripCharges(), SiStripClusterInfo::stripCharges(), StripSubClusterShapeFilterBase::testLastHit(), reco::InverseCrosstalkMatrix::unfold(), and SiStripClusterInfo::width().

◆ barycenter()

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 23 of file SiStripCluster.cc.

                                        {
   int sumx = 0;
   int suma = 0;
   auto asize = amplitudes_.size();
   for (auto i = 0U; i < asize; ++i) {
     sumx += 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.
   // Need to mask off the high bit of firstStrip_, which contains the merged status.
   return float((firstStrip_ & stripIndexMask)) + float(sumx) / float(suma) + 0.5f;
 }