AxialTrackCreatorHitLegendre.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/findlets/base/Findlet.h>
 
#include <tracking/trackFindingCDC/legendre/quadtree/AxialHitQuadTreeProcessor.h>
 
namespace Belle2 {
 
  namespace TrackingUtilities {
    class CDCTrack;
    class CDCWireHit;
  }
 
  namespace TrackFindingCDC {

    class AxialTrackCreatorHitLegendre : public
      TrackingUtilities::Findlet<const TrackingUtilities::CDCWireHit* const, TrackingUtilities::CDCTrack> {
 
    private:
      using Super = TrackingUtilities::Findlet<const TrackingUtilities::CDCWireHit* const, TrackingUtilities::CDCTrack>;
 
    public:
      enum class EPass {
        NonCurlers,
        NonCurlersWithIncreasingThreshold,
        FullRange,
        Straight
      };
 
    public:
      AxialTrackCreatorHitLegendre();

      explicit AxialTrackCreatorHitLegendre(EPass pass);

      std::string getDescription() final;

      void exposeParameters(ModuleParamList* moduleParamList, const std::string& prefix) final;

      void initialize() final;

      void apply(const std::vector<const TrackingUtilities::CDCWireHit*>& axialWireHits,
                 std::vector<TrackingUtilities::CDCTrack>& tracks) final;
 
    private:
      std::unique_ptr<AxialHitQuadTreeProcessor> constructQTProcessor(EPass pass);

      using CandidateReceiver = AxialHitQuadTreeProcessor::CandidateReceiver;
      void executeRelaxation(const CandidateReceiver& candidateReceiver,
                             AxialHitQuadTreeProcessor& qtProcessor);
 
    private: // Parameters
      EPass m_pass = EPass::NonCurlers;

      const double m_param_stepScale = 0.75;

      int m_param_minNHits = 10;

      double m_param_precision = 0.00000001;
    };
  }
}