MuonDetLayerMeasurements Class Reference

#include <MuonDetLayerMeasurements.h>

Public Types
typedef MuonTransientTrackingRecHit::MuonRecHitContainer	MuonRecHitContainer

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...
	MuonDetLayerMeasurements (edm::InputTag dtlabel, edm::InputTag csclabel, edm::InputTag rpclabel, edm::InputTag gemlabel, edm::InputTag me0label, edm::ConsumesCollector &iC, bool enableDT=true, bool enableCSC=true, bool enableRPC=true, bool enableGEM=true, bool enableME0=true)
MuonRecHitContainer	recHits (const DetLayer *layer, const edm::Event &iEvent)
	returns the rechits which are on the layer More...
MuonRecHitContainer	recHits (const DetLayer *layer)
	returns the rechits which are on the layer More...
void	setEvent (const edm::Event &)
	set event More...
virtual	~MuonDetLayerMeasurements ()

Private Member Functions
void	checkCSCRecHits ()
void	checkDTRecHits ()
void	checkEvent () const
	check that the event is set, and throw otherwise More...
void	checkGEMRecHits ()
void	checkME0RecHits ()
void	checkRPCRecHits ()
MuonRecHitContainer	recHits (const GeomDet *, const edm::Event &iEvent)
	obtain TrackingRecHits from a DetLayer More...

Private Attributes
edm::EDGetTokenT< CSCSegmentCollection >	cscToken_
edm::EDGetTokenT< DTRecSegment4DCollection >	dtToken_
bool	enableCSCMeasurement
bool	enableDTMeasurement
bool	enableGEMMeasurement
bool	enableME0Measurement
bool	enableRPCMeasurement
edm::EDGetTokenT< GEMRecHitCollection >	gemToken_
edm::EDGetTokenT< ME0SegmentCollection >	me0Token_
edm::EDGetTokenT< RPCRecHitCollection >	rpcToken_
edm::Event::CacheIdentifier_t	theCSCEventCacheID
edm::Handle< CSCSegmentCollection >	theCSCRecHits
edm::Event::CacheIdentifier_t	theDTEventCacheID
edm::Handle< DTRecSegment4DCollection >	theDTRecHits
const edm::Event *	theEvent
edm::Event::CacheIdentifier_t	theGEMEventCacheID
edm::Handle< GEMRecHitCollection >	theGEMRecHits
edm::Event::CacheIdentifier_t	theME0EventCacheID
edm::Handle< ME0SegmentCollection >	theME0RecHits
edm::Event::CacheIdentifier_t	theRPCEventCacheID
edm::Handle< RPCRecHitCollection >	theRPCRecHits

Detailed Description

The class to access recHits and TrajectoryMeasurements from DetLayer.

Author

C. Liu, R. Bellan, N. Amapane by C. Calabria to include GEMs by D. Nash to include ME0s

by C. Calabria & A. Sharma to include GEMs

The class to access recHits and TrajectoryMeasurements from DetLayer.

Author

C. Liu, R. Bellan, N. Amapane by C. Calabria to include GEMs by D. Nash to include ME0s

by C. Calabria & A.Sharma to include GEMs

Definition at line 41 of file MuonDetLayerMeasurements.h.

Member Typedef Documentation

typedef MuonTransientTrackingRecHit::MuonRecHitContainer MuonDetLayerMeasurements::MuonRecHitContainer

Definition at line 43 of file MuonDetLayerMeasurements.h.

Constructor & Destructor Documentation

MuonDetLayerMeasurements::MuonDetLayerMeasurements	(	edm::InputTag	dtlabel,
		edm::InputTag	csclabel,
		edm::InputTag	rpclabel,
		edm::InputTag	gemlabel,
		edm::InputTag	me0label,
		edm::ConsumesCollector &	iC,
		bool	enableDT = true,
		bool	enableCSC = true,
		bool	enableRPC = true,
		bool	enableGEM = true,
		bool	enableME0 = true
	)

Definition at line 25 of file MuonDetLayerMeasurements.cc.

