MuonR4::SpacePointPerLayerSorter Class Reference

The SpacePointPerLayerSorter groups the space points by their layer Identifier. More...

#include <SpacePointPerLayerSorter.h>

Collaboration diagram for MuonR4::SpacePointPerLayerSorter:

Public Types
using	HitVec = std::vector< const SpacePoint * >
using	HitLayVec = std::vector< HitVec >

Public Member Functions
	SpacePointPerLayerSorter (const SpacePointBucket &bucket)
	Constructor taking a complete bucket More...
	SpacePointPerLayerSorter (HitVec vec)
	Constructor taking a subset of SpacePoints. More...
const HitLayVec &	mdtHits () const
	Returns the sorted Mdt hits. More...
unsigned int	nMdtHits () const
	Returns the number of all Mdt hits in the seed. More...
const HitLayVec &	stripHits () const
	Returns the sorted strip hits. More...
unsigned int	nStripHits () const
	Returns the number of all strip hits in the seed. More...
unsigned int	firstLayerFrom2ndMl () const
	Returns the layer index with hits from the second multilayer More...

Private Attributes
HitLayVec	m_mdtLayers {}
	Sorted Mdt hits per tube layer. More...
HitLayVec	m_stripLayers {}
	Sorted Strip hits per gasGap strip More...
unsigned int	m_nMdtHits {0}
	Number of all Mdt tube hits More...
unsigned int	m_nStripHits {0}
	Number of all strip hits. More...
unsigned int	m_tubeLaySwitch {0}
	Index of the first tube-layer from the second multilayer. More...

Detailed Description

The SpacePointPerLayerSorter groups the space points by their layer Identifier.

It is defined as the Identifier of the first tube in layer for the Mdts or as the Identifier of the first strip in a gasGap expressed in an eta view. First, all hits are sorted by increasing chamber z - i.e. going outwards the detector, and then grouped into two sets of vectors. One for the Mdts and the other for the remaining strip detectors.

Definition at line 14 of file SpacePointPerLayerSorter.h.

Member Typedef Documentation

HitLayVec

using MuonR4::SpacePointPerLayerSorter::HitLayVec = std::vector<HitVec>

Definition at line 17 of file SpacePointPerLayerSorter.h.

HitVec

using MuonR4::SpacePointPerLayerSorter::HitVec = std::vector<const SpacePoint*>

Definition at line 16 of file SpacePointPerLayerSorter.h.

Constructor & Destructor Documentation

SpacePointPerLayerSorter() [1/2]

MuonR4::SpacePointPerLayerSorter::SpacePointPerLayerSorter

(

const SpacePointBucket &

bucket

)

Constructor taking a complete bucket

Definition at line 17 of file SpacePointPerLayerSorter.cxx.

                                                                                    :
        SpacePointPerLayerSorter(stripSmartPtr(bucket)){}

SpacePointPerLayerSorter() [2/2]

MuonR4::SpacePointPerLayerSorter::SpacePointPerLayerSorter

(

HitVec

vec

)

Constructor taking a subset of SpacePoints.

Sort space points by z

The hits are radially sorted from low local-z to high local z. Build the gasGap Identifier to find out to which layer the hit belongs to and then use the layer counting map as auxillary object fetch the indices for the sorted measurements

