MuonHoughTransformer_CurvedAtACylinder.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef MUONHOUGHPATTERNEVENT_MUONHOUGHTRANSFORMER_CURVEDATACYLINDER_H
#define MUONHOUGHPATTERNEVENT_MUONHOUGHTRANSFORMER_CURVEDATACYLINDER_H
 
#include "MuonHoughPatternEvent/MuonHoughTransformer.h"
 
class MuonHoughTransformer_CurvedAtACylinder : public MuonHoughTransformer {
public:
    MuonHoughTransformer_CurvedAtACylinder(int nbins, int nbins_angle, double detectorsize, double detectorsize_angle,
                                           double threshold_histo, int number_of_sectors = 1);
    ~MuonHoughTransformer_CurvedAtACylinder() = default;

    void fillHit(const std::shared_ptr<MuonHoughHit>& hit, double weight) override final;
    int fillHisto(double xbin, double theta, double weight, int sector)  override final;

    std::unique_ptr<MuonHoughPattern> hookAssociateHitsToMaximum(const MuonHoughHitContainer& event, std::pair<double, double> coordsmaximum,
                                                         double residu_mm, double residu_grad, int sector) const override final;

    int sector(const std::shared_ptr<MuonHoughHit>& hit) const override final;
    float weightHoughTransform(double r0) const  override final;
 
private:
    std::unique_ptr<double[]> m_invcurvature;
    std::unique_ptr<double[]> m_weightcurvature;
   
};
 
inline int MuonHoughTransformer_CurvedAtACylinder::sector(const std::shared_ptr<MuonHoughHit>& hit) const { return hit->phiSector(); }
 
#endif  // MUONHOUGHPATTERNEVENT_MUONHOUGHTRANSFORMER_CURVEDATACYLINDER_H