MuonDetIdAssociator Class Reference

#include <MuonDetIdAssociator.h>

Inheritance diagram for MuonDetIdAssociator:

Public Member Functions
virtual const GeomDet *	getGeomDet (const DetId &id) const override
	MuonDetIdAssociator ()
	MuonDetIdAssociator (const int nPhi, const int nEta, const double etaBinSize)
	MuonDetIdAssociator (const edm::ParameterSet &pSet)
virtual const char *	name () const override
virtual void	setConditions (const DetIdAssociatorRecord &iRecord) override
virtual void	setCSCBadChambers (const CSCBadChambers *ptr)
virtual void	setGeometry (const GlobalTrackingGeometry *ptr)
virtual void	setGeometry (const DetIdAssociatorRecord &iRecord) override
Public Member Functions inherited from DetIdAssociator
virtual void	buildMap ()
	make the look-up map More...
	DetIdAssociator (const int nPhi, const int nEta, const double etaBinSize)
double	etaBinSize () const
	look-up map bin size in eta dimension More...
virtual std::vector< DetId >	getCrossedDetIds (const std::set< DetId > &, const std::vector< GlobalPoint > &trajectory) const
virtual std::vector< DetId >	getCrossedDetIds (const std::set< DetId > &, const std::vector< SteppingHelixStateInfo > &trajectory, const double toleranceInSigmas=-1) const
virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &, const int iN=0) const
virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &direction, const unsigned int iNEtaPlus, const unsigned int iNEtaMinus, const unsigned int iNPhiPlus, const unsigned int iNPhiMinus) const
virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &direction, const MapRange &mapRange) const
virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &point, const double d=0) const
virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &point, const double dThetaPlus, const double dThetaMinus, const double dPhiPlus, const double dPhiMinus) const
virtual std::set< DetId >	getDetIdsInACone (const std::set< DetId > &, const std::vector< GlobalPoint > &trajectory, const double dR) const
virtual int	iEta (const GlobalPoint &) const
	look-up map eta index More...
virtual int	iPhi (const GlobalPoint &) const
	look-up map phi index More...
int	nEtaBins () const
	number of bins of the look-up map in eta dimension More...
int	nPhiBins () const
	number of bins of the look-up map in phi dimension More...
virtual bool	selectAllInACone (const double dR) const
	helper to see if getDetIdsInACone is useful More...
const FiducialVolume &	volume () const
	get active detector volume More...
virtual	~DetIdAssociator ()

Protected Member Functions
virtual void	check_setup () const override
virtual std::pair< const_iterator, const_iterator >	getDetIdPoints (const DetId &id, std::vector< GlobalPoint > &points) const override
virtual GlobalPoint	getPosition (const DetId &id) const override
virtual void	getValidDetIds (unsigned int, std::vector< DetId > &) const override
virtual bool	insideElement (const GlobalPoint &point, const DetId &id) const override
Protected Member Functions inherited from DetIdAssociator
virtual bool	crossedElement (const GlobalPoint &, const GlobalPoint &, const DetId &, const double toleranceInSigmas=-1, const SteppingHelixStateInfo *=0) const
virtual void	dumpMapContent (int, int) const
virtual void	dumpMapContent (int, int, int, int) const
void	fillSet (std::set< DetId > &set, unsigned int iEta, unsigned int iPhi) const
virtual const unsigned int	getNumberOfSubdetectors () const
unsigned int	index (unsigned int iEta, unsigned int iPhi) const
virtual bool	nearElement (const GlobalPoint &point, const DetId &id, const double distance) const

Protected Attributes
const CSCBadChambers *	cscbadchambers_
const GlobalTrackingGeometry *	geometry_
bool	includeBadChambers_
bool	includeGEM_
bool	includeME0_
Protected Attributes inherited from DetIdAssociator
std::vector< DetId >	container_
const double	etaBinSize_
std::vector< std::pair< unsigned int, unsigned int > >	lookupMap_
double	maxEta_
double	minTheta_
const int	nEta_
const int	nPhi_
bool	theMapIsValid_
FiducialVolume	volume_

Additional Inherited Members
Public Types inherited from DetIdAssociator
typedef std::vector< GlobalPoint >::const_iterator	const_iterator
enum	PropagationTarget { Barrel, ForwardEndcap, BackwardEndcap }

Detailed Description

Definition at line 33 of file MuonDetIdAssociator.h.

