LayerMeasurements Class Reference

#include <LayerMeasurements.h>

Public Types
using	SimpleHitContainer = std::vector< BaseTrackerRecHit * >

Public Member Functions
void	addInvalidMeas (std::vector< TrajectoryMeasurement > &measVec, const DetGroup &group, const DetLayer &layer) const
std::vector < TrajectoryMeasurementGroup >	groupedMeasurements (const DetLayer &layer, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
MeasurementDetWithData	idToDet (const DetId &id) const
	LayerMeasurements ()
	LayerMeasurements (const MeasurementDetSystem &detSystem, const MeasurementTrackerEvent &data)
std::vector < TrajectoryMeasurement >	measurements (const DetLayer &layer, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
bool	recHits (SimpleHitContainer &result, const DetLayer &layer, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const

Private Attributes
const MeasurementTrackerEvent *	theData
const MeasurementDetSystem *	theDetSystem

Detailed Description

Definition at line 21 of file LayerMeasurements.h.

Member Typedef Documentation

using LayerMeasurements::SimpleHitContainer = std::vector<BaseTrackerRecHit *>

Definition at line 24 of file LayerMeasurements.h.

Constructor & Destructor Documentation

LayerMeasurements::LayerMeasurements ( )

inline

Definition at line 28 of file LayerMeasurements.h.

: theDetSystem(0), theData(0) {}

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

inline

Definition at line 31 of file LayerMeasurements.h.

                                                                                                 :
    theDetSystem(&detSystem), theData(&data) {}

Member Function Documentation

void LayerMeasurements::addInvalidMeas	(	std::vector< TrajectoryMeasurement > &	measVec,
		const DetGroup &	group,
		const DetLayer &	layer
	)		const

Definition at line 182 of file LayerMeasurements.cc.

References TrackingRecHit::missing.

{
  if (!measVec.empty()) {
    // invalidMeas on Det of most compatible hit
    measVec.emplace_back(measVec.front().predictedState(), 
             std::make_shared<InvalidTrackingRecHit>(*measVec.front().recHit()->det(), TrackingRecHit::missing),
             0.,&layer);
  }
  else if (!group.empty()) {
    // invalid state on first compatible Det
    measVec.emplace_back(group.front().trajectoryState(), 
             std::make_shared<InvalidTrackingRecHit>(*group.front().det(), TrackingRecHit::missing), 0.,&layer);
  }
}

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

Definition at line 138 of file LayerMeasurements.cc.

References tracking::TempMeasurements::clear(), GeometricSearchDet::compatible(), tracking::TempMeasurements::distances, F(), GeomDet::geographicalId(), GeometricSearchDet::groupedCompatibleDets(), tracking::TempMeasurements::hits, i, TrackingRecHit::inactive, MeasurementDetWithData::isNull(), MeasurementDetWithData::measurements(), query::result, tracking::TempMeasurements::size(), python.multivaluedict::sort(), and GeometricSearchDet::surface().

                                                           {
  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 = theDetSystem->idToDet(det.det()->geographicalId(), *theData);
      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
    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;
}

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

inline

Definition at line 59 of file LayerMeasurements.h.

References MeasurementDetSystem::idToDet(), theData, and theDetSystem.

                                                        {
     return theDetSystem->idToDet(id, *theData);
  }

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

Definition at line 110 of file LayerMeasurements.cc.

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

Referenced by SiPixelHitEfficiencySource::analyze(), HitEff::analyze(), OutInConversionSeedFinder::completeSeed(), InOutConversionSeedFinder::completeSeed(), TSGFromPropagation::findMeasurements(), FastTSGFromPropagation::findMeasurements(), NuclearInteractionFinder::findMeasurementsFromTSOS(), SiStripElectronSeedGenerator::findSeedsFromCluster(), CosmicTrackingRegion::hits_(), PixelMatchNextLayers::PixelMatchNextLayers(), and OutInConversionSeedFinder::startSeed().

                                                        {

  typedef DetLayer::DetWithState   DetWithState;

  vector<DetWithState>  const & compatDets = layer.compatibleDets( startingState, prop, est);
  
  if (!compatDets.empty())  return get(theDetSystem, theData, 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;  
  
}

bool LayerMeasurements::recHits	(	SimpleHitContainer &	result,
		const DetLayer &	layer,
		const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est
	)		const

Definition at line 92 of file LayerMeasurements.cc.

References GeometricSearchDet::compatibleDets(), MeasurementDetWithData::recHits(), and run_regression::ret.

Referenced by RectangularEtaPhiTrackingRegion::hits().

                                                        {

  auto  const & compatDets = layer.compatibleDets( startingState, prop, est);  
  if (compatDets.empty()) return false;
  bool ret=false;
  for ( auto const & ds : compatDets) {
    auto mdet = theDetSystem->idToDet(ds.first->geographicalId(), *theData);
    ret |=mdet.recHits(result,ds.second,est);
  }
  return ret;
}

Member Data Documentation

const MeasurementTrackerEvent* LayerMeasurements::theData

private

Definition at line 66 of file LayerMeasurements.h.

Referenced by idToDet().

const MeasurementDetSystem* LayerMeasurements::theDetSystem

private

Definition at line 65 of file LayerMeasurements.h.

Referenced by idToDet().