ClusterCreator.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/filters/wireHitRelation/WholeWireHitRelationFilter.h>
 
#include <tracking/trackingUtilities/eventdata/segments/CDCWireHitCluster.h>
#include <tracking/trackingUtilities/eventdata/hits/CDCWireHit.h>
 
#include <tracking/trackingUtilities/ca/Clusterizer.h>
 
#include <tracking/trackingUtilities/filters/base/RelationFilterUtil.h>
 
#include <tracking/trackingUtilities/utilities/Algorithms.h>
 
#include <vector>
#include <string>
 
namespace Belle2 {
  namespace TrackingUtilities {
    class CDCWireHit;
  }
 
  namespace TrackFindingCDC {

    template <class AWireHitRelationFilter = WholeWireHitRelationFilter>
    class ClusterCreator : public TrackingUtilities::Findlet<TrackingUtilities::CDCWireHit&, TrackingUtilities::CDCWireHitCluster> {
 
    private:
      using Super = TrackingUtilities::Findlet<TrackingUtilities::CDCWireHit, TrackingUtilities::CDCWireHitCluster>;
 
    public:
      ClusterCreator()
      {
        this->addProcessingSignalListener(&m_wireHitRelationFilter);
      }

      std::string getDescription() final
      {
        return "Groups the wire hits into super by expanding the given wire hit relation";
      }
 
    public:
      void apply(std::vector<TrackingUtilities::CDCWireHit>& inputWireHits,
                 std::vector<TrackingUtilities::CDCWireHitCluster>& outputClusters) final
      {
 
        // Obtain the wire hits as pointers.
        const std::vector<TrackingUtilities::CDCWireHit*> wireHitPtrs = TrackingUtilities::as_pointers<TrackingUtilities::CDCWireHit>
            (inputWireHits);
 
        // create the neighborhood
        m_wireHitRelations.clear();
        TrackingUtilities::RelationFilterUtil::appendUsing(m_wireHitRelationFilter, wireHitPtrs, m_wireHitRelations);
 
        B2ASSERT("Expect wire hit neighborhood to be symmetric ",
                 TrackingUtilities::WeightedRelationUtil<TrackingUtilities::CDCWireHit>::areSymmetric(m_wireHitRelations));
 
        // Compose the clusters
        m_wirehitClusterizer.apply(wireHitPtrs, m_wireHitRelations, outputClusters);
      }
 
    private:
      TrackingUtilities::Clusterizer<TrackingUtilities::CDCWireHit, TrackingUtilities::CDCWireHitCluster> m_wirehitClusterizer;

      std::vector<TrackingUtilities::WeightedRelation<TrackingUtilities::CDCWireHit>> m_wireHitRelations;

      AWireHitRelationFilter m_wireHitRelationFilter;
    };
  }
}