#include <TkPhase2OTMeasurementDet.h>

Inheritance diagram for TkPhase2OTMeasurementDet:

Public Types
typedef detset::const_iterator	const_iterator

typedef edmNew::DetSet< Phase2TrackerCluster1D >	detset

typedef ClusterParameterEstimator< Phase2TrackerCluster1D >::LocalValues	LocalValues

typedef edm::Ref< edmNew::DetSetVector< Phase2TrackerCluster1D >, Phase2TrackerCluster1D >	Phase2TrackerCluster1DRef

Public Types inherited from MeasurementDet
typedef TrackingRecHit::ConstRecHitContainer	RecHitContainer

using	SimpleHitContainer = std::vector< BaseTrackerRecHit * >

typedef tracking::TempMeasurements	TempMeasurements

Public Member Functions
TrackingRecHit::RecHitPointer	buildRecHit (const Phase2TrackerCluster1DRef &cluster, const LocalTrajectoryParameters &ltp) const

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

int	index () const
	Sets the list of bad ROCs, identified by the positions of their centers in the local coordinate frame. More...

bool	isActive (const MeasurementTrackerEvent &data) const override
	Is this module active in reconstruction? It must be both 'setActiveThisEvent' and 'setActive'. More...

bool	isEmpty (const Phase2OTMeasurementDetSet &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 (const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &, const MeasurementTrackerEvent &data, RecHitContainer &result, std::vector< float > &) 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 (Phase2OTMeasurementDetSet &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	setEmpty (Phase2OTMeasurementDetSet &data)

void	setIndex (int i)

const PixelGeomDetUnit &	specificGeomDet () const

	TkPhase2OTMeasurementDet (const GeomDet *gdet, Phase2OTMeasurementConditionSet &conditionSet)

void	update (Phase2OTMeasurementDetSet &data, const detset &detSet)

	~TkPhase2OTMeasurementDet () 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 Member Functions
Phase2OTMeasurementConditionSet &	conditionSet ()

const Phase2OTMeasurementConditionSet &	conditionSet () const

const ClusterParameterEstimator< Phase2TrackerCluster1D > *	cpe () const

Private Attributes
unsigned int	id_

int	index_

Phase2OTMeasurementConditionSet *	theDetConditions

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

TrackingRecHit::ConstRecHitPointer	theMissingHit

Detailed Description

Definition at line 16 of file TkPhase2OTMeasurementDet.h.

Member Typedef Documentation

◆ const_iterator

typedef detset::const_iterator TkPhase2OTMeasurementDet::const_iterator

Definition at line 21 of file TkPhase2OTMeasurementDet.h.

◆ detset

typedef edmNew::DetSet<Phase2TrackerCluster1D> TkPhase2OTMeasurementDet::detset

Definition at line 20 of file TkPhase2OTMeasurementDet.h.

◆ LocalValues

typedef ClusterParameterEstimator<Phase2TrackerCluster1D>::LocalValues TkPhase2OTMeasurementDet::LocalValues

Definition at line 22 of file TkPhase2OTMeasurementDet.h.

◆ Phase2TrackerCluster1DRef

typedef edm::Ref<edmNew::DetSetVector<Phase2TrackerCluster1D>, Phase2TrackerCluster1D> TkPhase2OTMeasurementDet::Phase2TrackerCluster1DRef

Definition at line 18 of file TkPhase2OTMeasurementDet.h.

Constructor & Destructor Documentation

◆ TkPhase2OTMeasurementDet()

TkPhase2OTMeasurementDet::TkPhase2OTMeasurementDet	(	const GeomDet *	gdet,
		Phase2OTMeasurementConditionSet &	conditionSet
	)

Definition at line 10 of file TkPhase2OTMeasurementDet.cc.

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

◆ ~TkPhase2OTMeasurementDet()

TkPhase2OTMeasurementDet::~TkPhase2OTMeasurementDet ( )

inlineoverride

Definition at line 34 of file TkPhase2OTMeasurementDet.h.

34 {}

Member Function Documentation

◆ buildRecHit()

TrackingRecHit::RecHitPointer TkPhase2OTMeasurementDet::buildRecHit	(	const Phase2TrackerCluster1DRef &	cluster,
		const LocalTrajectoryParameters &	ltp
	)		const

Definition at line 121 of file TkPhase2OTMeasurementDet.cc.

                                                                                                                 {
   const PixelGeomDetUnit& gdu(specificGeomDet());
   auto&& params = cpe()->localParameters(*cluster, gdu);
  
   return std::make_shared<Phase2TrackerRecHit1D>(params.first, params.second, fastGeomDet(), cluster);
 }

References cpe(), MeasurementDet::fastGeomDet(), ClusterParameterEstimator< T >::localParameters(), CalibrationSummaryClient_cfi::params, and specificGeomDet().

Referenced by recHits().

◆ conditionSet() [1/2]

Phase2OTMeasurementConditionSet& TkPhase2OTMeasurementDet::conditionSet ( )

inlineprivate

Definition at line 88 of file TkPhase2OTMeasurementDet.h.

88 { return *theDetConditions; }

◆ conditionSet() [2/2]

const Phase2OTMeasurementConditionSet& TkPhase2OTMeasurementDet::conditionSet ( ) const

inlineprivate

Definition at line 89 of file TkPhase2OTMeasurementDet.h.

89 { return *theDetConditions; }

◆ cpe()

const ClusterParameterEstimator<Phase2TrackerCluster1D>* TkPhase2OTMeasurementDet::cpe ( ) const

inlineprivate

Definition at line 91 of file TkPhase2OTMeasurementDet.h.

91 { return conditionSet().cpe(); }

Referenced by buildRecHit().

◆ hasBadComponents()

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

overridevirtual

Implements MeasurementDet.

Definition at line 165 of file TkPhase2OTMeasurementDet.cc.

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

Referenced by measurements().

◆ index()

int TkPhase2OTMeasurementDet::index ( ) const

inline

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

Clear the list of bad ROCs

Definition at line 79 of file TkPhase2OTMeasurementDet.h.

79 { return index_; }

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

◆ isActive()

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

inlineoverridevirtual

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

Implements MeasurementDet.

Definition at line 69 of file TkPhase2OTMeasurementDet.h.

69 { return data.phase2OTData().isActive(index()); }

References data.

Referenced by TkStackMeasurementDet::isActive(), measurements(), and recHits().

◆ isEmpty()

bool TkPhase2OTMeasurementDet::isEmpty ( const Phase2OTMeasurementDetSet & data ) const

inline

Definition at line 32 of file TkPhase2OTMeasurementDet.h.

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

References data.

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

◆ measurements()

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

overridevirtual

Implements MeasurementDet.

Definition at line 18 of file TkPhase2OTMeasurementDet.cc.

                                                                             {
   if (!isActive(data)) {
     result.add(theInactiveHit, 0.F);
     return true;
   }
  
   if (recHits(stateOnThisDet, est, data, result.hits, result.distances))
     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 data, F(), hasBadComponents(), isActive(), recHits(), mps_fire::result, MeasurementDet::theInactiveHit, and MeasurementDet::theMissingHit.

◆ recHits() [1/3]

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

overridevirtual

Implements MeasurementDet.

Definition at line 129 of file TkPhase2OTMeasurementDet.cc.

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

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

Referenced by measurements().

◆ recHits() [2/3]

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

overridevirtual

use the usual 5 sigma cut from the Traj to identify the column....

Reimplemented from MeasurementDet.

Definition at line 36 of file TkPhase2OTMeasurementDet.cc.

                                                                       {
   if UNLIKELY ((!isActive(data)) || isEmpty(data.phase2OTData()))
     return false;
  
   auto oldSize = result.size();
  
   const detset& detSet = data.phase2OTData().detSet(index());
   auto begin = &(data.phase2OTData().handle()->data().front());
   auto reject = [&](auto ci) -> bool {
     return (!data.phase2OTClustersToSkip().empty()) && data.phase2OTClustersToSkip()[ci - begin];
   };
  
   auto firstCluster = detSet.begin();
   auto lastCluster = detSet.end();
  
   // 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();
   auto ye = le.yy();
   if (lape.valid()) {
     ye += lape.yy();
   }
   // 5 sigma to be on the safe side
   ye = 5.f * std::sqrt(ye);
   LocalVector maxD(0, ye, 0);
   // pixel topology is rectangular: x and y are independent
   auto ymin = specificGeomDet().specificTopology().measurementPosition(stateOnThisDet.localPosition() - maxD);
   auto ymax = specificGeomDet().specificTopology().measurementPosition(stateOnThisDet.localPosition() + maxD);
   int utraj = ymin.x();
   // do not apply for iteration not cutting on propagation
   if (est.maxSagitta() >= 0) {
     int colMin = ymin.y();
     int colMax = ymax.y();
     firstCluster = std::find_if(firstCluster, detSet.end(), [colMin](const Phase2TrackerCluster1D& hit) {
       return int(hit.column()) >= colMin;
     });
     lastCluster = std::find_if(
         firstCluster, detSet.end(), [colMax](const Phase2TrackerCluster1D& hit) { return int(hit.column()) > colMax; });
   }
  
   while (firstCluster != lastCluster) {  // loop on each column
     auto const col = firstCluster->column();
     auto endCluster = std::find_if(
         firstCluster, detSet.end(), [col](const Phase2TrackerCluster1D& hit) { return hit.column() != col; });
     // find trajectory position in this column
     auto rightCluster = std::find_if(
         firstCluster, endCluster, [utraj](const Phase2TrackerCluster1D& hit) { return int(hit.firstStrip()) > utraj; });
     // search for compatible clusters...
     if (rightCluster != firstCluster) {
       // there are hits on the left of the utraj
       auto leftCluster = rightCluster;
       while (--leftCluster >= firstCluster) {
         if (reject(leftCluster))
           continue;
         Phase2TrackerCluster1DRef cluster = detSet.makeRefTo(data.phase2OTData().handle(), leftCluster);
         auto hit = buildRecHit(cluster, stateOnThisDet.localParameters());
         auto diffEst = est.estimate(stateOnThisDet, *hit);
         if (!diffEst.first)
           break;  // exit loop on first incompatible hit
         result.push_back(hit);
         diffs.push_back(diffEst.second);
       }
     }
     for (; rightCluster != endCluster; rightCluster++) {
       if (reject(rightCluster))
         continue;
       Phase2TrackerCluster1DRef cluster = detSet.makeRefTo(data.phase2OTData().handle(), rightCluster);
       auto hit = buildRecHit(cluster, stateOnThisDet.localParameters());
       auto diffEst = est.estimate(stateOnThisDet, *hit);
       if (!diffEst.first)
         break;  // exit loop on first incompatible hit
       result.push_back(hit);
       diffs.push_back(diffEst.second);
     }
     firstCluster = endCluster;
   }  // loop over columns
   return result.size() > oldSize;
 }

References edmNew::DetSet< T >::begin(), buildRecHit(), cuy::col, data, edmNew::DetSet< T >::end(), MeasurementEstimator::estimate(), MeasurementDet::fastGeomDet(), index(), createfilelist::int, isActive(), isEmpty(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), TrajectoryStateOnSurface::localPosition(), edmNew::DetSet< T >::makeRefTo(), MeasurementEstimator::maxSagitta(), Topology::measurementPosition(), LocalTrajectoryError::positionError(), reject, mps_fire::result, specificGeomDet(), PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), UNLIKELY, LocalError::valid(), L1TOccupancyClient_cfi::ymax, L1TOccupancyClient_cfi::ymin, and LocalError::yy().

◆ recHits() [3/3]

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

inlineoverridevirtual

Reimplemented from MeasurementDet.

Definition at line 44 of file TkPhase2OTMeasurementDet.h.

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

References cms::cuda::assert().

◆ setActive()

void TkPhase2OTMeasurementDet::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 62 of file TkPhase2OTMeasurementDet.h.

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

◆ setActiveThisEvent()

void TkPhase2OTMeasurementDet::setActiveThisEvent	(	Phase2OTMeasurementDetSet &	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 65 of file TkPhase2OTMeasurementDet.h.

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

References data.

◆ setEmpty()

void TkPhase2OTMeasurementDet::setEmpty ( Phase2OTMeasurementDetSet & data )

inline

Definition at line 31 of file TkPhase2OTMeasurementDet.h.

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

References data.

◆ setIndex()

void TkPhase2OTMeasurementDet::setIndex ( int i )

inline

Definition at line 80 of file TkPhase2OTMeasurementDet.h.

80 { index_ = i; }

References mps_fire::i.

◆ specificGeomDet()

const PixelGeomDetUnit& TkPhase2OTMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 56 of file TkPhase2OTMeasurementDet.h.

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

References MeasurementDet::fastGeomDet().

Referenced by buildRecHit(), and recHits().

◆ update()

void TkPhase2OTMeasurementDet::update	(	Phase2OTMeasurementDetSet &	data,
		const detset &	detSet
	)

inline

Definition at line 26 of file TkPhase2OTMeasurementDet.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

◆ id_

unsigned int TkPhase2OTMeasurementDet::id_

private

Definition at line 83 of file TkPhase2OTMeasurementDet.h.

◆ index_

int TkPhase2OTMeasurementDet::index_

private

Definition at line 86 of file TkPhase2OTMeasurementDet.h.

◆ theDetConditions

Phase2OTMeasurementConditionSet* TkPhase2OTMeasurementDet::theDetConditions

private

Definition at line 87 of file TkPhase2OTMeasurementDet.h.

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ const_iterator

◆ detset

◆ LocalValues

◆ Phase2TrackerCluster1DRef

Constructor & Destructor Documentation

◆ TkPhase2OTMeasurementDet()

◆ ~TkPhase2OTMeasurementDet()

Member Function Documentation

◆ buildRecHit()

◆ conditionSet() [1/2]

◆ conditionSet() [2/2]

◆ cpe()

◆ hasBadComponents()

◆ index()

◆ isActive()

◆ isEmpty()

◆ measurements()

◆ recHits() [1/3]

◆ recHits() [2/3]

◆ recHits() [3/3]

◆ setActive()

◆ setActiveThisEvent()

◆ setEmpty()

◆ setIndex()

◆ specificGeomDet()

◆ update()

Member Data Documentation

◆ id_

◆ index_

◆ theDetConditions