#include <TkStripMeasurementDet.h>

Inheritance diagram for TkStripMeasurementDet:

Public Types
typedef StMeasurementConditionSet::BadStripBlock	BadStripBlock

typedef StMeasurementConditionSet::BadStripCuts	BadStripCuts

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

typedef edmNew::DetSet< SiStripCluster >	detset

typedef StripClusterParameterEstimator::LocalValues	LocalValues

typedef detset::const_iterator	new_const_iterator

typedef SiStripRecHit2D::ClusterRef	SiStripClusterRef

typedef StripClusterParameterEstimator::VLocalValues	VLocalValues

Public Types inherited from MeasurementDet
typedef TrackingRecHit::ConstRecHitContainer	RecHitContainer

using	SimpleHitContainer = std::vector< BaseTrackerRecHit * >

typedef tracking::TempMeasurements	TempMeasurements

Public Member Functions
bool	accept (SiStripClusterRef const &r, const std::vector< bool > &skipClusters) const

bool	accept (unsigned int key, const std::vector< bool > &skipClusters) const

void	advance (TkStripRecHitIter &hi) const

std::vector< BadStripBlock > const &	badStripBlocks () const

template<class ClusterRefT >
TrackingRecHit::RecHitPointer	buildRecHit (const ClusterRefT &cluster, const TrajectoryStateOnSurface &ltp) const

template<class ClusterRefT >
void	buildRecHits (const ClusterRefT &cluster, const TrajectoryStateOnSurface &ltp, const RecHitContainer &_res) const

const detset &	detSet (const StMeasurementDetSet &theDets) const

bool	empty (const MeasurementTrackerEvent &data) const

template<class ClusterRefT >
bool	filteredRecHits (const ClusterRefT &cluster, StripCPE::AlgoParam const &cpepar, const TrajectoryStateOnSurface &ltp, const MeasurementEstimator &est, const std::vector< bool > &skipClusters, RecHitContainer &result, std::vector< float > &diffs) const

template<class ClusterRefT >
bool	filteredRecHits (const ClusterRefT &cluster, StripCPE::AlgoParam const &cpepar, const TrajectoryStateOnSurface &ltp, const MeasurementEstimator &est, const std::vector< bool > &skipClusters, std::vector< SiStripRecHit2D > &result) const

std::vector< BadStripBlock > &	getBadStripBlocks ()

bool	hasAllGoodChannels () const
	does this module have at least one bad strip, APV or channel? More...

bool	hasBadComponents (const TrajectoryStateOnSurface &tsos, const MeasurementTrackerEvent &data) const override

SiStripRecHit2D	hit (TkStripRecHitIter const &hi) const

std::tuple< TkStripRecHitIter, TkStripRecHitIter >	hitRange (const TrajectoryStateOnSurface &, const MeasurementTrackerEvent &data) const

int	index () const

bool	isActive (const MeasurementTrackerEvent &data) const override
	Is this module active in reconstruction? It must be both 'setActiveThisEvent' and 'setActive'. More...

bool	isEmpty (const StMeasurementDetSet &theDets) const

bool	maskBad128StripBlocks () const

bool	measurements (const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &est, const MeasurementTrackerEvent &data, TempMeasurements &result) const override

unsigned int	rawId () const

RecHitContainer	recHits (const TrajectoryStateOnSurface &, const MeasurementTrackerEvent &data) const override

bool	recHits (const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &est, const MeasurementTrackerEvent &data, RecHitContainer &result, std::vector< float > &diffs) const override

bool	recHits (SimpleHitContainer &result, const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &, const MeasurementTrackerEvent &data) const override

void	set128StripStatus (bool good, int idx=-1)
	Sets the status of a block of 128 strips (or all blocks if idx=-1) More...

void	setActiveThisEvent (StMeasurementDetSet &theDets, bool active) const
	Turn on/off the module for reconstruction for one events. This per-event flag is cleared by any call to 'update' or 'setEmpty' More...

void	setActiveThisPeriod (StMeasurementDetSet &theDets, bool active)
	Turn on/off the module for reconstruction, for the full run or lumi (using info from DB, usually). More...

void	setEmpty (StMeasurementDetSet &theDets) const

void	setIndex (int i)

bool	simpleRecHits (const TrajectoryStateOnSurface &ts, const MeasurementEstimator &est, const MeasurementTrackerEvent &data, std::vector< SiStripRecHit2D > &result) const

void	simpleRecHits (const TrajectoryStateOnSurface &ts, const MeasurementTrackerEvent &data, std::vector< SiStripRecHit2D > &result) const

const StripGeomDetUnit &	specificGeomDet () const

unsigned char	subId () const

bool	testStrips (float utraj, float uerr) const
	return true if there are 'enough' good strips in the utraj +/- 3 uerr range. More...

const detset &	theSet (const StMeasurementDetSet &theDets) const

	TkStripMeasurementDet (const GeomDet *gdet, StMeasurementConditionSet &conditionSet)

	~TkStripMeasurementDet () override

Public Member Functions inherited from MeasurementDet
const GeomDet &	fastGeomDet () const

std::vector< TrajectoryMeasurement >	fastMeasurements (const TrajectoryStateOnSurface &stateOnThisDet, const TrajectoryStateOnSurface &, const Propagator &, const MeasurementEstimator &est, const MeasurementTrackerEvent &data) const

virtual const GeomDet &	geomDet () const

	MeasurementDet (const GeomDet *gdet)

const Surface::PositionType &	position () const

const Surface &	surface () const

virtual	~MeasurementDet ()=default

Private Types
using	AClusters = StripClusterParameterEstimator::AClusters

using	ALocalValues = StripClusterParameterEstimator::ALocalValues

