The muon space point is the combination of two uncalibrated measurements one of them measures the eta and the other the phi coordinate. More...

#include <SpacePoint.h>

Collaboration diagram for MuonR4::SpacePoint:

Public Member Functions
	SpacePoint (const ActsGeometryContext &gctx, const xAOD::UncalibratedMeasurement primMeas, const xAOD::UncalibratedMeasurement secondMeas=nullptr)

const xAOD::UncalibratedMeasurement *	primaryMeasurement () const

const xAOD::UncalibratedMeasurement *	secondaryMeasurement () const

const MuonGMR4::SpectrometerSector *	msSector () const

const MuonGMR4::Chamber *	chamber () const
	: Pointer to the associated chamber More...

const Amg::Vector3D &	positionInChamber () const

const Amg::Vector3D &	directionInChamber () const

const Amg::Vector3D &	normalInChamber () const

Amg::Vector3D	planeNormal () const
	Returns the vector pointing out of the measurement plane. More...

xAOD::UncalibMeasType	type () const

bool	measuresPhi () const
	: Does the space point contain a phi measurement More...

bool	measuresEta () const
	: Does the space point contain an eta measurement More...

const Identifier &	identify () const
	: Identifier of the primary measurement More...

double	driftRadius () const
	: Returns the size of the drift radius More...

Amg::Vector2D	uncertainty () const
	: Returns the uncertainties on the space point More...

const	AmgSymMatrix (2) &covariance() const

bool	operator== (const SpacePoint &other) const
	: Equality check by checking the prd pointers More...

void	setInstanceCounts (std::shared_ptr< unsigned > etaCounts, std::shared_ptr< unsigned > phiCounts)
	Set the number of space points built with the same eta / phi prd. More...

unsigned int	nEtaInstanceCounts () const
	How many space points have been built in total with the same eta prd. More...

unsigned int	nPhiInstanceCounts () const
	How many space points have been built in total with the same phi prd More...

unsigned int	dimension () const
	Is the space point a 1D or combined 2D measurement. More...

Private Member Functions
	AmgSymMatrix (2) m_measCovariance
	: Measurement covariance If the spacePoint represents an 1D measurement the second coordinate is the length of the channel (e.g halfLength of the wire or of the associated strip) the uncertainty of the other coordinate, otherwise More...

Private Attributes
const xAOD::UncalibratedMeasurement *	m_primaryMeas {nullptr}

const xAOD::UncalibratedMeasurement *	m_secondaryMeas {nullptr}

Identifier	m_id {xAOD::identify(m_primaryMeas)}

const MuonGMR4::Chamber *	m_chamber {xAOD::readoutElement(m_primaryMeas)->chamber()}

const MuonGMR4::SpectrometerSector *	m_msSector {m_chamber->parent()}

Amg::Vector3D	m_pos {Amg::Vector3D::Zero()}

Amg::Vector3D	m_dir {Amg::Vector3D::Zero()}

Amg::Vector3D	m_normal {Amg::Vector3D::Zero()}

std::shared_ptr< const unsigned >	m_etaInstances {std::make_shared<unsigned>(1)}
	In how many space points is the eta measurement used. More...

std::shared_ptr< const unsigned >	m_phiInstances {std::make_shared<unsigned>(1)}
	In how many space points is the phi measurement used. More...

Detailed Description

The muon space point is the combination of two uncalibrated measurements one of them measures the eta and the other the phi coordinate.

In cases, without a complementary measurment the spacepoint just represents the single measurement and hence has a uncertainty into the other direction corresponding to the half-length of the measurement channel

Definition at line 18 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/MuonSpacePoint/SpacePoint.h.

Constructor & Destructor Documentation

◆ SpacePoint()

SpacePoint::SpacePoint	(	const ActsGeometryContext &	gctx,
		const xAOD::UncalibratedMeasurement *	primMeas,
		const xAOD::UncalibratedMeasurement *	secondMeas = `nullptr`
	)

Position of the measurements expressed in the sector frame

Direction along which the measurement strips point to

Intersect the two channels to define the space point

Do nothing for the pads

In case of 2D measurements like sTgc-pads or BI-RPC strips we can directly take the covariance from the measurement itself. To indicate that the space point measures both, eta & phi coordinate set the secondary measurement to be the primary one

