MTDDetLayerMeasurements Class Reference

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

{}

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(), phase1PixelTopology::layer, 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 (const 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, phase1PixelTopology::layer, recHits(), mps_fire::result, and jetUpdater_cfi::sort.

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(), phase1PixelTopology::layer, 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 iEvent, phase1PixelTopology::layer, measurements(), mps_fire::result, and jetUpdater_cfi::sort.

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, phase1PixelTopology::layer, LogDebug, mtdRecHits_cfi::mtdRecHits, recHits(), mps_fire::result, and jetUpdater_cfi::sort.

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(), phase1PixelTopology::layer, 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 iEvent, phase1PixelTopology::layer, LogDebug, measurements(), mps_fire::result, and jetUpdater_cfi::sort.

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(), phase1PixelTopology::layer, 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 iEvent, and phase1PixelTopology::layer.

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.

{ theEvent = &event; }

References edmPickEvents::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().