TkStripMeasurementDet Class Reference

#include <TkStripMeasurementDet.h>

Inheritance diagram for TkStripMeasurementDet:

Classes
struct	BadStripBlock
struct	BadStripCuts

Public Types
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 edm::LazyGetter < SiStripCluster >::value_ref	SiStripRegionalClusterRef
typedef StripClusterParameterEstimator::VLocalValues	VLocalValues
Public Types inherited from MeasurementDet
typedef TransientTrackingRecHit::ConstRecHitContainer	RecHitContainer

Public Member Functions
bool	accept (SiStripClusterRef &r) const
bool	accept (SiStripRegionalClusterRef &r) const
TransientTrackingRecHit::RecHitPointer	buildRecHit (const SiStripClusterRef &, const TrajectoryStateOnSurface &ltp) const
TransientTrackingRecHit::RecHitPointer	buildRecHit (const SiStripRegionalClusterRef &, const TrajectoryStateOnSurface &ltp) const
TkStripMeasurementDet::RecHitContainer	buildRecHits (const SiStripClusterRef &, const TrajectoryStateOnSurface &ltp) const
TkStripMeasurementDet::RecHitContainer	buildRecHits (const SiStripRegionalClusterRef &, const TrajectoryStateOnSurface &ltp) const
virtual std::vector < TrajectoryMeasurement >	fastMeasurements (const TrajectoryStateOnSurface &stateOnThisDet, const TrajectoryStateOnSurface &startingState, const Propagator &, const MeasurementEstimator &) 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
bool	isActive () const
	Is this module active in reconstruction? It must be both 'setActiveThisEvent' and 'setActive'. More...
bool	isEmpty ()
unsigned int	rawId () const
virtual RecHitContainer	recHits (const TrajectoryStateOnSurface &) const
void	set128StripStatus (bool good, int idx=-1)
	Sets the status of a block of 128 strips (or all blocks if idx=-1) More...
void	setActive (bool active)
	Turn on/off the module for reconstruction, for the full run or lumi (using info from DB, usually). This also resets the 'setActiveThisEvent' to true. More...
void	setActiveThisEvent (bool active)
	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	setBadStripCuts (BadStripCuts cuts)
template<typename IT >
void	setClusterToSkip (IT begin, IT end)
void	setEmpty ()
void	setMaskBad128StripBlocks (bool maskThem)
template<typename IT >
void	setRegionalClustersToSkip (IT begin, IT end)
void	simpleRecHits (const TrajectoryStateOnSurface &ts, std::vector< SiStripRecHit2D > &result) const
int	size () const
const StripGeomDetUnit &	specificGeomDet () 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 ()
	TkStripMeasurementDet (const GeomDet *gdet, const StripClusterParameterEstimator *cpe, bool regional)
void	unset ()
void	update (const detset &detSet, const edm::Handle< edmNew::DetSetVector< SiStripCluster > > h, unsigned int id)
void	update (std::vector< SiStripCluster >::const_iterator begin, std::vector< SiStripCluster >::const_iterator end, const edm::Handle< edm::LazyGetter< SiStripCluster > > h, unsigned int id)
virtual	~TkStripMeasurementDet ()
Public Member Functions inherited from MeasurementDet
virtual const GeomDet &	geomDet () const
	MeasurementDet (const GeomDet *gdet)
const Surface::PositionType &	position () const
const BoundSurface &	surface () const

Private Member Functions
template<class ClusterRefT >
void	buildSimpleRecHit (const ClusterRefT &cluster, const TrajectoryStateOnSurface &ltp, std::vector< SiStripRecHit2D > &res) const
bool	isMasked (const SiStripCluster &cluster) const

Private Attributes
bool	activeThisEvent_
bool	activeThisPeriod_
bool	bad128Strip_ [6]
std::vector< BadStripBlock >	badStripBlocks_
BadStripCuts	badStripCuts_
uint	beginClusterI_
detset	detSet_
bool	empty
uint	endClusterI_
edm::Handle < edmNew::DetSetVector < SiStripCluster > >	handle_
bool	hasAny128StripBad_
unsigned int	id_
bool	isRegional
bool	maskBad128StripBlocks_
edm::Handle< edm::LazyGetter < SiStripCluster > >	regionalHandle_
std::set< SiStripClusterRef >	skipClusters_
std::set < SiStripRegionalClusterRef >	skipRegClusters_
const StripClusterParameterEstimator *	theCPE
const StripGeomDetUnit *	theStripGDU
int	totalStrips_

Detailed Description

Definition at line 22 of file TkStripMeasurementDet.h.

Member Typedef Documentation

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

Definition at line 35 of file TkStripMeasurementDet.h.

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

Definition at line 32 of file TkStripMeasurementDet.h.

typedef StripClusterParameterEstimator::LocalValues TkStripMeasurementDet::LocalValues

Definition at line 25 of file TkStripMeasurementDet.h.