Private Member Functions
BadStripCuts const &	badStripCuts () const

void	buildSimpleRecHits (AClusters const &clusters, const MeasurementTrackerEvent &data, const detset &detSet, const TrajectoryStateOnSurface &ltp, std::vector< SiStripRecHit2D > &res) const

StMeasurementConditionSet &	conditionSet ()

const StMeasurementConditionSet &	conditionSet () const

const StripCPE *	cpe () const

bool	hasAny128StripBad () const

bool	isMasked (const SiStripCluster &cluster) const

int	totalStrips () const

Private Attributes
int	index_

StMeasurementConditionSet *	theDetConditions

Additional Inherited Members
Protected Attributes inherited from MeasurementDet
TrackingRecHit::ConstRecHitPointer	theInactiveHit

TrackingRecHit::ConstRecHitPointer	theMissingHit

Detailed Description

Definition at line 74 of file TkStripMeasurementDet.h.

Member Typedef Documentation

◆ AClusters

using TkStripMeasurementDet::AClusters = StripClusterParameterEstimator::AClusters

private

Definition at line 245 of file TkStripMeasurementDet.h.

◆ ALocalValues

using TkStripMeasurementDet::ALocalValues = StripClusterParameterEstimator::ALocalValues

private

Definition at line 246 of file TkStripMeasurementDet.h.

◆ BadStripBlock

typedef StMeasurementConditionSet::BadStripBlock TkStripMeasurementDet::BadStripBlock

Definition at line 237 of file TkStripMeasurementDet.h.

◆ BadStripCuts

typedef StMeasurementConditionSet::BadStripCuts TkStripMeasurementDet::BadStripCuts

Definition at line 232 of file TkStripMeasurementDet.h.

◆ const_iterator

typedef std::vector<SiStripCluster>::const_iterator TkStripMeasurementDet::const_iterator

Definition at line 84 of file TkStripMeasurementDet.h.

◆ detset

typedef edmNew::DetSet<SiStripCluster> TkStripMeasurementDet::detset

Definition at line 81 of file TkStripMeasurementDet.h.

◆ LocalValues

typedef StripClusterParameterEstimator::LocalValues TkStripMeasurementDet::LocalValues

Definition at line 76 of file TkStripMeasurementDet.h.

◆ new_const_iterator

typedef detset::const_iterator TkStripMeasurementDet::new_const_iterator

Definition at line 82 of file TkStripMeasurementDet.h.

◆ SiStripClusterRef

typedef SiStripRecHit2D::ClusterRef TkStripMeasurementDet::SiStripClusterRef

Definition at line 79 of file TkStripMeasurementDet.h.

◆ VLocalValues

typedef StripClusterParameterEstimator::VLocalValues TkStripMeasurementDet::VLocalValues

Definition at line 77 of file TkStripMeasurementDet.h.

Constructor & Destructor Documentation

◆ ~TkStripMeasurementDet()

TkStripMeasurementDet::~TkStripMeasurementDet ( )

inlineoverride

Definition at line 86 of file TkStripMeasurementDet.h.

86 {}

◆ TkStripMeasurementDet()

TkStripMeasurementDet::TkStripMeasurementDet	(	const GeomDet *	gdet,
		StMeasurementConditionSet &	conditionSet
	)

Definition at line 16 of file TkStripMeasurementDet.cc.

     : MeasurementDet(gdet), index_(-1), theDetConditions(&conditions) {
   if (dynamic_cast<const StripGeomDetUnit*>(gdet) == nullptr) {
     throw MeasurementDetException("TkStripMeasurementDet constructed with a GeomDet which is not a StripGeomDetUnit");
   }
 }

Member Function Documentation

◆ accept() [1/2]

bool TkStripMeasurementDet::accept	(	SiStripClusterRef const &	r,
		const std::vector< bool > &	skipClusters
	)		const

inline

Definition at line 277 of file TkStripMeasurementDet.h.

                                                                                             {
     return accept(r.key(), skipClusters);
   }

References accept(), alignCSCRings::r, and HLTSiStripMonitoring_cff::skipClusters.

Referenced by advance(), empty(), esMonitoring.FDJsonServer::handle_accept(), recHits(), and simpleRecHits().

◆ accept() [2/2]

bool TkStripMeasurementDet::accept	(	unsigned int	key,
		const std::vector< bool > &	skipClusters
	)		const

inline

Definition at line 281 of file TkStripMeasurementDet.h.

                                                                                   {
     if (skipClusters.empty())
       return true;
     if (key >= skipClusters.size()) {
       LogDebug("TkStripMeasurementDet")
           << key << " is larger than: " << skipClusters.size()
           << "\n This must be a new cluster, and therefore should not be skiped most likely.";
       return true;
     }
     return (not(skipClusters[key]));
   }

References crabWrapper::key, LogDebug, and HLTSiStripMonitoring_cff::skipClusters.

Referenced by esMonitoring.FDJsonServer::handle_accept().

◆ advance()

void TkStripMeasurementDet::advance ( TkStripRecHitIter & hi ) const

Definition at line 250 of file TkStripMeasurementDet.cc.

                                                                {
   while (!hi.empty()) {
     auto ci = hi.clusterI;
     auto const& data = *hi.data;
     if (isMasked(*ci))
       continue;
     SiStripClusterRef cluster = edmNew::makeRefTo(data.stripData().handle(), ci);
     if (accept(cluster, data.stripClustersToSkip()))
       return;
     ++hi.clusterI;
   }
 }

References accept(), data, isMasked(), and edmNew::makeRefTo().

Referenced by TkStripRecHitIter::advance().

◆ badStripBlocks()

std::vector<BadStripBlock> const& TkStripMeasurementDet::badStripBlocks ( ) const

