TkPixelMeasurementDet Class Reference

#include <TkPixelMeasurementDet.h>

Inheritance diagram for TkPixelMeasurementDet:

Public Types
typedef detset::const_iterator	const_iterator
typedef edmNew::DetSet< SiPixelCluster >	detset
typedef PixelClusterParameterEstimator::LocalValues	LocalValues
typedef edm::Ref< edmNew::DetSetVector< SiPixelCluster >, SiPixelCluster >	SiPixelClusterRef
Public Types inherited from MeasurementDet
typedef TrackingRecHit::ConstRecHitContainer	RecHitContainer
using	SimpleHitContainer = std::vector< BaseTrackerRecHit * >
typedef tracking::TempMeasurements	TempMeasurements

Public Member Functions
bool	accept (SiPixelClusterRefNew &r, const std::vector< bool > skipClusters) const
TrackingRecHit::RecHitPointer	buildRecHit (const SiPixelClusterRef &cluster, const LocalTrajectoryParameters &ltp) const
void	clearBadRocPositions ()
	Clear the list of bad ROCs. More...
bool	hasBadComponents (const TrajectoryStateOnSurface &tsos, const MeasurementTrackerEvent &dat) const
int	index () const
bool	isActive (const MeasurementTrackerEvent &data) const
	Is this module active in reconstruction? It must be both 'setActiveThisEvent' and 'setActive'. More...
bool	isEmpty (const PxMeasurementDetSet &data) const
virtual bool	measurements (const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &est, const MeasurementTrackerEvent &dat, TempMeasurements &result) const
virtual RecHitContainer	recHits (const TrajectoryStateOnSurface &, const MeasurementTrackerEvent &dat) const
virtual bool	recHits (SimpleHitContainer &result, const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &, const MeasurementTrackerEvent &data) const
void	setActive (bool active)
	Turn on/off the module for reconstruction, for the full run or lumi (using info from DB, usually). More...
