#include <TkStackMeasurementDet.h>

Inheritance diagram for TkStackMeasurementDet:

Public Types
typedef detset::const_iterator	const_iterator

typedef edmNew::DetSet< Phase2TrackerCluster1D >	detset

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
bool	hasBadComponents (const TrajectoryStateOnSurface &tsos, const MeasurementTrackerEvent &data) const override
	return TRUE if at least one of the lower and upper components has badChannels More...

void	init (const MeasurementDet lowerDet, const MeasurementDet upperDet)

bool	isActive (const MeasurementTrackerEvent &data) const override
	return TRUE if both lower and upper components are active More...

bool	isEmpty (const Phase2OTMeasurementDetSet &data) const

const TkPhase2OTMeasurementDet *	lowerDet () const

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

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

const StackGeomDet &	specificGeomDet () const

	TkStackMeasurementDet (const StackGeomDet gdet, const PixelClusterParameterEstimator cpe)

const TkPhase2OTMeasurementDet *	upperDet () const

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

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

const Surface &	surface () const

virtual	~MeasurementDet ()=default

Private Attributes
const TkPhase2OTMeasurementDet *	theLowerDet

const PixelClusterParameterEstimator *	thePixelCPE

const TkPhase2OTMeasurementDet *	theUpperDet

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

TrackingRecHit::ConstRecHitPointer	theMissingHit

Detailed Description

Definition at line 14 of file TkStackMeasurementDet.h.

Member Typedef Documentation

◆ const_iterator

typedef detset::const_iterator TkStackMeasurementDet::const_iterator

Definition at line 25 of file TkStackMeasurementDet.h.

◆ detset

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

Definition at line 24 of file TkStackMeasurementDet.h.

◆ Phase2TrackerCluster1DRef

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

Definition at line 22 of file TkStackMeasurementDet.h.

Constructor & Destructor Documentation

◆ TkStackMeasurementDet()

TkStackMeasurementDet::TkStackMeasurementDet	(	const StackGeomDet *	gdet,
		const PixelClusterParameterEstimator *	cpe
	)

Definition at line 7 of file TkStackMeasurementDet.cc.

8 : MeasurementDet(gdet), thePixelCPE(cpe), theLowerDet(nullptr), theUpperDet(nullptr) {}

TkStackMeasurementDet::theLowerDet

const TkPhase2OTMeasurementDet * theLowerDet

Definition: TkStackMeasurementDet.h:50

MeasurementDet::MeasurementDet

MeasurementDet(const GeomDet *gdet)

Definition: MeasurementDet.h:23

TkStackMeasurementDet::thePixelCPE

const PixelClusterParameterEstimator * thePixelCPE

Definition: TkStackMeasurementDet.h:49

TkStackMeasurementDet::theUpperDet

const TkPhase2OTMeasurementDet * theUpperDet

Definition: TkStackMeasurementDet.h:51

Member Function Documentation

◆ hasBadComponents()

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

inlineoverridevirtual

return TRUE if at least one of the lower and upper components has badChannels

Implements MeasurementDet.

Definition at line 44 of file TkStackMeasurementDet.h.

References data, lowerDet(), and upperDet().

                                                                                                                   {
     return (lowerDet()->hasBadComponents(tsos, data) || upperDet()->hasBadComponents(tsos, data));
   }

◆ init()

void TkStackMeasurementDet::init	(	const MeasurementDet *	lowerDet,
		const MeasurementDet *	upperDet
	)

Definition at line 10 of file TkStackMeasurementDet.cc.

References lowerDet(), theLowerDet, theUpperDet, and upperDet().

Referenced by MeasurementTrackerImpl::initStackDet().

                                                                                                {
   theLowerDet = dynamic_cast<const TkPhase2OTMeasurementDet*>(lowerDet);
   theUpperDet = dynamic_cast<const TkPhase2OTMeasurementDet*>(upperDet);
 
   if ((theLowerDet == nullptr) || (theUpperDet == nullptr)) {
     throw MeasurementDetException(
         "TkStackMeasurementDet ERROR: Trying to glue a det which is not a TkPhase2OTMeasurementDet");
   }
 }

◆ isActive()

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

inlineoverridevirtual

return TRUE if both lower and upper components are active

Implements MeasurementDet.

Definition at line 36 of file TkStackMeasurementDet.h.

References data, TkPhase2OTMeasurementDet::isActive(), lowerDet(), and upperDet().

Referenced by measurements(), and recHits().

                                                                     {
     return lowerDet()->isActive(data) && upperDet()->isActive(data);
   }

◆ isEmpty()

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

inline

Definition at line 39 of file TkStackMeasurementDet.h.

References data, lowerDet(), and upperDet().

Referenced by plotting.Plot::clone().

                                                             {
     return data.empty(lowerDet()->index()) || data.empty(upperDet()->index());
   }

◆ lowerDet()

const TkPhase2OTMeasurementDet* TkStackMeasurementDet::lowerDet ( ) const

inline

Definition at line 32 of file TkStackMeasurementDet.h.

References theLowerDet.

Referenced by hasBadComponents(), init(), isActive(), isEmpty(), and recHits().

32 { return theLowerDet; }

TkStackMeasurementDet::theLowerDet

const TkPhase2OTMeasurementDet * theLowerDet

Definition: TkStackMeasurementDet.h:50