inline

Definition at line 240 of file TkStripMeasurementDet.h.

240 { return conditionSet().badStripBlocks(index()); }

Referenced by testStrips().

◆ badStripCuts()

BadStripCuts const& TkStripMeasurementDet::badStripCuts ( ) const

inlineprivate

Definition at line 257 of file TkStripMeasurementDet.h.

257 { return conditionSet().badStripCuts(index()); }

Referenced by testStrips().

◆ buildRecHit()

template<class ClusterRefT >

TrackingRecHit::RecHitPointer TkStripMeasurementDet::buildRecHit	(	const ClusterRefT &	cluster,
		const TrajectoryStateOnSurface &	ltp
	)		const

inline

Definition at line 150 of file TkStripMeasurementDet.h.

                                                                                                                  {
     const GeomDetUnit& gdu(specificGeomDet());
     LocalValues lv = cpe()->localParameters(*cluster, gdu, ltp);
     return std::make_shared<SiStripRecHit2D>(lv.first, lv.second, fastGeomDet(), cluster);
   }

Referenced by recHits().

◆ buildRecHits()

template<class ClusterRefT >

void TkStripMeasurementDet::buildRecHits	(	const ClusterRefT &	cluster,
		const TrajectoryStateOnSurface &	ltp,
		const RecHitContainer &	_res
	)		const

inline

Definition at line 157 of file TkStripMeasurementDet.h.

                                                        {
     RecHitContainer res = _res;
     const GeomDetUnit& gdu(specificGeomDet());
     VLocalValues vlv = cpe()->localParametersV(*cluster, gdu, ltp);
     for (VLocalValues::const_iterator it = vlv.begin(); it != vlv.end(); ++it)
       res.push_back(std::make_shared<SiStripRecHit2D>(it->first, it->second, fastGeomDet(), cluster));
   }

◆ buildSimpleRecHits()

void TkStripMeasurementDet::buildSimpleRecHits	(	AClusters const &	clusters,
		const MeasurementTrackerEvent &	data,
		const detset &	detSet,
		const TrajectoryStateOnSurface &	ltp,
		std::vector< SiStripRecHit2D > &	res
	)		const

inlineprivate

Definition at line 263 of file TkStripMeasurementDet.h.

                                                                  {
     const GeomDetUnit& gdu(specificGeomDet());
     declareDynArray(LocalValues, clusters.size(), alv);
     cpe()->localParameters(clusters, alv, gdu, ltp.localParameters());
     res.reserve(alv.size());
     for (unsigned int i = 0; i < clusters.size(); ++i)
       res.emplace_back(alv[i].first, alv[i].second, gdu, detSet.makeRefTo(data.stripData().handle(), clusters[i]));
   }

References bsc_activity_cfg::clusters, data, declareDynArray, mps_fire::i, TrajectoryStateOnSurface::localParameters(), and edmNew::DetSet< T >::makeRefTo().

Referenced by simpleRecHits().

◆ conditionSet() [1/2]

StMeasurementConditionSet& TkStripMeasurementDet::conditionSet ( )

inlineprivate

Definition at line 250 of file TkStripMeasurementDet.h.

250 { return *theDetConditions; }

◆ conditionSet() [2/2]

const StMeasurementConditionSet& TkStripMeasurementDet::conditionSet ( ) const

inlineprivate

Definition at line 251 of file TkStripMeasurementDet.h.

251 { return *theDetConditions; }

◆ cpe()

const StripCPE* TkStripMeasurementDet::cpe ( ) const

inlineprivate

Definition at line 253 of file TkStripMeasurementDet.h.

253 { return static_cast<const StripCPE*>(conditionSet().stripCPE()); }

Referenced by hit(), recHits(), and simpleRecHits().

◆ detSet()

const detset& TkStripMeasurementDet::detSet ( const StMeasurementDetSet & theDets ) const

inline

Definition at line 102 of file TkStripMeasurementDet.h.

102 { return theDets.detSet(index()); }

References StMeasurementDetSet::detSet().

Referenced by empty(), hitRange(), recHits(), and simpleRecHits().

◆ empty()

bool TkStripMeasurementDet::empty ( const MeasurementTrackerEvent & data ) const

Definition at line 24 of file TkStripMeasurementDet.cc.

                                                                            {
   if UNLIKELY ((!isActive(data)) || isEmpty(data.stripData()))
     return true;
  
   const detset& detSet = data.stripData().detSet(index());
   for (auto ci = detSet.begin(); ci != detSet.end(); ++ci) {
     if (isMasked(*ci))
       continue;
     SiStripClusterRef cluster = detSet.makeRefTo(data.stripData().handle(), ci);
     if (accept(cluster, data.stripClustersToSkip()))
       return false;
   }
   return true;
 }

References accept(), edmNew::DetSet< T >::begin(), data, detSet(), edmNew::DetSet< T >::end(), index(), isActive(), isEmpty(), isMasked(), edmNew::DetSet< T >::makeRefTo(), and UNLIKELY.

◆ filteredRecHits() [1/2]

template<class ClusterRefT >

bool TkStripMeasurementDet::filteredRecHits	(	const ClusterRefT &	cluster,
		StripCPE::AlgoParam const &	cpepar,
		const TrajectoryStateOnSurface &	ltp,
		const MeasurementEstimator &	est,
		const std::vector< bool > &	skipClusters,
		RecHitContainer &	result,
		std::vector< float > &	diffs
	)		const

inline

