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 86 of file TkStripMeasurementDet.h.

Member Typedef Documentation

using TkStripMeasurementDet::AClusters = StripClusterParameterEstimator::AClusters

private

Definition at line 239 of file TkStripMeasurementDet.h.

using TkStripMeasurementDet::ALocalValues = StripClusterParameterEstimator::ALocalValues

private

Definition at line 240 of file TkStripMeasurementDet.h.

typedef StMeasurementConditionSet::BadStripBlock TkStripMeasurementDet::BadStripBlock

Definition at line 231 of file TkStripMeasurementDet.h.

typedef StMeasurementConditionSet::BadStripCuts TkStripMeasurementDet::BadStripCuts

Definition at line 226 of file TkStripMeasurementDet.h.

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

Definition at line 97 of file TkStripMeasurementDet.h.

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

Definition at line 94 of file TkStripMeasurementDet.h.

typedef StripClusterParameterEstimator::LocalValues TkStripMeasurementDet::LocalValues

Definition at line 89 of file TkStripMeasurementDet.h.

typedef detset::const_iterator TkStripMeasurementDet::new_const_iterator

Definition at line 95 of file TkStripMeasurementDet.h.

typedef SiStripRecHit2D::ClusterRef TkStripMeasurementDet::SiStripClusterRef

Definition at line 92 of file TkStripMeasurementDet.h.

typedef StripClusterParameterEstimator::VLocalValues TkStripMeasurementDet::VLocalValues

Definition at line 90 of file TkStripMeasurementDet.h.

Constructor & Destructor Documentation

TkStripMeasurementDet::~TkStripMeasurementDet ( )

inlineoverride

Definition at line 99 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 282 of file TkStripMeasurementDet.h.

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

Referenced by advance(), Vispa.Gui.BoxContentDialog.BoxContentDialog::apply(), Vispa.Plugins.ConfigEditor.ToolDialog.ToolDialog::apply(), 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 286 of file TkStripMeasurementDet.h.

References LogDebug.

Referenced by Vispa.Gui.BoxContentDialog.BoxContentDialog::apply(), Vispa.Plugins.ConfigEditor.ToolDialog.ToolDialog::apply(), and 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 242 of file TkStripMeasurementDet.cc.

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

                                                                 {
  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 234 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::badStripBlocks().

Referenced by testStrips().

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

BadStripCuts const& TkStripMeasurementDet::badStripCuts ( ) const

inlineprivate

Definition at line 252 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::badStripCuts().

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 156 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 165 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 262 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 245 of file TkStripMeasurementDet.h.

{ return *theDetConditions; }

const StMeasurementConditionSet& TkStripMeasurementDet::conditionSet ( ) const

inlineprivate

Definition at line 246 of file TkStripMeasurementDet.h.

{ return *theDetConditions; }

const StripCPE* TkStripMeasurementDet::cpe ( ) const

inlineprivate

Definition at line 248 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::stripCPE().

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 116 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 26 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 175 of file TkStripMeasurementDet.h.

References accept(), MeasurementEstimator::estimate(), LogDebug, eostools::move(), MeasurementEstimator::preFilter(), and rpcPointValidation_cfi::recHit.

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::move(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 194 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 233 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::getBadStripBlocks().

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

bool TkStripMeasurementDet::hasAllGoodChannels ( ) const

inline

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

Definition at line 219 of file TkStripMeasurementDet.h.

Referenced by TkGluedMeasurementDet::measurements().

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

bool TkStripMeasurementDet::hasAny128StripBad ( ) const

inlineprivate

Definition at line 254 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::hasAny128StripBad().

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

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

inlineoverridevirtual

Implements MeasurementDet.

Definition at line 123 of file TkStripMeasurementDet.h.

References data, relativeConstraints::empty, MeasurementDet::measurements(), MeasurementDet::recHits(), and mps_fire::result.

{return false;}

SiStripRecHit2D TkStripMeasurementDet::hit

(

TkStripRecHitIter const &

hi

)

const

Definition at line 253 of file TkStripMeasurementDet.cc.

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

Referenced by 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 233 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 109 of file TkStripMeasurementDet.h.

Referenced by empty(), hitRange(), BeautifulSoup.PageElement::insert(), 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 120 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 107 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 257 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::isMasked().

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

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

bool TkStripMeasurementDet::maskBad128StripBlocks ( ) const

inline

Definition at line 236 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::maskBad128StripBlocks().

{ return conditionSet().maskBad128StripBlocks();}

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

overridevirtual

Implements MeasurementDet.

Definition at line 167 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(), TrackInfoProducer_cfi::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 111 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::id().

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 41 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(), StMeasurementDetSet::handle(), hit(), index(), isActive(), isEmpty(), TrajectoryStateOnSurface::localParameters(), TrajectoryStateOnSurface::localPosition(), edmNew::DetSet< T >::makeRefTo(), Topology::measurementPosition(), eostools::move(), specificGeomDet(), StripGeomDetUnit::specificTopology(), MeasurementTrackerEvent::stripClustersToSkip(), MeasurementTrackerEvent::stripData(), 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(std::move(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(std::move(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 136 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 222 of file TkStripMeasurementDet.h.

References training_settings::idx, and StMeasurementConditionSet::set128StripStatus().

                                                {
    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 216 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 212 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::setActive().

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

void TkStripMeasurementDet::setEmpty ( StMeasurementDetSet &  theDets ) const

inline

Definition at line 105 of file TkStripMeasurementDet.h.

References StMeasurementDetSet::setEmpty().

{ theDets.setEmpty(index()); }

void TkStripMeasurementDet::setIndex ( int  i )

inline

Definition at line 103 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 213 of file TkStripMeasurementDet.cc.

References accept(), buildSimpleRecHits(), fastPrimaryVertexProducer_cfi::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 101 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 151 of file TkStripMeasurementDet.h.

References 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 112 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::subId().

{ 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 265 of file TkStripMeasurementDet.cc.

References badStripBlocks(), badStripCuts(), end, 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 115 of file TkStripMeasurementDet.h.

References StMeasurementDetSet::detSet().

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

int TkStripMeasurementDet::totalStrips ( ) const

inlineprivate

Definition at line 251 of file TkStripMeasurementDet.h.

References StMeasurementConditionSet::totalStrips().

Referenced by testStrips().

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

Member Data Documentation

int TkStripMeasurementDet::index_

private

Definition at line 243 of file TkStripMeasurementDet.h.

StMeasurementConditionSet* TkStripMeasurementDet::theDetConditions

private

Definition at line 244 of file TkStripMeasurementDet.h.