◆ measurements()

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

overridevirtual

Implements MeasurementDet.

Definition at line 100 of file TkStackMeasurementDet.cc.

References data, MeasurementEstimator::estimate(), F(), isActive(), LogDebug, LogTrace, eostools::move(), recHits(), mps_fire::result, MeasurementDet::theInactiveHit, and MeasurementDet::theMissingHit.

                                                                          {
   LogDebug("MeasurementTracker") << "TkStackMeasurementDet::measurements";
 
   if (!isActive(data)) {
     result.add(theInactiveHit, 0.F);
     return true;
   }
 
   LogTrace("MeasurementTracker") << " is active";
 
   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) {
       LogDebug("MeasurementTracker") << "New vh added with chi2: " << diffEst.second;
       result.add(std::move(hit), diffEst.second);
     }
   }
 
   if (result.size() > oldSize)
     return true;
 
   // create a TrajectoryMeasurement with an invalid RecHit and zero estimate
   result.add(theMissingHit, 0.F);
   LogDebug("MeasurementTracker") << "adding missing hit";
   return false;
 }

◆ recHits()

TkStackMeasurementDet::RecHitContainer TkStackMeasurementDet::recHits	(	const TrajectoryStateOnSurface &	ts,
		const MeasurementTrackerEvent &	data
	)		const

overridevirtual

Implements MeasurementDet.

Definition at line 20 of file TkStackMeasurementDet.cc.

References edmNew::DetSet< T >::begin(), data, HLT_2023v12_cff::distance, edmNew::DetSet< T >::empty(), edmNew::DetSet< T >::end(), GeomDet::geographicalId(), isActive(), LogTrace, lowerDet(), edmNew::DetSet< T >::offset(), mps_fire::result, edmNew::DetSet< T >::size(), specificGeomDet(), and upperDet().

Referenced by measurements().

                                                                                                                  {
   RecHitContainer result;
 
   if (data.phase2OTVectorHits().empty())
     return result;
   LogTrace("MeasurementTracker") << " is not empty";
   if (!isActive(data))
     return result;
   LogTrace("MeasurementTracker") << " and is active";
 
   //find clusters to skip
   const detset& lowerDetSet = data.phase2OTData().detSet(lowerDet()->index());
   const detset& upperDetSet = data.phase2OTData().detSet(upperDet()->index());
   std::vector<bool> skipClustersUpper(data.phase2OTClustersToSkip().empty() ? 0 : upperDetSet.size(), false);
   std::vector<bool> skipClustersLower(data.phase2OTClustersToSkip().empty() ? 0 : lowerDetSet.size(), false);
 
   const Phase2TrackerCluster1D* begin = nullptr;
   if (!data.phase2OTData().handle()->data().empty()) {
     begin = &(data.phase2OTData().handle()->data().front());
   }
   if (!data.phase2OTClustersToSkip().empty()) {
     if (!lowerDetSet.empty()) {
       for (const_iterator cil = lowerDetSet.begin(); cil != lowerDetSet.end(); ++cil) {
         if (cil < begin) {
           edm::LogError("IndexMisMatch") << "TkStackMeasurementDet cannot create hit because of index mismatch.";
           return result;
         }
         unsigned int indexl = cil - begin;
         if (data.phase2OTClustersToSkip()[indexl]) {
           int iLocalL = std::distance(lowerDetSet.begin(), cil);
           skipClustersLower[iLocalL] = true;
         }
       }
     }
     if (!upperDetSet.empty()) {
       for (const_iterator ciu = upperDetSet.begin(); ciu != upperDetSet.end(); ++ciu) {
         if (ciu < begin) {
           edm::LogError("IndexMisMatch") << "TkStackMeasurementDet cannot create hit because of index mismatch.";
           return result;
         }
         unsigned int indexu = ciu - begin;
         if (data.phase2OTClustersToSkip()[indexu]) {
           int iLocalU = std::distance(upperDetSet.begin(), ciu);
           skipClustersUpper[iLocalU] = true;
         }
       }
     }
   }
   DetId detIdStack = specificGeomDet().geographicalId();
 
   auto iterator = data.phase2OTVectorHits().find(detIdStack);
   if (iterator == data.phase2OTVectorHits().end())
     return result;
   for (const auto& vecHit : data.phase2OTVectorHits()[detIdStack]) {
     if (!data.phase2OTClustersToSkip().empty()) {
       if (skipClustersLower[vecHit.lowerCluster().key() - lowerDetSet.offset()])
         continue;
       if (skipClustersUpper[vecHit.upperCluster().key() - upperDetSet.offset()])
         continue;
     }
     result.push_back(std::make_shared<VectorHit>(vecHit));
   }
 
   iterator = data.phase2OTVectorHitsRej().find(detIdStack);
   if (iterator == data.phase2OTVectorHitsRej().end())
     return result;
   for (const auto& vecHit : data.phase2OTVectorHitsRej()[detIdStack]) {
     if (!data.phase2OTClustersToSkip().empty()) {
       if (skipClustersLower[vecHit.lowerCluster().key() - lowerDetSet.offset()])
         continue;
       if (skipClustersUpper[vecHit.upperCluster().key() - upperDetSet.offset()])
         continue;
     }
     result.push_back(std::make_shared<VectorHit>(vecHit));
   }
 
   return result;
 }