Definition at line 168 of file TkStripMeasurementDet.h.

                                                       {
     if (isMasked(*cluster))
       return true;
     if (!accept(cluster, skipClusters))
       return true;
     if (!est.preFilter(ltp, ClusterFilterPayload(rawId(), &*cluster)))
       return true;  // avoids shadow; consistent with previous statement...
     auto const& vl = cpe()->localParameters(*cluster, cpepar);
     SiStripRecHit2D recHit(vl.first,
                            vl.second,
                            fastGeomDet(),
                            cluster);  // FIXME add cluster count in OmniRef (and move again to multiple sub-clusters..)
     std::pair<bool, double> diffEst = est.estimate(ltp, recHit);
     LogDebug("TkStripMeasurementDet") << " chi2=" << diffEst.second;
     if (diffEst.first) {
       result.push_back(std::make_shared<SiStripRecHit2D>(recHit));
       diffs.push_back(diffEst.second);
     }
     return diffEst.first;
   }

References accept(), MeasurementEstimator::estimate(), LogDebug, MeasurementEstimator::preFilter(), rpcPointValidation_cfi::recHit, mps_fire::result, and HLTSiStripMonitoring_cff::skipClusters.

Referenced by recHits(), and simpleRecHits().

◆ filteredRecHits() [2/2]

template<class ClusterRefT >

bool TkStripMeasurementDet::filteredRecHits	(	const ClusterRefT &	cluster,
		StripCPE::AlgoParam const &	cpepar,
		const TrajectoryStateOnSurface &	ltp,
		const MeasurementEstimator &	est,
		const std::vector< bool > &	skipClusters,
		std::vector< SiStripRecHit2D > &	result
	)		const

inline

Definition at line 196 of file TkStripMeasurementDet.h.

                                                                  {
     if (isMasked(*cluster))
       return true;
     if (!accept(cluster, skipClusters))
       return true;
     if (!est.preFilter(ltp, ClusterFilterPayload(rawId(), &*cluster)))
       return true;  // avoids shadow; consistent with previous statement...
     auto const& vl = cpe()->localParameters(*cluster, cpepar);
     result.emplace_back(vl.first, vl.second, fastGeomDet(), cluster);  // FIXME add cluster count in OmniRef
     std::pair<bool, double> diffEst = est.estimate(ltp, result.back());
     LogDebug("TkStripMeasurementDet") << " chi2=" << diffEst.second;
     if (!diffEst.first)
       result.pop_back();
     return diffEst.first;
   }

References accept(), MeasurementEstimator::estimate(), LogDebug, MeasurementEstimator::preFilter(), mps_fire::result, and HLTSiStripMonitoring_cff::skipClusters.

◆ getBadStripBlocks()

std::vector<BadStripBlock>& TkStripMeasurementDet::getBadStripBlocks ( )

inline

Definition at line 239 of file TkStripMeasurementDet.h.

239 { return conditionSet().getBadStripBlocks(index()); }

◆ hasAllGoodChannels()

bool TkStripMeasurementDet::hasAllGoodChannels ( ) const

inline

does this module have at least one bad strip, APV or channel?

Definition at line 227 of file TkStripMeasurementDet.h.

227 { return (!hasAny128StripBad()) && badStripBlocks().empty(); }

Referenced by TkGluedMeasurementDet::measurements().

◆ hasAny128StripBad()

bool TkStripMeasurementDet::hasAny128StripBad ( ) const

inlineprivate

Definition at line 259 of file TkStripMeasurementDet.h.

259 { return conditionSet().hasAny128StripBad(index()); }

◆ hasBadComponents()

bool TkStripMeasurementDet::hasBadComponents	(	const TrajectoryStateOnSurface &	tsos,
		const MeasurementTrackerEvent &	data
	)		const

inlineoverridevirtual

Implements MeasurementDet.

Definition at line 108 of file TkStripMeasurementDet.h.

                                                                                                                   {
     return false;
   }

◆ hit()

SiStripRecHit2D TkStripMeasurementDet::hit ( TkStripRecHitIter const & hi ) const

Definition at line 263 of file TkStripMeasurementDet.cc.

                                                                             {
   const GeomDetUnit& gdu(specificGeomDet());
   auto ci = hi.clusterI;
   auto const& data = *hi.data;
   auto const& ltp = *hi.tsos;
  
   SiStripClusterRef cluster = edmNew::makeRefTo(data.stripData().handle(), ci);
   LocalValues lv = cpe()->localParameters(*cluster, gdu, ltp);
   return SiStripRecHit2D(lv.first, lv.second, gdu, cluster);
 }

References cpe(), data, StripCPE::localParameters(), edmNew::makeRefTo(), and specificGeomDet().

Referenced by TkStripRecHitIter::buildHit(), recHits(), and simpleRecHits().

◆ hitRange()

std::tuple< TkStripRecHitIter, TkStripRecHitIter > TkStripMeasurementDet::hitRange	(	const TrajectoryStateOnSurface &	ts,
		const MeasurementTrackerEvent &	data
	)		const

Definition at line 241 of file TkStripMeasurementDet.cc.

                                                                                    {
   if (isEmpty(data.stripData()) || !isActive(data))
     return std::tuple<TkStripRecHitIter, TkStripRecHitIter>();
   const detset& detSet = data.stripData().detSet(index());
   return std::make_tuple(TkStripRecHitIter(detSet.begin(), detSet.end(), *this, ts, data),
                          TkStripRecHitIter(detSet.end(), detSet.end(), *this, ts, data));
 }

References edmNew::DetSet< T >::begin(), data, detSet(), edmNew::DetSet< T >::end(), index(), isActive(), and isEmpty().

◆ index()

int TkStripMeasurementDet::index ( ) const

inline

Definition at line 96 of file TkStripMeasurementDet.h.

96 { return index_; }

Referenced by empty(), hitRange(), BeautifulSoup.PageElement::insert(), recHits(), and simpleRecHits().