typedef detset::const_iterator TkStripMeasurementDet::new_const_iterator

Definition at line 33 of file TkStripMeasurementDet.h.

typedef SiStripRecHit2D::ClusterRef TkStripMeasurementDet::SiStripClusterRef

Definition at line 28 of file TkStripMeasurementDet.h.

typedef edm::LazyGetter<SiStripCluster>::value_ref TkStripMeasurementDet::SiStripRegionalClusterRef

Definition at line 30 of file TkStripMeasurementDet.h.

typedef StripClusterParameterEstimator::VLocalValues TkStripMeasurementDet::VLocalValues

Definition at line 26 of file TkStripMeasurementDet.h.

Constructor & Destructor Documentation

virtual TkStripMeasurementDet::~TkStripMeasurementDet ( )

inlinevirtual

Definition at line 37 of file TkStripMeasurementDet.h.

{}

TkStripMeasurementDet::TkStripMeasurementDet	(	const GeomDet *	gdet,
		const StripClusterParameterEstimator *	cpe,
		bool	regional
	)

Definition at line 17 of file TkStripMeasurementDet.cc.

References GeomDet::geographicalId(), id_, StripTopology::nstrips(), DetId::rawId(), specificGeomDet(), StripGeomDetUnit::specificTopology(), theStripGDU, and totalStrips_.

                                         : 
    MeasurementDet (gdet),
    isRegional(regional),
    empty(true),
    activeThisEvent_(true), activeThisPeriod_(true),
    theCPE(cpe)
  {
    theStripGDU = dynamic_cast<const StripGeomDetUnit*>(gdet);
    if (theStripGDU == 0) {
      throw MeasurementDetException( "TkStripMeasurementDet constructed with a GeomDet which is not a StripGeomDetUnit");
    }

    //intialize the detId !
    id_ = gdet->geographicalId().rawId();
    //initalize the total number of strips
    totalStrips_ =  specificGeomDet().specificTopology().nstrips();
  }

Member Function Documentation

bool TkStripMeasurementDet::accept ( SiStripClusterRef &  r ) const

inline

Definition at line 199 of file TkStripMeasurementDet.h.

References skipClusters_.

Referenced by python.Vispa.Gui.BoxContentDialog.BoxContentDialog::apply(), python.Vispa.Plugins.ConfigEditor.ToolDialog.ToolDialog::apply(), fastMeasurements(), recHits(), and simpleRecHits().

                                                  {
    return (skipClusters_.find(r) == skipClusters_.end());
  }

bool TkStripMeasurementDet::accept ( SiStripRegionalClusterRef &  r ) const

inline

Definition at line 202 of file TkStripMeasurementDet.h.

References skipRegClusters_.

Referenced by python.Vispa.Gui.BoxContentDialog.BoxContentDialog::apply(), and python.Vispa.Plugins.ConfigEditor.ToolDialog.ToolDialog::apply().

                                                        {
    return (skipRegClusters_.find(r) == skipRegClusters_.end());
  }

TransientTrackingRecHit::RecHitPointer TkStripMeasurementDet::buildRecHit	(	const SiStripClusterRef &	cluster,
		const TrajectoryStateOnSurface &	ltp
	)		const

Definition at line 195 of file TkStripMeasurementDet.cc.

References TSiStripRecHit2DLocalPos::build(), MeasurementDet::geomDet(), ClusterParameterEstimator< T >::localParameters(), specificGeomDet(), and theCPE.

Referenced by recHits().

{
  const GeomDetUnit& gdu( specificGeomDet());
  LocalValues lv = theCPE->localParameters( *cluster, gdu, ltp);
  return TSiStripRecHit2DLocalPos::build( lv.first, lv.second, &geomDet(), cluster, theCPE);
}

TransientTrackingRecHit::RecHitPointer TkStripMeasurementDet::buildRecHit	(	const SiStripRegionalClusterRef &	cluster,
		const TrajectoryStateOnSurface &	ltp
	)		const

Definition at line 204 of file TkStripMeasurementDet.cc.

References TSiStripRecHit2DLocalPos::build(), MeasurementDet::geomDet(), ClusterParameterEstimator< T >::localParameters(), specificGeomDet(), and theCPE.

{
  const GeomDetUnit& gdu( specificGeomDet());
  LocalValues lv = theCPE->localParameters( *cluster, gdu, ltp);
  return TSiStripRecHit2DLocalPos::build( lv.first, lv.second, &geomDet(), cluster, theCPE);
}

TkStripMeasurementDet::RecHitContainer TkStripMeasurementDet::buildRecHits	(	const SiStripClusterRef &	cluster,
		const TrajectoryStateOnSurface &	ltp
	)		const

Definition at line 215 of file TkStripMeasurementDet.cc.

References TSiStripRecHit2DLocalPos::build(), MeasurementDet::geomDet(), ClusterParameterEstimator< T >::localParametersV(), specificGeomDet(), and theCPE.