Constructor & Destructor Documentation

MuonDetIdAssociator::MuonDetIdAssociator ( )

inline

Definition at line 35 of file MuonDetIdAssociator.h.

:DetIdAssociator(48, 48 , 0.125),geometry_(0),cscbadchambers_(0),includeBadChambers_(0){};

MuonDetIdAssociator::MuonDetIdAssociator ( const int  nPhi, const int  nEta, const double  etaBinSize  )

inline

Definition at line 36 of file MuonDetIdAssociator.h.

     :DetIdAssociator(nPhi, nEta, etaBinSize),geometry_(0),cscbadchambers_(0),includeBadChambers_(0){};

MuonDetIdAssociator::MuonDetIdAssociator ( const edm::ParameterSet &  pSet )

inline

Definition at line 39 of file MuonDetIdAssociator.h.

     :DetIdAssociator(pSet.getParameter<int>("nPhi"),pSet.getParameter<int>("nEta"),pSet.getParameter<double>("etaBinSize")),geometry_(0),cscbadchambers_(0),includeBadChambers_(pSet.getParameter<bool>("includeBadChambers")),includeGEM_(pSet.getParameter<bool>("includeGEM")),includeME0_(pSet.getParameter<bool>("includeME0")){};

Member Function Documentation

void MuonDetIdAssociator::check_setup ( ) const

overrideprotectedvirtual

Reimplemented from DetIdAssociator.

Definition at line 32 of file MuonDetIdAssociator.cc.

References DetIdAssociator::check_setup(), cscbadchambers_, Exception, and geometry_.

Referenced by name().

                                            {
   if (geometry_==0) throw cms::Exception("ConfigurationProblem") << "GlobalTrackingGeomtry is not set\n";
   if (cscbadchambers_==0) throw cms::Exception("ConfigurationProblem") << "CSCBadChambers is not set\n";
   DetIdAssociator::check_setup();
}

std::pair< DetIdAssociator::const_iterator, DetIdAssociator::const_iterator > MuonDetIdAssociator::getDetIdPoints ( const DetId &  id, std::vector< GlobalPoint > &  points  ) const

overrideprotectedvirtual

Implements DetIdAssociator.

Definition at line 113 of file MuonDetIdAssociator.cc.

References Surface::bounds(), geometry_, getGeomDet(), GlobalTrackingGeometry::idToDet(), Bounds::length(), GeomDet::surface(), Bounds::thickness(), GeomDet::toGlobal(), and Bounds::width().

Referenced by name().

{
   points.clear();
   points.reserve(8);
   const GeomDet* geomDet = getGeomDet( id );
   
   // the coners of muon detector elements are not stored and can be only calculated
   // based on methods defined in the interface class Bounds:
   //   width() - x
   //   length() - y 
   //   thinkness() - z
   // NOTE: this convention is implementation specific and can fail. Both
   //       RectangularPlaneBounds and TrapezoidalPlaneBounds use it.
   // Even though the CSC geomtry is more complicated (trapezoid),  it's enough 
   // to estimate which bins should contain this element. For the distance
   // calculation from the edge, we will use exact geometry to get it right.
        
   const Bounds* bounds = &(geometry_->idToDet(id)->surface().bounds());
   points.push_back(geomDet->toGlobal(LocalPoint(+bounds->width()/2,+bounds->length()/2,+bounds->thickness()/2)));
   points.push_back(geomDet->toGlobal(LocalPoint(-bounds->width()/2,+bounds->length()/2,+bounds->thickness()/2)));
   points.push_back(geomDet->toGlobal(LocalPoint(+bounds->width()/2,-bounds->length()/2,+bounds->thickness()/2)));
   points.push_back(geomDet->toGlobal(LocalPoint(-bounds->width()/2,-bounds->length()/2,+bounds->thickness()/2)));
   points.push_back(geomDet->toGlobal(LocalPoint(+bounds->width()/2,+bounds->length()/2,-bounds->thickness()/2)));
   points.push_back(geomDet->toGlobal(LocalPoint(-bounds->width()/2,+bounds->length()/2,-bounds->thickness()/2)));
   points.push_back(geomDet->toGlobal(LocalPoint(+bounds->width()/2,-bounds->length()/2,-bounds->thickness()/2)));
   points.push_back(geomDet->toGlobal(LocalPoint(-bounds->width()/2,-bounds->length()/2,-bounds->thickness()/2)));
   
   return std::pair<const_iterator,const_iterator>(points.begin(),points.end());
}

