#include <MTDDetLayerMeasurements.h>

Public Types
typedef std::pair< const GeomDet *, TrajectoryStateOnSurface >	DetWithState

typedef std::vector< TrajectoryMeasurement >	MeasurementContainer

typedef std::vector< GenericTransientTrackingRecHit::RecHitPointer >	MTDRecHitContainer

Public Member Functions
MeasurementContainer	fastMeasurements (const DetLayer *layer, const TrajectoryStateOnSurface &theStateOnDet, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est)
	faster version in case the TrajectoryState on the surface of the GeomDet is already available More...

MeasurementContainer	fastMeasurements (const DetLayer *layer, const TrajectoryStateOnSurface &theStateOnDet, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, const edm::Event &iEvent)
	faster version in case the TrajectoryState on the surface of the GeomDet is already available More...

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

std::vector< TrajectoryMeasurementGroup >	groupedMeasurements (const DetLayer *layer, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, const edm::Event &iEvent)

MeasurementContainer	measurements (const DetLayer layer, const GeomDet det, const TrajectoryStateOnSurface &stateOnDet, const MeasurementEstimator &est, const edm::Event &iEvent)

MeasurementContainer	measurements (const DetLayer *layer, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est)
	returns TMeasurements in a DetLayer compatible with the TSOS. More...

MeasurementContainer	measurements (const DetLayer *layer, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, const edm::Event &iEvent)
	returns TMeasurements in a DetLayer compatible with the TSOS. More...

	MTDDetLayerMeasurements (edm::InputTag mtdlabel, edm::ConsumesCollector &iC)

MTDRecHitContainer	recHits (const DetLayer *layer)
	returns the rechits which are on the layer More...

MTDRecHitContainer	recHits (const DetLayer *layer, const edm::Event &iEvent)
	returns the rechits which are on the layer More...

void	setEvent (const edm::Event &)
	set event More...

virtual	~MTDDetLayerMeasurements ()

Private Member Functions
void	checkEvent () const
	check that the event is set, and throw otherwise More...

void	checkMTDRecHits ()

MTDRecHitContainer	recHits (const GeomDet *, const edm::Event &iEvent)
	obtain TrackingRecHits from a DetLayer More...

Private Attributes
edm::EDGetTokenT< MTDTrackingRecHit >	mtdToken_

const edm::Event *	theEvent

edm::Event::CacheIdentifier_t	theMTDEventCacheID

edm::Handle< edmNew::DetSetVector< MTDTrackingRecHit > >	theMTDRecHits

Detailed Description

The class to access recHits and TrajectoryMeasurements from DetLayer.

Author: B. Tannenwald Adapted from RecoMuon version.

The class to access recHits and TrajectoryMeasurements from DetLayer.

Author: B. Tannenwald Adapted from RecoMuon version.

Definition at line 30 of file MTDDetLayerMeasurements.h.

Member Typedef Documentation

◆ DetWithState

typedef std::pair<const GeomDet*, TrajectoryStateOnSurface> MTDDetLayerMeasurements::DetWithState

Definition at line 33 of file MTDDetLayerMeasurements.h.

◆ MeasurementContainer

typedef std::vector<TrajectoryMeasurement> MTDDetLayerMeasurements::MeasurementContainer

Definition at line 32 of file MTDDetLayerMeasurements.h.

◆ MTDRecHitContainer

typedef std::vector<GenericTransientTrackingRecHit::RecHitPointer> MTDDetLayerMeasurements::MTDRecHitContainer

Definition at line 34 of file MTDDetLayerMeasurements.h.

Constructor & Destructor Documentation

◆ MTDDetLayerMeasurements()

MTDDetLayerMeasurements::MTDDetLayerMeasurements	(	edm::InputTag	mtdlabel,
		edm::ConsumesCollector &	iC
	)

Definition at line 22 of file MTDDetLayerMeasurements.cc.

     : theMTDRecHits(), theMTDEventCacheID(0), theEvent(nullptr) {
   mtdToken_ = iC.consumes<MTDTrackingRecHit>(mtdlabel);
 }

References edm::ConsumesCollector::consumes(), and mtdToken_.

◆ ~MTDDetLayerMeasurements()

MTDDetLayerMeasurements::~MTDDetLayerMeasurements ( )

virtual

