AuthorHierachy.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONTRACKEVENT_AuthorHierachy_H
#define MUONTRACKEVENT_AuthorHierachy_H
 
#include "xAODMuon/Muon.h"
 
#include "Acts/Utilities/PointerTraits.hpp"
#include "AthContainers/DataVector.h"
 
#include <concepts>

namespace MuonR4 {
 
    namespace detail {
        template <typename ObjType> concept hasAuthor = requires(const ObjType& obj) {
            { obj.author() } -> std::same_as<xAOD::Muon::Author>;
        };

        static_assert(hasAuthor<xAOD::Muon>);
    }

    struct AuthorHierachy {
        using Author = xAOD::Muon::Author;
        unsigned authorRank(const Author author) const;
        bool operator()(const Author a, const Author b) const;
        template <detail::hasAuthor ObjType>
        bool operator()(const ObjType&a ,const  ObjType& b) const {
            return (*this)(a.author(), b.author());
        }

        template <Acts::PointerConcept ObjPtr>
        bool operator()(const ObjPtr& a, const ObjPtr& b) const  {
            assert(a != nullptr);
            assert(b != nullptr);
            return (*this)(*a, *b);
        }

        template <detail::hasAuthor ObjType>
        bool operator()(const DataVector<ObjType>&a, const DataVector<ObjType>& b) const {
            if (a.empty() || b.empty()) {
                return !a.empty();
            }
            // Check in the debug build that there is only one author in each vector
            assert(std::all_of(a.begin(), a.end(), [&](const auto obj){
                              return obj->author() == a.at(0)->author();
                            }));
            assert(std::all_of(b.begin(), b.end(), [&](const auto obj){
                              return obj->author() == b.at(0)->author();
                            }));
 
            return (*this)(a.front(), b.front());
        }

        template <detail::hasAuthor ObjType>
        bool operator()(const DataVector<ObjType>* a, const DataVector<ObjType>* b) const {
            assert( a != nullptr);
            assert( b != nullptr);
            return (*this)(*a, *b);
        }
    };
 
}
 
#endif