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

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.

{}

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.

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

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

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

accept() [2/2]

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

inline

Definition at line 281 of file TkStripMeasurementDet.h.

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

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

advance()

void TkStripMeasurementDet::advance

(

TkStripRecHitIter &

hi

)

const

Definition at line 250 of file TkStripMeasurementDet.cc.

References accept(), data, 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;
  }
}

badStripBlocks()

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

inline

Definition at line 240 of file TkStripMeasurementDet.h.

Referenced by testStrips().

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

badStripCuts()

BadStripCuts const& TkStripMeasurementDet::badStripCuts ( ) const

inlineprivate

Definition at line 257 of file TkStripMeasurementDet.h.

Referenced by testStrips().

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

buildRecHit()

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

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.

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

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

conditionSet() [1/2]

StMeasurementConditionSet& TkStripMeasurementDet::conditionSet ( )

inlineprivate

Definition at line 250 of file TkStripMeasurementDet.h.

{ return *theDetConditions; }

conditionSet() [2/2]

const StMeasurementConditionSet& TkStripMeasurementDet::conditionSet ( ) const

inlineprivate

Definition at line 251 of file TkStripMeasurementDet.h.

{ return *theDetConditions; }

cpe()

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

detSet()

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

empty()

bool TkStripMeasurementDet::empty

(

const MeasurementTrackerEvent &

data

)

const

Definition at line 24 of file TkStripMeasurementDet.cc.

References accept(), edmNew::DetSet< T >::begin(), data, detSet(), edmNew::DetSet< T >::end(), index(), isActive(), isEmpty(), isMasked(), edmNew::DetSet< T >::makeRefTo(), 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;
}

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.

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

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

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.

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

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

getBadStripBlocks()

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

inline

Definition at line 239 of file TkStripMeasurementDet.h.

{ 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.

Referenced by TkGluedMeasurementDet::measurements().

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

hasAny128StripBad()

bool TkStripMeasurementDet::hasAny128StripBad ( ) const

inlineprivate

Definition at line 259 of file TkStripMeasurementDet.h.

{ 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.

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

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

hitRange()

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(), edmNew::DetSet< T >::end(), index(), isActive(), and isEmpty().

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

index()

int TkStripMeasurementDet::index ( ) const

inline

Definition at line 96 of file TkStripMeasurementDet.h.

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

{ return index_; }

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.

References data.

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

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

isEmpty()

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

isMasked()

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

maskBad128StripBlocks()

bool TkStripMeasurementDet::maskBad128StripBlocks ( ) const

inline

Definition at line 242 of file TkStripMeasurementDet.h.

{ 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.

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.

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

rawId()

unsigned int TkStripMeasurementDet::rawId ( ) const

inline

Definition at line 98 of file TkStripMeasurementDet.h.

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

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

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.

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().

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

recHits() [2/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.

References edmNew::DetSet< T >::begin(), cpe(), data, detSet(), edmNew::DetSet< T >::end(), filteredRecHits(), StripCPE::getAlgoParam(), h, 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.

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

recHits() [3/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.

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.

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

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.

References caHitNtupletGeneratorKernels::good, and heavyIonCSV_trainingSettings::idx.

{ conditionSet().set128StripStatus(index(), good, 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.

References StMeasurementDetSet::setActiveThisEvent().

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

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.

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

setEmpty()

void TkStripMeasurementDet::setEmpty ( StMeasurementDetSet &  theDets ) const

inline

Definition at line 92 of file TkStripMeasurementDet.h.

References StMeasurementDetSet::setEmpty().

{ theDets.setEmpty(index()); }

setIndex()

void TkStripMeasurementDet::setIndex ( int  i )

inline

Definition at line 90 of file TkStripMeasurementDet.h.

References mps_fire::i.

{ index_ = i; }

simpleRecHits() [1/2]

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(), 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().

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

simpleRecHits() [2/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.

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.

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

specificGeomDet()

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

subId()

unsigned char TkStripMeasurementDet::subId ( ) const

inline

Definition at line 99 of file TkStripMeasurementDet.h.

{ 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.

References badStripBlocks(), badStripCuts(), mps_fire::end, f, SiStripPI::max, StMeasurementConditionSet::BadStripCuts::maxBad, StMeasurementConditionSet::BadStripCuts::maxConsecutiveBad, SiStripPI::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;
}

theSet()

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

inline

Definition at line 101 of file TkStripMeasurementDet.h.

References StMeasurementDetSet::detSet().

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

totalStrips()

int TkStripMeasurementDet::totalStrips ( ) const

inlineprivate

Definition at line 256 of file TkStripMeasurementDet.h.

Referenced by testStrips().

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

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.