Referenced by fastMeasurements().

{
  const GeomDetUnit& gdu( specificGeomDet());
  VLocalValues vlv = theCPE->localParametersV( *cluster, gdu, ltp);
  RecHitContainer res;
  for(VLocalValues::const_iterator it=vlv.begin();it!=vlv.end();++it){
    res.push_back(TSiStripRecHit2DLocalPos::build( it->first, it->second, &geomDet(), cluster, theCPE));
  }
  return res; 
}

TkStripMeasurementDet::RecHitContainer TkStripMeasurementDet::buildRecHits	(	const SiStripRegionalClusterRef &	cluster,
		const TrajectoryStateOnSurface &	ltp
	)		const

Definition at line 229 of file TkStripMeasurementDet.cc.

References TSiStripRecHit2DLocalPos::build(), MeasurementDet::geomDet(), ClusterParameterEstimator< T >::localParametersV(), specificGeomDet(), and theCPE.

{
  const GeomDetUnit& gdu( specificGeomDet());
  VLocalValues vlv = theCPE->localParametersV( *cluster, gdu, ltp);
  RecHitContainer res;
  for(VLocalValues::const_iterator it=vlv.begin();it!=vlv.end();++it){
    res.push_back(TSiStripRecHit2DLocalPos::build( it->first, it->second, &geomDet(), cluster, theCPE));
  }
  return res; 
}

template<class ClusterRefT >

void TkStripMeasurementDet::buildSimpleRecHit ( const ClusterRefT &  cluster, const TrajectoryStateOnSurface &  ltp, std::vector< SiStripRecHit2D > &  res  ) const

private

Definition at line 276 of file TkStripMeasurementDet.cc.

References MeasurementDet::geomDet(), ClusterParameterEstimator< T >::localParametersV(), specificGeomDet(), and theCPE.

Referenced by simpleRecHits().

{
  const GeomDetUnit& gdu( specificGeomDet());
  VLocalValues vlv = theCPE->localParametersV( *cluster, gdu, ltp);
  for(VLocalValues::const_iterator it=vlv.begin();it!=vlv.end();++it){
    res.push_back(SiStripRecHit2D( it->first, it->second, geomDet().geographicalId(), cluster));
  }
}

std::vector< TrajectoryMeasurement > TkStripMeasurementDet::fastMeasurements ( const TrajectoryStateOnSurface &  stateOnThisDet, const TrajectoryStateOnSurface &  startingState, const Propagator &  , const MeasurementEstimator &    ) const

virtual

faster version in case the TrajectoryState on the surface of the Det is already available. The first TrajectoryStateOnSurface is on the surface of this Det, and the second TrajectoryStateOnSurface is the statrting state, usually not on the surface of this Det. The stateOnThisDet should the result of
prop.propagate( startingState, this->surface())

Implements MeasurementDet.

Definition at line 39 of file TkStripMeasurementDet.cc.

References accept(), edmNew::DetSet< T >::begin(), beginClusterI_, InvalidTransientRecHit::build(), buildRecHits(), detSet_, empty, edmNew::DetSet< T >::end(), endClusterI_, MeasurementEstimator::estimate(), GeomDet::geographicalId(), MeasurementDet::geomDet(), handle_, TrajectoryStateOnSurface::hasError(), TrackingRecHit::inactive, isActive(), isMasked(), isRegional, edm::Ref< C, T, F >::key(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), LogDebug, edmNew::makeRefTo(), edm::makeRefToLazyGetter(), Topology::measurementError(), Topology::measurementPosition(), TrackingRecHit::missing, LocalTrajectoryError::positionError(), DetId::rawId(), recHits(), regionalHandle_, query::result, size(), python.multivaluedict::sort(), StripGeomDetUnit::specificTopology(), mathSSE::sqrt(), testStrips(), theStripGDU, and x.

