ClusterCreator.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/findlets/base/Findlet.h>
 
#include <tracking/trackFindingCDC/filters/wireHitRelation/WholeWireHitRelationFilter.h>
 
#include <tracking/trackFindingCDC/eventdata/segments/CDCWireHitCluster.h>
#include <tracking/trackFindingCDC/eventdata/hits/CDCWireHit.h>
 
#include <tracking/trackFindingCDC/ca/Clusterizer.h>
 
#include <tracking/trackFindingCDC/filters/base/RelationFilterUtil.h>
 
#include <tracking/trackFindingCDC/utilities/Algorithms.h>
 
#include <vector>
#include <string>
 
namespace Belle2 {
  namespace TrackFindingCDC {
 
    template <class AWireHitRelationFilter = WholeWireHitRelationFilter>
    class ClusterCreator : public Findlet<CDCWireHit&, CDCWireHitCluster> {
 
    private:
      using Super = Findlet<CDCWireHit, 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<CDCWireHit>& inputWireHits,
                 std::vector<CDCWireHitCluster>& outputClusters) final {
 
        // Obtain the wire hits as pointers.
        const std::vector<CDCWireHit*> wireHitPtrs = as_pointers<CDCWireHit>(inputWireHits);
 
        // create the neighborhood
        m_wireHitRelations.clear();
        RelationFilterUtil::appendUsing(m_wireHitRelationFilter, wireHitPtrs, m_wireHitRelations);
 
        B2ASSERT("Expect wire hit neighborhood to be symmetric ",
        WeightedRelationUtil<CDCWireHit>::areSymmetric(m_wireHitRelations));
 
        // Compose the clusters
        m_wirehitClusterizer.apply(wireHitPtrs, m_wireHitRelations, outputClusters);
      }
 
    private:
      Clusterizer<CDCWireHit, CDCWireHitCluster> m_wirehitClusterizer;
 
      std::vector<WeightedRelation<CDCWireHit>> m_wireHitRelations;
 
      AWireHitRelationFilter m_wireHitRelationFilter;
    };
  }
}