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

RecHitContainer	compHits (const TrajectoryStateOnSurface &, const MeasurementTrackerEvent &dat, float xl, float yl) const

const PxMeasurementDetSet::BadFEDChannelPositions *	getBadFEDChannelPositions (const MeasurementTrackerEvent &data) const

bool	hasBadComponents (const TrajectoryStateOnSurface &tsos, const MeasurementTrackerEvent &dat) const override

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 PxMeasurementDetSet &data) const

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

RecHitContainer	recHits (const TrajectoryStateOnSurface &, const MeasurementTrackerEvent &dat) const override

bool	recHits (SimpleHitContainer &result, const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &, const MeasurementTrackerEvent &data) const override

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)

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

virtual bool	recHits (const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &, const MeasurementTrackerEvent &data, RecHitContainer &result, std::vector< float > &) const

const Surface &	surface () const

virtual	~MeasurementDet ()=default

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

◆ const_iterator

typedef detset::const_iterator TkPixelMeasurementDet::const_iterator

Definition at line 22 of file TkPixelMeasurementDet.h.

◆ detset

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

Definition at line 21 of file TkPixelMeasurementDet.h.

◆ LocalValues

typedef PixelClusterParameterEstimator::LocalValues TkPixelMeasurementDet::LocalValues

Definition at line 23 of file TkPixelMeasurementDet.h.

◆ SiPixelClusterRef

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

Definition at line 19 of file TkPixelMeasurementDet.h.

Constructor & Destructor Documentation

◆ TkPixelMeasurementDet()

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) == nullptr) {
     throw MeasurementDetException("TkPixelMeasurementDet constructed with a GeomDet which is not a PixelGeomDetUnit");
   }
 }

◆ ~TkPixelMeasurementDet()

TkPixelMeasurementDet::~TkPixelMeasurementDet ( )

inlineoverride

Definition at line 35 of file TkPixelMeasurementDet.h.

35 {}

Member Function Documentation

◆ accept()

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

inline

Definition at line 99 of file TkPixelMeasurementDet.h.

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

References alignCSCRings::r, and HLTSiStripMonitoring_cff::skipClusters.

Referenced by esMonitoring.FDJsonServer::handle_accept().

◆ buildRecHit()

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

Definition at line 69 of file TkPixelMeasurementDet.cc.

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

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

Referenced by compHits().

◆ clearBadRocPositions()

void TkPixelMeasurementDet::clearBadRocPositions ( )

inline

Clear the list of bad ROCs.

Definition at line 77 of file TkPixelMeasurementDet.h.

77 { badRocPositions_.clear(); }

◆ compHits()

TkPixelMeasurementDet::RecHitContainer TkPixelMeasurementDet::compHits	(	const TrajectoryStateOnSurface &	ts,
		const MeasurementTrackerEvent &	dat,
		float	xl,
		float	yl
	)		const

