#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

float	getSplitClusterError () const

void	setSplitClusterError (float errx)

	SiStripCluster ()

	SiStripCluster (const SiStripDigiRange &range)

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

Private Attributes
std::vector< uint8_t >	amplitudes_

float	error_x

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.

18 : error_x( -99999.9 ) {}

SiStripCluster::error_x

float error_x

Definition: SiStripCluster.h:73

SiStripCluster::SiStripCluster ( const SiStripDigiRange & range )

explicit

check if digis consecutive

Definition at line 4 of file SiStripCluster.cc.

References amplitudes_, i, and j.

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

                                :
      amplitudes_(begin,end), firstStrip_(firstStrip), 
   // ggiurgiu@fnal.gov, 01/05/12
   // Initialize the split cluster errors to un-physical values.
   // The CPE will check these errors and if they are not un-physical,
   // it will recognize the clusters as split and assign these (increased)
   // errors to the corresponding rechit.
   error_x(-99999.9){}

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

References amplitudes_.

Referenced by SiStripGainFromCalibTree::algoAnalyze(), TrackerHitAssociator::associateSimpleRecHitCluster(), SimpleCosmicBONSeeder::checkCharge(), ClusterShapeHitFilter::checkClusterCharge(), DeDxDiscriminatorProducer::ClusterSaturatingStrip(), ClusterVariables::clusterSize(), SiStripGainFromData::ComputeChargeOverPath(), DeDxEstimatorProducer::getCharge(), DeDxDiscriminatorProducer::GetProbability(), ClusterShapeHitFilter::getSizes(), SiStripRecHitConverterAlgorithm::isMasked(), StMeasurementConditionSet::isMasked(), CalibrationTrackSelector::isOkCharge(), DeDxDiscriminatorLearner::Learn(), StripCPEfromTrackAngle::localParameters(), StripCPEgeometric::localParameters(), StripCPEfromTemplate::localParameters(), PrintRecoObjects::print(), ClusterTPAssociationProducer::produce(), ShallowGainCalibration::produce(), TrackClusterSplitter::splitCluster(), ClusterVariables::stripCharges(), SiStripClusterInfo::stripCharges(), and SiStripClusterInfo::width().

48 {return amplitudes_;}

SiStripCluster::amplitudes_

std::vector< uint8_t > amplitudes_

Definition: SiStripCluster.h:61

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_, f, firstStrip_, and i.

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

                                       {
   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.5f;
 }