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

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.

Definition at line 32 of file MTDDetLayerMeasurements.h.

Member Typedef Documentation

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

Definition at line 35 of file MTDDetLayerMeasurements.h.

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

Definition at line 34 of file MTDDetLayerMeasurements.h.

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

Definition at line 36 of file MTDDetLayerMeasurements.h.

Constructor & Destructor Documentation

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

Definition at line 23 of file MTDDetLayerMeasurements.cc.

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

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

MTDDetLayerMeasurements::~MTDDetLayerMeasurements ( )

virtual

Definition at line 32 of file MTDDetLayerMeasurements.cc.

32 {}

Member Function Documentation

void MTDDetLayerMeasurements::checkEvent ( ) const

private

check that the event is set, and throw otherwise

Definition at line 226 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";
 }

void MTDDetLayerMeasurements::checkMTDRecHits ( )

private

Definition at line 57 of file MTDDetLayerMeasurements.cc.

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

Referenced by recHits().

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

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

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

Referenced by fastMeasurements().

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

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

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

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

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

Definition at line 190 of file MTDDetLayerMeasurements.cc.

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

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());
     }
     
     if (!groupMeasurements.empty()) 
       std::sort( groupMeasurements.begin(), groupMeasurements.end(), TrajMeasLessEstim());  
     
     result.push_back(TrajectoryMeasurementGroup(groupMeasurements, grp));
   }
 
   return result;
 }

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

Definition at line 180 of file MTDDetLayerMeasurements.cc.

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

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

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

Definition at line 112 of file MTDDetLayerMeasurements.cc.

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

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

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

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

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

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

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

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

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

returns the rechits which are on the layer

Definition at line 231 of file MTDDetLayerMeasurements.cc.

References GeometricSearchDet::basicComponents().

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

MTDRecHitContainer MTDDetLayerMeasurements::recHits ( const DetLayer * layer )

returns the rechits which are on the layer

Definition at line 244 of file MTDDetLayerMeasurements.cc.

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

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

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

private

obtain TrackingRecHits from a DetLayer

Definition at line 34 of file MTDDetLayerMeasurements.cc.

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

 {
   DetId geoId = geomDet->geographicalId();
   theEvent = &iEvent;
   MTDRecHitContainer result;
 
   checkMTDRecHits();
   
   // Create the ChamberId
   DetId detId(geoId.rawId());
   LogDebug("Track|RecoMTD|MTDDetLayerMeasurements") << "(MTD): "<<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;
 }

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

set event

Definition at line 221 of file MTDDetLayerMeasurements.cc.

References event(), and theEvent.

                                                             {
   theEvent = &event;
 }

Member Data Documentation

edm::EDGetTokenT<MTDTrackingRecHit> MTDDetLayerMeasurements::mtdToken_

private

Definition at line 115 of file MTDDetLayerMeasurements.h.

Referenced by checkMTDRecHits(), and MTDDetLayerMeasurements().

const edm::Event* MTDDetLayerMeasurements::theEvent

private

Definition at line 125 of file MTDDetLayerMeasurements.h.

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

edm::Event::CacheIdentifier_t MTDDetLayerMeasurements::theMTDEventCacheID

private

Definition at line 123 of file MTDDetLayerMeasurements.h.

Referenced by checkMTDRecHits().

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

private

Definition at line 118 of file MTDDetLayerMeasurements.h.

Referenced by checkMTDRecHits().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation