TkPhase2OTMeasurementDet Class Reference

#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

typedef detset::const_iterator TkPhase2OTMeasurementDet::const_iterator

Definition at line 21 of file TkPhase2OTMeasurementDet.h.

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

Definition at line 20 of file TkPhase2OTMeasurementDet.h.

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

Definition at line 22 of file TkPhase2OTMeasurementDet.h.

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

Definition at line 18 of file TkPhase2OTMeasurementDet.h.

Constructor & Destructor Documentation

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

inlineoverride

Definition at line 34 of file TkPhase2OTMeasurementDet.h.

{}

Member Function Documentation

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

Definition at line 121 of file TkPhase2OTMeasurementDet.cc.

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

Referenced by recHits().

                                                                                                                {
  const PixelGeomDetUnit& gdu(specificGeomDet());
  auto&& params = cpe()->localParameters(*cluster, gdu);

  return std::make_shared<Phase2TrackerRecHit1D>(params.first, params.second, fastGeomDet(), cluster);
}

Phase2OTMeasurementConditionSet& TkPhase2OTMeasurementDet::conditionSet ( )

inlineprivate

Definition at line 88 of file TkPhase2OTMeasurementDet.h.

{ return *theDetConditions; }

const Phase2OTMeasurementConditionSet& TkPhase2OTMeasurementDet::conditionSet ( ) const

inlineprivate

Definition at line 89 of file TkPhase2OTMeasurementDet.h.

{ return *theDetConditions; }

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

inlineprivate

Definition at line 91 of file TkPhase2OTMeasurementDet.h.

Referenced by buildRecHit().

{ return conditionSet().cpe(); }

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

overridevirtual

Implements MeasurementDet.

Definition at line 165 of file TkPhase2OTMeasurementDet.cc.

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

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

{ return index_; }

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.

References Phase2OTMeasurementDetSet::isActive(), and MeasurementTrackerEvent::phase2OTData().

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

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

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

inline

Definition at line 32 of file TkPhase2OTMeasurementDet.h.

References Phase2OTMeasurementDetSet::empty().

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

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

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.

References tracking::TempMeasurements::add(), tracking::TempMeasurements::distances, F(), hasBadComponents(), tracking::TempMeasurements::hits, isActive(), recHits(), MeasurementDet::theInactiveHit, and MeasurementDet::theMissingHit.

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

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

overridevirtual

Implements MeasurementDet.

Definition at line 129 of file TkPhase2OTMeasurementDet.cc.

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

Referenced by measurements().

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

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.

References SplitLinear::begin, edmNew::DetSet< T >::begin(), buildRecHit(), cuy::col, runTheMatrix::const, Phase2OTMeasurementDetSet::detSet(), edmNew::DetSet< T >::end(), MeasurementEstimator::estimate(), MeasurementDet::fastGeomDet(), Phase2TrackerCluster1D::firstStrip(), Phase2OTMeasurementDetSet::handle(), index(), isActive(), isEmpty(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), TrajectoryStateOnSurface::localPosition(), edmNew::DetSet< T >::makeRefTo(), MeasurementEstimator::maxSagitta(), Topology::measurementPosition(), MeasurementTrackerEvent::phase2OTClustersToSkip(), MeasurementTrackerEvent::phase2OTData(), LocalTrajectoryError::positionError(), reject, specificGeomDet(), PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), UNLIKELY, LocalError::valid(), SiStrip_OfflineMonitoring_cff::ymax, SiStrip_OfflineMonitoring_cff::ymin, and LocalError::yy().

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

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.

References cms::cuda::assert().

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

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.

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

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.

References Phase2OTMeasurementDetSet::setActiveThisEvent().

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

void TkPhase2OTMeasurementDet::setEmpty ( Phase2OTMeasurementDetSet &  data )

inline

Definition at line 31 of file TkPhase2OTMeasurementDet.h.

References Phase2OTMeasurementDetSet::setEmpty().

{ data.setEmpty(index()); }

void TkPhase2OTMeasurementDet::setIndex ( int  i )

inline

Definition at line 80 of file TkPhase2OTMeasurementDet.h.

References mps_fire::i.

{ index_ = i; }

const PixelGeomDetUnit& TkPhase2OTMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 56 of file TkPhase2OTMeasurementDet.h.

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

Referenced by buildRecHit(), and recHits().

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

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

inline

Definition at line 26 of file TkPhase2OTMeasurementDet.h.

References Phase2OTMeasurementDetSet::setActiveThisEvent(), and Phase2OTMeasurementDetSet::update().

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

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

Member Data Documentation

unsigned int TkPhase2OTMeasurementDet::id_

private

Definition at line 83 of file TkPhase2OTMeasurementDet.h.

int TkPhase2OTMeasurementDet::index_

private

Definition at line 86 of file TkPhase2OTMeasurementDet.h.

Phase2OTMeasurementConditionSet* TkPhase2OTMeasurementDet::theDetConditions

private

Definition at line 87 of file TkPhase2OTMeasurementDet.h.