SpacePointLoaderAndPreparer.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 <framework/datastore/StoreArray.h>
#include <framework/core/ModuleParamList.h>
#include <framework/geometry/B2Vector3.h>
#include <framework/database/DBObjPtr.h>
#include <mdst/dbobjects/BeamSpot.h>
#include <tracking/trackFindingCDC/findlets/base/Findlet.h>
#include <tracking/trackFindingCDC/utilities/StringManipulation.h>
#include <tracking/vxdHoughTracking/entities/VXDHoughState.h>
#include <tracking/spacePointCreation/SpacePoint.h>
#include <vxd/geometry/GeoCache.h>
 
#include <string>
#include <vector>
 
namespace Belle2 {
  namespace vxdHoughTracking {
    class SpacePointLoaderAndPreparer : public TrackFindingCDC::Findlet<const SpacePoint*, VXDHoughState> {
      using Super = TrackFindingCDC::Findlet<const SpacePoint*, VXDHoughState>;
 
    public:
      SpacePointLoaderAndPreparer() {};

      void exposeParameters(ModuleParamList* moduleParamList, const std::string& prefix) override
      {
        Super::exposeParameters(moduleParamList, prefix);
 
        moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "SVDSpacePointStoreArrayName"), m_SVDSpacePointStoreArrayName,
                                      "Name of the SVDSpacePoints Store Array.", m_SVDSpacePointStoreArrayName);
 
        moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "useAllSpacePoints"), m_useAllSpacePoints,
                                      "Use all SVDSpacePoints for track finding or only unassigned ones?", m_useAllSpacePoints);
      };

      void initialize() override
      {
        Super::initialize();
        m_storeSpacePoints.isRequired(m_SVDSpacePointStoreArrayName);
      };

      void beginRun() override
      {
        Super::beginRun();
 
        DBObjPtr<BeamSpot> beamSpotDB;
        if (beamSpotDB.isValid()) {
          m_BeamSpotPosition = (*beamSpotDB).getIPPosition();
        } else {
          m_BeamSpotPosition.SetXYZ(0., 0., 0.);
        }
      }

      void apply(std::vector<const SpacePoint*>& spacePoints, std::vector<VXDHoughState>& hits) override
      {
        if (m_storeSpacePoints.getEntries() == 0) return;
 
        hits.reserve(m_storeSpacePoints.getEntries());
        spacePoints.reserve(m_storeSpacePoints.getEntries());
        for (auto& spacePoint : m_storeSpacePoints) {
          if (not m_useAllSpacePoints and spacePoint.getAssignmentState() == true) {
            continue;
          }
          hits.emplace_back(VXDHoughState(&spacePoint, m_BeamSpotPosition));
          spacePoints.emplace_back(&spacePoint);
        }
      };
 
    private:
      std::string m_SVDSpacePointStoreArrayName = "SVDSpacePoints";

      StoreArray<SpacePoint> m_storeSpacePoints;

      B2Vector3D m_BeamSpotPosition;

      bool m_useAllSpacePoints = false;
    };
 
  }
}