{ 
  std::vector<TrajectoryMeasurement> result;

  if (isActive() == false) {
    result.push_back( TrajectoryMeasurement( stateOnThisDet, 
            InvalidTransientRecHit::build(&geomDet(), TrackingRecHit::inactive), 
        0.F));
    //    LogDebug("TkStripMeasurementDet") << " DetID " << id_ << " inactive";
    return result;
  }
 
  float utraj =  theStripGDU->specificTopology().measurementPosition( stateOnThisDet.localPosition()).x();
  float uerr;
  //  if (theClusterRange.first == theClusterRange.second) { // empty
  if (empty  == true){
    //LogDebug("TkStripMeasurementDet") << " DetID " << id_ << " empty ";
    if (stateOnThisDet.hasError()){
    uerr= sqrt(theStripGDU->specificTopology().measurementError(stateOnThisDet.localPosition(),stateOnThisDet.localError().positionError()).uu());
     if (testStrips(utraj,uerr)) {
        result.push_back( TrajectoryMeasurement( stateOnThisDet, InvalidTransientRecHit::build(&geomDet(), TrackingRecHit::missing), 0.F));
     } else { 
        result.push_back( TrajectoryMeasurement( stateOnThisDet, InvalidTransientRecHit::build(&geomDet(), TrackingRecHit::inactive), 0.F));
     }
    }else{
      result.push_back( TrajectoryMeasurement( stateOnThisDet, InvalidTransientRecHit::build(&geomDet(), TrackingRecHit::missing), 0.F));
    }
    return result;
  }
  
  if(!isRegional){//old implemetation with DetSet
    new_const_iterator rightCluster = 
      std::find_if( detSet_.begin(), detSet_.end(), StripClusterAboveU( utraj)); //FIXME

    if ( rightCluster != detSet_.begin()) {
      // there are hits on the left of the utraj
      new_const_iterator leftCluster = rightCluster;
      while ( --leftCluster >=  detSet_.begin()) {
        if (isMasked(*leftCluster)) continue;
    SiStripClusterRef clusterref = edmNew::makeRefTo( handle_, leftCluster ); 
    if (accept(clusterref)){
    RecHitContainer recHits = buildRecHits(clusterref,stateOnThisDet); 
    bool isCompatible(false);
    for(RecHitContainer::const_iterator recHit=recHits.begin();recHit!=recHits.end();++recHit){   
      std::pair<bool,double> diffEst = est.estimate(stateOnThisDet, **recHit);
      if ( diffEst.first ) {
        result.push_back( TrajectoryMeasurement( stateOnThisDet, *recHit, 
                             diffEst.second));
        isCompatible = true;
      }
    }
    if(!isCompatible) break; // exit loop on first incompatible hit
    }
    else LogDebug("TkStripMeasurementDet")<<"skipping this str from last iteration on"<<geomDet().geographicalId().rawId()<<" key: "<<clusterref.key();
      }
    }
    
    for ( ; rightCluster != detSet_.end(); rightCluster++) {
      if (isMasked(*rightCluster)) continue;
      SiStripClusterRef clusterref = edmNew::makeRefTo( handle_, rightCluster ); 
      if (accept(clusterref)){
      RecHitContainer recHits = buildRecHits(clusterref,stateOnThisDet); 
      bool isCompatible(false);
      for(RecHitContainer::const_iterator recHit=recHits.begin();recHit!=recHits.end();++recHit){     
    std::pair<bool,double> diffEst = est.estimate(stateOnThisDet, **recHit);
    if ( diffEst.first ) {
      result.push_back( TrajectoryMeasurement( stateOnThisDet, *recHit, 
                           diffEst.second));
      isCompatible = true;
    }
      }
      if(!isCompatible) break; // exit loop on first incompatible hit
      }
      else LogDebug("TkStripMeasurementDet")<<"skipping this str from last iteration on"<<geomDet().geographicalId().rawId()<<" key: "<<clusterref.key();
    }
  }// end block with DetSet
  else{
    result.reserve(size());
    uint rightCluster = beginClusterI_;
    for (; rightCluster!= endClusterI_;++rightCluster){
      SiStripRegionalClusterRef clusterref = edm::makeRefToLazyGetter(regionalHandle_,rightCluster);
      if (clusterref->barycenter() > utraj) break;
    }

    uint leftCluster = 1;
    for (uint iReadBackWard=1; iReadBackWard<=(rightCluster-beginClusterI_) ; ++iReadBackWard){
        leftCluster=rightCluster-iReadBackWard;
    SiStripRegionalClusterRef clusterref = edm::makeRefToLazyGetter(regionalHandle_,leftCluster);
        if (isMasked(*clusterref)) continue;
    if (accept(clusterref)){
    RecHitContainer recHits = buildRecHits(clusterref,stateOnThisDet); 
    bool isCompatible(false);
    for(RecHitContainer::const_iterator recHit=recHits.begin();recHit!=recHits.end();++recHit){   
      std::pair<bool,double> diffEst = est.estimate(stateOnThisDet, **recHit);
      if ( diffEst.first ) {
        result.push_back( TrajectoryMeasurement( stateOnThisDet, *recHit, 
                             diffEst.second));
        isCompatible = true;
      }
    }
    if(!isCompatible) break; // exit loop on first incompatible hit
    }
    else LogDebug("TkStripMeasurementDet")<<"skipping this reg str from last iteration on"<<geomDet().geographicalId().rawId()<<" key: "<<clusterref.key();
    }
    
    for ( ; rightCluster != endClusterI_; ++rightCluster) {
      SiStripRegionalClusterRef clusterref = edm::makeRefToLazyGetter(regionalHandle_,rightCluster);
      if (isMasked(*clusterref)) continue;
      if (accept(clusterref)){
      RecHitContainer recHits = buildRecHits(clusterref,stateOnThisDet); 
      bool isCompatible(false);
      for(RecHitContainer::const_iterator recHit=recHits.begin();recHit!=recHits.end();++recHit){     
    std::pair<bool,double> diffEst = est.estimate(stateOnThisDet, **recHit);
    if ( diffEst.first ) {
      result.push_back( TrajectoryMeasurement( stateOnThisDet, *recHit, 
                           diffEst.second));
      isCompatible = true;
    }
      }
      if(!isCompatible) break; // exit loop on first incompatible hit
      }
      else LogDebug("TkStripMeasurementDet")<<"skipping this reg str from last iteration on"<<geomDet().geographicalId().rawId()<<" key: "<<clusterref.key();
    }
  }


  if ( result.empty()) {
    // create a TrajectoryMeasurement with an invalid RecHit and zero estimate
    if (stateOnThisDet.hasError()){
    uerr= sqrt(theStripGDU->specificTopology().measurementError(stateOnThisDet.localPosition(),stateOnThisDet.localError().positionError()).uu());
     if (testStrips(utraj,uerr)) {
       //LogDebug("TkStripMeasurementDet") << " DetID " << id_ << " empty after search, but active ";
       result.push_back( TrajectoryMeasurement( stateOnThisDet, InvalidTransientRecHit::build(&geomDet(), TrackingRecHit::missing), 0.F));
     } else { 
       //LogDebug("TkStripMeasurementDet") << " DetID " << id_ << " empty after search, and inactive ";
       result.push_back( TrajectoryMeasurement( stateOnThisDet, InvalidTransientRecHit::build(&geomDet(), TrackingRecHit::inactive), 0.F));
     }
    }else{
      result.push_back( TrajectoryMeasurement( stateOnThisDet, InvalidTransientRecHit::build(&geomDet(), TrackingRecHit::missing), 0.F));
    }
  }
  else {
    //LogDebug("TkStripMeasurementDet") << " DetID " << id_ << " full: " << (result.size()) << " compatible hits";
    // sort results according to estimator value
    if ( result.size() > 1) {
      sort( result.begin(), result.end(), TrajMeasLessEstim());
    }
  }
  return result;
}

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