Definition at line 85 of file TkPixelMeasurementDet.cc.

                                                                                        {
   RecHitContainer result;
   if (isEmpty(data.pixelData()) == true)
     return result;
   if (isActive(data) == false)
     return result;
   const SiPixelCluster* begin = nullptr;
   if (!data.pixelData().handle()->data().empty()) {
     begin = &(data.pixelData().handle()->data().front());
   }
   const detset& detSet = data.pixelData().detSet(index());
   result.reserve(detSet.size());
  
   // pixel topology is rectangular, all positions are independent
   LocalVector maxD(xl, yl, 0);
   auto PMinus = specificGeomDet().specificTopology().measurementPosition(ts.localPosition() - maxD);
   auto PPlus = specificGeomDet().specificTopology().measurementPosition(ts.localPosition() + maxD);
  
   int xminus = PMinus.x();
   int yminus = PMinus.y();
   int xplus = PPlus.x() + 0.5f;
   int yplus = PPlus.y() + 0.5f;
  
   // rechits are sorted in x...
   auto rightCluster = std::find_if(
       detSet.begin(), detSet.end(), [xplus](const SiPixelCluster& cl) { return cl.minPixelRow() > xplus; });
  
   // std::cout << "px xlim " << xl << ' ' << xminus << '/' << xplus << ' ' << rightCluster-detSet.begin() << ',' << detSet.end()-rightCluster << std::endl;
  
   // consider only compatible clusters
   for (auto ci = detSet.begin(); ci != rightCluster; ++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 (ci->maxPixelRow() < xminus)
       continue;
     // also check compatibility in y... (does not add much)
     if (ci->minPixelCol() > yplus)
       continue;
     if (ci->maxPixelCol() < yminus)
       continue;
  
     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;
 }

References edmNew::DetSet< T >::begin(), buildRecHit(), GetRecoTauVFromDQM_MC_cff::cl, data, edmNew::DetSet< T >::end(), MeasurementDet::fastGeomDet(), GeomDet::geographicalId(), index(), isActive(), isEmpty(), TrajectoryStateOnSurface::localParameters(), TrajectoryStateOnSurface::localPosition(), LogDebug, edmNew::DetSet< T >::makeRefTo(), Topology::measurementPosition(), or, DetId::rawId(), mps_fire::result, edmNew::DetSet< T >::size(), specificGeomDet(), PixelGeomDetUnit::specificTopology(), and PV2DBase< T, PVType, FrameType >::x().

Referenced by measurements(), and recHits().

◆ conditionSet() [1/2]

PxMeasurementConditionSet& TkPixelMeasurementDet::conditionSet ( )

inlineprivate

Definition at line 93 of file TkPixelMeasurementDet.h.

93 { return *theDetConditions; }

◆ conditionSet() [2/2]

const PxMeasurementConditionSet& TkPixelMeasurementDet::conditionSet ( ) const

inlineprivate

Definition at line 94 of file TkPixelMeasurementDet.h.

94 { return *theDetConditions; }

◆ cpe()

const PixelClusterParameterEstimator* TkPixelMeasurementDet::cpe ( ) const

inlineprivate

Definition at line 96 of file TkPixelMeasurementDet.h.

96 { return conditionSet().pixelCPE(); }

Referenced by buildRecHit().

◆ getBadFEDChannelPositions()

const PxMeasurementDetSet::BadFEDChannelPositions* TkPixelMeasurementDet::getBadFEDChannelPositions ( const MeasurementTrackerEvent & data ) const

inline

Definition at line 79 of file TkPixelMeasurementDet.h.

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

References data.

Referenced by hasBadComponents().

◆ hasBadComponents()

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

overridevirtual

Implements MeasurementDet.

Definition at line 149 of file TkPixelMeasurementDet.cc.

                                                                                         {
   auto badFEDChannelPositions = getBadFEDChannelPositions(data);
   if (badRocPositions_.empty() && badFEDChannelPositions == nullptr)
     return false;
  
   auto lp = tsos.localPosition();
   auto le = tsos.localError().positionError();
   for (auto const& broc : badRocPositions_) {
     auto dx = std::abs(broc.x() - lp.x()) - theRocWidth;
     auto dy = std::abs(broc.y() - lp.y()) - theRocHeight;
     if ((dx <= 0.f) & (dy <= 0.f))
       return true;
     if ((dx * dx < 9.f * le.xx()) && (dy * dy < 9.f * le.yy()))
       return true;
   }
  
   if (badFEDChannelPositions == nullptr)
     return false;
   float dx = 3.f * std::sqrt(le.xx()) + theRocWidth, dy = 3.f * std::sqrt(le.yy()) + theRocHeight;
   for (auto const& p : *badFEDChannelPositions) {
     if (lp.x() > (p.first.x() - dx) && lp.x() < (p.second.x() + dx) && lp.y() > (p.first.y() - dy) &&
         lp.y() < (p.second.y() + dy)) {
       return true;
     }
   }
  
   return false;
 }

References funct::abs(), badRocPositions_, data, PVValHelper::dx, PVValHelper::dy, f, getBadFEDChannelPositions(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), AlCaHLTBitMon_ParallelJobs::p, LocalTrajectoryError::positionError(), and mathSSE::sqrt().

Referenced by measurements().

◆ index()

int TkPixelMeasurementDet::index ( ) const

inline

Definition at line 84 of file TkPixelMeasurementDet.h.

84 { return index_; }

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

◆ isActive()

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

inlineoverridevirtual

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

Implements MeasurementDet.

Definition at line 70 of file TkPixelMeasurementDet.h.

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

References data.

Referenced by compHits(), and measurements().

◆ isEmpty()

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

inline

Definition at line 33 of file TkPixelMeasurementDet.h.

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

References data.

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

◆ measurements()

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

overridevirtual

Implements MeasurementDet.

Definition at line 24 of file TkPixelMeasurementDet.cc.

                                                                          {
   if (!isActive(data)) {
     result.add(theInactiveHit, 0.F);
     return true;
   }
  
   auto xl = 100.f;
   auto yl = 100.f;
   // do not apply for iteration not cutting on propagation
   if (est.maxSagitta() >= 0) {
     // do not use this as it does not account for APE...
     // auto xyLimits = est.maximalLocalDisplacement(stateOnThisDet,fastGeomDet().specificSurface());
     auto le = stateOnThisDet.localError().positionError();
     LocalError lape = static_cast<TrackerGeomDet const&>(fastGeomDet()).localAlignmentError();
     xl = le.xx();
     yl = le.yy();
     if (lape.valid()) {
       xl += lape.xx();
       yl += lape.yy();
     }
     // 5 sigma to be on the safe side
     xl = 5.f * std::sqrt(xl);
     yl = 5.f * std::sqrt(yl);
   }
  
   auto oldSize = result.size();
   MeasurementDet::RecHitContainer&& allHits = compHits(stateOnThisDet, data, xl, yl);
   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;
 }

References compHits(), data, MeasurementEstimator::estimate(), F(), MeasurementDet::fastGeomDet(), hasBadComponents(), isActive(), TrajectoryStateOnSurface::localError(), MeasurementEstimator::maxSagitta(), eostools::move(), LocalTrajectoryError::positionError(), mps_fire::result, mathSSE::sqrt(), MeasurementDet::theInactiveHit, MeasurementDet::theMissingHit, LocalError::valid(), LocalError::xx(), and LocalError::yy().

◆ recHits() [1/2]

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

overridevirtual

Implements MeasurementDet.

Definition at line 78 of file TkPixelMeasurementDet.cc.

                                                                                                                  {
   float xl = 100.f;  // larger than any detector
   float yl = 100.f;
   return compHits(ts, data, xl, yl);
 }

References compHits(), and data.

◆ recHits() [2/2]

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

inlineoverridevirtual

Reimplemented from MeasurementDet.

Definition at line 47 of file TkPixelMeasurementDet.h.

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

References cms::cuda::assert().

◆ setActive()

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 65 of file TkPixelMeasurementDet.h.

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

◆ setActiveThisEvent()

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 68 of file TkPixelMeasurementDet.h.

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

References data.

◆ setBadRocPositions()

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 75 of file TkPixelMeasurementDet.h.

75 { badRocPositions_.swap(positions); }

◆ setEmpty()

void TkPixelMeasurementDet::setEmpty ( PxMeasurementDetSet & data )

inline

Definition at line 32 of file TkPixelMeasurementDet.h.

32 { data.setEmpty(index()); }

References data.

◆ setIndex()

void TkPixelMeasurementDet::setIndex ( int i )

inline

Definition at line 85 of file TkPixelMeasurementDet.h.

85 { index_ = i; }

References mps_fire::i.

◆ specificGeomDet()

const PixelGeomDetUnit& TkPixelMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 59 of file TkPixelMeasurementDet.h.

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

References MeasurementDet::fastGeomDet().

Referenced by buildRecHit(), and compHits().

◆ update()

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

inline

Definition at line 27 of file TkPixelMeasurementDet.h.

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

References data.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

Member Data Documentation

◆ badRocPositions_

std::vector<LocalPoint> TkPixelMeasurementDet::badRocPositions_

private

Definition at line 89 of file TkPixelMeasurementDet.h.

Referenced by hasBadComponents().

◆ id_

unsigned int TkPixelMeasurementDet::id_

private

Definition at line 88 of file TkPixelMeasurementDet.h.

◆ index_

int TkPixelMeasurementDet::index_

private

Definition at line 91 of file TkPixelMeasurementDet.h.

◆ theDetConditions

PxMeasurementConditionSet* TkPixelMeasurementDet::theDetConditions

private

Definition at line 92 of file TkPixelMeasurementDet.h.

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ const_iterator

◆ detset

◆ LocalValues

◆ SiPixelClusterRef

Constructor & Destructor Documentation

◆ TkPixelMeasurementDet()

◆ ~TkPixelMeasurementDet()

Member Function Documentation

◆ accept()

◆ buildRecHit()

◆ clearBadRocPositions()

◆ compHits()

◆ conditionSet() [1/2]

◆ conditionSet() [2/2]

◆ cpe()

◆ getBadFEDChannelPositions()

◆ hasBadComponents()

◆ index()

◆ isActive()

◆ isEmpty()

◆ measurements()

◆ recHits() [1/2]

◆ recHits() [2/2]

◆ setActive()

◆ setActiveThisEvent()

◆ setBadRocPositions()

◆ setEmpty()

◆ setIndex()

◆ specificGeomDet()

◆ update()

Member Data Documentation

◆ badRocPositions_

◆ id_

◆ index_

◆ theDetConditions