Definition at line 27 of file MTDDetLayerMeasurements.cc.

27 {}

Member Function Documentation

◆ checkEvent()

void MTDDetLayerMeasurements::checkEvent ( ) const

private

check that the event is set, and throw otherwise

Definition at line 194 of file MTDDetLayerMeasurements.cc.

                                                {
   if (!theEvent)
     throw cms::Exception("MTDDetLayerMeasurements") << "The event has not been set";
 }

References Exception, and theEvent.

Referenced by checkMTDRecHits(), fastMeasurements(), groupedMeasurements(), measurements(), and recHits().

◆ checkMTDRecHits()

void MTDDetLayerMeasurements::checkMTDRecHits ( )

private

Definition at line 50 of file MTDDetLayerMeasurements.cc.

                                               {
   LogDebug("Track|RecoMTD|MTDDetLayerMeasurements") << "Checking MTD RecHits";
   checkEvent();
   auto cacheID = theEvent->cacheIdentifier();
   if (cacheID == theMTDEventCacheID)
     return;
  
   {
     theEvent->getByToken(mtdToken_, theMTDRecHits);
     theMTDEventCacheID = cacheID;
   }
   if (!theMTDRecHits.isValid()) {
     throw cms::Exception("MTDDetLayerMeasurements") << "Cannot get MTD RecHits";
     LogDebug("Track|RecoMTD|MTDDetLayerMeasurements") << "Cannot get MTD RecHits";
   }
 }

References edm::Event::cacheIdentifier(), checkEvent(), Exception, edm::Event::getByToken(), edm::HandleBase::isValid(), LogDebug, mtdToken_, theEvent, theMTDEventCacheID, and theMTDRecHits.

Referenced by recHits().

◆ fastMeasurements() [1/2]

MeasurementContainer MTDDetLayerMeasurements::fastMeasurements	(	const DetLayer *	layer,
		const TrajectoryStateOnSurface &	theStateOnDet,
		const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est
	)

faster version in case the TrajectoryState on the surface of the GeomDet is already available

fastMeasurements method if already got the Event

Definition at line 145 of file MTDDetLayerMeasurements.cc.

                                                                                                 {
   checkEvent();
   return fastMeasurements(layer, theStateOnDet, startingState, prop, est, *theEvent);
 }

References checkEvent(), fastMeasurements(), and theEvent.

◆ fastMeasurements() [2/2]

MeasurementContainer MTDDetLayerMeasurements::fastMeasurements	(	const DetLayer *	layer,
		const TrajectoryStateOnSurface &	theStateOnDet,
		const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est,
		const edm::Event &	iEvent
	)

faster version in case the TrajectoryState on the surface of the GeomDet is already available

Definition at line 122 of file MTDDetLayerMeasurements.cc.

                                                                                        {
   MeasurementContainer result;
   MTDRecHitContainer rhs = recHits(layer, iEvent);
   for (MTDRecHitContainer::value_type irh : rhs) {
     MeasurementEstimator::HitReturnType estimate = est.estimate(theStateOnDet, (*irh));
     if (estimate.first) {
       result.push_back(TrajectoryMeasurement(theStateOnDet, irh, estimate.second, layer));
     }
   }
  
   if (!result.empty()) {
     sort(result.begin(), result.end(), TrajMeasLessEstim());
   }
  
   return result;
 }

References MeasurementEstimator::estimate(), iEvent, recHits(), and mps_fire::result.

Referenced by fastMeasurements().

◆ groupedMeasurements() [1/2]

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

Definition at line 154 of file MTDDetLayerMeasurements.cc.

                                      {
   checkEvent();
   return groupedMeasurements(layer, startingState, prop, est, *theEvent);
 }

References checkEvent(), groupedMeasurements(), and theEvent.

◆ groupedMeasurements() [2/2]

std::vector< TrajectoryMeasurementGroup > MTDDetLayerMeasurements::groupedMeasurements	(	const DetLayer *	layer,
		const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est,
		const edm::Event &	iEvent
	)

