LayerMeasurements Class Reference

#include <LayerMeasurements.h>

Public Member Functions
std::vector< TrajectoryMeasurementGroup >	groupedMeasurements (const DetLayer &layer, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
MeasurementDetWithData	idToDet (const DetId &id) const
	LayerMeasurements (const MeasurementDetSystem &detSystem, const MeasurementTrackerEvent &data)
std::vector< TrajectoryMeasurement >	measurements (const DetLayer &layer, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
std::vector< BaseTrackerRecHit * >	recHits (const DetLayer &layer, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const

Private Attributes
MeasurementTrackerEvent const &	data_
MeasurementDetSystem const &	detSystem_

Detailed Description

Definition at line 18 of file LayerMeasurements.h.

Constructor & Destructor Documentation

LayerMeasurements()

LayerMeasurements::LayerMeasurements ( const MeasurementDetSystem &  detSystem, const MeasurementTrackerEvent &  data  )

inline

Definition at line 20 of file LayerMeasurements.h.

      : detSystem_(detSystem), data_(data) {}

Member Function Documentation

groupedMeasurements()

vector< TrajectoryMeasurementGroup > LayerMeasurements::groupedMeasurements	(	const DetLayer &	layer,
		const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est
	)		const

Definition at line 150 of file LayerMeasurements.cc.

References tracking::TempMeasurements::clear(), tracking::TempMeasurements::distances, F(), GeomDet::geographicalId(), tracking::TempMeasurements::hits, mps_fire::i, TrackingRecHit::inactive, MeasurementDetWithData::isNull(), phase1PixelTopology::layer, LogDebug, MeasurementDetWithData::measurements(), eostools::move(), mps_fire::result, tracking::TempMeasurements::size(), and jetUpdater_cfi::sort.

                                                                                                                 {
  vector<TrajectoryMeasurementGroup> result;

  vector<DetGroup>&& groups = layer.groupedCompatibleDets(startingState, prop, est);
  result.reserve(groups.size());

  tracking::TempMeasurements tmps;
  for (auto& grp : groups) {
    if (grp.empty())
      continue;

    vector<TrajectoryMeasurement> tmpVec;
    for (auto const& det : grp) {
      MeasurementDetWithData mdet = detSystem_.idToDet(det.det()->geographicalId(), data_);
      if (mdet.isNull()) {
        throw MeasurementDetException("MeasurementDet not found");
      }
      if (mdet.measurements(det.trajectoryState(), est, tmps))
        for (std::size_t i = 0; i != tmps.size(); ++i)
          tmpVec.emplace_back(det.trajectoryState(), std::move(tmps.hits[i]), tmps.distances[i], &layer);
      tmps.clear();
    }

    // sort the final result
    LogDebug("LayerMeasurements") << "Sorting " << tmpVec.size() << " measurements in this grp.";
    sort(tmpVec.begin(), tmpVec.end(), TrajMeasLessEstim());
    addInvalidMeas(tmpVec, grp, layer);
    result.emplace_back(std::move(tmpVec), std::move(grp));
  }

  // if the result is empty check if the layer is compatible (for invalid measurement)
  if (result.empty()) {
    pair<bool, TrajectoryStateOnSurface> compat = layer.compatible(startingState, prop, est);
    if (compat.first) {
      vector<TrajectoryMeasurement> tmVec;
      tmVec.emplace_back(compat.second,
                         std::make_shared<InvalidTrackingRecHitNoDet>(layer.surface(), TrackingRecHit::inactive),
                         0.F,
                         &layer);
      result.emplace_back(std::move(tmVec), DetGroup());
    }
  }
  return result;
}

idToDet()

MeasurementDetWithData LayerMeasurements::idToDet ( const DetId &  id ) const

inline

Definition at line 39 of file LayerMeasurements.h.

References data_, detSystem_, and MeasurementDetSystem::idToDet().

{ return detSystem_.idToDet(id, data_); }

measurements()

vector< TrajectoryMeasurement > LayerMeasurements::measurements	(	const DetLayer &	layer,
		const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est
	)		const

Definition at line 123 of file LayerMeasurements.cc.

References F(), TrackingRecHit::inactive, phase1PixelTopology::layer, LogDebug, and mps_fire::result.

Referenced by SiPixelHitEfficiencySource::analyze(), OutInConversionSeedFinder::completeSeed(), InOutConversionSeedFinder::completeSeed(), NuclearInteractionFinder::findMeasurementsFromTSOS(), CosmicTrackingRegion::hits_(), and OutInConversionSeedFinder::startSeed().

                                                                                                     {
  typedef DetLayer::DetWithState DetWithState;

  vector<DetWithState> const& compatDets = layer.compatibleDets(startingState, prop, est);

  if (!compatDets.empty())
    return get(detSystem_, data_, layer, compatDets, startingState, prop, est);

  vector<TrajectoryMeasurement> result;
  pair<bool, TrajectoryStateOnSurface> compat = layer.compatible(startingState, prop, est);

  if (compat.first) {
    result.push_back(
        TrajectoryMeasurement(compat.second,
                              std::make_shared<InvalidTrackingRecHitNoDet>(layer.surface(), TrackingRecHit::inactive),
                              0.F,
                              &layer));
    LogDebug("LayerMeasurements") << "adding a missing hit.";
  } else
    LogDebug("LayerMeasurements") << "adding not measurement.";

  return result;
}

recHits()

std::vector< BaseTrackerRecHit * > LayerMeasurements::recHits	(	const DetLayer &	layer,
		const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est
	)		const

Definition at line 108 of file LayerMeasurements.cc.

References phase1PixelTopology::layer, MeasurementDetWithData::recHits(), and mps_fire::result.

                                                                                                  {
  std::vector<BaseTrackerRecHit*> result;
  auto const& compatDets = layer.compatibleDets(startingState, prop, est);
  if (compatDets.empty())
    return result;
  for (auto const& ds : compatDets) {
    auto mdet = detSystem_.idToDet(ds.first->geographicalId(), data_);
    mdet.recHits(result, ds.second, est);
  }
  return result;
}

Member Data Documentation

data_

MeasurementTrackerEvent const& LayerMeasurements::data_

private

Definition at line 43 of file LayerMeasurements.h.

Referenced by idToDet().

detSystem_

MeasurementDetSystem const& LayerMeasurements::detSystem_

private

Definition at line 42 of file LayerMeasurements.h.

Referenced by idToDet().