trackAlgoPriorityOrder.h

Go to the documentation of this file.

#ifndef RecoTracker_FinalTrackSelectors_trackAlgoPriorityOrder_h
#define RecoTracker_FinalTrackSelectors_trackAlgoPriorityOrder_h

#include "DataFormats/TrackReco/interface/TrackBase.h"

#include <array>

namespace impl {

  constexpr reco::TrackBase::TrackAlgorithm algoPriorityOrder[] ={
    reco::TrackBase::undefAlgorithm,
    reco::TrackBase::ctf,
    reco::TrackBase::cosmics,
    reco::TrackBase::duplicateMerge,
    reco::TrackBase::initialStep,
    reco::TrackBase::highPtTripletStep,
    reco::TrackBase::detachedQuadStep,
    reco::TrackBase::detachedTripletStep,
    reco::TrackBase::lowPtQuadStep,
    reco::TrackBase::lowPtTripletStep,
    reco::TrackBase::pixelPairStep,
    reco::TrackBase::mixedTripletStep,
    reco::TrackBase::pixelLessStep,
    reco::TrackBase::tobTecStep,
    reco::TrackBase::jetCoreRegionalStep,
    reco::TrackBase::conversionStep,
    reco::TrackBase::muonSeededStepInOut,
    reco::TrackBase::muonSeededStepOutIn,
    reco::TrackBase::outInEcalSeededConv,
    reco::TrackBase::inOutEcalSeededConv,
    reco::TrackBase::nuclInter,
    reco::TrackBase::standAloneMuon,
    reco::TrackBase::globalMuon,
    reco::TrackBase::cosmicStandAloneMuon,
    reco::TrackBase::cosmicGlobalMuon,
    reco::TrackBase::reservedForUpgrades1,
    reco::TrackBase::reservedForUpgrades2,
    reco::TrackBase::bTagGhostTracks,
    reco::TrackBase::beamhalo,
    reco::TrackBase::gsf,
    reco::TrackBase::hltPixel,
    reco::TrackBase::hltIter0,
    reco::TrackBase::hltIter1,
    reco::TrackBase::hltIter2,
    reco::TrackBase::hltIter3,
    reco::TrackBase::hltIter4,
    reco::TrackBase::hltIterX,
    reco::TrackBase::hiRegitMuInitialStep,
    reco::TrackBase::hiRegitMuPixelPairStep,
    reco::TrackBase::hiRegitMuMixedTripletStep,
    reco::TrackBase::hiRegitMuPixelLessStep,
    reco::TrackBase::hiRegitMuDetachedTripletStep,
    reco::TrackBase::hiRegitMuMuonSeededStepInOut,
    reco::TrackBase::hiRegitMuMuonSeededStepOutIn,
    reco::TrackBase::hiRegitMuLowPtTripletStep,
    reco::TrackBase::hiRegitMuTobTecStep
  };

  static_assert(reco::TrackBase::algoSize == sizeof(algoPriorityOrder)/sizeof(unsigned int), "Please update me too after adding new enumerators to reco::TrackBase::TrackAlgorithm");

  constexpr unsigned int findIndex(const reco::TrackBase::TrackAlgorithm algo, const unsigned int index) {
    return index < sizeof(algoPriorityOrder)/sizeof(unsigned int) ?
      (algo == algoPriorityOrder[index] ? index : findIndex(algo, index+1)) :
      throw "Index out of bounds, this means that some reco::TrackBase::TrackAlgorithm enumerator is missing from impl::algoPriorityOrder array.";
  }

  constexpr unsigned int priorityForAlgo(const reco::TrackBase::TrackAlgorithm algo) {
    return findIndex(algo, 0);
  }


  template <typename T, size_t N, size_t I>
  struct MakeArray {
    template <typename ...Args>
    constexpr static
    std::array<T, N> value(Args&&... args) {
      return MakeArray<T, N, I-1>::value(priorityForAlgo(static_cast<reco::TrackBase::TrackAlgorithm>(I-1)), std::forward<Args>(args)...);
    }
  };

  template <typename T, size_t N>
  struct MakeArray<T, N, 0> {
    template <typename ...Args>
    constexpr static
    std::array<T, N> value(Args&&... args) {
      return std::array<T, N>{{std::forward<Args>(args)...}};
    }
  };


  template <typename T, size_t N>
  constexpr
  std::array<T, N> makeArray() {
    return MakeArray<T, N, N>::value();
  }

}

constexpr std::array<unsigned int, reco::TrackBase::algoSize> trackAlgoPriorityOrder = impl::makeArray<unsigned int, reco::TrackBase::algoSize>();


#endif // DataFormats_TrackReco_trackAlgoPriorityOrder_h