const GeomDet * MuonDetIdAssociator::getGeomDet ( const DetId &  id ) const

overridevirtual

Implements DetIdAssociator.

Definition at line 38 of file MuonDetIdAssociator.cc.

References Exception, geometry_, and GlobalTrackingGeometry::idToDet().

Referenced by getDetIdPoints(), getPosition(), and setConditions().

{
   if (geometry_==0) throw cms::Exception("ConfigurationProblem") << "GlobalTrackingGeomtry is not set\n";
   const GeomDet* gd = geometry_->idToDet(id);
   if (gd == 0) throw cms::Exception("NoGeometry") << "Cannot find GeomDet for DetID: " << id.rawId() <<"\n";
   return gd;
}

GlobalPoint MuonDetIdAssociator::getPosition ( const DetId &  id ) const

overrideprotectedvirtual

Implements DetIdAssociator.

Definition at line 47 of file MuonDetIdAssociator.cc.

References getGeomDet(), point, position, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by name().

                                                                  {
   Surface::PositionType point(getGeomDet(id)->surface().position());
   return GlobalPoint(point.x(),point.y(),point.z());
}

void MuonDetIdAssociator::getValidDetIds ( unsigned int  subDectorIndex, std::vector< DetId > &  validIds  ) const

overrideprotectedvirtual

Implements DetIdAssociator.

Definition at line 52 of file MuonDetIdAssociator.cc.

References cscbadchambers_, TrackingGeometry::dets(), dt, Exception, geometry_, includeBadChambers_, includeGEM_, includeME0_, CSCBadChambers::isInBadChamber(), alignCSCRings::r, and GlobalTrackingGeometry::slaveGeometry().

Referenced by name().

                                                                                                      {
  validIds.clear();
  if (geometry_==0) throw cms::Exception("ConfigurationProblem") << "GlobalTrackingGeomtry is not set\n";
  if (subDectorIndex!=0) throw cms::Exception("FatalError") << 
    "Muon sub-dectors are all handle as one sub-system, but subDetectorIndex is not zero.\n";

  // CSC 
  if (! geometry_->slaveGeometry(CSCDetId()) ) throw cms::Exception("FatalError") << "Cannnot CSCGeometry\n";
  auto const & geomDetsCSC = geometry_->slaveGeometry(CSCDetId())->dets();
  for(auto it = geomDetsCSC.begin(); it != geomDetsCSC.end(); ++it)
    if (auto csc = dynamic_cast<const CSCChamber*>(*it)) {
      if ((! includeBadChambers_) && (cscbadchambers_->isInBadChamber(CSCDetId(csc->id())))) continue;
      validIds.push_back(csc->id());
    }
  
  // DT
  if (! geometry_->slaveGeometry(DTChamberId()) ) throw cms::Exception("FatalError") << "Cannnot DTGeometry\n";
  auto const & geomDetsDT = geometry_->slaveGeometry(DTChamberId())->dets();
  for(auto it = geomDetsDT.begin(); it != geomDetsDT.end(); ++it)
    if (auto dt = dynamic_cast<const DTChamber*>(*it)) validIds.push_back(dt->id());

  // RPC
  if (! geometry_->slaveGeometry(RPCDetId()) ) throw cms::Exception("FatalError") << "Cannnot RPCGeometry\n";
  auto const & geomDetsRPC = geometry_->slaveGeometry(RPCDetId())->dets();
  for(auto it = geomDetsRPC.begin(); it != geomDetsRPC.end(); ++it)
    if (auto rpc = dynamic_cast<const RPCChamber*>(*it)) {
      std::vector< const RPCRoll*> rolls = (rpc->rolls());
      for(std::vector<const RPCRoll*>::iterator r = rolls.begin(); r != rolls.end(); ++r)
    validIds.push_back((*r)->id().rawId());
    }

  // GEM
  if (includeGEM_){
    if (! geometry_->slaveGeometry(GEMDetId()) ) throw cms::Exception("FatalError") << "Cannnot GEMGeometry\n";
    auto const & geomDetsGEM = geometry_->slaveGeometry(GEMDetId())->dets();
    for(auto it = geomDetsGEM.begin(); it != geomDetsGEM.end(); ++it){
      if (auto gem = dynamic_cast<const GEMSuperChamber*>(*it)) {
    validIds.push_back(gem->id());
      }
    }
  }
  // ME0
  if (includeME0_){
    if (! geometry_->slaveGeometry(ME0DetId()) ) throw cms::Exception("FatalError") << "Cannnot ME0Geometry\n";
    auto const & geomDetsME0 = geometry_->slaveGeometry(ME0DetId())->dets();
    for(auto it = geomDetsME0.begin(); it != geomDetsME0.end(); ++it){
      if (auto me0 = dynamic_cast<const ME0Chamber*>(*it)) {
    validIds.push_back(me0->id());
      }
    }
  }
  
}

