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::ConsumesCollector &iC, bool enableDT=true, bool enableCSC=true, bool enableRPC=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	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	enableRPCMeasurement
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	theRPCEventCacheID
edm::Handle< RPCRecHitCollection >	theRPCRecHits

Detailed Description

The class to access recHits and TrajectoryMeasurements from DetLayer.

Author

C. Liu, R. Bellan, N. Amapane

Definition at line 39 of file MuonDetLayerMeasurements.h.

Member Typedef Documentation

typedef MuonTransientTrackingRecHit::MuonRecHitContainer MuonDetLayerMeasurements::MuonRecHitContainer

Definition at line 41 of file MuonDetLayerMeasurements.h.

Constructor & Destructor Documentation

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

Definition at line 24 of file MuonDetLayerMeasurements.cc.

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

                                                                         : 
  enableDTMeasurement(enableDT),
  enableCSCMeasurement(enableCSC),
  enableRPCMeasurement(enableRPC),
  theDTRecHits(),
  theCSCRecHits(),
  theRPCRecHits(),
  theDTEventCacheID(0),
  theCSCEventCacheID(0),
  theRPCEventCacheID(0),
  theEvent(0)
{

  dtToken_ = iC.consumes<DTRecSegment4DCollection>(dtlabel);
  cscToken_ = iC.consumes<CSCSegmentCollection>(csclabel);
  rpcToken_ = iC.consumes<RPCRecHitCollection>(rpclabel);


  static 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();
  procInstance++;
}

MuonDetLayerMeasurements::~MuonDetLayerMeasurements ( )

virtual

Definition at line 59 of file MuonDetLayerMeasurements.cc.

{}

Member Function Documentation

void MuonDetLayerMeasurements::checkCSCRecHits ( )

private

Definition at line 148 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 132 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 340 of file MuonDetLayerMeasurements.cc.

References Exception, and theEvent.

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

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

void MuonDetLayerMeasurements::checkRPCRecHits ( )

private

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

References MeasurementEstimator::estimate(), recHits(), query::result, and python.multivaluedict::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 278 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 299 of file MuonDetLayerMeasurements.cc.

References GeometricSearchDet::groupedCompatibleDets(), iEvent, measurements(), query::result, and python.multivaluedict::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 289 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 221 of file MuonDetLayerMeasurements.cc.

References MeasurementEstimator::estimate(), LogTrace, recHits(), query::result, and python.multivaluedict::sort().

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

                                            {
  MeasurementContainer result;
  
  // 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);
    LogTrace("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 194 of file MuonDetLayerMeasurements.cc.

References GeometricSearchDet::compatibleDets(), LogTrace, measurements(), query::result, and python.multivaluedict::sort().

                                               {
  
  MeasurementContainer result;
  
  std::vector<DetWithState> dss = layer->compatibleDets(startingState, prop, est);
  LogTrace("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 184 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 345 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 359 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 61 of file MuonDetLayerMeasurements.cc.

References checkCSCRecHits(), checkDTRecHits(), checkRPCRecHits(), MuonSubdetId::CSC, DetId::det(), MuonSubdetId::DT, enableCSCMeasurement, enableDTMeasurement, enableRPCMeasurement, Exception, GeomDet::geographicalId(), iEvent, DetId::rawId(), query::result, MuonSubdetId::RPC, MuonTransientTrackingRecHit::specificBuild(), DetId::subdetId(), theCSCRecHits, theDTRecHits, theEvent, 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());
      // LogTrace("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());
      //    LogTrace("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());
      // LogTrace("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 {
    // wrong type
    throw cms::Exception("MuonDetLayerMeasurements") << "The DetLayer with det " << geoId.det() << " subdet " << geoId.subdetId() << " is not a valid Muon DetLayer. ";
  }
  return result;
}

void MuonDetLayerMeasurements::setEvent

(

const edm::Event &

event

)

set event

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

Referenced by checkCSCRecHits(), and MuonDetLayerMeasurements().

edm::EDGetTokenT<DTRecSegment4DCollection> MuonDetLayerMeasurements::dtToken_

private

Definition at line 126 of file MuonDetLayerMeasurements.h.

Referenced by checkDTRecHits(), and MuonDetLayerMeasurements().

bool MuonDetLayerMeasurements::enableCSCMeasurement

private

Definition at line 132 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

bool MuonDetLayerMeasurements::enableDTMeasurement

private

Definition at line 131 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

bool MuonDetLayerMeasurements::enableRPCMeasurement

private

Definition at line 133 of file MuonDetLayerMeasurements.h.

Referenced by recHits().

edm::EDGetTokenT<RPCRecHitCollection> MuonDetLayerMeasurements::rpcToken_

private

Definition at line 128 of file MuonDetLayerMeasurements.h.

Referenced by checkRPCRecHits(), and MuonDetLayerMeasurements().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theCSCEventCacheID

private

Definition at line 146 of file MuonDetLayerMeasurements.h.

Referenced by checkCSCRecHits().

edm::Handle<CSCSegmentCollection> MuonDetLayerMeasurements::theCSCRecHits

private

Definition at line 137 of file MuonDetLayerMeasurements.h.

Referenced by checkCSCRecHits(), and recHits().

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theDTEventCacheID

private

Definition at line 145 of file MuonDetLayerMeasurements.h.

Referenced by checkDTRecHits().

edm::Handle<DTRecSegment4DCollection> MuonDetLayerMeasurements::theDTRecHits

private

Definition at line 136 of file MuonDetLayerMeasurements.h.

Referenced by checkDTRecHits(), and recHits().

const edm::Event* MuonDetLayerMeasurements::theEvent

private

Definition at line 149 of file MuonDetLayerMeasurements.h.

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

edm::Event::CacheIdentifier_t MuonDetLayerMeasurements::theRPCEventCacheID

private

Definition at line 147 of file MuonDetLayerMeasurements.h.

Referenced by checkRPCRecHits().

edm::Handle<RPCRecHitCollection> MuonDetLayerMeasurements::theRPCRecHits

private

Definition at line 138 of file MuonDetLayerMeasurements.h.

Referenced by checkRPCRecHits(), and recHits().