inline

Definition at line 138 of file TkStripMeasurementDet.h.

References badStripBlocks_.

{ return badStripBlocks_; }

bool TkStripMeasurementDet::hasAllGoodChannels ( ) const

inline

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

Definition at line 115 of file TkStripMeasurementDet.h.

References badStripBlocks_, and hasAny128StripBad_.

Referenced by TkGluedMeasurementDet::fastMeasurements().

{ return !hasAny128StripBad_ && badStripBlocks_.empty(); }

bool TkStripMeasurementDet::hasBadComponents ( const TrajectoryStateOnSurface &  tsos ) const

inlinevirtual

Implements MeasurementDet.

Definition at line 69 of file TkStripMeasurementDet.h.

{return false;}

bool TkStripMeasurementDet::isActive ( ) const

inlinevirtual

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

Implements MeasurementDet.

Definition at line 66 of file TkStripMeasurementDet.h.

References activeThisEvent_, and activeThisPeriod_.

Referenced by TkGluedMeasurementDet::fastMeasurements(), fastMeasurements(), TkGluedMeasurementDet::isActive(), recHits(), and simpleRecHits().

{ return activeThisEvent_ && activeThisPeriod_; }

bool TkStripMeasurementDet::isEmpty ( void  )

inline

Definition at line 99 of file TkStripMeasurementDet.h.

References empty.

{return empty;}

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

inlineprivate

Definition at line 174 of file TkStripMeasurementDet.h.

References SiStripCluster::amplitudes(), bad128Strip_, SiStripCluster::barycenter(), and SiStripCluster::firstStrip().

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

                                                            {
      if ( bad128Strip_[cluster.firstStrip() >> 7] ) {
          if ( bad128Strip_[(cluster.firstStrip()+cluster.amplitudes().size())  >> 7] ||
               bad128Strip_[static_cast<int32_t>(cluster.barycenter()-0.499999) >> 7] ) {
              return true;
          }
      } else {
          if ( bad128Strip_[(cluster.firstStrip()+cluster.amplitudes().size())  >> 7] &&
               bad128Strip_[static_cast<int32_t>(cluster.barycenter()-0.499999) >> 7] ) {
              return true;
          }
      }
      return false;
  }

unsigned int TkStripMeasurementDet::rawId ( ) const

inline

Definition at line 101 of file TkStripMeasurementDet.h.

References id_.

{ return id_; }

TkStripMeasurementDet::RecHitContainer TkStripMeasurementDet::recHits ( const TrajectoryStateOnSurface &  ts ) const

virtual

Implements MeasurementDet.

Definition at line 244 of file TkStripMeasurementDet.cc.

