TkPixelMeasurementDet Class Reference

#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_
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), index_(0), 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.

{}

Member Function Documentation

accept()

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

inline

Definition at line 98 of file TkPixelMeasurementDet.h.

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

Referenced by esMonitoring.FDJsonServer::handle_accept().

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

buildRecHit()

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

Definition at line 69 of file TkPixelMeasurementDet.cc.

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

Referenced by compHits().

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

clearBadRocPositions()

void TkPixelMeasurementDet::clearBadRocPositions ( )

inline

Clear the list of bad ROCs.

Definition at line 77 of file TkPixelMeasurementDet.h.

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

References edmNew::DetSet< T >::begin(), buildRecHit(), haddnano::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().

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

conditionSet() [1/2]

PxMeasurementConditionSet& TkPixelMeasurementDet::conditionSet ( )

inlineprivate

Definition at line 92 of file TkPixelMeasurementDet.h.

{ return *theDetConditions; }

conditionSet() [2/2]

const PxMeasurementConditionSet& TkPixelMeasurementDet::conditionSet ( ) const

inlineprivate

Definition at line 93 of file TkPixelMeasurementDet.h.

{ return *theDetConditions; }

cpe()

const PixelClusterParameterEstimator* TkPixelMeasurementDet::cpe ( ) const

inlineprivate

Definition at line 95 of file TkPixelMeasurementDet.h.

Referenced by buildRecHit().

{ return conditionSet().pixelCPE(); }

getBadFEDChannelPositions()

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

inline

Definition at line 79 of file TkPixelMeasurementDet.h.

References data.

Referenced by hasBadComponents().

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

hasBadComponents()

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

overridevirtual

Implements MeasurementDet.

Definition at line 149 of file TkPixelMeasurementDet.cc.

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

Referenced by measurements().

                                                                                        {
  auto badFEDChannelPositions = getBadFEDChannelPositions(data);
  if (badRocPositions_.empty() && badFEDChannelPositions == nullptr)
    return false;

  auto lp = tsos.localPosition();
  float lexx = tsos.hasError() ? tsos.localError().positionError().xx() : 0.0f;
  float leyy = tsos.hasError() ? tsos.localError().positionError().yy() : 0.0f;
  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 * lexx) && (dy * dy < 9.f * leyy))
      return true;
  }

  if (badFEDChannelPositions == nullptr)
    return false;
  float dx = 3.f * std::sqrt(lexx) + theRocWidth, dy = 3.f * std::sqrt(leyy) + 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;
}

index()

int TkPixelMeasurementDet::index ( ) const

inline

Definition at line 84 of file TkPixelMeasurementDet.h.

Referenced by compHits().

{ return index_; }

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.

References data.

Referenced by compHits(), and measurements().

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

isEmpty()

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

inline

Definition at line 33 of file TkPixelMeasurementDet.h.

References data.

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

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

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.

References compHits(), runTheMatrix::const, 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().

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

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.

References compHits(), and data.

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

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.

References cms::cuda::assert().

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

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.

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

References data.

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

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.

{ badRocPositions_.swap(positions); }

setEmpty()

void TkPixelMeasurementDet::setEmpty ( PxMeasurementDetSet &  data )

inline

Definition at line 32 of file TkPixelMeasurementDet.h.

References data.

{ data.setEmpty(index()); }

setIndex()

void TkPixelMeasurementDet::setIndex ( int  i )

inline

Definition at line 85 of file TkPixelMeasurementDet.h.

References mps_fire::i.

{ index_ = i; }

specificGeomDet()

const PixelGeomDetUnit& TkPixelMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 59 of file TkPixelMeasurementDet.h.

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

Referenced by buildRecHit(), and compHits().

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

update()

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

inline

Definition at line 27 of file TkPixelMeasurementDet.h.

References data.

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

badRocPositions_

std::vector<LocalPoint> TkPixelMeasurementDet::badRocPositions_

private

Definition at line 88 of file TkPixelMeasurementDet.h.

Referenced by hasBadComponents().

index_

int TkPixelMeasurementDet::index_

private

Definition at line 90 of file TkPixelMeasurementDet.h.

theDetConditions

PxMeasurementConditionSet* TkPixelMeasurementDet::theDetConditions

private

Definition at line 91 of file TkPixelMeasurementDet.h.