Definition at line 20 of file SpacePointPerLayerSorter.cxx.

                                                                  {
        if (hits.empty()) return;
 
        std::ranges::stable_sort(hits, [](const SpacePoint* a, const SpacePoint*b){
                        const double dZ = a->positionInChamber().z() - b->positionInChamber().z(); 
                        if (std::abs(dZ) > 50.*Gaudi::Units::micrometer) {
                            return dZ < 0.;
                        }
                        return a->positionInChamber().y() < b->positionInChamber().y();
        });
        const Muon::IMuonIdHelperSvc* idHelperSvc{hits.front()->msSector()->idHelperSvc()};
        m_mdtLayers.reserve(8);
 
        using LayerCounting = std::unordered_map<Identifier, unsigned int>;
        LayerCounting mdtLayerCounting{}, stripLayerCounting{};
        for (const SpacePoint* hit : hits) {
            const Identifier& id {hit->identify()};
            if (hit->type() == xAOD::UncalibMeasType::MdtDriftCircleType) {
                const MdtIdHelper& idHelper{idHelperSvc->mdtIdHelper()};
                const Identifier layId = idHelper.channelID(idHelper.stationName(id), 1, idHelper.stationPhi(id), 
                                                            idHelper.multilayer(id), idHelper.tubeLayer(id), 1);
 
                const unsigned int layer = mdtLayerCounting.insert(std::make_pair(layId, mdtLayerCounting.size())).first->second;
                if (layer >= m_mdtLayers.size()) {
                    m_mdtLayers.resize(layer + 1);
                }
                HitVec& pushTo{m_mdtLayers[layer]};
                if (pushTo.capacity() == pushTo.size()){
                    pushTo.reserve(pushTo.size() + 5);
                }
                pushTo.push_back(hit);
                ++m_nMdtHits;
                if (!m_tubeLaySwitch && layer &&  
                    m_mdtLayers[layer -1].front()->primaryMeasurement()->identifierHash() != hit->primaryMeasurement()->identifierHash()) {
 
                }
            } else {
                const Identifier layId = idHelperSvc->gasGapId(id);
                const unsigned int layer = stripLayerCounting.insert(std::make_pair(layId, stripLayerCounting.size())).first->second;
            
                if (layer >= m_stripLayers.size()) {
                    m_stripLayers.resize(layer + 1);
                }
                HitVec& pushTo{m_stripLayers[layer]};
                if (pushTo.capacity() == pushTo.size()){
                    pushTo.reserve(pushTo.size() + 5);
                }
                pushTo.push_back(hit);
                ++m_nStripHits;
            }
        }
    }

Member Function Documentation

firstLayerFrom2ndMl()

unsigned int MuonR4::SpacePointPerLayerSorter::firstLayerFrom2ndMl ( ) const

inline

Returns the layer index with hits from the second multilayer

Definition at line 39 of file SpacePointPerLayerSorter.h.

                                                     {
              return m_tubeLaySwitch;
            }

mdtHits()

const HitLayVec& MuonR4::SpacePointPerLayerSorter::mdtHits ( ) const

inline

Returns the sorted Mdt hits.

Definition at line 23 of file SpacePointPerLayerSorter.h.

                                             {
              return m_mdtLayers;
            }

nMdtHits()

unsigned int MuonR4::SpacePointPerLayerSorter::nMdtHits ( ) const

inline

Returns the number of all Mdt hits in the seed.

Definition at line 27 of file SpacePointPerLayerSorter.h.

                                          {
              return m_nMdtHits;
            }

nStripHits()

unsigned int MuonR4::SpacePointPerLayerSorter::nStripHits ( ) const

inline

Returns the number of all strip hits in the seed.

Definition at line 35 of file SpacePointPerLayerSorter.h.

                                            {
              return m_nStripHits;
            }

stripHits()

const HitLayVec& MuonR4::SpacePointPerLayerSorter::stripHits ( ) const

inline

Returns the sorted strip hits.

Definition at line 31 of file SpacePointPerLayerSorter.h.

                                               {
              return m_stripLayers;
            }

Member Data Documentation

m_mdtLayers

HitLayVec MuonR4::SpacePointPerLayerSorter::m_mdtLayers {}

private

Sorted Mdt hits per tube layer.

Definition at line 44 of file SpacePointPerLayerSorter.h.

m_nMdtHits

unsigned int MuonR4::SpacePointPerLayerSorter::m_nMdtHits {0}

private

Number of all Mdt tube hits

Definition at line 48 of file SpacePointPerLayerSorter.h.

m_nStripHits

unsigned int MuonR4::SpacePointPerLayerSorter::m_nStripHits {0}

private

Number of all strip hits.

Definition at line 50 of file SpacePointPerLayerSorter.h.

m_stripLayers

HitLayVec MuonR4::SpacePointPerLayerSorter::m_stripLayers {}

private

Sorted Strip hits per gasGap strip

Definition at line 46 of file SpacePointPerLayerSorter.h.

m_tubeLaySwitch

unsigned int MuonR4::SpacePointPerLayerSorter::m_tubeLaySwitch {0}

private

Index of the first tube-layer from the second multilayer.

Definition at line 52 of file SpacePointPerLayerSorter.h.

---

The documentation for this class was generated from the following files:

• SpacePointPerLayerSorter.h
• SpacePointPerLayerSorter.cxx