References accept(), edmNew::DetSet< T >::begin(), beginClusterI_, buildRecHit(), detSet_, empty, edmNew::DetSet< T >::end(), endClusterI_, GeomDet::geographicalId(), MeasurementDet::geomDet(), handle_, isActive(), isMasked(), isRegional, edm::Ref< C, T, F >::key(), LogDebug, edmNew::makeRefTo(), edm::makeRefToLazyGetter(), DetId::rawId(), regionalHandle_, query::result, edmNew::DetSet< T >::size(), and size().

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

{
  RecHitContainer result;
  if (empty == true) return result;
  if (isActive() == false) return result; // GIO

  if(!isRegional){//old implemetation with DetSet
    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 = edmNew::makeRefTo( handle_, ci ); 
      if (accept(cluster))
    result.push_back( buildRecHit( cluster, ts));
      else LogDebug("TkStripMeasurementDet")<<"skipping this str from last iteration on"<<geomDet().geographicalId().rawId()<<" key: "<<cluster.key();
    }
  }else{
    result.reserve(size());
    for (uint ci = beginClusterI_ ; ci!= endClusterI_;++ci){
      SiStripRegionalClusterRef clusterRef = edm::makeRefToLazyGetter(regionalHandle_,ci);
      if (isMasked(*clusterRef)) continue;
      if (accept(clusterRef))
    result.push_back( buildRecHit( clusterRef, ts));
      else LogDebug("TkStripMeasurementDet")<<"skipping this reg str from last iteration on"<<geomDet().geographicalId().rawId()<<" key: "<<clusterRef.key();
      }
  }
  return result;

}

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

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

Definition at line 317 of file TkStripMeasurementDet.cc.

References bad128Strip_, hasAny128StripBad_, and i.

                                                           { 
   if (idx == -1) {
       std::fill(bad128Strip_, bad128Strip_+6, !good);
       hasAny128StripBad_ = !good;
   } else {
       bad128Strip_[idx] = !good;
       if (good == false) {
            hasAny128StripBad_ = false;
       } else { // this should not happen, as usually you turn on all fibers
                // and then turn off the bad ones, and not vice-versa,
                // so I don't care if it's not optimized
            hasAny128StripBad_ = true;
            for (int i = 0; i < (totalStrips_ >> 7); i++) {
                if (bad128Strip_[i] == false) hasAny128StripBad_ = false;
            }
       }    
   }
    
}

void TkStripMeasurementDet::setActive ( bool  active )

inline

Turn on/off the module for reconstruction, for the full run or lumi (using info from DB, usually). This also resets the 'setActiveThisEvent' to true.

Definition at line 109 of file TkStripMeasurementDet.h.

References activeThisEvent_, activeThisPeriod_, and empty.

{ activeThisPeriod_ = active; activeThisEvent_ = true; if (!active) empty = true; }

void TkStripMeasurementDet::setActiveThisEvent ( bool  active )

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

References activeThisEvent_, and empty.

Referenced by OnDemandMeasurementTracker::assign().

{ activeThisEvent_ = active;  if (!active) empty = true; }

void TkStripMeasurementDet::setBadStripCuts ( BadStripCuts  cuts )

inline

Definition at line 131 of file TkStripMeasurementDet.h.

References badStripCuts_.

{ badStripCuts_ = cuts; }

template<typename IT >

void TkStripMeasurementDet::setClusterToSkip ( IT  begin, IT  end  )

inline

Definition at line 211 of file TkStripMeasurementDet.h.

References skipClusters_.

                                           {
    skipClusters_.clear();
    skipClusters_.insert(begin,end);
  }

void TkStripMeasurementDet::setEmpty ( )

inline

Definition at line 72 of file TkStripMeasurementDet.h.

References activeThisEvent_, and empty.

Referenced by OnDemandMeasurementTracker::assign().

{empty = true; activeThisEvent_ = true; }

void TkStripMeasurementDet::setMaskBad128StripBlocks ( bool  maskThem )

inline

Definition at line 140 of file TkStripMeasurementDet.h.

References maskBad128StripBlocks_.

{ maskBad128StripBlocks_ = maskThem; }

template<typename IT >

void TkStripMeasurementDet::setRegionalClustersToSkip ( IT  begin, IT  end  )

inline

Definition at line 216 of file TkStripMeasurementDet.h.

References skipRegClusters_.

Referenced by OnDemandMeasurementTracker::assign(), and MeasurementTracker::updateStrips().

                                                    {
    skipRegClusters_.clear();
    skipRegClusters_.insert(begin,end);
  }

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

Definition at line 289 of file TkStripMeasurementDet.cc.

References accept(), edmNew::DetSet< T >::begin(), beginClusterI_, buildSimpleRecHit(), detSet_, empty, edmNew::DetSet< T >::end(), endClusterI_, GeomDet::geographicalId(), MeasurementDet::geomDet(), handle_, isActive(), isMasked(), isRegional, edm::Ref< C, T, F >::key(), LogDebug, edmNew::makeRefTo(), edm::makeRefToLazyGetter(), DetId::rawId(), regionalHandle_, edmNew::DetSet< T >::size(), and size().