References edm::ConsumesCollector::consumes(), cscToken_, dtToken_, gemToken_, me0Token_, and rpcToken_.

    : enableDTMeasurement(enableDT),
      enableCSCMeasurement(enableCSC),
      enableRPCMeasurement(enableRPC),
      enableGEMMeasurement(enableGEM),
      enableME0Measurement(enableME0),
      theDTRecHits(),
      theCSCRecHits(),
      theRPCRecHits(),
      theGEMRecHits(),
      theDTEventCacheID(0),
      theCSCEventCacheID(0),
      theRPCEventCacheID(0),
      theGEMEventCacheID(0),
      theME0EventCacheID(0),
      theEvent(nullptr) {
  dtToken_ = iC.consumes<DTRecSegment4DCollection>(dtlabel);
  cscToken_ = iC.consumes<CSCSegmentCollection>(csclabel);
  rpcToken_ = iC.consumes<RPCRecHitCollection>(rpclabel);
  gemToken_ = iC.consumes<GEMRecHitCollection>(gemlabel);
  me0Token_ = iC.consumes<ME0SegmentCollection>(me0label);

  static std::atomic<int> procInstance{0};
  std::ostringstream sDT;
  sDT << "MuonDetLayerMeasurements::checkDTRecHits::" << procInstance;
  //  theDTCheckName = sDT.str();
  std::ostringstream sRPC;
  sRPC << "MuonDetLayerMeasurements::checkRPCRecHits::" << procInstance;
  //theRPCCheckName = sRPC.str();
  std::ostringstream sCSC;
  sCSC << "MuonDetLayerMeasurements::checkCSCRecHits::" << procInstance;
  //theCSCCheckName = sCSC.str();
  std::ostringstream sGEM;
  sGEM << "MuonDetLayerMeasurements::checkGEMRecHits::" << procInstance;
  //theGEMCheckName = sGEM.str();
  std::ostringstream sME0;
  sME0 << "MuonDetLayerMeasurements::checkME0RecHits::" << procInstance;
  //theME0CheckName = sME0.str();
  procInstance++;
}

MuonDetLayerMeasurements::~MuonDetLayerMeasurements ( )

virtual

Definition at line 76 of file MuonDetLayerMeasurements.cc.

{}

Member Function Documentation

void MuonDetLayerMeasurements::checkCSCRecHits ( )

private

Definition at line 207 of file MuonDetLayerMeasurements.cc.

References edm::Event::cacheIdentifier(), checkEvent(), cscToken_, Exception, edm::Event::getByToken(), edm::HandleBase::isValid(), theCSCEventCacheID, theCSCRecHits, and theEvent.

Referenced by recHits().

                                               {
  checkEvent();
  auto cacheID = theEvent->cacheIdentifier();
  if (cacheID == theCSCEventCacheID)
    return;

  {
    theEvent->getByToken(cscToken_, theCSCRecHits);
    theCSCEventCacheID = cacheID;
  }
  if (!theCSCRecHits.isValid()) {
    throw cms::Exception("MuonDetLayerMeasurements") << "Cannot get CSC RecHits";
  }
}

void MuonDetLayerMeasurements::checkDTRecHits ( )

private

Definition at line 195 of file MuonDetLayerMeasurements.cc.

References edm::Event::cacheIdentifier(), checkEvent(), dtToken_, Exception, edm::Event::getByToken(), edm::HandleBase::isValid(), theDTEventCacheID, theDTRecHits, and theEvent.

Referenced by recHits().

                                              {
  checkEvent();
  auto const cacheID = theEvent->cacheIdentifier();
  if (cacheID == theDTEventCacheID)
    return;

  { theEvent->getByToken(dtToken_, theDTRecHits); }
  if (!theDTRecHits.isValid()) {
    throw cms::Exception("MuonDetLayerMeasurements") << "Cannot get DT RecHits";
  }
}

void MuonDetLayerMeasurements::checkEvent ( ) const

private

check that the event is set, and throw otherwise

Definition at line 412 of file MuonDetLayerMeasurements.cc.

References Exception, and theEvent.

Referenced by checkCSCRecHits(), checkDTRecHits(), checkGEMRecHits(), checkME0RecHits(), checkRPCRecHits(), fastMeasurements(), groupedMeasurements(), measurements(), and recHits().

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

void MuonDetLayerMeasurements::checkGEMRecHits ( )

private

Definition at line 237 of file MuonDetLayerMeasurements.cc.

References edm::Event::cacheIdentifier(), checkEvent(), Exception, gemToken_, edm::Event::getByToken(), edm::HandleBase::isValid(), theEvent, theGEMEventCacheID, and theGEMRecHits.

