#include <SortMeasurementsByPosition.h>

Collaboration diagram for Muon::SortTSOSs:

Public Member Functions
bool	operator() (const Trk::TrackStateOnSurface tsos1, const Trk::TrackStateOnSurface tsos2) const
bool	operator() (const std::unique_ptr< const Trk::TrackStateOnSurface > &tsos1, const std::unique_ptr< const Trk::TrackStateOnSurface > &tsos2) const
	SortTSOSs (const IMuonEDMHelperSvc h, const IMuonIdHelperSvc idh)

Public Attributes
const IMuonEDMHelperSvc *	m_helperSvc
const IMuonIdHelperSvc *	m_idHelperSvc

Detailed Description

Definition at line 80 of file SortMeasurementsByPosition.h.

Constructor & Destructor Documentation

◆ SortTSOSs()

Muon::SortTSOSs::SortTSOSs	(	const IMuonEDMHelperSvc *	h,
		const IMuonIdHelperSvc *	idh )

inline

Definition at line 146 of file SortMeasurementsByPosition.h.

146: m_helperSvc(h), m_idHelperSvc(idh) {}

Muon::SortTSOSs::m_idHelperSvc

const IMuonIdHelperSvc * m_idHelperSvc

Definition SortMeasurementsByPosition.h:149

Muon::SortTSOSs::m_helperSvc

const IMuonEDMHelperSvc * m_helperSvc

Definition SortMeasurementsByPosition.h:148

Member Function Documentation

◆ operator()() [1/2]

bool Muon::SortTSOSs::operator()	(	const std::unique_ptr< const Trk::TrackStateOnSurface > &	tsos1,
		const std::unique_ptr< const Trk::TrackStateOnSurface > &	tsos2 ) const

inline

Definition at line 141 of file SortMeasurementsByPosition.h.

                                                                                        {
            return this->operator()(tsos1.get(), tsos2.get());
        }

◆ operator()() [2/2]

bool Muon::SortTSOSs::operator()	(	const Trk::TrackStateOnSurface *	tsos1,
		const Trk::TrackStateOnSurface *	tsos2 ) const

inline

Definition at line 82 of file SortMeasurementsByPosition.h.

                                                                                                        {
            if (!tsos1->trackParameters()) {
                std::cout << "Muon::SortTSOSs: state 1 without parameters " << std::endl;
                return false;
            }
            if (!tsos2->trackParameters()) {
                std::cout << "Muon::SortTSOSs: state 2 without parameters " << std::endl;
                return true;
            }
            // first, address rare problems with funky track directions leading to strange orderings by checking for comparisons of calo
            // deposits and MS hits
            const Trk::MeasurementBase* meas1 = tsos1->measurementOnTrack();
            Identifier id1 = meas1 ? m_helperSvc->getIdentifier(*meas1) : Identifier();
 
            const Trk::MeasurementBase* meas2 = tsos2->measurementOnTrack();
            Identifier id2 = meas2 ? m_helperSvc->getIdentifier(*meas2) : Identifier();
 
            bool okId1 = id1.is_valid() && m_idHelperSvc->isMuon(id1) ? true : false;
            bool okId2 = id2.is_valid() && m_idHelperSvc->isMuon(id2) ? true : false;
 
            if (tsos1->type(Trk::TrackStateOnSurface::TrackStateOnSurfaceType::CaloDeposit) ||
                tsos2->type(Trk::TrackStateOnSurface::TrackStateOnSurfaceType::CaloDeposit)) {
                if (okId1 && tsos2->type(Trk::TrackStateOnSurface::TrackStateOnSurfaceType::CaloDeposit)) return false;
                if (tsos1->type(Trk::TrackStateOnSurface::TrackStateOnSurfaceType::CaloDeposit) && okId2) return true;
            }
 
            // get average direction of the 2 TSOSs
            Amg::Vector3D trackDir = tsos1->trackParameters()->momentum().unit();
            trackDir += tsos2->trackParameters()->momentum().unit();
            const Amg::Vector3D& pos1 = tsos1->trackParameters()->position();
            const Amg::Vector3D& pos2 = tsos2->trackParameters()->position();
            double dist = (pos2 - pos1).dot(trackDir);
 
            if (std::abs(dist) < 1e-5) {
                // one is a good muon hit and one is not: good muon hit comes after
                if (okId1 && !okId2) return true;
                if (!okId1 && okId2) return false;
                // both invalid or non-muon: consider them equal
                if (!okId1 && !okId2) return false;
                // now we have 2 valid muon Ids
                bool measPhi1 = m_idHelperSvc->measuresPhi(id1);
                bool measPhi2 = m_idHelperSvc->measuresPhi(id2);
                // put phi measurement in front of eta measurements
                if (measPhi1 && !measPhi2) return true;
                if (!measPhi1 && measPhi2) return false;
                // now both are phi or both are eta
                // decide on identifier (increasing away from IP within same technology)
                double inOrOut = (pos1 + pos2).dot(trackDir);
                if (inOrOut >= 0.0) {
                    return id1 < id2;
                } else {
                    return id1 > id2;
                }
 
                // if we get here, there was no decision
                return false;  // consider them equal
            }
            return dist > 0.;
        }