#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 30 of file MTDDetLayerMeasurements.h.

Member Typedef Documentation

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

Definition at line 33 of file MTDDetLayerMeasurements.h.

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

Definition at line 32 of file MTDDetLayerMeasurements.h.

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

Definition at line 34 of file MTDDetLayerMeasurements.h.

Constructor & Destructor Documentation

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

Definition at line 22 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 27 of file MTDDetLayerMeasurements.cc.

27 {}

Member Function Documentation

void MTDDetLayerMeasurements::checkEvent ( ) const

private

check that the event is set, and throw otherwise

Definition at line 194 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 50 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 122 of file MTDDetLayerMeasurements.cc.

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

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

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

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

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

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

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

                                                                                    {
   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 68 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 199 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 211 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 29 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 192 of file MTDDetLayerMeasurements.cc.

References event(), and theEvent.

192 { theEvent = &event; }

event

How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past

MTDDetLayerMeasurements::theEvent

const edm::Event * theEvent

Definition: MTDDetLayerMeasurements.h:112