CDCSegment.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/trackingUtilities/eventdata/trajectories/CDCTrajectory2D.h>
 
#include <cdc/topology/ISuperLayer.h>
#include <cdc/topology/EStereoKind.h>
 
#include <vector>
 
namespace Belle2 {
  namespace TrackingUtilities {

    template <class T>
    class CDCSegment : public std::vector<T> {
 
    public:
      CDCSegment() = default;

      CDC::EStereoKind getStereoKind() const
      {
        return  CDC::ISuperLayerUtil::getStereoKind(getISuperLayer());
      }

      bool isAxial() const
      {
        return getStereoKind() ==  CDC::EStereoKind::c_Axial;
      }

      CDC::ISuperLayer getISuperLayer() const
      {
        return  CDC::ISuperLayerUtil::getFrom(this->front());
      }

      const std::vector<T>& items() const
      {
        return *this;
      }

      CDCTrajectory2D& getTrajectory2D() const
      {
        return m_trajectory2D;
      }

      void setTrajectory2D(const CDCTrajectory2D& trajectory2D) const
      {
        m_trajectory2D = trajectory2D;
      }

      void sortByArcLength2D()
      {
        std::stable_sort(this->begin(), this->end(), [](const T & recoHit, const T & otherRecoHit) {
          return recoHit.getArcLength2D() < otherRecoHit.getArcLength2D();
        });
      }

      double getAliasScore() const
      {
        return m_aliasScore;
      }

      void setAliasScore(double aliasScore)
      {
        m_aliasScore = aliasScore;
      }
 
    private:
      mutable CDCTrajectory2D m_trajectory2D;

      // Stefano S: this was not initialized, I was not sure if to use 0 was safe
      // then I set the default value to NaN to solve a warning.
      // Not clear to me why the variable is "double" and not boolean
      double m_aliasScore = std::numeric_limits<double>::quiet_NaN();
 
    };
 
  }
}