◆ isActive()

bool TkStripMeasurementDet::isActive ( const MeasurementTrackerEvent & data ) const

inlineoverridevirtual

Is this module active in reconstruction? It must be both 'setActiveThisEvent' and 'setActive'.

Implements MeasurementDet.

Definition at line 105 of file TkStripMeasurementDet.h.

105 { return data.stripData().isActive(index()); }

References data.

Referenced by empty(), hitRange(), TkGluedMeasurementDet::measurements(), measurements(), TkGluedMeasurementDet::recHits(), recHits(), and simpleRecHits().

◆ isEmpty()

bool TkStripMeasurementDet::isEmpty ( const StMeasurementDetSet & theDets ) const

inline

Definition at line 94 of file TkStripMeasurementDet.h.

94 { return theDets.empty(index()); }

References StMeasurementDetSet::empty().

Referenced by plotting.Plot::clone(), empty(), hitRange(), measurements(), recHits(), and simpleRecHits().

◆ isMasked()

bool TkStripMeasurementDet::isMasked ( const SiStripCluster & cluster ) const

inlineprivate

Definition at line 261 of file TkStripMeasurementDet.h.

261 { return conditionSet().isMasked(index(), cluster); }

Referenced by advance(), empty(), recHits(), and simpleRecHits().

◆ maskBad128StripBlocks()

bool TkStripMeasurementDet::maskBad128StripBlocks ( ) const

inline

Definition at line 242 of file TkStripMeasurementDet.h.

242 { return conditionSet().maskBad128StripBlocks(); }

◆ measurements()

bool TkStripMeasurementDet::measurements	(	const TrajectoryStateOnSurface &	stateOnThisDet,
		const MeasurementEstimator &	est,
		const MeasurementTrackerEvent &	data,
		TempMeasurements &	result
	)		const

overridevirtual

Implements MeasurementDet.

Definition at line 178 of file TkStripMeasurementDet.cc.

                                                                          {
   if (!isActive(data)) {
     LogDebug("TkStripMeasurementDet") << " found an inactive module " << rawId();
     result.add(theInactiveHit, 0.F);
     return true;
   }
  
   if (!isEmpty(data.stripData())) {
     LogDebug("TkStripMeasurementDet") << " found hit on this module " << rawId();
     RecHitContainer rechits;
     std::vector<float> diffs;
     if (recHits(stateOnThisDet, est, data, result.hits, result.distances))
       return true;
   }
  
   // create a TrajectoryMeasurement with an invalid RecHit and zero estimate
  
   if (!stateOnThisDet.hasError()) {
     result.add(theMissingHit, 0.F);
     return false;
   }
  
   float utraj = specificGeomDet().specificTopology().measurementPosition(stateOnThisDet.localPosition()).x();
   float uerr = sqrt(specificGeomDet()
                         .specificTopology()
                         .measurementError(stateOnThisDet.localPosition(), stateOnThisDet.localError().positionError())
                         .uu());
   if (testStrips(utraj, uerr)) {
     //LogDebug("TkStripMeasurementDet") << " DetID " << id_ << " empty after search, but active ";
     result.add(theMissingHit, 0.F);
     return false;
   }
  
   //LogDebug("TkStripMeasurementDet") << " DetID " << id_ << " empty after search, and inactive ";
   result.add(theInactiveHit, 0.F);
   return true;
 }

References data, F(), TrajectoryStateOnSurface::hasError(), isActive(), isEmpty(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), LogDebug, Topology::measurementPosition(), LocalTrajectoryError::positionError(), rawId(), HI_PhotonSkim_cff::rechits, recHits(), mps_fire::result, specificGeomDet(), StripGeomDetUnit::specificTopology(), mathSSE::sqrt(), testStrips(), MeasurementDet::theInactiveHit, MeasurementDet::theMissingHit, and x.

◆ rawId()

unsigned int TkStripMeasurementDet::rawId ( ) const

inline

Definition at line 98 of file TkStripMeasurementDet.h.

98 { return conditionSet().id(index()); }

Referenced by measurements(), recHits(), and simpleRecHits().

◆ recHits() [1/3]

TkStripMeasurementDet::RecHitContainer TkStripMeasurementDet::recHits	(	const TrajectoryStateOnSurface &	ts,
		const MeasurementTrackerEvent &	data
	)		const

overridevirtual

Implements MeasurementDet.

Definition at line 39 of file TkStripMeasurementDet.cc.

                                                                                                                  {
   RecHitContainer result;
   if UNLIKELY ((!isActive(data)) || isEmpty(data.stripData()))
     return result;
   const detset& detSet = data.stripData().detSet(index());
   result.reserve(detSet.size());
   for (new_const_iterator ci = detSet.begin(); ci != detSet.end(); ++ci) {
     if (isMasked(*ci))
       continue;
     // for ( ClusterIterator ci=theClusterRange.first; ci != theClusterRange.second; ci++) {
     SiStripClusterRef cluster = detSet.makeRefTo(data.stripData().handle(), ci);
     if (accept(cluster, data.stripClustersToSkip()))
       result.push_back(buildRecHit(cluster, ts));
     else
       LogDebug("TkStripMeasurementDet") << "skipping this str from last iteration on" << rawId()
                                         << " key: " << cluster.key();
   }
   return result;
 }

References accept(), edmNew::DetSet< T >::begin(), buildRecHit(), data, detSet(), edmNew::DetSet< T >::end(), index(), isActive(), isEmpty(), isMasked(), edm::Ref< C, T, F >::key(), LogDebug, edmNew::DetSet< T >::makeRefTo(), rawId(), mps_fire::result, edmNew::DetSet< T >::size(), and UNLIKELY.