Definition at line 163 of file MTDDetLayerMeasurements.cc.

                             {
   std::vector<TrajectoryMeasurementGroup> result;
   // if we want to use the concept of InvalidRecHits,
   // we can reuse LayerMeasurements from TrackingTools/MeasurementDet
   std::vector<DetGroup> groups(layer->groupedCompatibleDets(startingState, prop, est));
  
   for (const auto& grp : groups) {
     std::vector<TrajectoryMeasurement> groupMeasurements;
     for (const auto& detAndStateItr : grp) {
       std::vector<TrajectoryMeasurement> detMeasurements =
           measurements(layer, detAndStateItr.det(), detAndStateItr.trajectoryState(), est, iEvent);
       groupMeasurements.insert(groupMeasurements.end(), detMeasurements.begin(), detMeasurements.end());
     }
  
     if (!groupMeasurements.empty())
       std::sort(groupMeasurements.begin(), groupMeasurements.end(), TrajMeasLessEstim());
  
     result.push_back(TrajectoryMeasurementGroup(groupMeasurements, grp));
   }
  
   return result;
 }

References GeometricSearchDet::groupedCompatibleDets(), iEvent, measurements(), and mps_fire::result.

Referenced by groupedMeasurements().

◆ measurements() [1/3]

MeasurementContainer MTDDetLayerMeasurements::measurements	(	const DetLayer *	layer,
		const GeomDet *	det,
		const TrajectoryStateOnSurface &	stateOnDet,
		const MeasurementEstimator &	est,
		const edm::Event &	iEvent
	)

Definition at line 97 of file MTDDetLayerMeasurements.cc.

                                                                                    {
   MeasurementContainer result;
  
   // Get the Segments which relies on the GeomDet given by compatibleDets
   MTDRecHitContainer mtdRecHits = recHits(det, iEvent);
  
   // Create the Trajectory Measurement
   for (auto rechit = mtdRecHits.begin(); rechit != mtdRecHits.end(); ++rechit) {
     MeasurementEstimator::HitReturnType estimate = est.estimate(stateOnDet, **rechit);
     LogDebug("RecoMTD") << "Dimension: " << (*rechit)->dimension() << " Chi2: " << estimate.second << std::endl;
     if (estimate.first) {
       result.push_back(TrajectoryMeasurement(stateOnDet, *rechit, estimate.second, layer));
     }
   }
  
   if (!result.empty())
     sort(result.begin(), result.end(), TrajMeasLessEstim());
  
   return result;
 }

References MeasurementEstimator::estimate(), iEvent, LogDebug, mtdRecHits_cfi::mtdRecHits, recHits(), and mps_fire::result.

Referenced by groupedMeasurements(), and measurements().

◆ measurements() [2/3]

MeasurementContainer MTDDetLayerMeasurements::measurements	(	const DetLayer *	layer,
		const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est
	)

returns TMeasurements in a DetLayer compatible with the TSOS.

measurements method if already got the Event

Definition at line 68 of file MTDDetLayerMeasurements.cc.

                                                                                             {
   checkEvent();
   return measurements(layer, startingState, prop, est, *theEvent);
 }

References checkEvent(), measurements(), and theEvent.

◆ measurements() [3/3]

MeasurementContainer MTDDetLayerMeasurements::measurements	(	const DetLayer *	layer,
		const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est,
		const edm::Event &	iEvent
	)

returns TMeasurements in a DetLayer compatible with the TSOS.

Definition at line 76 of file MTDDetLayerMeasurements.cc.

                                                                                    {
   MeasurementContainer result;
  
   std::vector<DetWithState> dss = layer->compatibleDets(startingState, prop, est);
   LogDebug("RecoMTD") << "compatibleDets: " << dss.size() << std::endl;
  
   for (const auto& dws : dss) {
     MeasurementContainer detMeasurements = measurements(layer, dws.first, dws.second, est, iEvent);
     result.insert(result.end(), detMeasurements.begin(), detMeasurements.end());
   }
  
   if (!result.empty())
     sort(result.begin(), result.end(), TrajMeasLessEstim());
  
   return result;
 }

References GeometricSearchDet::compatibleDets(), iEvent, LogDebug, measurements(), and mps_fire::result.

◆ recHits() [1/3]

MTDRecHitContainer MTDDetLayerMeasurements::recHits ( const DetLayer * layer )

returns the rechits which are on the layer

Definition at line 211 of file MTDDetLayerMeasurements.cc.

                                                                          {
   checkEvent();
   return recHits(layer, *theEvent);
 }

References checkEvent(), recHits(), and theEvent.