void	setActiveThisEvent (PxMeasurementDetSet &data, 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	setBadRocPositions (std::vector< LocalPoint > &positions)
	Sets the list of bad ROCs, identified by the positions of their centers in the local coordinate frame. More...
void	setEmpty (PxMeasurementDetSet &data)
void	setIndex (int i)
const PixelGeomDetUnit &	specificGeomDet () const
	TkPixelMeasurementDet (const GeomDet *gdet, PxMeasurementConditionSet &conditionSet)
void	update (PxMeasurementDetSet &data, const detset &detSet)
virtual	~TkPixelMeasurementDet ()
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
virtual bool	recHits (const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &, const MeasurementTrackerEvent &data, RecHitContainer &result, std::vector< float > &) const
const Surface &	surface () const

Private Member Functions
PxMeasurementConditionSet &	conditionSet ()
const PxMeasurementConditionSet &	conditionSet () const
const PixelClusterParameterEstimator *	cpe () const

Private Attributes
std::vector< LocalPoint >	badRocPositions_
unsigned int	id_
int	index_
PxMeasurementConditionSet *	theDetConditions

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

Detailed Description

Definition at line 17 of file TkPixelMeasurementDet.h.

Member Typedef Documentation

typedef detset::const_iterator TkPixelMeasurementDet::const_iterator

Definition at line 23 of file TkPixelMeasurementDet.h.

typedef edmNew::DetSet<SiPixelCluster> TkPixelMeasurementDet::detset

Definition at line 22 of file TkPixelMeasurementDet.h.

typedef PixelClusterParameterEstimator::LocalValues TkPixelMeasurementDet::LocalValues

Definition at line 24 of file TkPixelMeasurementDet.h.

typedef edm::Ref<edmNew::DetSetVector<SiPixelCluster>, SiPixelCluster> TkPixelMeasurementDet::SiPixelClusterRef

Definition at line 20 of file TkPixelMeasurementDet.h.

Constructor & Destructor Documentation

TkPixelMeasurementDet::TkPixelMeasurementDet	(	const GeomDet *	gdet,
		PxMeasurementConditionSet &	conditionSet
	)

Definition at line 17 of file TkPixelMeasurementDet.cc.

                                                                   : 
    MeasurementDet (gdet),
    theDetConditions(&conditions)
  {
    if ( dynamic_cast<const PixelGeomDetUnit*>(gdet) == 0) {
      throw MeasurementDetException( "TkPixelMeasurementDet constructed with a GeomDet which is not a PixelGeomDetUnit");
    }
  }

virtual TkPixelMeasurementDet::~TkPixelMeasurementDet ( )

inlinevirtual

Definition at line 37 of file TkPixelMeasurementDet.h.

References MeasurementDet::recHits().

{ }

Member Function Documentation

bool TkPixelMeasurementDet::accept ( SiPixelClusterRefNew &  r, const std::vector< bool >  skipClusters  ) const

inline

Definition at line 91 of file TkPixelMeasurementDet.h.

References edm::Ref< C, T, F >::key().

Referenced by Vispa.Gui.BoxContentDialog.BoxContentDialog::apply(), Vispa.Plugins.ConfigEditor.ToolDialog.ToolDialog::apply(), and esMonitoring.FDJsonServer::handle_accept().

                                                                                         {
    
    if(skipClusters.empty()) return true;
    if (r.key()>=skipClusters.size()){
      edm::LogError("IndexMisMatch")<<r.key()<<" is larger than: "<<skipClusters.size()<<" no skipping done";
      return true;
    }
    return not skipClusters[r.key()];
  }

TrackingRecHit::RecHitPointer TkPixelMeasurementDet::buildRecHit	(	const SiPixelClusterRef &	cluster,
		const LocalTrajectoryParameters &	ltp
	)		const

Definition at line 55 of file TkPixelMeasurementDet.cc.

References cpe(), MeasurementDet::fastGeomDet(), PixelClusterParameterEstimator::getParameters(), and specificGeomDet().

Referenced by recHits().

{
  const GeomDetUnit& gdu( specificGeomDet());

  auto && params = cpe()->getParameters( * cluster, gdu, ltp );
  return std::make_shared<SiPixelRecHit>( std::get<0>(params), std::get<1>(params), std::get<2>(params), fastGeomDet(), cluster);
}

void TkPixelMeasurementDet::clearBadRocPositions ( )

inline

Clear the list of bad ROCs.

Definition at line 73 of file TkPixelMeasurementDet.h.

{ badRocPositions_.clear(); }

PxMeasurementConditionSet& TkPixelMeasurementDet::conditionSet ( )

inlineprivate

Definition at line 84 of file TkPixelMeasurementDet.h.

{ return *theDetConditions; }

const PxMeasurementConditionSet& TkPixelMeasurementDet::conditionSet ( ) const

inlineprivate

Definition at line 85 of file TkPixelMeasurementDet.h.

{ return *theDetConditions; }

const PixelClusterParameterEstimator* TkPixelMeasurementDet::cpe ( ) const

inlineprivate

Definition at line 87 of file TkPixelMeasurementDet.h.

Referenced by buildRecHit().

{ return conditionSet().pixelCPE(); }

bool TkPixelMeasurementDet::hasBadComponents ( const TrajectoryStateOnSurface &  tsos, const MeasurementTrackerEvent &  dat  ) const

virtual

Implements MeasurementDet.

Definition at line 98 of file TkPixelMeasurementDet.cc.

References funct::abs(), badRocPositions_, TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), LocalTrajectoryError::positionError(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by measurements().

                                                                                                                          {
    if (badRocPositions_.empty()) return false;
    LocalPoint lp = tsos.localPosition();
    LocalError le = tsos.localError().positionError();
    double dx = 3*std::sqrt(le.xx()) + theRocWidth, dy = 3*std::sqrt(le.yy()) + theRocHeight;
    for (std::vector<LocalPoint>::const_iterator it = badRocPositions_.begin(), ed = badRocPositions_.end(); it != ed; ++it) {
        if ( (std::abs(it->x() - lp.x()) < dx) &&
             (std::abs(it->y() - lp.y()) < dy) ) return true;
    } 
    return false;
}

int TkPixelMeasurementDet::index ( ) const

inline

Definition at line 75 of file TkPixelMeasurementDet.h.

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

{ return index_; }

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

inlinevirtual

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

Implements MeasurementDet.

Definition at line 66 of file TkPixelMeasurementDet.h.

References MeasurementDet::hasBadComponents(), diffTreeTool::index, PxMeasurementDetSet::isActive(), and MeasurementTrackerEvent::pixelData().

Referenced by measurements(), and recHits().

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

bool TkPixelMeasurementDet::isEmpty ( const PxMeasurementDetSet &  data ) const

inline

Definition at line 35 of file TkPixelMeasurementDet.h.

References PxMeasurementDetSet::empty(), and diffTreeTool::index.

Referenced by plotting.Plot::clone(), and recHits().

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

bool TkPixelMeasurementDet::measurements ( const TrajectoryStateOnSurface &  stateOnThisDet, const MeasurementEstimator &  est, const MeasurementTrackerEvent &  dat, TempMeasurements &  result  ) const

virtual

Implements MeasurementDet.

Definition at line 27 of file TkPixelMeasurementDet.cc.

References tracking::TempMeasurements::add(), MeasurementEstimator::estimate(), F(), hasBadComponents(), isActive(), eostools::move(), recHits(), tracking::TempMeasurements::size(), MeasurementDet::theInactiveHit, and MeasurementDet::theMissingHit.

                                                       {

  if (!isActive(data)) {
    result.add(theInactiveHit, 0.F);
    return true;
  }
  
  auto oldSize = result.size();
  MeasurementDet::RecHitContainer && allHits = recHits(stateOnThisDet, data);
  for (auto && hit : allHits) {
    std::pair<bool,double> diffEst = est.estimate( stateOnThisDet, *hit);
    if ( diffEst.first)
      result.add(std::move(hit), diffEst.second);
  }

  if (result.size()>oldSize) return true;

  // create a TrajectoryMeasurement with an invalid RecHit and zero estimate
  bool inac = hasBadComponents(stateOnThisDet, data);
  result.add(inac ? theInactiveHit : theMissingHit, 0.F);
  return inac;

}

TkPixelMeasurementDet::RecHitContainer TkPixelMeasurementDet::recHits ( const TrajectoryStateOnSurface &  ts, const MeasurementTrackerEvent &  dat  ) const

virtual

Implements MeasurementDet.

Definition at line 65 of file TkPixelMeasurementDet.cc.

References begin, edmNew::DetSet< T >::begin(), buildRecHit(), edmNew::DetSetVector< T >::data(), PxMeasurementDetSet::detSet(), edmNew::DetSet< T >::end(), MeasurementDet::fastGeomDet(), GeomDet::geographicalId(), PxMeasurementDetSet::handle(), index(), isActive(), isEmpty(), TrajectoryStateOnSurface::localParameters(), LogDebug, edmNew::DetSet< T >::makeRefTo(), or, MeasurementTrackerEvent::pixelClustersToSkip(), MeasurementTrackerEvent::pixelData(), DetId::rawId(), mps_fire::result, edmNew::DetSet< T >::size(), and SiPixelCluster::size().

Referenced by measurements().

{
  RecHitContainer result;
  if (isEmpty(data.pixelData())== true ) return result;
  if (isActive(data) == false) return result;
  const SiPixelCluster* begin=0;
  if (0 != data.pixelData().handle()->data().size()) {
     begin = &(data.pixelData().handle()->data().front());
  }
  const detset & detSet = data.pixelData().detSet(index());
  result.reserve(detSet.size());
  for ( const_iterator ci = detSet.begin(); ci != detSet.end(); ++ ci ) {
    
    if (ci < begin){
      edm::LogError("IndexMisMatch")<<"TkPixelMeasurementDet cannot create hit because of index mismatch.";
      return result;
    }
     unsigned int index = ci-begin;
     if (!data.pixelClustersToSkip().empty() &&  index>=data.pixelClustersToSkip().size()){
       edm::LogError("IndexMisMatch")<<"TkPixelMeasurementDet cannot create hit because of index mismatch. i.e "<<index<<" >= "<<data.pixelClustersToSkip().size();
       return result;
     }
     if(data.pixelClustersToSkip().empty() or (not data.pixelClustersToSkip()[index]) ) {
       SiPixelClusterRef cluster = detSet.makeRefTo( data.pixelData().handle(), ci );
       result.push_back( buildRecHit( cluster, ts.localParameters() ) );
     }else{   
       LogDebug("TkPixelMeasurementDet")<<"skipping this cluster from last iteration on "<<fastGeomDet().geographicalId().rawId()<<" key: "<<index;
     }
  }
  return result;
}

virtual bool TkPixelMeasurementDet::recHits ( SimpleHitContainer &  result, const TrajectoryStateOnSurface &  stateOnThisDet, const MeasurementEstimator &  , const MeasurementTrackerEvent &  data  ) const

inlinevirtual

Reimplemented from MeasurementDet.

Definition at line 42 of file TkPixelMeasurementDet.h.

References MeasurementDet::measurements(), and mps_fire::result.

                                                                                                                                        {
    assert("not implemented for Pixel yet"==nullptr);
  }

void TkPixelMeasurementDet::setActive ( bool  active )

inline

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

Definition at line 61 of file TkPixelMeasurementDet.h.

References diffTreeTool::index.

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

void TkPixelMeasurementDet::setActiveThisEvent ( PxMeasurementDetSet &  data, 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 64 of file TkPixelMeasurementDet.h.

References diffTreeTool::index, and PxMeasurementDetSet::setActiveThisEvent().

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

void TkPixelMeasurementDet::setBadRocPositions ( std::vector< LocalPoint > &  positions )

inline

Sets the list of bad ROCs, identified by the positions of their centers in the local coordinate frame.

Definition at line 71 of file TkPixelMeasurementDet.h.

{ badRocPositions_.swap(positions); }

void TkPixelMeasurementDet::setEmpty ( PxMeasurementDetSet &  data )

inline

Definition at line 34 of file TkPixelMeasurementDet.h.

References diffTreeTool::index, and PxMeasurementDetSet::setEmpty().

{ data.setEmpty(index());  }

void TkPixelMeasurementDet::setIndex ( int  i )

inline

Definition at line 76 of file TkPixelMeasurementDet.h.

References i.

{ index_ = i; }

const PixelGeomDetUnit& TkPixelMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 54 of file TkPixelMeasurementDet.h.

References MeasurementDet::fastGeomDet().

Referenced by buildRecHit().

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

void TkPixelMeasurementDet::update ( PxMeasurementDetSet &  data, const detset &  detSet  )

inline

Definition at line 29 of file TkPixelMeasurementDet.h.

References diffTreeTool::index, PxMeasurementDetSet::setActiveThisEvent(), and PxMeasurementDetSet::update().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                                 { 
    data.update(index(), detSet);
    data.setActiveThisEvent(index(), true);
  }

Member Data Documentation

std::vector< LocalPoint > TkPixelMeasurementDet::badRocPositions_

private

Definition at line 80 of file TkPixelMeasurementDet.h.

Referenced by hasBadComponents().

unsigned int TkPixelMeasurementDet::id_

private

Definition at line 79 of file TkPixelMeasurementDet.h.

int TkPixelMeasurementDet::index_

private

Definition at line 82 of file TkPixelMeasurementDet.h.

PxMeasurementConditionSet* TkPixelMeasurementDet::theDetConditions

private

Definition at line 83 of file TkPixelMeasurementDet.h.