WireNeighborKind.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Oliver Frost                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include <tracking/trackFindingCDC/numerics/Modulo.h>
#include <tracking/trackFindingCDC/numerics/ESign.h>
 
namespace Belle2 {
  namespace TrackFindingCDC {
 
    class WireNeighborKind {
 
    public:
      WireNeighborKind() = default;
 
      WireNeighborKind(short cellDistance, short oClockDirection)
        : m_cellDistance(cellDistance)
        , m_oClockDirection(moduloFast(oClockDirection, 12))
      {}
 
      short getILayerDifference() const
      {
        return getILayerDifferenceSign() * getAbsILayerDifference();
      }
 
      short getAbsILayerDifference() const
      {
        short slope = abs(symmetricModuloFast(getOClockDirection() + 3, 6));
        short iAbsLayerDifference = slope * getCellDistance() / 2;
        if (iAbsLayerDifference > 2) iAbsLayerDifference = 2;
        return iAbsLayerDifference;
      }
 
      short getILayerDifferenceSign() const
      {
        return sign(symmetricModuloFast(getOClockDirection() + 3, 12));
      }
 
      short getCellDistance() const
      { return m_cellDistance; }
 
      short getOClockDirection() const
      { return m_oClockDirection; }
 
      bool isValid() const
      {
        return
          m_cellDistance >= 1 and
          m_cellDistance <= 2 and
          m_oClockDirection >= 0 and
          m_oClockDirection < 12;
      }
 
    private:
      short m_cellDistance = 0;
 
      short m_oClockDirection = 0;
 
    };
 
  }
}