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 45 of file MuonDetLayerMeasurements.h.

Member Typedef Documentation

typedef MuonTransientTrackingRecHit::MuonRecHitContainer MuonDetLayerMeasurements::MuonRecHitContainer

Definition at line 47 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 29 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 78 of file MuonDetLayerMeasurements.cc.

{}

Member Function Documentation

void MuonDetLayerMeasurements::checkCSCRecHits ( )

private

Definition at line 224 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 208 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 460 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 257 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 273 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 241 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 372 of file MuonDetLayerMeasurements.cc.

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

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 398 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 419 of file MuonDetLayerMeasurements.cc.

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

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 409 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 331 of file MuonDetLayerMeasurements.cc.

References enableME0Measurement, MeasurementEstimator::estimate(), GeomDet::geographicalId(), LogDebug, MuonSubdetId::ME0, DetId::rawId(), recHits(), mps_fire::result, jetUpdater_cfi::sort, 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 303 of file MuonDetLayerMeasurements.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("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 293 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 465 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 479 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 80 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, 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 455 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 136 of file MuonDetLayerMeasurements.h.

Referenced by checkCSCRecHits(), and MuonDetLayerMeasurements().

edm::EDGetTokenT<DTRecSegment4DCollection> MuonDetLayerMeasurements::dtToken_

private

Definition at line 135 of file MuonDetLayerMeasurements.h.

Referenced by checkDTRecHits(), and MuonDetLayerMeasurements().

bool MuonDetLayerMeasurements::enableCSCMeasurement

private

Definition at line 143 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

bool MuonDetLayerMeasurements::enableDTMeasurement

private

Definition at line 142 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

bool MuonDetLayerMeasurements::enableGEMMeasurement

private

Definition at line 145 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

bool MuonDetLayerMeasurements::enableME0Measurement

private

Definition at line 146 of file MuonDetLayerMeasurements.h.

Referenced by measurements(), and recHits().

bool MuonDetLayerMeasurements::enableRPCMeasurement

private

Definition at line 144 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

edm::EDGetTokenT<GEMRecHitCollection> MuonDetLayerMeasurements::gemToken_

private

Definition at line 138 of file MuonDetLayerMeasurements.h.

Referenced by checkGEMRecHits(), and MuonDetLayerMeasurements().

edm::EDGetTokenT<ME0SegmentCollection> MuonDetLayerMeasurements::me0Token_

private

Definition at line 139 of file MuonDetLayerMeasurements.h.

Referenced by checkME0RecHits(), and MuonDetLayerMeasurements().

edm::EDGetTokenT<RPCRecHitCollection> MuonDetLayerMeasurements::rpcToken_

private

Definition at line 137 of file MuonDetLayerMeasurements.h.

Referenced by checkRPCRecHits(), and MuonDetLayerMeasurements().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theCSCEventCacheID

private

Definition at line 163 of file MuonDetLayerMeasurements.h.

Referenced by checkCSCRecHits().

edm::Handle<CSCSegmentCollection> MuonDetLayerMeasurements::theCSCRecHits

private

Definition at line 150 of file MuonDetLayerMeasurements.h.

Referenced by checkCSCRecHits(), and recHits().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theDTEventCacheID

private

Definition at line 162 of file MuonDetLayerMeasurements.h.

Referenced by checkDTRecHits().

edm::Handle<DTRecSegment4DCollection> MuonDetLayerMeasurements::theDTRecHits

private

Definition at line 149 of file MuonDetLayerMeasurements.h.

Referenced by checkDTRecHits(), and recHits().

const edm::Event* MuonDetLayerMeasurements::theEvent

private

Definition at line 168 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 165 of file MuonDetLayerMeasurements.h.

Referenced by checkGEMRecHits().

edm::Handle<GEMRecHitCollection> MuonDetLayerMeasurements::theGEMRecHits

private

Definition at line 152 of file MuonDetLayerMeasurements.h.

Referenced by checkGEMRecHits(), and recHits().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theME0EventCacheID

private

Definition at line 166 of file MuonDetLayerMeasurements.h.

Referenced by checkME0RecHits().

edm::Handle<ME0SegmentCollection> MuonDetLayerMeasurements::theME0RecHits

private

Definition at line 153 of file MuonDetLayerMeasurements.h.

Referenced by checkME0RecHits(), and recHits().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theRPCEventCacheID

private

Definition at line 164 of file MuonDetLayerMeasurements.h.

Referenced by checkRPCRecHits().

edm::Handle<RPCRecHitCollection> MuonDetLayerMeasurements::theRPCRecHits

private

Definition at line 151 of file MuonDetLayerMeasurements.h.

Referenced by checkRPCRecHits(), and recHits().