#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, const edm::Event &iEvent)
	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)
	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, const edm::Event &iEvent)

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

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, const edm::Event &iEvent)
	returns TMeasurements in a DetLayer compatible with the TSOS. More...

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

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

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

MTDRecHitContainer	recHits (const DetLayer *layer)
	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...

template<class T >
T	sortResult (T &)

Private Attributes
const edm::Event *	theEvent

edm::Event::CacheIdentifier_t	theMTDEventCacheID

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

edm::EDGetTokenT< MTDTrackingRecHit >	theMTDToken

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	(	const edm::InputTag &	mtdlabel,
		edm::ConsumesCollector &	iC
	)

Definition at line 22 of file MTDDetLayerMeasurements.cc.

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

◆ ~MTDDetLayerMeasurements()

MTDDetLayerMeasurements::~MTDDetLayerMeasurements ( )

virtual

Definition at line 28 of file MTDDetLayerMeasurements.cc.

28 {}

Member Function Documentation

◆ checkEvent()

void MTDDetLayerMeasurements::checkEvent ( ) const

private

check that the event is set, and throw otherwise

Definition at line 189 of file MTDDetLayerMeasurements.cc.

References Exception, and theEvent.

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

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

◆ checkMTDRecHits()

void MTDDetLayerMeasurements::checkMTDRecHits ( )

private

Definition at line 50 of file MTDDetLayerMeasurements.cc.

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

Referenced by recHits().

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

◆ fastMeasurements() [1/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 124 of file MTDDetLayerMeasurements.cc.

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

Referenced by fastMeasurements().

                                                                                        {
   MeasurementContainer result;
   MTDRecHitContainer rhs = recHits(layer, iEvent);
   for (const auto& irh : rhs) {
     MeasurementEstimator::HitReturnType estimate = est.estimate(theStateOnDet, (*irh));
     if (estimate.first) {
       result.push_back(TrajectoryMeasurement(theStateOnDet, irh, estimate.second, layer));
     }
   }
 
   return sortResult(result);
 }

◆ fastMeasurements() [2/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 143 of file MTDDetLayerMeasurements.cc.

References checkEvent(), fastMeasurements(), pixelTopology::layer, and theEvent.

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

◆ groupedMeasurements() [1/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 161 of file MTDDetLayerMeasurements.cc.

References iEvent, pixelTopology::layer, measurements(), mps_fire::result, and sortResult().

Referenced by groupedMeasurements().

                             {
   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());
     }
 
     result.push_back(TrajectoryMeasurementGroup(sortResult(groupMeasurements), grp));
   }
 
   return result;
 }

◆ groupedMeasurements() [2/2]

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

Definition at line 152 of file MTDDetLayerMeasurements.cc.

References checkEvent(), groupedMeasurements(), pixelTopology::layer, and theEvent.

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

◆ 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 102 of file MTDDetLayerMeasurements.cc.

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

Referenced by groupedMeasurements(), and measurements().

                                                                                    {
   MeasurementContainer result;
 
   // Get the Segments which relies on the GeomDet given by compatibleDets
   MTDRecHitContainer mtdRecHits = recHits(det, iEvent);
 
   // Create the Trajectory Measurement
   for (const auto& rechit : mtdRecHits) {
     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));
     }
   }
 
   return sortResult(result);
 }

◆ measurements() [2/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 84 of file MTDDetLayerMeasurements.cc.

References iEvent, pixelTopology::layer, LogDebug, measurements(), mps_fire::result, and sortResult().

                                                                                    {
   MeasurementContainer result;
 
   const auto& dss = layer->compatibleDets(startingState, prop, est);
   LogDebug("MTDDetLayerMeasurements") << "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());
   }
 
   return sortResult(result);
 }

◆ measurements() [3/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 76 of file MTDDetLayerMeasurements.cc.

References checkEvent(), pixelTopology::layer, measurements(), and theEvent.

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

◆ recHits() [1/3]

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

returns the rechits which are on the layer

Definition at line 194 of file MTDDetLayerMeasurements.cc.

References iEvent, and pixelTopology::layer.

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

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

◆ recHits() [2/3]

MTDRecHitContainer MTDDetLayerMeasurements::recHits ( const DetLayer * layer )

returns the rechits which are on the layer

Definition at line 206 of file MTDDetLayerMeasurements.cc.

References checkEvent(), pixelTopology::layer, recHits(), and theEvent.

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

◆ recHits() [3/3]

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

private

obtain TrackingRecHits from a DetLayer

Definition at line 30 of file MTDDetLayerMeasurements.cc.

References GenericTransientTrackingRecHit::build(), checkMTDRecHits(), GeomDet::geographicalId(), totem::nt2::vfat::geoId(), iEvent, LogDebug, mps_fire::result, and theEvent.

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

◆ setEvent()

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

set event

Definition at line 187 of file MTDDetLayerMeasurements.cc.

References edmPickEvents::event, and theEvent.

187 { theEvent = &event; }

edmPickEvents.event

event

Definition: edmPickEvents.py:280

MTDDetLayerMeasurements::theEvent

const edm::Event * theEvent

Definition: MTDDetLayerMeasurements.h:116

◆ sortResult()

template<class T >

T MTDDetLayerMeasurements::sortResult ( T & result )

private

Definition at line 67 of file MTDDetLayerMeasurements.cc.

References mps_fire::result, and jetUpdater_cfi::sort.

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

                                                {
   if (!result.empty()) {
     sort(result.begin(), result.end(), TrajMeasLessEstim());
   }
 
   return result;
 }

Member Data Documentation

◆ theEvent

const edm::Event* MTDDetLayerMeasurements::theEvent

private

Definition at line 116 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 114 of file MTDDetLayerMeasurements.h.

Referenced by checkMTDRecHits().

◆ theMTDRecHits

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

private

Definition at line 109 of file MTDDetLayerMeasurements.h.

Referenced by checkMTDRecHits().

◆ theMTDToken

edm::EDGetTokenT<MTDTrackingRecHit> MTDDetLayerMeasurements::theMTDToken

private

Definition at line 106 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()

◆ sortResult()

Member Data Documentation

◆ theEvent

◆ theMTDEventCacheID

◆ theMTDRecHits

◆ theMTDToken