Referenced by TkGluedMeasurementDet::collectRecHits(), and measurements().

◆ recHits() [2/3]

bool TkStripMeasurementDet::recHits	(	const TrajectoryStateOnSurface &	stateOnThisDet,
		const MeasurementEstimator &	est,
		const MeasurementTrackerEvent &	data,
		RecHitContainer &	result,
		std::vector< float > &	diffs
	)		const

overridevirtual

Reimplemented from MeasurementDet.

Definition at line 139 of file TkStripMeasurementDet.cc.

                                                                    {
   if UNLIKELY ((!isActive(data)) || isEmpty(data.stripData()))
     return false;
  
   auto oldSize = result.size();
   auto const& cpepar = cpe()->getAlgoParam(specificGeomDet(), stateOnThisDet.localParameters());
  
   int utraj = specificGeomDet().specificTopology().measurementPosition(stateOnThisDet.localPosition()).x();
  
   const detset& detSet = data.stripData().detSet(index());
   auto rightCluster = std::find_if(
       detSet.begin(), detSet.end(), [utraj](const SiStripCluster& hit) { return hit.firstStrip() > utraj; });
  
   if (rightCluster != detSet.begin()) {
     // there are hits on the left of the utraj
     auto leftCluster = rightCluster;
     while (--leftCluster >= detSet.begin()) {
       SiStripClusterRef clusterref = detSet.makeRefTo(data.stripData().handle(), leftCluster);
       bool isCompatible =
           filteredRecHits(clusterref, cpepar, stateOnThisDet, est, data.stripClustersToSkip(), result, diffs);
       if (!isCompatible)
         break;  // exit loop on first incompatible hit
     }
   }
   for (; rightCluster != detSet.end(); rightCluster++) {
     SiStripClusterRef clusterref = detSet.makeRefTo(data.stripData().handle(), rightCluster);
     bool isCompatible =
         filteredRecHits(clusterref, cpepar, stateOnThisDet, est, data.stripClustersToSkip(), result, diffs);
     if (!isCompatible)
       break;  // exit loop on first incompatible hit
   }
  
   return result.size() > oldSize;
 }

References edmNew::DetSet< T >::begin(), cpe(), data, detSet(), edmNew::DetSet< T >::end(), filteredRecHits(), StripCPE::getAlgoParam(), hit(), index(), isActive(), isEmpty(), TrajectoryStateOnSurface::localParameters(), TrajectoryStateOnSurface::localPosition(), edmNew::DetSet< T >::makeRefTo(), Topology::measurementPosition(), mps_fire::result, specificGeomDet(), StripGeomDetUnit::specificTopology(), UNLIKELY, and x.

◆ recHits() [3/3]

bool TkStripMeasurementDet::recHits	(	SimpleHitContainer &	result,
		const TrajectoryStateOnSurface &	stateOnThisDet,
		const MeasurementEstimator &	est,
		const MeasurementTrackerEvent &	data
	)		const

overridevirtual

Reimplemented from MeasurementDet.

Definition at line 61 of file TkStripMeasurementDet.cc.

                                                                                {
   if UNLIKELY ((!isActive(data)) || isEmpty(data.stripData()))
     return false;
   auto oldSize = result.size();
  
   int utraj = specificGeomDet().specificTopology().measurementPosition(stateOnThisDet.localPosition()).x();
   const detset& detSet = data.stripData().detSet(index());
   auto const& cpepar = cpe()->getAlgoParam(specificGeomDet(), stateOnThisDet.localParameters());
  
   auto rightCluster = std::find_if(
       detSet.begin(), detSet.end(), [utraj](const SiStripCluster& hit) { return hit.firstStrip() > utraj; });
  
   std::vector<SiStripRecHit2D> tmp;
   if (rightCluster != detSet.begin()) {
     // there are hits on the left of the utraj
     auto leftCluster = rightCluster;
     while (--leftCluster >= detSet.begin()) {
       SiStripClusterRef clusterref = detSet.makeRefTo(data.stripData().handle(), leftCluster);
       bool isCompatible = filteredRecHits(clusterref, cpepar, stateOnThisDet, est, data.stripClustersToSkip(), tmp);
       if (!isCompatible)
         break;  // exit loop on first incompatible hit
       for (auto&& h : tmp)
         result.push_back(new SiStripRecHit2D(h));
       tmp.clear();
     }
   }
   for (; rightCluster != detSet.end(); rightCluster++) {
     SiStripClusterRef clusterref = detSet.makeRefTo(data.stripData().handle(), rightCluster);
     bool isCompatible = filteredRecHits(clusterref, cpepar, stateOnThisDet, est, data.stripClustersToSkip(), tmp);
     if (!isCompatible)
       break;  // exit loop on first incompatible hit
     for (auto&& h : tmp)
       result.push_back(new SiStripRecHit2D(h));
     tmp.clear();
   }
  
   return result.size() > oldSize;
 }

References edmNew::DetSet< T >::begin(), cpe(), data, detSet(), edmNew::DetSet< T >::end(), filteredRecHits(), StripCPE::getAlgoParam(), hit(), index(), isActive(), isEmpty(), TrajectoryStateOnSurface::localParameters(), TrajectoryStateOnSurface::localPosition(), edmNew::DetSet< T >::makeRefTo(), Topology::measurementPosition(), mps_fire::result, specificGeomDet(), StripGeomDetUnit::specificTopology(), createJobs::tmp, UNLIKELY, and x.

◆ set128StripStatus()

void TkStripMeasurementDet::set128StripStatus	(	bool	good,
		int	idx = `-1`
	)

inline

Sets the status of a block of 128 strips (or all blocks if idx=-1)

Definition at line 230 of file TkStripMeasurementDet.h.

