MdtAnalyticRtCalibAlg.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef MUONCALIBR4_MDTANALYTICCALIBALG_H
#define MUONCALIBR4_MDTANALYTICCALIBALG_H
 
#include <AthenaBaseComps/AthHistogramAlgorithm.h>
 
#include "MuonIdHelpers/IMuonIdHelperSvc.h"
#include "MdtCalibData/MdtCalibDataContainer.h"
 
#include "StoreGate/ReadCondHandleKey.h"
#include "StoreGate/WriteCondHandleKey.h"
 
class TGraph;
namespace MuonCalib{
  class SamplePoint;
}
namespace MuonCalibR4{
    class MdtAnalyticRtCalibAlg : public AthHistogramAlgorithm {
      public:
        using AthHistogramAlgorithm::AthHistogramAlgorithm;
        virtual ~MdtAnalyticRtCalibAlg() = default;
 
        virtual StatusCode initialize() override final;
        virtual StatusCode execute() override final;
 
        enum class PolyType{
            ChebyChev,
            Legendre,
            Simple
        };
        using RtRelationPtr = MuonCalib::MdtFullCalibData::RtRelationPtr;
 
      private:
 
          void drawRt(const EventContext& ctx, 
                      const Identifier& detId,
                      const std::vector<MuonCalib::SamplePoint>& rtPoints,
                      const MuonCalib::MdtRtRelation& inRel) const;
 
          void drawResoFunc(const EventContext& ctx,
                            const Identifier& detId,
                            const std::vector<MuonCalib::SamplePoint>& resoPoints,
                            const MuonCalib::IRtResolution& inReso) const;
 
          RtRelationPtr translateRt(const EventContext& ctx,
                                    const Identifier& detId,
                                    const MuonCalib::MdtRtRelation& inRel) const;
          
          
          std::vector<Identifier> tubeIds(const Identifier& chId) const;
          
          void saveGraph(const std::string& path, std::unique_ptr<TGraph>&& graph) const;
          ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc{this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};
          SG::ReadCondHandleKey<MuonCalib::MdtCalibDataContainer> m_readKey{this, "ReadKey", "MdtCalibConstantsR4"};
          SG::WriteCondHandleKey<MuonCalib::MdtCalibDataContainer> m_writeKey{this, "WriteKey", "MdtCalibConstants",
                                                                                            "Conditions object containing the calibrations"};   
 
          Gaudi::Property<unsigned> m_maxOrder{this, "maxOrder", 12};
          Gaudi::Property<int> m_polyTypeRt{this, "PolyTypeRt", static_cast<int>(PolyType::ChebyChev)};
          Gaudi::Property<int> m_polyTypeTr{this, "PolyTypeTr", static_cast<int>(PolyType::Legendre)};
          Gaudi::Property<bool> m_fitRtReso{this, "FitRtReso", true};
          Gaudi::Property<double> m_relUncReso{this, "RelUncertOnReso", 0.01};
          Gaudi::Property<unsigned> m_maxOrderReso{this, "maxOrderReso", 25};
          Gaudi::Property<float> m_chiCutOff{this,"chi2CutOff", 0.05};
 
          Gaudi::Property<bool> m_saveDiagnostic{this, "saveDiagnosticHist", true};
          Gaudi::Property<std::string> m_outStream{this, "OutStream", "MDTANALYTICRTS"};
          Gaudi::Property<unsigned> m_precCutOff{this,"precCutOff", 6};
          Gaudi::Property<bool> m_fillMissingCh{this, "fillMissingCh", true};
          Gaudi::Property<float> m_missingT0{this, "missingT0", 5.4597301};
    };
}
 
 
 
#endif