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_param_SVDSpacePointStoreArrayName,
                                      "Name of the SVDSpacePoints Store Array.", m_param_SVDSpacePointStoreArrayName);
 
        moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "minimumUClusterTime"), m_param_minimumUClusterTime,
                                      "Minimum time of the u cluster (in ns).", m_param_minimumUClusterTime);
 
        moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "minimumVClusterTime"), m_param_minimumVClusterTime,
                                      "Minimum time of the v cluster (in ns).", m_param_minimumVClusterTime);
 
        moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "maximumUClusterTime"), m_param_maximumUClusterTime,
                                      "Maximum time of the u cluster (in ns).", m_param_maximumUClusterTime);
 
        moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "maximumVClusterTime"), m_param_maximumVClusterTime,
                                      "Maximum time of the v cluster (in ns).", m_param_maximumVClusterTime);
      };
 
      void initialize() override
      {
        Super::initialize();
        m_storeSpacePoints.isRequired(m_param_SVDSpacePointStoreArrayName);
      };
 
      void beginRun() override
      {
        Super::beginRun();
 
        if (m_BeamSpotDB.isValid()) {
          m_BeamSpot = *m_BeamSpotDB;
          const TVector3& BeamSpotPosition = m_BeamSpot.getIPPosition();
          m_BeamSpotPosition.SetXYZ(BeamSpotPosition.X(), BeamSpotPosition.Y(), BeamSpotPosition.Z());
        } else {
          m_BeamSpotPosition.SetXYZ(0., 0., 0.);
        }
        B2INFO("DATCON uses following BeamSpot: " <<
               m_BeamSpotPosition.X() << ", " << m_BeamSpotPosition.Y() << ", " << m_BeamSpotPosition.Z());
      }
 
      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 (spacePoint.TimeU() >= m_param_minimumUClusterTime and
              spacePoint.TimeV() >= m_param_minimumVClusterTime and
              spacePoint.TimeU() <= m_param_maximumUClusterTime and
              spacePoint.TimeV() <= m_param_maximumVClusterTime) {
            hits.emplace_back(VXDHoughState(&spacePoint, m_BeamSpotPosition));
            spacePoints.emplace_back(&spacePoint);
          }
        }
      };
 
    private:
      std::string m_param_SVDSpacePointStoreArrayName = "SVDSpacePoints";
 
      double m_param_minimumUClusterTime = -25;
      double m_param_minimumVClusterTime = -25;
      double m_param_maximumUClusterTime = 40;
      double m_param_maximumVClusterTime = 40;
 
      StoreArray<SpacePoint> m_storeSpacePoints;
 
      DBObjPtr<BeamSpot> m_BeamSpotDB;
      BeamSpot m_BeamSpot;
      B2Vector3D m_BeamSpotPosition;
    };
 
  }
}