230 { conditionSet().set128StripStatus(index(), good, idx); }

References heavyIonCSV_trainingSettings::idx.

◆ setActiveThisEvent()

void TkStripMeasurementDet::setActiveThisEvent	(	StMeasurementDetSet &	theDets,
		bool	active
	)		const

inline

Turn on/off the module for reconstruction for one events. This per-event flag is cleared by any call to 'update' or 'setEmpty'

Definition at line 222 of file TkStripMeasurementDet.h.

                                                                            {
     theDets.setActiveThisEvent(index(), active);
   }

References StMeasurementDetSet::setActiveThisEvent().

◆ setActiveThisPeriod()

void TkStripMeasurementDet::setActiveThisPeriod	(	StMeasurementDetSet &	theDets,
		bool	active
	)

inline

Turn on/off the module for reconstruction, for the full run or lumi (using info from DB, usually).

Definition at line 218 of file TkStripMeasurementDet.h.

218 { conditionSet().setActive(index(), active); }

◆ setEmpty()

void TkStripMeasurementDet::setEmpty ( StMeasurementDetSet & theDets ) const

inline

Definition at line 92 of file TkStripMeasurementDet.h.

92 { theDets.setEmpty(index()); }

References StMeasurementDetSet::setEmpty().

◆ setIndex()

void TkStripMeasurementDet::setIndex ( int i )

inline

Definition at line 90 of file TkStripMeasurementDet.h.

90 { index_ = i; }

References mps_fire::i.

◆ simpleRecHits() [1/2]

bool TkStripMeasurementDet::simpleRecHits	(	const TrajectoryStateOnSurface &	ts,
		const MeasurementEstimator &	est,
		const MeasurementTrackerEvent &	data,
		std::vector< SiStripRecHit2D > &	result
	)		const

Definition at line 103 of file TkStripMeasurementDet.cc.

                                                                                     {
   if UNLIKELY ((!isActive(data)) || isEmpty(data.stripData()))
     return false;
  
   auto oldSize = result.size();
  
   int utraj = specificGeomDet().specificTopology().measurementPosition(stateOnThisDet.localPosition()).x();
   const detset& detSet = data.stripData().detSet(index());
   auto const& cpepar = cpe()->getAlgoParam(specificGeomDet(), stateOnThisDet.localParameters());
  
   auto rightCluster = std::find_if(
       detSet.begin(), detSet.end(), [utraj](const SiStripCluster& hit) { return hit.firstStrip() > utraj; });
  
   if (rightCluster != detSet.begin()) {
     // there are hits on the left of the utraj
     auto leftCluster = rightCluster;
     while (--leftCluster >= detSet.begin()) {
       SiStripClusterRef clusterref = detSet.makeRefTo(data.stripData().handle(), leftCluster);
       bool isCompatible = filteredRecHits(clusterref, cpepar, stateOnThisDet, est, data.stripClustersToSkip(), result);
       if (!isCompatible)
         break;  // exit loop on first incompatible hit
     }
   }
   for (; rightCluster != detSet.end(); rightCluster++) {
     SiStripClusterRef clusterref = detSet.makeRefTo(data.stripData().handle(), rightCluster);
     bool isCompatible = filteredRecHits(clusterref, cpepar, stateOnThisDet, est, data.stripClustersToSkip(), result);
     if (!isCompatible)
       break;  // exit loop on first incompatible hit
   }
  
   return result.size() > oldSize;
 }

References edmNew::DetSet< T >::begin(), cpe(), data, detSet(), edmNew::DetSet< T >::end(), filteredRecHits(), StripCPE::getAlgoParam(), hit(), index(), isActive(), isEmpty(), TrajectoryStateOnSurface::localParameters(), TrajectoryStateOnSurface::localPosition(), edmNew::DetSet< T >::makeRefTo(), Topology::measurementPosition(), mps_fire::result, specificGeomDet(), StripGeomDetUnit::specificTopology(), UNLIKELY, and x.

◆ simpleRecHits() [2/2]

void TkStripMeasurementDet::simpleRecHits	(	const TrajectoryStateOnSurface &	ts,
		const MeasurementTrackerEvent &	data,
		std::vector< SiStripRecHit2D > &	result
	)		const

Definition at line 219 of file TkStripMeasurementDet.cc.

                                                                                     {
   if (isEmpty(data.stripData()) || !isActive(data))
     return;
  
   const detset& detSet = data.stripData().detSet(index());
   unInitDynArray(AClusters::value_type, detSet.size(), clusters);
   assert(clusters.empty());
   for (auto const& ci : detSet) {
     if (isMasked(ci))
       continue;
     if (accept(detSet.makeKeyOf(&ci), data.stripClustersToSkip()))
       clusters.push_back(&ci);
     else
       LogDebug("TkStripMeasurementDet") << "skipping this str from last iteration on" << rawId()
                                         << " key: " << detSet.makeKeyOf(&ci);
   }
   if (!clusters.empty())
     buildSimpleRecHits(clusters, data, detSet, ts, result);
 }

References accept(), cms::cuda::assert(), buildSimpleRecHits(), bsc_activity_cfg::clusters, data, detSet(), index(), isActive(), isEmpty(), isMasked(), LogDebug, edmNew::DetSet< T >::makeKeyOf(), rawId(), mps_fire::result, edmNew::DetSet< T >::size(), and unInitDynArray.

Referenced by TkGluedMeasurementDet::collectRecHits().

◆ specificGeomDet()

const StripGeomDetUnit& TkStripMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 147 of file TkStripMeasurementDet.h.

147 { return static_cast<StripGeomDetUnit const&>(fastGeomDet()); }

References MeasurementDet::fastGeomDet().