Referenced by recHits().

                                               {
  checkEvent();
  auto cacheID = theEvent->cacheIdentifier();
  if (cacheID == theGEMEventCacheID)
    return;

  {
    theEvent->getByToken(gemToken_, theGEMRecHits);
    theGEMEventCacheID = cacheID;
  }
  if (!theGEMRecHits.isValid()) {
    throw cms::Exception("MuonDetLayerMeasurements") << "Cannot get GEM RecHits";
  }
}

void MuonDetLayerMeasurements::checkME0RecHits ( )

private

Definition at line 252 of file MuonDetLayerMeasurements.cc.

References edm::Event::cacheIdentifier(), checkEvent(), Exception, edm::Event::getByToken(), edm::HandleBase::isValid(), LogDebug, me0Token_, theEvent, theME0EventCacheID, and theME0RecHits.

Referenced by recHits().

                                               {
  LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "Checking ME0 RecHits";
  checkEvent();
  auto cacheID = theEvent->cacheIdentifier();
  if (cacheID == theME0EventCacheID)
    return;

  {
    theEvent->getByToken(me0Token_, theME0RecHits);
    theME0EventCacheID = cacheID;
  }
  if (!theME0RecHits.isValid()) {
    throw cms::Exception("MuonDetLayerMeasurements") << "Cannot get ME0 RecHits";
    LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "Cannot get ME0 RecHits";
  }
}

void MuonDetLayerMeasurements::checkRPCRecHits ( )

private

Definition at line 222 of file MuonDetLayerMeasurements.cc.

References edm::Event::cacheIdentifier(), checkEvent(), Exception, edm::Event::getByToken(), edm::HandleBase::isValid(), rpcToken_, theEvent, theRPCEventCacheID, and theRPCRecHits.

Referenced by recHits().

                                               {
  checkEvent();
  auto cacheID = theEvent->cacheIdentifier();
  if (cacheID == theRPCEventCacheID)
    return;

  {
    theEvent->getByToken(rpcToken_, theRPCRecHits);
    theRPCEventCacheID = cacheID;
  }
  if (!theRPCRecHits.isValid()) {
    throw cms::Exception("MuonDetLayerMeasurements") << "Cannot get RPC RecHits";
  }
}

MeasurementContainer MuonDetLayerMeasurements::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 336 of file MuonDetLayerMeasurements.cc.

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