◆ recHits() [2/3]

MTDRecHitContainer MTDDetLayerMeasurements::recHits	(	const DetLayer *	layer,
		const edm::Event &	iEvent
	)

returns the rechits which are on the layer

Definition at line 199 of file MTDDetLayerMeasurements.cc.

                                                                                                  {
   MTDRecHitContainer rhs;
  
   std::vector<const GeomDet*> gds = layer->basicComponents();
  
   for (const GeomDet* igd : gds) {
     MTDRecHitContainer detHits = recHits(igd, iEvent);
     rhs.insert(rhs.end(), detHits.begin(), detHits.end());
   }
   return rhs;
 }

References GeometricSearchDet::basicComponents(), and iEvent.

Referenced by fastMeasurements(), measurements(), and recHits().

◆ recHits() [3/3]

MTDRecHitContainer MTDDetLayerMeasurements::recHits	(	const GeomDet *	geomDet,
		const edm::Event &	iEvent
	)

private

obtain TrackingRecHits from a DetLayer

Definition at line 29 of file MTDDetLayerMeasurements.cc.

                                                                                                   {
   DetId geoId = geomDet->geographicalId();
   theEvent = &iEvent;
   MTDRecHitContainer result;
  
   checkMTDRecHits();
  
   // Create the ChamberId
   DetId detId(geoId.rawId());
   LogDebug("Track|RecoMTD|MTDDetLayerMeasurements") << "(MTD): " << static_cast<MTDDetId>(detId) << std::endl;
  
   // Get the MTD-Segment which relies on this chamber
   //auto cmp = [](const unsigned one, const unsigned two) -> bool { return one < two; };
   auto detset = (*theMTDRecHits)[detId];
  
   for (const auto& rechit : detset)
     result.push_back(GenericTransientTrackingRecHit::build(geomDet, &rechit));
  
   return result;
 }

References GenericTransientTrackingRecHit::build(), checkMTDRecHits(), GeomDet::geographicalId(), iEvent, LogDebug, DetId::rawId(), mps_fire::result, and theEvent.

◆ setEvent()

void MTDDetLayerMeasurements::setEvent ( const edm::Event & event )

set event

Definition at line 192 of file MTDDetLayerMeasurements.cc.

192 { theEvent = &event; }

References event(), and theEvent.

Member Data Documentation

◆ mtdToken_

edm::EDGetTokenT<MTDTrackingRecHit> MTDDetLayerMeasurements::mtdToken_

private

Definition at line 102 of file MTDDetLayerMeasurements.h.

Referenced by checkMTDRecHits(), and MTDDetLayerMeasurements().

◆ theEvent

const edm::Event* MTDDetLayerMeasurements::theEvent

private

Definition at line 112 of file MTDDetLayerMeasurements.h.

Referenced by checkEvent(), checkMTDRecHits(), fastMeasurements(), groupedMeasurements(), measurements(), recHits(), and setEvent().

◆ theMTDEventCacheID

edm::Event::CacheIdentifier_t MTDDetLayerMeasurements::theMTDEventCacheID

private

Definition at line 110 of file MTDDetLayerMeasurements.h.

Referenced by checkMTDRecHits().

◆ theMTDRecHits

edm::Handle<edmNew::DetSetVector<MTDTrackingRecHit> > MTDDetLayerMeasurements::theMTDRecHits

private

Definition at line 105 of file MTDDetLayerMeasurements.h.

Referenced by checkMTDRecHits().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ DetWithState

◆ MeasurementContainer

◆ MTDRecHitContainer

Constructor & Destructor Documentation

◆ MTDDetLayerMeasurements()

◆ ~MTDDetLayerMeasurements()

Member Function Documentation

◆ checkEvent()

◆ checkMTDRecHits()

◆ fastMeasurements() [1/2]

◆ fastMeasurements() [2/2]

◆ groupedMeasurements() [1/2]

◆ groupedMeasurements() [2/2]

◆ measurements() [1/3]

◆ measurements() [2/3]

◆ measurements() [3/3]

◆ recHits() [1/3]

◆ recHits() [2/3]

◆ recHits() [3/3]

◆ setEvent()

Member Data Documentation

◆ mtdToken_

◆ theEvent

◆ theMTDEventCacheID

◆ theMTDRecHits