MuonHoughEventData.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef MUONR4__MUONHOUGHEVENTDATA__H 
#define MUONR4__MUONHOUGHEVENTDATA__H
 
#include "xAODMeasurementBase/UncalibratedMeasurement.h"
#include "xAODMuonPrepData/MdtDriftCircleContainer.h"
#include "xAODMuonPrepData/RpcStripContainer.h"
#include "xAODMuonPrepData/TgcStripContainer.h"
#include "xAODMuonPrepData/MMClusterContainer.h"
#include "MuonSpacePoint/MuonSpacePointContainer.h"
#include <map> 
#include <memory> 
 
namespace MuonR4{ 
 
  inline void updateSearchWindow(std::pair<double,double> & searchWindow, double value){
    searchWindow.first = std::min(searchWindow.first, value); 
    searchWindow.second = std::max(searchWindow.second, value); 
  }
 
  template <typename maximum_t, typename houghPlane_t, typename peakFinder_t, typename peakFinderConfig_t> 
  struct MuonHoughEventData_impl{
            MuonHoughEventData_impl(const ActsGeometryContext& _gctx):
                gctx{_gctx} {}
            // cached ACTS geometry context 
            const ActsGeometryContext& gctx; 
            // Hough accumulator 
            std::unique_ptr<houghPlane_t> houghPlane{nullptr}; 
            // peak finder instance
            std::unique_ptr<peakFinder_t> peakFinder{nullptr}; 
            struct HoughSetupForBucket{
                // Bucket containing the input space points 
                const MuonSpacePointBucket* bucket{nullptr};
                // placeholder boundary for the search space  
                static constexpr double dummyBoundary{1.e9};
                // search window for the intercept, optimised for the bucket
                std::pair<double, double> searchWindowIntercept{dummyBoundary, -dummyBoundary};
                // search window for the angle, optimised for the bucket
                std::pair<double, double> searchWindowTanAngle{dummyBoundary, -dummyBoundary};
 
                void updateHoughWindow(double interceptMeas, double TanAngleMeas) {
                    updateSearchWindow(searchWindowIntercept, interceptMeas); 
                    updateSearchWindow(searchWindowTanAngle, TanAngleMeas); 
                }
            };
            // the maxima found by the hough transform
            std::vector<maximum_t> maxima{};
            // the hough setup for each logical muon chamber
            std::map<const MuonGMR4::MuonChamber*, std::vector<HoughSetupForBucket>> houghSetups{}; 
            // the axis ranges currently mapped to the accumulator
            Acts::HoughTransformUtils::HoughAxisRanges currAxisRanges; 
            // current search window for the intercept
            std::pair<double,double> searchSpaceIntercept{10000000,-100000000.}; 
            // current search window for the angle
            std::pair<double,double> searchSpaceTanAngle{100000000.,-100000.}; 
            // hit counter for phi-measurements
            size_t phiHitsOnMax{0};
    };
}
 
#endif