bool MuonDetIdAssociator::insideElement ( const GlobalPoint &  point, const DetId &  id  ) const

overrideprotectedvirtual

Implements DetIdAssociator.

Definition at line 106 of file MuonDetIdAssociator.cc.

References Surface::bounds(), Exception, geometry_, GlobalTrackingGeometry::idToDet(), Bounds::inside(), GeomDet::surface(), and GeomDet::toLocal().

Referenced by name().

                                                                                       {
   if (geometry_==0) throw cms::Exception("ConfigurationProblem") << "GlobalTrackingGeomtry is not set\n";
   LocalPoint lp = geometry_->idToDet(id)->toLocal(point);
   return geometry_->idToDet(id)->surface().bounds().inside(lp);
}

virtual const char* MuonDetIdAssociator::name ( void  ) const

inlineoverridevirtual

Implements DetIdAssociator.

Definition at line 56 of file MuonDetIdAssociator.h.

References check_setup(), getDetIdPoints(), getPosition(), getValidDetIds(), insideElement(), and point.

Referenced by config.CFG::__str__(), and validation.Sample::digest().

{ return "AllMuonDetectors"; }

virtual void MuonDetIdAssociator::setConditions ( const DetIdAssociatorRecord &  iRecord )

inlineoverridevirtual

Reimplemented from DetIdAssociator.

Definition at line 48 of file MuonDetIdAssociator.h.

References getGeomDet(), edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::getRecord(), edm::ESHandle< T >::product(), and setCSCBadChambers().

                                                                            {
      edm::ESHandle<CSCBadChambers> cscbadchambersH;
      iRecord.getRecord<CSCBadChambersRcd>().get(cscbadchambersH);
      setCSCBadChambers(cscbadchambersH.product());
   };

virtual void MuonDetIdAssociator::setCSCBadChambers ( const CSCBadChambers *  ptr )

inlinevirtual

Definition at line 46 of file MuonDetIdAssociator.h.

References cscbadchambers_.

Referenced by setConditions().

{ cscbadchambers_ = ptr; }

virtual void MuonDetIdAssociator::setGeometry ( const GlobalTrackingGeometry *  ptr )

inlinevirtual

Definition at line 42 of file MuonDetIdAssociator.h.

References geometry_.

Referenced by setGeometry().

{ geometry_ = ptr; }

void MuonDetIdAssociator::setGeometry ( const DetIdAssociatorRecord &  iRecord )

overridevirtual

Implements DetIdAssociator.

Definition at line 143 of file MuonDetIdAssociator.cc.

References edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::getRecord(), edm::ESHandle< T >::product(), and setGeometry().

{
  edm::ESHandle<GlobalTrackingGeometry> geometryH;
  iRecord.getRecord<GlobalTrackingGeometryRecord>().get(geometryH);
  setGeometry(geometryH.product());
}

Member Data Documentation

const CSCBadChambers* MuonDetIdAssociator::cscbadchambers_

protected

Definition at line 72 of file MuonDetIdAssociator.h.

Referenced by check_setup(), getValidDetIds(), and setCSCBadChambers().

const GlobalTrackingGeometry* MuonDetIdAssociator::geometry_

protected

Definition at line 70 of file MuonDetIdAssociator.h.

Referenced by check_setup(), getDetIdPoints(), getGeomDet(), getValidDetIds(), insideElement(), and setGeometry().

bool MuonDetIdAssociator::includeBadChambers_

protected

Definition at line 73 of file MuonDetIdAssociator.h.

Referenced by getValidDetIds().

bool MuonDetIdAssociator::includeGEM_

protected

Definition at line 74 of file MuonDetIdAssociator.h.

Referenced by getValidDetIds().

bool MuonDetIdAssociator::includeME0_

protected

Definition at line 75 of file MuonDetIdAssociator.h.

Referenced by getValidDetIds().