Definition at line 33 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                                                             :
         m_primaryMeas{primaryMeas},
         m_secondaryMeas{secondaryMeas},
         m_dir{xAOD::channelDirInChamber(gctx, primaryMeas)},
         m_normal{xAOD::channelNormalInChamber(gctx, primaryMeas)} {
         
         AmgSymMatrix(2) Jac{AmgSymMatrix(2)::Identity()}, uvcov {AmgSymMatrix(2)::Identity()};
         
  
         if (primaryMeas->numDimensions() == 1) {
             uvcov(0,0) = primaryMeas->localCovariance<1>()[0];
         }
         Jac.col(0)  = m_normal.block<2,1>(0,0).unit();
         if (secondaryMeas) {
             const Amg::Vector3D pos1{xAOD::positionInChamber(gctx, primaryMeas)};
             const Amg::Vector3D pos2{xAOD::positionInChamber(gctx, secondaryMeas)};
             const Amg::Vector3D dir2{xAOD::channelDirInChamber(gctx, secondaryMeas)};
             m_pos = pos1 + Amg::intersect<3>(pos2,dir2, pos1, m_dir).value_or(0) * m_dir;
             Jac.col(1)  = xAOD::channelNormalInChamber(gctx, secondaryMeas).block<2,1>(0,0).unit();
             uvcov(1,1) = secondaryMeas->localCovariance<1>()[0]; 
         } else { 
             m_pos = xAOD::positionInChamber(gctx, primaryMeas);
             Jac.col(1) = m_dir.block<2,1>(0,0).unit();
             if (primaryMeas->type() == xAOD::UncalibMeasType::MdtDriftCircleType) {
                 const xAOD::MdtDriftCircle* dc = static_cast<const xAOD::MdtDriftCircle*>(primaryMeas);
                 uvcov(1,1) = 0.5* dc->readoutElement()->activeTubeLength(dc->measurementHash());
             } else if (primaryMeas->type() == xAOD::UncalibMeasType::RpcStripType) {
                 if (primaryMeas->numDimensions() == 1) {
                     const xAOD::RpcStrip* strip = static_cast<const xAOD::RpcStrip*>(primaryMeas);
                     uvcov(1,1) = strip->measuresPhi() ? 0.5* strip->readoutElement()->stripPhiLength():
                                                         0.5* strip->readoutElement()->stripEtaLength();
                 }
             } else if (primaryMeas->type() == xAOD::UncalibMeasType::TgcStripType) {
                 const xAOD::TgcStrip* strip = static_cast<const xAOD::TgcStrip*>(primaryMeas);
                 if (strip->measuresPhi()) {
                     uvcov(1,1) = 0.5 * strip->readoutElement()->stripLayout(strip->gasGap()).stripLength(strip->channelNumber());
                 } else {
                     uvcov(1,1) = 0.5 * strip->readoutElement()->wireGangLayout(strip->gasGap()).stripLength(strip->channelNumber());
                 }
             } else if (primaryMeas->type() == xAOD::UncalibMeasType::MMClusterType) {
                 const xAOD::MMCluster* clust = static_cast<const xAOD::MMCluster*>(primaryMeas);
                 uvcov(1,1) = 0.5 * clust->readoutElement()->stripLayer(clust->measurementHash()).design().stripLength(clust->channelNumber());
             } else if (primaryMeas->type() == xAOD::UncalibMeasType::sTgcStripType) {
                 const xAOD::sTgcMeasurement* meas = static_cast<const xAOD::sTgcMeasurement*>(primaryMeas);
                 switch (meas->channelType()) {
                     case sTgcIdHelper::sTgcChannelTypes::Strip:
                         uvcov(1,1) = 0.5 *  meas->readoutElement()->stripDesign(meas->measurementHash()).stripLength(meas->channelNumber());
                         break;
                     case sTgcIdHelper::sTgcChannelTypes::Wire:
                         uvcov(1,1) = 0.5 *  meas->readoutElement()->wireDesign(meas->measurementHash()).stripLength(meas->channelNumber());
                         break;
                     case sTgcIdHelper::sTgcChannelTypes::Pad:
                         break;
                 }
             }
             uvcov(1,1) = std::pow(uvcov(1,1), 2);
         }
  
         if (primaryMeas->numDimensions() == 2) {
             uvcov = xAOD::toEigen(primaryMeas->localCovariance<2>());
             m_secondaryMeas = m_primaryMeas;
         }
         m_measCovariance = Jac.inverse() * uvcov * Jac;
     }

Member Function Documentation

◆ AmgSymMatrix() [1/2]

const MuonR4::SpacePoint::AmgSymMatrix ( 2 ) const &

◆ AmgSymMatrix() [2/2]

MuonR4::SpacePoint::AmgSymMatrix ( 2 )

inlineprivate

: Measurement covariance If the spacePoint represents an 1D measurement the second coordinate is the length of the channel (e.g halfLength of the wire or of the associated strip) the uncertainty of the other coordinate, otherwise