Referenced by TkGluedMeasurementDet::collectRecHits().

{
  if (empty || !isActive()) return;

  if(!isRegional){//old implemetation with DetSet
    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 = edmNew::makeRefTo( handle_, ci ); 
      if (accept(cluster))
    buildSimpleRecHit( cluster, ts,result);
      else LogDebug("TkStripMeasurementDet")<<"skipping this str from last iteration on"<<geomDet().geographicalId().rawId()<<" key: "<<cluster.key();
    }
  }else{
    result.reserve(size());
    for (uint ci = beginClusterI_ ; ci!= endClusterI_;++ci){
      SiStripRegionalClusterRef clusterRef = edm::makeRefToLazyGetter(regionalHandle_,ci);
      if (isMasked(*clusterRef)) continue;
      if (accept(clusterRef))
    buildSimpleRecHit( clusterRef, ts,result);
      else LogDebug("TkStripMeasurementDet")<<"skipping this reg str from last iteration on"<<geomDet().geographicalId().rawId()<<" key: "<<clusterRef.key();
    }
  }
}

int TkStripMeasurementDet::size ( void  ) const

inline

Definition at line 105 of file TkStripMeasurementDet.h.

References beginClusterI_, and endClusterI_.

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

{return endClusterI_ - beginClusterI_ ; }

const StripGeomDetUnit& TkStripMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 83 of file TkStripMeasurementDet.h.

References theStripGDU.

Referenced by buildRecHit(), buildRecHits(), buildSimpleRecHit(), TkGluedMeasurementDet::testStrips(), and TkStripMeasurementDet().

{return *theStripGDU;}

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

References cond::ecalcond::bad(), badStripBlocks_, badStripCuts_, end, max(), TkStripMeasurementDet::BadStripCuts::maxBad, TkStripMeasurementDet::BadStripCuts::maxConsecutiveBad, min, convertSQLiteXML::ok, and totalStrips_.

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

                                                               {
    int16_t start = (int16_t) std::max<float>(utraj - 3*uerr, 0);
    int16_t end   = (int16_t) std::min<float>(utraj + 3*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;
    BSBIT bsbc = badStripBlocks_.begin(), bsbe = badStripBlocks_.end();

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

inline

Definition at line 104 of file TkStripMeasurementDet.h.

References detSet_.

{return detSet_;}

void TkStripMeasurementDet::unset ( )

inline

Definition at line 206 of file TkStripMeasurementDet.h.

References skipClusters_, and skipRegClusters_.

Referenced by OnDemandMeasurementTracker::assign().

              {
    skipClusters_.clear();
    skipRegClusters_.clear();
  }

void TkStripMeasurementDet::update ( const detset &  detSet, const edm::Handle< edmNew::DetSetVector< SiStripCluster > >  h, unsigned int  id  )

inline

Definition at line 43 of file TkStripMeasurementDet.h.

References detSet_, empty, h, handle_, id_, and isRegional.

Referenced by OnDemandMeasurementTracker::assign(), python.Vispa.Gui.VispaWidget.VispaWidget::autosize(), python.Vispa.Views.LineDecayView.LineDecayContainer::createObject(), python.Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), python.Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), python.Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), python.Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), python.Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), python.Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), python.Vispa.Gui.FindDialog.FindDialog::reset(), python.Vispa.Gui.PortConnection.PointToPointConnection::select(), python.Vispa.Gui.VispaWidget.VispaWidget::select(), python.Vispa.Views.LineDecayView.LineDecayContainer::select(), python.Vispa.Gui.VispaWidget.VispaWidget::setText(), python.Vispa.Gui.VispaWidget.VispaWidget::setTitle(), python.Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), python.Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), python.Vispa.Gui.PortConnection.PointToPointConnection::updateConnection(), and MeasurementTracker::updateStrips().

                             { 
    detSet_ = detSet; 
    handle_ = h;
    id_ = id;
    empty = false;
    isRegional = false;
  }

void TkStripMeasurementDet::update ( std::vector< SiStripCluster >::const_iterator  begin, std::vector< SiStripCluster >::const_iterator  end, const edm::Handle< edm::LazyGetter< SiStripCluster > >  h, unsigned int  id  )

inline

Definition at line 53 of file TkStripMeasurementDet.h.

References activeThisEvent_, beginClusterI_, empty, endClusterI_, h, id_, isRegional, and regionalHandle_.

Referenced by python.Vispa.Gui.VispaWidget.VispaWidget::autosize(), python.Vispa.Views.LineDecayView.LineDecayContainer::createObject(), python.Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), python.Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), python.Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), python.Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), python.Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), python.Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), python.Vispa.Gui.FindDialog.FindDialog::reset(), python.Vispa.Gui.PortConnection.PointToPointConnection::select(), python.Vispa.Gui.VispaWidget.VispaWidget::select(), python.Vispa.Views.LineDecayView.LineDecayContainer::select(), python.Vispa.Gui.VispaWidget.VispaWidget::setText(), python.Vispa.Gui.VispaWidget.VispaWidget::setTitle(), python.Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), python.Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and python.Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                             { 
    regionalHandle_ = h;
    beginClusterI_ = begin - regionalHandle_->begin_record();
    endClusterI_ = end - regionalHandle_->begin_record();
    id_ = id;
    empty = false;
    activeThisEvent_ = true;
    isRegional = true;
  }

