TkStripMeasurementDet Class Reference

#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 (SimpleHitContainer &result, const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &, const MeasurementTrackerEvent &data) const override
bool	recHits (const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &est, const MeasurementTrackerEvent &data, RecHitContainer &result, std::vector< float > &diffs) 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)
void	simpleRecHits (const TrajectoryStateOnSurface &ts, const MeasurementTrackerEvent &data, std::vector< SiStripRecHit2D > &result) const
bool	simpleRecHits (const TrajectoryStateOnSurface &ts, const MeasurementEstimator &est, 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

using TkStripMeasurementDet::AClusters = StripClusterParameterEstimator::AClusters

private

Definition at line 245 of file TkStripMeasurementDet.h.

using TkStripMeasurementDet::ALocalValues = StripClusterParameterEstimator::ALocalValues

private

Definition at line 246 of file TkStripMeasurementDet.h.

typedef StMeasurementConditionSet::BadStripBlock TkStripMeasurementDet::BadStripBlock

Definition at line 237 of file TkStripMeasurementDet.h.

typedef StMeasurementConditionSet::BadStripCuts TkStripMeasurementDet::BadStripCuts

Definition at line 232 of file TkStripMeasurementDet.h.

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

Definition at line 84 of file TkStripMeasurementDet.h.

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

Definition at line 81 of file TkStripMeasurementDet.h.

typedef StripClusterParameterEstimator::LocalValues TkStripMeasurementDet::LocalValues

Definition at line 76 of file TkStripMeasurementDet.h.

typedef detset::const_iterator TkStripMeasurementDet::new_const_iterator

Definition at line 82 of file TkStripMeasurementDet.h.

typedef SiStripRecHit2D::ClusterRef TkStripMeasurementDet::SiStripClusterRef

Definition at line 79 of file TkStripMeasurementDet.h.

typedef StripClusterParameterEstimator::VLocalValues TkStripMeasurementDet::VLocalValues

Definition at line 77 of file TkStripMeasurementDet.h.

Constructor & Destructor Documentation

TkStripMeasurementDet::~TkStripMeasurementDet ( )

inlineoverride

Definition at line 86 of file TkStripMeasurementDet.h.

{}

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

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

inline

Definition at line 277 of file TkStripMeasurementDet.h.

References accept(), edm::Ref< C, T, F >::key(), and HLT_FULL_cff::skipClusters.

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

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

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

inline

Definition at line 281 of file TkStripMeasurementDet.h.

References LogDebug.

Referenced by esMonitoring.FDJsonServer::handle_accept().

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

void TkStripMeasurementDet::advance

(

TkStripRecHitIter &

hi

)

const

Definition at line 250 of file TkStripMeasurementDet.cc.

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

Referenced by TkStripRecHitIter::advance().

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

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

inline

Definition at line 240 of file TkStripMeasurementDet.h.

Referenced by testStrips().

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

BadStripCuts const& TkStripMeasurementDet::badStripCuts ( ) const

inlineprivate

Definition at line 257 of file TkStripMeasurementDet.h.

Referenced by testStrips().

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

template<class ClusterRefT >

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

inline

Definition at line 150 of file TkStripMeasurementDet.h.

References MeasurementDet::fastGeomDet().

Referenced by recHits().

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

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

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.

References declareDynArray, StMeasurementDetSet::handle(), mps_fire::i, TrajectoryStateOnSurface::localParameters(), edmNew::DetSet< T >::makeRefTo(), DynArray< T >::size(), and MeasurementTrackerEvent::stripData().

Referenced by simpleRecHits().

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

StMeasurementConditionSet& TkStripMeasurementDet::conditionSet ( )

inlineprivate

Definition at line 250 of file TkStripMeasurementDet.h.

{ return *theDetConditions; }

const StMeasurementConditionSet& TkStripMeasurementDet::conditionSet ( ) const

inlineprivate

Definition at line 251 of file TkStripMeasurementDet.h.

{ return *theDetConditions; }

const StripCPE* TkStripMeasurementDet::cpe ( ) const

inlineprivate

Definition at line 253 of file TkStripMeasurementDet.h.

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

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

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

inline

Definition at line 102 of file TkStripMeasurementDet.h.

References StMeasurementDetSet::detSet().

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

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

bool TkStripMeasurementDet::empty

(

const MeasurementTrackerEvent &

data

)

const

Definition at line 24 of file TkStripMeasurementDet.cc.

References accept(), edmNew::DetSet< T >::begin(), detSet(), StMeasurementDetSet::detSet(), edmNew::DetSet< T >::end(), StMeasurementDetSet::handle(), index(), isActive(), isEmpty(), isMasked(), edmNew::DetSet< T >::makeRefTo(), MeasurementTrackerEvent::stripClustersToSkip(), MeasurementTrackerEvent::stripData(), and UNLIKELY.

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

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.

References accept(), MeasurementEstimator::estimate(), LogDebug, and MeasurementEstimator::preFilter().

Referenced by recHits(), and simpleRecHits().

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

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.

References accept(), MeasurementEstimator::estimate(), LogDebug, and MeasurementEstimator::preFilter().

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

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

inline

Definition at line 239 of file TkStripMeasurementDet.h.

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

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.

Referenced by TkGluedMeasurementDet::measurements().

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

bool TkStripMeasurementDet::hasAny128StripBad ( ) const

inlineprivate

Definition at line 259 of file TkStripMeasurementDet.h.

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

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

inlineoverridevirtual

Implements MeasurementDet.

Definition at line 108 of file TkStripMeasurementDet.h.

                                                                                                                  {
    return false;
  }

SiStripRecHit2D TkStripMeasurementDet::hit

(

TkStripRecHitIter const &

hi

)

const

Definition at line 263 of file TkStripMeasurementDet.cc.

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

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

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

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

Definition at line 241 of file TkStripMeasurementDet.cc.

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

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

int TkStripMeasurementDet::index ( ) const

inline

Definition at line 96 of file TkStripMeasurementDet.h.

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

{ return index_; }

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.

References StMeasurementDetSet::isActive(), and MeasurementTrackerEvent::stripData().

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

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

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

inline

Definition at line 94 of file TkStripMeasurementDet.h.

References StMeasurementDetSet::empty().

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

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

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

inlineprivate

Definition at line 261 of file TkStripMeasurementDet.h.

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

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

bool TkStripMeasurementDet::maskBad128StripBlocks ( ) const

inline

Definition at line 242 of file TkStripMeasurementDet.h.

{ return conditionSet().maskBad128StripBlocks(); }

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.

References tracking::TempMeasurements::add(), tracking::TempMeasurements::distances, F(), TrajectoryStateOnSurface::hasError(), tracking::TempMeasurements::hits, isActive(), isEmpty(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), LogDebug, Topology::measurementPosition(), LocalTrajectoryError::positionError(), rawId(), HI_PhotonSkim_cff::rechits, recHits(), specificGeomDet(), StripGeomDetUnit::specificTopology(), mathSSE::sqrt(), MeasurementTrackerEvent::stripData(), testStrips(), MeasurementDet::theInactiveHit, MeasurementDet::theMissingHit, and x.

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

unsigned int TkStripMeasurementDet::rawId ( ) const

inline

Definition at line 98 of file TkStripMeasurementDet.h.

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

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

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

overridevirtual

Implements MeasurementDet.

Definition at line 39 of file TkStripMeasurementDet.cc.

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

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

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

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.

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

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

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.

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

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

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.

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

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.

References StMeasurementDetSet::setActiveThisEvent().

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

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.

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

void TkStripMeasurementDet::setEmpty ( StMeasurementDetSet &  theDets ) const

inline

Definition at line 92 of file TkStripMeasurementDet.h.

References StMeasurementDetSet::setEmpty().

{ theDets.setEmpty(index()); }

void TkStripMeasurementDet::setIndex ( int  i )

inline

Definition at line 90 of file TkStripMeasurementDet.h.

References mps_fire::i.

{ index_ = i; }

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

Definition at line 219 of file TkStripMeasurementDet.cc.

References accept(), cms::cuda::assert(), buildSimpleRecHits(), HLT_FULL_cff::clusters, detSet(), StMeasurementDetSet::detSet(), index(), isActive(), isEmpty(), isMasked(), LogDebug, rawId(), edmNew::DetSet< T >::size(), MeasurementTrackerEvent::stripClustersToSkip(), MeasurementTrackerEvent::stripData(), and unInitDynArray.

Referenced by TkGluedMeasurementDet::collectRecHits().

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

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.

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

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

const StripGeomDetUnit& TkStripMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 147 of file TkStripMeasurementDet.h.

References runTheMatrix::const, and MeasurementDet::fastGeomDet().

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

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

unsigned char TkStripMeasurementDet::subId ( ) const

inline

Definition at line 99 of file TkStripMeasurementDet.h.

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

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.

References badStripBlocks(), badStripCuts(), dataset::end, validate-o2o-wbm::f, SiStripPI::max, StMeasurementConditionSet::BadStripCuts::maxBad, StMeasurementConditionSet::BadStripCuts::maxConsecutiveBad, min(), convertSQLiteXML::ok, command_line::start, and totalStrips().

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

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

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

inline

Definition at line 101 of file TkStripMeasurementDet.h.

References StMeasurementDetSet::detSet().

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

int TkStripMeasurementDet::totalStrips ( ) const

inlineprivate

Definition at line 256 of file TkStripMeasurementDet.h.

Referenced by testStrips().

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

Member Data Documentation

int TkStripMeasurementDet::index_

private

Definition at line 248 of file TkStripMeasurementDet.h.

StMeasurementConditionSet* TkStripMeasurementDet::theDetConditions

private

Definition at line 249 of file TkStripMeasurementDet.h.