Definition at line 89 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/MuonSpacePoint/SpacePoint.h.

89 {AmgSymMatrix(2)::Identity()};

◆ chamber()

const MuonGMR4::Chamber * SpacePoint::chamber ( ) const

: Pointer to the associated chamber

Definition at line 113 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                                      {
         return m_chamber;
     }

◆ dimension()

unsigned int SpacePoint::dimension ( ) const

Is the space point a 1D or combined 2D measurement.

Definition at line 160 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                              { 
         return (secondaryMeasurement() != nullptr) + 1;
     }

◆ directionInChamber()

const Amg::Vector3D & SpacePoint::directionInChamber ( ) const

Definition at line 122 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                                             {
         return m_dir;
     } 

◆ driftRadius()

double SpacePoint::driftRadius ( ) const

: Returns the size of the drift radius

Definition at line 143 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                          { 
         return m_primaryMeas->type() == xAOD::UncalibMeasType::MdtDriftCircleType ?  
                static_cast<const xAOD::MdtDriftCircle*>(m_primaryMeas)->driftRadius() : 0.;
     }

◆ identify()

const Identifier & SpacePoint::identify ( ) const

: Identifier of the primary measurement

Definition at line 140 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                                  {
         return m_id;
     }

◆ measuresEta()

bool SpacePoint::measuresEta ( ) const

: Does the space point contain an eta measurement

Definition at line 137 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                        {
         return secondaryMeasurement() ||  !chamber()->idHelperSvc()->measuresPhi(identify());
     }

◆ measuresPhi()

bool SpacePoint::measuresPhi ( ) const

: Does the space point contain a phi measurement

Definition at line 134 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                        {
         return secondaryMeasurement() ||  chamber()->idHelperSvc()->measuresPhi(identify());
     }

◆ msSector()

const MuonGMR4::SpectrometerSector * SpacePoint::msSector ( ) const

Definition at line 116 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                                                  {
         return m_msSector;
     }

◆ nEtaInstanceCounts()

unsigned int SpacePoint::nEtaInstanceCounts ( ) const

How many space points have been built in total with the same eta prd.

Definition at line 158 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

158 { return (*m_etaInstances); }

◆ normalInChamber()

const Amg::Vector3D & SpacePoint::normalInChamber ( ) const

Definition at line 125 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                                          {
         return m_normal;
     }

◆ nPhiInstanceCounts()

unsigned int SpacePoint::nPhiInstanceCounts ( ) const

How many space points have been built in total with the same phi prd

Definition at line 159 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

159 { return (*m_phiInstances); }

◆ operator==()

bool MuonR4::SpacePoint::operator== ( const SpacePoint & other ) const

inline

: Equality check by checking the prd pointers

Definition at line 60 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/MuonSpacePoint/SpacePoint.h.

                                                            {
                 return primaryMeasurement() == other.primaryMeasurement() &&
                        secondaryMeasurement() == other.secondaryMeasurement();
             }

◆ planeNormal()

Amg::Vector3D SpacePoint::planeNormal ( ) const

Returns the vector pointing out of the measurement plane.

Definition at line 128 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                               {
         return directionInChamber().cross(normalInChamber()).unit();
     }

◆ positionInChamber()

const Amg::Vector3D & SpacePoint::positionInChamber ( ) const

Definition at line 119 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                                            {
         return m_pos;
     }

◆ primaryMeasurement()

const xAOD::UncalibratedMeasurement * SpacePoint::primaryMeasurement ( ) const

Definition at line 107 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                                                             {
        return m_primaryMeas;
     }

◆ secondaryMeasurement()

const xAOD::UncalibratedMeasurement * SpacePoint::secondaryMeasurement ( ) const

Definition at line 110 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                                                               {
        return m_secondaryMeas;
     }

◆ setInstanceCounts()

void MuonR4::SpacePoint::setInstanceCounts	(	std::shared_ptr< unsigned >	etaCounts,
		std::shared_ptr< unsigned >	phiCounts
	)

Set the number of space points built with the same eta / phi prd.

◆ type()

xAOD::UncalibMeasType SpacePoint::type ( ) const

Definition at line 131 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                                {
         return primaryMeasurement()->type();
     }

◆ uncertainty()

Amg::Vector2D SpacePoint::uncertainty ( ) const

: Returns the uncertainties on the space point

Definition at line 147 of file MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/SpacePoint.cxx.

                                               {
         return Amg::Vector2D{ Amg::error(m_measCovariance,0),Amg::error(m_measCovariance,1)  };
     }