Referenced by hit(), measurements(), recHits(), simpleRecHits(), and TkGluedMeasurementDet::testStrips().

◆ subId()

unsigned char TkStripMeasurementDet::subId ( ) const

inline

Definition at line 99 of file TkStripMeasurementDet.h.

99 { return conditionSet().subId(index()); }

◆ testStrips()

bool TkStripMeasurementDet::testStrips	(	float	utraj,
		float	uerr
	)		const

return true if there are 'enough' good strips in the utraj +/- 3 uerr range.

Definition at line 274 of file TkStripMeasurementDet.cc.

                                                                     {
   int16_t start = (int16_t)std::max<float>(utraj - 3.f * uerr, 0);
   int16_t end = (int16_t)std::min<float>(utraj + 3.f * uerr, totalStrips());
  
   if (start >= end) {  // which means either end <=0 or start >= totalStrips_
     /* LogDebug("TkStripMeasurementDet") << "Testing module " << id_ <<","<<
             " U = " << utraj << " +/- " << uerr << 
             "; Range [" << (utraj - 3*uerr) << ", " << (utraj + 3*uerr) << "] " << 
             ": YOU'RE COMPLETELY OFF THE MODULE."; */
     //return false;
     return true;  // Wolfgang thinks this way is better
                   // and solves some problems with grouped ckf
   }
  
   typedef std::vector<BadStripBlock>::const_iterator BSBIT;
  
   int16_t bad = 0, largestBadBlock = 0;
   for (BSBIT bsbc = badStripBlocks().begin(), bsbe = badStripBlocks().end(); bsbc != bsbe; ++bsbc) {
     if (bsbc->last < start)
       continue;
     if (bsbc->first > end)
       break;
     int16_t thisBad = std::min(bsbc->last, end) - std::max(bsbc->first, start);
     if (thisBad > largestBadBlock)
       largestBadBlock = thisBad;
     bad += thisBad;
   }
  
   bool ok = (bad < (end - start)) && (uint16_t(bad) <= badStripCuts().maxBad) &&
             (uint16_t(largestBadBlock) <= badStripCuts().maxConsecutiveBad);
  
   //    if (bad) {
   //       edm::LogWarning("TkStripMeasurementDet") << "Testing module " << id_ <<" (subdet: "<< SiStripDetId(id_).subdetId() << ")" <<
   //            " U = " << utraj << " +/- " << uerr <<
   //            "; Range [" << (utraj - 3*uerr) << ", " << (utraj + 3*uerr) << "] " <<
   //            "= [" << start << "," << end << "]" <<
   //            " total strips:" << (end-start) << ", good:" << (end-start-bad) << ", bad:" << bad << ", largestBadBlock: " << largestBadBlock <<
   //            ". " << (ok ? "OK" : "NO");
   //    }
   return ok;
 }

References badStripBlocks(), badStripCuts(), mps_fire::end, f, SiStripPI::max, StMeasurementConditionSet::BadStripCuts::maxBad, StMeasurementConditionSet::BadStripCuts::maxConsecutiveBad, min(), convertSQLiteXML::ok, command_line::start, and totalStrips().

Referenced by measurements(), and TkGluedMeasurementDet::testStrips().

◆ theSet()

const detset& TkStripMeasurementDet::theSet ( const StMeasurementDetSet & theDets ) const

inline

Definition at line 101 of file TkStripMeasurementDet.h.

101 { return theDets.detSet(index()); }

References StMeasurementDetSet::detSet().

◆ totalStrips()

int TkStripMeasurementDet::totalStrips ( ) const

inlineprivate

Definition at line 256 of file TkStripMeasurementDet.h.

256 { return conditionSet().totalStrips(index()); }

Referenced by testStrips().

Member Data Documentation

◆ index_

int TkStripMeasurementDet::index_

private

Definition at line 248 of file TkStripMeasurementDet.h.

◆ theDetConditions

StMeasurementConditionSet* TkStripMeasurementDet::theDetConditions

private

Definition at line 249 of file TkStripMeasurementDet.h.

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ AClusters

◆ ALocalValues

◆ BadStripBlock

◆ BadStripCuts

◆ const_iterator

◆ detset

◆ LocalValues

◆ new_const_iterator

◆ SiStripClusterRef

◆ VLocalValues

Constructor & Destructor Documentation

◆ ~TkStripMeasurementDet()

◆ TkStripMeasurementDet()

Member Function Documentation

◆ accept() [1/2]

◆ accept() [2/2]

◆ advance()

◆ badStripBlocks()

◆ badStripCuts()

◆ buildRecHit()

◆ buildRecHits()

◆ buildSimpleRecHits()

◆ conditionSet() [1/2]

◆ conditionSet() [2/2]

◆ cpe()

◆ detSet()

◆ empty()

◆ filteredRecHits() [1/2]

◆ filteredRecHits() [2/2]

◆ getBadStripBlocks()

◆ hasAllGoodChannels()

◆ hasAny128StripBad()

◆ hasBadComponents()

◆ hit()

◆ hitRange()

◆ index()

◆ isActive()

◆ isEmpty()

◆ isMasked()

◆ maskBad128StripBlocks()

◆ measurements()

◆ rawId()

◆ recHits() [1/3]

◆ recHits() [2/3]

◆ recHits() [3/3]

◆ set128StripStatus()

◆ setActiveThisEvent()

◆ setActiveThisPeriod()

◆ setEmpty()

◆ setIndex()

◆ simpleRecHits() [1/2]

◆ simpleRecHits() [2/2]

◆ specificGeomDet()

◆ subId()

◆ testStrips()

◆ theSet()

◆ totalStrips()

Member Data Documentation

◆ index_

◆ theDetConditions