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

21 : 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 149 of file LayerMeasurements.cc.

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

References tracking::TempMeasurements::clear(), GeometricSearchDet::compatible(), tracking::TempMeasurements::distances, F(), GeomDet::geographicalId(), GeometricSearchDet::groupedCompatibleDets(), tracking::TempMeasurements::hits, mps_fire::i, TrackingRecHit::inactive, MeasurementDetWithData::isNull(), LogDebug, MeasurementDetWithData::measurements(), eostools::move(), mps_fire::result, tracking::TempMeasurements::size(), and GeometricSearchDet::surface().

◆ idToDet()

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

inline

Definition at line 39 of file LayerMeasurements.h.

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

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

◆ measurements()

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

Definition at line 122 of file LayerMeasurements.cc.

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

References GeometricSearchDet::compatible(), GeometricSearchDet::compatibleDets(), F(), get, TrackingRecHit::inactive, LogDebug, mps_fire::result, and GeometricSearchDet::surface().

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

◆ recHits()

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

Definition at line 107 of file LayerMeasurements.cc.

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

References GeometricSearchDet::compatibleDets(), MeasurementDetWithData::recHits(), and mps_fire::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_