GlobalPattern.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonFastRecoEvent/GlobalPattern.h"
#include "MuonDetDescrUtils/MuonSectorMapping.h"
 
namespace MuonR4 {
 
GlobalPattern::GlobalPattern(HitCollection&& hitPerStation)
    : m_hitsInStation(std::move(hitPerStation)) {};
 
std::vector<GlobalPattern::StIndex> GlobalPattern::getStations() const {
    std::vector<StIndex> out{};
    out.reserve(m_hitsInStation.size());
    std::ranges::transform(m_hitsInStation, std::back_inserter(out), 
        [](const auto& pair){ return pair.first;});
    return out;
}
 
const std::vector<GlobalPattern::HitType>& GlobalPattern::hitsInStation(StIndex station) const {
    const auto& it = m_hitsInStation.find(station);
    if (it != m_hitsInStation.end()) {
        return it->second;
    }
    static const std::vector<HitType> empty{};
    return empty;
}
 
double GlobalPattern::sectorPhi() const { 
    static const Muon::MuonSectorMapping sectorMap{};
    return sectorMap.sectorOverlapPhi(sector(), secondarySector()); 
}
 
void GlobalPattern::print(std::ostream& ostr) const {
    ostr<<"SpacePoint Pattern, Sector: "<< sector() << "  & " <<  (isSectorOverlap() ? std::to_string(secondarySector()) : "-")<< " Sector Phi: "<<sectorPhi()<<", theta: "<<theta()
        <<", nPrecisionHits: "<<nPrecisionHits()<<", nEtaNonPrecisionHits: "<<nEtaNonPrecisionHits()<<", nPhiHits: "<<nPhiHits()
        <<", total residual: "<<totalResidual()<<", total normalized residual: "<<totalNormalizedResidual();    
    ostr<<", Hit per station: \n";
    for (const auto& [station, hits] : m_hitsInStation) {
        ostr<<"  Station "<<stName(station)<<": "<<hits.size()<<" hits\n";
        for (const auto& hit : hits) {
            ostr<<"    "<<*hit<<"\n";
        }
    }
}
 
}