#include <MuonDetIdAssociator.h>

Inheritance diagram for MuonDetIdAssociator:

Public Member Functions
const GeomDet *	getGeomDet (const DetId &id) const override

	MuonDetIdAssociator ()

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

	MuonDetIdAssociator (int nPhi, int nEta, double etaBinSize, const GlobalTrackingGeometry geom, const CSCBadChambers badChambers, bool includeBadChambers, bool includeGEM, bool includeME0)

const char *	name () const override

virtual void	setCSCBadChambers (const CSCBadChambers *ptr)

virtual void	setGeometry (const GlobalTrackingGeometry *ptr)

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
void	check_setup () const override

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

GlobalPoint	getPosition (const DetId &id) const override

void	getValidDetIds (unsigned int, std::vector< DetId > &) const override

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 *=nullptr) 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 29 of file MuonDetIdAssociator.h.

Constructor & Destructor Documentation

◆ MuonDetIdAssociator() [1/3]

MuonDetIdAssociator::MuonDetIdAssociator ( )

inline

Definition at line 31 of file MuonDetIdAssociator.h.

32 : DetIdAssociator(48, 48, 0.125), geometry_(nullptr), cscbadchambers_(nullptr), includeBadChambers_(false){};

MuonDetIdAssociator::cscbadchambers_

const CSCBadChambers * cscbadchambers_

Definition: MuonDetIdAssociator.h:76

MuonDetIdAssociator::includeBadChambers_

bool includeBadChambers_

Definition: MuonDetIdAssociator.h:77

DetIdAssociator::DetIdAssociator

DetIdAssociator(const int nPhi, const int nEta, const double etaBinSize)

Definition: DetIdAssociator.cc:24

MuonDetIdAssociator::geometry_

const GlobalTrackingGeometry * geometry_

Definition: MuonDetIdAssociator.h:74

◆ MuonDetIdAssociator() [2/3]

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

inline

Definition at line 33 of file MuonDetIdAssociator.h.

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

◆ MuonDetIdAssociator() [3/3]

MuonDetIdAssociator::MuonDetIdAssociator	(	int	nPhi,
		int	nEta,
		double	etaBinSize,
		const GlobalTrackingGeometry *	geom,
		const CSCBadChambers *	badChambers,
		bool	includeBadChambers,
		bool	includeGEM,
		bool	includeME0
	)

inline

Definition at line 39 of file MuonDetIdAssociator.h.

       : DetIdAssociator(nPhi, nEta, etaBinSize),
         geometry_(geom),
         cscbadchambers_(badChambers),
         includeBadChambers_(includeBadChambers),
         includeGEM_(includeGEM),
         includeME0_(includeME0){};

Member Function Documentation

◆ check_setup()

void MuonDetIdAssociator::check_setup ( ) const

overrideprotectedvirtual

Reimplemented from DetIdAssociator.

Definition at line 31 of file MuonDetIdAssociator.cc.

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

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

◆ getDetIdPoints()

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

overrideprotectedvirtual

Implements DetIdAssociator.

Definition at line 126 of file MuonDetIdAssociator.cc.

References Surface::bounds(), geometry_, getGeomDet(), GlobalTrackingGeometry::idToDet(), hiPixelPairStep_cff::points, GeomDet::surface(), and GeomDet::toGlobal().

                                                            {
   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());
 }

◆ getGeomDet()

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

overridevirtual

Implements DetIdAssociator.

Definition at line 39 of file MuonDetIdAssociator.cc.

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

Referenced by getDetIdPoints(), and getPosition().

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

◆ getPosition()

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

overrideprotectedvirtual

Implements DetIdAssociator.

Definition at line 48 of file MuonDetIdAssociator.cc.

References getGeomDet(), point, and position.

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

◆ getValidDetIds()

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

overrideprotectedvirtual

Implements DetIdAssociator.

Definition at line 53 of file MuonDetIdAssociator.cc.

References cscbadchambers_, TrackingGeometry::dets(), dt, Exception, mixOne_premix_on_sim_cfi::gem, geometry_, includeBadChambers_, includeGEM_, includeME0_, CSCBadChambers::isInBadChamber(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, mixOne_premix_on_sim_cfi::me0, alignCSCRings::r, omtf::DataWord64::rpc, and GlobalTrackingGeometry::slaveGeometry().

                                                                                                       {
   validIds.clear();
   if (geometry_ == nullptr)
     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)) {
         if (gem->id().station() == 0)
           validIds.push_back(gem->id());
         else
           for (auto ch : gem->chambers())
             validIds.push_back(ch->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());
       }
     }
   }
 }

◆ insideElement()

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

overrideprotectedvirtual

Implements DetIdAssociator.

Definition at line 119 of file MuonDetIdAssociator.cc.

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

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