#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) {}

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.;
         }