IRtRelation.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef MUONCALIB_IRTRELATION_H
#define MUONCALIB_IRTRELATION_H
 
#include "MdtCalibData/CalibFunc.h"
#include "MuonCalibMath/UtilFunc.h"
#include "GeoModelUtilities/TransientConstSharedPtr.h"
 
#include <optional>
 
namespace MuonCalib {
 
    class IRtRelation;
    using IRtRelationPtr = GeoModel::TransientConstSharedPtr<IRtRelation>;
    class IRtRelation : public CalibFunc {
    public:
        using CalibFunc::CalibFunc;
        virtual ~IRtRelation() = default;
        virtual std::string typeName() const override final { return "IRtRelation"; }

        virtual double radius(double t) const = 0;
        virtual double driftVelocity(double t) const = 0;
        virtual double driftAcceleration(double t) const = 0;
        virtual double tLower() const = 0;
        virtual double tUpper() const = 0;
        virtual double tBinWidth() const =0;
        virtual unsigned nDoF() const = 0;
        double GetTmaxDiff() const { return m_tmax_diff.value_or(0.); }
 
        bool hasTmaxDiff() const { return m_tmax_diff.has_value(); }

        void SetTmaxDiff(const double d) { m_tmax_diff = d; }

        double getReducedTime(const double  t) const {
           return mapToUnitInterval(t, tLower(), tUpper());
        }
        /* @brief Returns the derivative term of the reduced time w.r.t. the time*/
        double dReducedTimeDt() const {
            return unitIntervalPrime(tLower(), tUpper());
        }
 
      private:
        std::optional<double> m_tmax_diff{std::nullopt};
      protected:
        static constexpr double s_tBinWidth = 1.e-3;
    };
 
}  // namespace MuonCalib
 
#endif