TkStackMeasurementDet Class Reference

#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

typedef detset::const_iterator TkStackMeasurementDet::const_iterator

Definition at line 25 of file TkStackMeasurementDet.h.

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

Definition at line 24 of file TkStackMeasurementDet.h.

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

Definition at line 22 of file TkStackMeasurementDet.h.

Constructor & Destructor Documentation

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

Definition at line 7 of file TkStackMeasurementDet.cc.

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

Member Function Documentation

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 lowerDet(), and upperDet().

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

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

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 TkPhase2OTMeasurementDet::isActive(), lowerDet(), and upperDet().

Referenced by measurements(), and recHits().

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

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

inline

Definition at line 39 of file TkStackMeasurementDet.h.

References Phase2OTMeasurementDetSet::empty(), lowerDet(), and upperDet().

Referenced by plotting.Plot::clone().

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

const TkPhase2OTMeasurementDet* TkStackMeasurementDet::lowerDet ( ) const

inline

Definition at line 32 of file TkStackMeasurementDet.h.

References theLowerDet.

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

{ return theLowerDet; }

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 tracking::TempMeasurements::add(), MeasurementEstimator::estimate(), F(), isActive(), LogDebug, LogTrace, eostools::move(), recHits(), tracking::TempMeasurements::size(), 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;
}

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

overridevirtual

Implements MeasurementDet.

Definition at line 20 of file TkStackMeasurementDet.cc.

References SplitLinear::begin, edmNew::DetSet< T >::begin(), HLT_FULL_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;
}

const StackGeomDet& TkStackMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 21 of file TkStackMeasurementDet.h.

References runTheMatrix::const, and MeasurementDet::fastGeomDet().

Referenced by MeasurementTrackerImpl::initStackDet(), and recHits().

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

const TkPhase2OTMeasurementDet* TkStackMeasurementDet::upperDet ( ) const

inline

Definition at line 33 of file TkStackMeasurementDet.h.

References theUpperDet.

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

{ return theUpperDet; }

Member Data Documentation

const TkPhase2OTMeasurementDet* TkStackMeasurementDet::theLowerDet

private

Definition at line 50 of file TkStackMeasurementDet.h.

Referenced by init(), and lowerDet().

const PixelClusterParameterEstimator* TkStackMeasurementDet::thePixelCPE

private

Definition at line 49 of file TkStackMeasurementDet.h.

const TkPhase2OTMeasurementDet* TkStackMeasurementDet::theUpperDet

private

Definition at line 51 of file TkStackMeasurementDet.h.

Referenced by init(), and upperDet().