ImpactSweeped.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2015 - 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/hough/perigee/ImpactRep.h>
#include <tracking/trackFindingCDC/hough/boxes/SweepBox.h>
#include <tracking/trackFindingCDC/hough/baseelements/SameSignChecker.h>
#include <tracking/trackFindingCDC/topology/ILayer.h>
 
namespace Belle2 {
  namespace TrackFindingCDC {
 
    template<class T, class AImpact = ContinuousImpact>
    class AImpactSweeped : public T {
 
    public:
      using T::T;
 
    public:
      using HoughBox = SweepBox<AImpact, typename T::HoughBox>;
 
      using SubordinaryHoughBox = typename T::HoughBox;
 
    public:
      ESign getDistanceSign(const HoughBox& houghBox,
                            float x,
                            float y,
                            float signedDriftLength,
                            float dxdz = 0,
                            float dydz = 0,
                            ILayer iCLayer = -1) const
      {
        float lowerImpact(houghBox.template getLowerBound<AImpact>());
        float upperImpact(houghBox.template getUpperBound<AImpact>());
        const SubordinaryHoughBox& subordinaryHoughBox = houghBox.getSubordinaryBox();
 
        const float lowerY = y - lowerImpact;
        const ESign lowerDistSign = T::getDistanceSign(subordinaryHoughBox,
                                                       x, lowerY,
                                                       signedDriftLength,
                                                       dxdz, dydz,
                                                       iCLayer);
 
        const float upperY = y - upperImpact;
        const ESign upperDistSign = T::getDistanceSign(subordinaryHoughBox,
                                                       x, upperY,
                                                       signedDriftLength,
                                                       dxdz, dydz,
                                                       iCLayer);
 
        return ESignUtil::common(lowerDistSign, upperDistSign);
      }
    };
  }
}