Referenced by fastMeasurements().

                                                                                        {
  MeasurementContainer result;
  MuonRecHitContainer rhs = recHits(layer, iEvent);
  for (MuonRecHitContainer::const_iterator irh = rhs.begin(); irh != rhs.end(); irh++) {
    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 MuonDetLayerMeasurements::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 359 of file MuonDetLayerMeasurements.cc.

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

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

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

Definition at line 377 of file MuonDetLayerMeasurements.cc.

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

Referenced by CosmicMuonTrajectoryBuilder::findBestMeasurements(), StandAloneMuonFilter::findBestMeasurements(), and 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));

  // this should be fixed either in RecoMuon/MeasurementDet/MuonDetLayerMeasurements or
  // RecoMuon/DetLayers/MuRingForwardDoubleLayer
  // and removed the reverse operation in StandAloneMuonFilter::findBestMeasurements

  for (std::vector<DetGroup>::const_iterator grp = groups.begin(); grp != groups.end(); ++grp) {
    std::vector<TrajectoryMeasurement> groupMeasurements;
    for (DetGroup::const_iterator detAndStateItr = grp->begin(); detAndStateItr != grp->end(); ++detAndStateItr) {
      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 > MuonDetLayerMeasurements::groupedMeasurements	(	const DetLayer *	layer,
		const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est
	)

Definition at line 368 of file MuonDetLayerMeasurements.cc.

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

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

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

Definition at line 301 of file MuonDetLayerMeasurements.cc.

References enableME0Measurement, MeasurementEstimator::estimate(), GeomDet::geographicalId(), LogDebug, MuonSubdetId::ME0, DetId::rawId(), recHits(), mps_fire::result, and DetId::subdetId().

Referenced by CosmicMuonTrajectoryBuilder::findBestMeasurements(), StandAloneMuonFilter::findBestMeasurements(), groupedMeasurements(), and measurements().

                                                                                    {
  MeasurementContainer result;

  DetId geoId = det->geographicalId();
  //no chamber here that is actually an me0....
  if (geoId.subdetId() == MuonSubdetId::ME0) {
    if (enableME0Measurement) {
      ME0DetId chamberId(geoId.rawId());
      LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements")
          << "ME0 Chamber ID in measurements: " << chamberId << std::endl;
    }
  }

  // Get the Segments which relies on the GeomDet given by compatibleDets
  MuonRecHitContainer muonRecHits = recHits(det, iEvent);

  // Create the Trajectory Measurement
  for (MuonRecHitContainer::const_iterator rechit = muonRecHits.begin(); rechit != muonRecHits.end(); ++rechit) {
    MeasurementEstimator::HitReturnType estimate = est.estimate(stateOnDet, **rechit);
    LogDebug("RecoMuon") << "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 MuonDetLayerMeasurements::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 278 of file MuonDetLayerMeasurements.cc.

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

                                                                                    {
  MeasurementContainer result;

  std::vector<DetWithState> dss = layer->compatibleDets(startingState, prop, est);
  LogDebug("RecoMuon") << "compatibleDets: " << dss.size() << std::endl;

  for (std::vector<DetWithState>::const_iterator detWithStateItr = dss.begin(); detWithStateItr != dss.end();
       ++detWithStateItr) {
    MeasurementContainer detMeasurements =
        measurements(layer, detWithStateItr->first, detWithStateItr->second, est, iEvent);
    result.insert(result.end(), detMeasurements.begin(), detMeasurements.end());
  }

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

  return result;
}

MeasurementContainer MuonDetLayerMeasurements::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 270 of file MuonDetLayerMeasurements.cc.

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

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

MuonRecHitContainer MuonDetLayerMeasurements::recHits	(	const DetLayer *	layer,
		const edm::Event &	iEvent
	)

returns the rechits which are on the layer

Definition at line 417 of file MuonDetLayerMeasurements.cc.

References GeometricSearchDet::basicComponents().

Referenced by MuonSeedBuilder::build(), fastMeasurements(), measurements(), MuonSeedOrcaPatternRecognition::produce(), recHits(), CosmicMuonTrajectoryBuilder::trajectories(), and CosmicMuonTrajectoryBuilder::unusedHits().

                                                                                                   {
  MuonRecHitContainer rhs;

  std::vector<const GeomDet*> gds = layer->basicComponents();

  for (std::vector<const GeomDet*>::const_iterator igd = gds.begin(); igd != gds.end(); igd++) {
    MuonRecHitContainer detHits = recHits(*igd, iEvent);
    rhs.insert(rhs.end(), detHits.begin(), detHits.end());
  }
  return rhs;
}

MuonRecHitContainer MuonDetLayerMeasurements::recHits

(

const DetLayer *

layer

)

returns the rechits which are on the layer

Definition at line 429 of file MuonDetLayerMeasurements.cc.

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

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

MuonRecHitContainer MuonDetLayerMeasurements::recHits ( const GeomDet *  geomDet, const edm::Event &  iEvent  )

private

obtain TrackingRecHits from a DetLayer

Definition at line 78 of file MuonDetLayerMeasurements.cc.

References checkCSCRecHits(), checkDTRecHits(), checkGEMRecHits(), checkME0RecHits(), checkRPCRecHits(), MuonSubdetId::CSC, DetId::det(), MuonSubdetId::DT, enableCSCMeasurement, enableDTMeasurement, enableGEMMeasurement, enableME0Measurement, enableRPCMeasurement, Exception, MuonSubdetId::GEM, GeomDet::geographicalId(), iEvent, LogDebug, MuonSubdetId::ME0, FastTimerService_cff::range, DetId::rawId(), mps_fire::result, MuonSubdetId::RPC, MuonTransientTrackingRecHit::specificBuild(), DetId::subdetId(), theCSCRecHits, theDTRecHits, theEvent, theGEMRecHits, theME0RecHits, and theRPCRecHits.

                                                                                                    {
  DetId geoId = geomDet->geographicalId();
  theEvent = &iEvent;
  MuonRecHitContainer result;

  if (geoId.subdetId() == MuonSubdetId::DT) {
    if (enableDTMeasurement) {
      checkDTRecHits();

      // Create the ChamberId
      DTChamberId chamberId(geoId.rawId());
      LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "(DT): " << chamberId << std::endl;

      // Get the DT-Segment which relies on this chamber
      DTRecSegment4DCollection::range range = theDTRecHits->get(chamberId);

      // Create the MuonTransientTrackingRechit
      for (DTRecSegment4DCollection::const_iterator rechit = range.first; rechit != range.second; ++rechit)
        result.push_back(MuonTransientTrackingRecHit::specificBuild(geomDet, &*rechit));
    }
  }

  else if (geoId.subdetId() == MuonSubdetId::CSC) {
    if (enableCSCMeasurement) {
      checkCSCRecHits();

      // Create the chamber Id
      CSCDetId chamberId(geoId.rawId());
      LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "(CSC): " << chamberId << std::endl;

      // Get the CSC-Segment which relies on this chamber
      CSCSegmentCollection::range range = theCSCRecHits->get(chamberId);

      // Create the MuonTransientTrackingRecHit
      for (CSCSegmentCollection::const_iterator rechit = range.first; rechit != range.second; ++rechit)
        result.push_back(MuonTransientTrackingRecHit::specificBuild(geomDet, &*rechit));
    }
  }

  else if (geoId.subdetId() == MuonSubdetId::RPC) {
    if (enableRPCMeasurement) {
      checkRPCRecHits();

      // Create the chamber Id
      RPCDetId chamberId(geoId.rawId());
      LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "(RPC): " << chamberId << std::endl;

      // Get the RPC-Segment which relies on this chamber
      RPCRecHitCollection::range range = theRPCRecHits->get(chamberId);

      // Create the MuonTransientTrackingRecHit
      for (RPCRecHitCollection::const_iterator rechit = range.first; rechit != range.second; ++rechit)
        result.push_back(MuonTransientTrackingRecHit::specificBuild(geomDet, &*rechit));
    }
  } else if (geoId.subdetId() == MuonSubdetId::GEM) {
    if (enableGEMMeasurement) {
      checkGEMRecHits();

      // Create the chamber Id
      GEMDetId chamberId(geoId.rawId());

      LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "(GEM): " << chamberId << std::endl;

      // Get the GEM-Segment which relies on this chamber
      GEMRecHitCollection::range range = theGEMRecHits->get(chamberId);

      LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements")
          << "Number of GEM rechits available =  " << theGEMRecHits->size() << ", from chamber: " << chamberId
          << std::endl;

      // Create the MuonTransientTrackingRecHit
      for (GEMRecHitCollection::const_iterator rechit = range.first; rechit != range.second; ++rechit)
        result.push_back(MuonTransientTrackingRecHit::specificBuild(geomDet, &*rechit));
      LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "Number of GEM rechits = " << result.size() << std::endl;
    }
  }

  else if (geoId.subdetId() == MuonSubdetId::ME0) {
    LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "(ME0): identified" << std::endl;
    if (enableME0Measurement) {
      checkME0RecHits();

      // Create the chamber Id
      ME0DetId chamberId(geoId.rawId());

      // Get the ME0-Segment which relies on this chamber
      // Getting rechits right now, not segments - maybe it should be segments?
      ME0SegmentCollection::range range = theME0RecHits->get(chamberId);

      LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements")
          << "Number of ME0 rechits available =  " << theME0RecHits->size() << ", from chamber: " << chamberId
          << std::endl;

      // Create the MuonTransientTrackingRecHit
      for (ME0SegmentCollection::const_iterator rechit = range.first; rechit != range.second; ++rechit) {
        LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "On ME0 iteration " << std::endl;
        result.push_back(MuonTransientTrackingRecHit::specificBuild(geomDet, &*rechit));
      }
      LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "Number of ME0 rechits = " << result.size() << std::endl;
    }
  } else {
    // wrong type
    throw cms::Exception("MuonDetLayerMeasurements") << "The DetLayer with det " << geoId.det() << " subdet "
                                                     << geoId.subdetId() << " is not a valid Muon DetLayer. ";
    LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "The DetLayer with det " << geoId.det() << " subdet "
                                                       << geoId.subdetId() << " is not a valid Muon DetLayer. ";
  }
  if (enableME0Measurement) {
    LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "(ME0): enabled" << std::endl;
  }

  if (enableGEMMeasurement) {
    LogDebug("Muon|RecoMuon|MuonDetLayerMeasurements") << "(GEM): enabled" << std::endl;
  }
  return result;
}

void MuonDetLayerMeasurements::setEvent

(

const edm::Event &

event

)

set event

Definition at line 410 of file MuonDetLayerMeasurements.cc.

References event(), and theEvent.

Referenced by CosmicMuonTrajectoryBuilder::setEvent(), and StandAloneMuonFilter::setEvent().

{ theEvent = &event; }

Member Data Documentation

edm::EDGetTokenT<CSCSegmentCollection> MuonDetLayerMeasurements::cscToken_

private

Definition at line 122 of file MuonDetLayerMeasurements.h.

Referenced by checkCSCRecHits(), and MuonDetLayerMeasurements().

edm::EDGetTokenT<DTRecSegment4DCollection> MuonDetLayerMeasurements::dtToken_

private

Definition at line 121 of file MuonDetLayerMeasurements.h.

Referenced by checkDTRecHits(), and MuonDetLayerMeasurements().

bool MuonDetLayerMeasurements::enableCSCMeasurement

private

Definition at line 128 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

bool MuonDetLayerMeasurements::enableDTMeasurement

private

Definition at line 127 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

bool MuonDetLayerMeasurements::enableGEMMeasurement

private

Definition at line 130 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

bool MuonDetLayerMeasurements::enableME0Measurement

private

Definition at line 131 of file MuonDetLayerMeasurements.h.

Referenced by measurements(), and recHits().

bool MuonDetLayerMeasurements::enableRPCMeasurement

private

Definition at line 129 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

edm::EDGetTokenT<GEMRecHitCollection> MuonDetLayerMeasurements::gemToken_

private

Definition at line 124 of file MuonDetLayerMeasurements.h.

Referenced by checkGEMRecHits(), and MuonDetLayerMeasurements().

edm::EDGetTokenT<ME0SegmentCollection> MuonDetLayerMeasurements::me0Token_

private

Definition at line 125 of file MuonDetLayerMeasurements.h.

Referenced by checkME0RecHits(), and MuonDetLayerMeasurements().

edm::EDGetTokenT<RPCRecHitCollection> MuonDetLayerMeasurements::rpcToken_

private

Definition at line 123 of file MuonDetLayerMeasurements.h.

Referenced by checkRPCRecHits(), and MuonDetLayerMeasurements().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theCSCEventCacheID

private

Definition at line 148 of file MuonDetLayerMeasurements.h.

Referenced by checkCSCRecHits().

edm::Handle<CSCSegmentCollection> MuonDetLayerMeasurements::theCSCRecHits

private

Definition at line 135 of file MuonDetLayerMeasurements.h.

Referenced by checkCSCRecHits(), and recHits().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theDTEventCacheID

private

Definition at line 147 of file MuonDetLayerMeasurements.h.

Referenced by checkDTRecHits().

edm::Handle<DTRecSegment4DCollection> MuonDetLayerMeasurements::theDTRecHits

private

Definition at line 134 of file MuonDetLayerMeasurements.h.

Referenced by checkDTRecHits(), and recHits().

const edm::Event* MuonDetLayerMeasurements::theEvent

private

Definition at line 153 of file MuonDetLayerMeasurements.h.

Referenced by checkCSCRecHits(), checkDTRecHits(), checkEvent(), checkGEMRecHits(), checkME0RecHits(), checkRPCRecHits(), fastMeasurements(), groupedMeasurements(), measurements(), recHits(), and setEvent().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theGEMEventCacheID

private

Definition at line 150 of file MuonDetLayerMeasurements.h.

Referenced by checkGEMRecHits().

edm::Handle<GEMRecHitCollection> MuonDetLayerMeasurements::theGEMRecHits

private

Definition at line 137 of file MuonDetLayerMeasurements.h.

Referenced by checkGEMRecHits(), and recHits().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theME0EventCacheID

private

Definition at line 151 of file MuonDetLayerMeasurements.h.

Referenced by checkME0RecHits().

edm::Handle<ME0SegmentCollection> MuonDetLayerMeasurements::theME0RecHits

private

Definition at line 138 of file MuonDetLayerMeasurements.h.

Referenced by checkME0RecHits(), and recHits().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theRPCEventCacheID

private

Definition at line 149 of file MuonDetLayerMeasurements.h.

Referenced by checkRPCRecHits().

edm::Handle<RPCRecHitCollection> MuonDetLayerMeasurements::theRPCRecHits

private

Definition at line 136 of file MuonDetLayerMeasurements.h.

Referenced by checkRPCRecHits(), and recHits().