StandardAxes.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#pragma once
 
#include <tracking/trackFindingCDC/hough/axes/DiscreteValue.h>
#include <tracking/trackFindingCDC/hough/axes/ContinuousValue.h>
 
#include <cdc/topology/ILayer.h>
 
#include <Math/Vector2D.h>
 
namespace Belle2 {
  namespace TrackFindingCDC {
    class ImpactTag;

    using DiscreteImpact = DiscreteValue<float, ImpactTag>;

    using ContinuousImpact = ContinuousValue<double, ImpactTag>;

    class CurvTag;

    using DiscreteCurv = DiscreteValue<float, CurvTag>;

    using ContinuousCurv = ContinuousValue<double, CurvTag>;

    class CurvWithArcLength2DCache {
    public:
      explicit CurvWithArcLength2DCache(float curv);

      explicit operator float() const
      {
        return m_curv;
      }

      float getArcLength2D(CDC::ILayer iCLayer, bool secondArm = false) const
      {
        return secondArm ? m_secondaryArcLength2DByICLayer[iCLayer] : m_arcLength2DByICLayer[iCLayer];
      }

      friend std::ostream& operator<<(std::ostream& output, const CurvWithArcLength2DCache& value)
      {
        return output << value.m_curv;
      }
 
    private:
      std::array<float, 56> m_arcLength2DByICLayer;

      std::array<float, 56> m_secondaryArcLength2DByICLayer;

      float m_curv;
    };

    using DiscreteCurvWithArcLength2DCache = DiscreteValue<CurvWithArcLength2DCache, CurvTag>;

    class TanLTag;

    using DiscreteTanL = DiscreteValue<float, TanLTag>;

    using ContinuousTanL = ContinuousValue<double, TanLTag>;

    class Z0Tag;

    using DiscreteZ0 = DiscreteValue<float, Z0Tag>;

    class Phi0Tag;

    using DiscretePhi0 = DiscreteValue<ROOT::Math::XYVector, Phi0Tag>;
 
    // NOTE The following is for quadratic and hyperbolic tracks
    class PTag;

    using DiscreteP = DiscreteValue<float, PTag>;

    using ContinuousP = ContinuousValue<double, PTag>;

    class QTag;

    using DiscreteQ = DiscreteValue<float, QTag>;

    using ContinuousQ = ContinuousValue<double, QTag>;
  }
}