Member Data Documentation

bool TkStripMeasurementDet::activeThisEvent_

private

Definition at line 148 of file TkStripMeasurementDet.h.

Referenced by isActive(), setActive(), setActiveThisEvent(), setEmpty(), and update().

bool TkStripMeasurementDet::activeThisPeriod_

private

Definition at line 148 of file TkStripMeasurementDet.h.

Referenced by isActive(), and setActive().

bool TkStripMeasurementDet::bad128Strip_[6]

private

Definition at line 161 of file TkStripMeasurementDet.h.

Referenced by isMasked(), and set128StripStatus().

std::vector<BadStripBlock> TkStripMeasurementDet::badStripBlocks_

private

Definition at line 164 of file TkStripMeasurementDet.h.

Referenced by getBadStripBlocks(), hasAllGoodChannels(), and testStrips().

BadStripCuts TkStripMeasurementDet::badStripCuts_

private

Definition at line 166 of file TkStripMeasurementDet.h.

Referenced by setBadStripCuts(), and testStrips().

uint TkStripMeasurementDet::beginClusterI_

private

Definition at line 171 of file TkStripMeasurementDet.h.

Referenced by fastMeasurements(), recHits(), simpleRecHits(), size(), and update().

detset TkStripMeasurementDet::detSet_

private

Definition at line 152 of file TkStripMeasurementDet.h.

Referenced by fastMeasurements(), recHits(), simpleRecHits(), theSet(), and update().

bool TkStripMeasurementDet::empty

private

Definition at line 146 of file TkStripMeasurementDet.h.

Referenced by fastMeasurements(), isEmpty(), recHits(), setActive(), setActiveThisEvent(), python.Vispa.Gui.VispaWidget.TextField::setAutosizeFont(), python.Vispa.Gui.VispaWidget.TextField::setAutotruncate(), setEmpty(), simpleRecHits(), and update().

uint TkStripMeasurementDet::endClusterI_

private

Definition at line 172 of file TkStripMeasurementDet.h.

Referenced by fastMeasurements(), recHits(), simpleRecHits(), size(), and update().

edm::Handle<edmNew::DetSetVector<SiStripCluster> > TkStripMeasurementDet::handle_

private

Definition at line 153 of file TkStripMeasurementDet.h.

Referenced by fastMeasurements(), recHits(), simpleRecHits(), and update().

bool TkStripMeasurementDet::hasAny128StripBad_

private

Definition at line 162 of file TkStripMeasurementDet.h.

Referenced by hasAllGoodChannels(), and set128StripStatus().

unsigned int TkStripMeasurementDet::id_

private

Definition at line 150 of file TkStripMeasurementDet.h.

Referenced by rawId(), TkStripMeasurementDet(), and update().

bool TkStripMeasurementDet::isRegional

private

Definition at line 144 of file TkStripMeasurementDet.h.

Referenced by fastMeasurements(), recHits(), simpleRecHits(), and update().

bool TkStripMeasurementDet::maskBad128StripBlocks_

private

Definition at line 162 of file TkStripMeasurementDet.h.

Referenced by setMaskBad128StripBlocks().

edm::Handle<edm::LazyGetter<SiStripCluster> > TkStripMeasurementDet::regionalHandle_

private

Definition at line 170 of file TkStripMeasurementDet.h.

Referenced by fastMeasurements(), recHits(), simpleRecHits(), and update().

std::set<SiStripClusterRef> TkStripMeasurementDet::skipClusters_

private

Definition at line 195 of file TkStripMeasurementDet.h.

Referenced by accept(), setClusterToSkip(), and unset().

std::set<SiStripRegionalClusterRef> TkStripMeasurementDet::skipRegClusters_

private

Definition at line 196 of file TkStripMeasurementDet.h.

Referenced by accept(), setRegionalClustersToSkip(), and unset().

const StripClusterParameterEstimator* TkStripMeasurementDet::theCPE

private

Definition at line 157 of file TkStripMeasurementDet.h.

Referenced by buildRecHit(), buildRecHits(), and buildSimpleRecHit().

const StripGeomDetUnit* TkStripMeasurementDet::theStripGDU

private

Definition at line 156 of file TkStripMeasurementDet.h.

Referenced by fastMeasurements(), specificGeomDet(), and TkStripMeasurementDet().

int TkStripMeasurementDet::totalStrips_

private

Definition at line 165 of file TkStripMeasurementDet.h.

Referenced by testStrips(), and TkStripMeasurementDet().