#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

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

Definition at line 12 of file SiStripCluster.h.

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

Definition at line 13 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 23 of file SiStripCluster.h.

23 {}

SiStripCluster::SiStripCluster ( const SiStripDigiRange & range )

explicit

check if digis consecutive

Definition at line 4 of file SiStripCluster.cc.

References amplitudes_, and mps_fire::i.

                                                             :
   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()); 
   }
 }

template<typename Iter >

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

inline

Definition at line 28 of file SiStripCluster.h.

29 :

30 amplitudes_(begin,end), firstStrip_(firstStrip) {}

SiStripCluster::firstStrip

uint16_t firstStrip() const

Definition: SiStripCluster.h:41

end

#define end

Definition: vmac.h:39

begin

#define begin

Definition: vmac.h:32

SiStripCluster::amplitudes_

std::vector< uint8_t > amplitudes_

Definition: SiStripCluster.h:78

SiStripCluster::firstStrip_

uint16_t firstStrip_

Definition: SiStripCluster.h:80

template<typename Iter >

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

inline

Definition at line 33 of file SiStripCluster.h.

References firstStrip_, and mergedValueMask.

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

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.

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

References amplitudes_, and barycenter().

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

54 {return amplitudes_;}

SiStripCluster::amplitudes_

std::vector< uint8_t > amplitudes_

Definition: SiStripCluster.h:78

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

References amplitudes_, firstStrip_, objects.autophobj::float, mps_fire::i, stripIndexMask, and mitigatedMETSequence_cff::U.

Referenced by amplitudes(), SiStripClusterInfo::baryStrip(), SiStripFineDelayHit::closestCluster(), SiStripRecHitConverterAlgorithm::isMasked(), StMeasurementConditionSet::isMasked(), StripCPEfromTemplate::localParameters(), StripCPE::localParameters(), StripCPEfromTrackAngle::localParameters(), PrintRecoObjects::print(), and TrackClusterSplitter::splitCluster().

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