VXDHoughState.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/spacePointCreation/SpacePoint.h>
#include <tracking/trackFindingCDC/ca/AutomatonCell.h>
#include <vxd/dataobjects/VxdID.h>
 
namespace Belle2 {
  class SpacePoint;
  namespace vxdHoughTracking {
 
 
    class VXDHoughState {
    public:
      VXDHoughState() = default;
      VXDHoughState(VXDHoughState&&) = default;
      VXDHoughState(const VXDHoughState&) = delete;
      VXDHoughState& operator=(const VXDHoughState&) = delete;
      VXDHoughState& operator=(VXDHoughState&&) = default;
      ~VXDHoughState() = default;
 
      explicit VXDHoughState(const SpacePoint* hit, const B2Vector3D& BeamSpotPosition) : m_hit(hit)
      {
        m_dataCache.sensorID = hit->getVxdID();
        m_dataCache.layer = hit->getVxdID().getLayerNumber();
        m_dataCache.ladder = hit->getVxdID().getLadderNumber();
        const VXD::SensorInfoBase& sensorInfo = VXD::GeoCache::getInstance().getSensorInfo(hit->getVxdID());
        m_dataCache.x = hit->X();
        m_dataCache.y = hit->Y();
        m_dataCache.z = hit->Z();
        const double conformalTransform = 2. / ((hit->X() - BeamSpotPosition.X()) * (hit->X() - BeamSpotPosition.X()) +
                                                (hit->Y() - BeamSpotPosition.Y()) * (hit->Y() - BeamSpotPosition.Y()));
        m_dataCache.xConformal = hit->X() * conformalTransform;
        m_dataCache.yConformal = hit->Y() * conformalTransform;
        m_dataCache.localNormalizedu = hit->getNormalizedLocalU();
        m_dataCache.localNormalizedv = hit->getNormalizedLocalV();
        m_dataCache.sensorCenterPhi = sensorInfo.pointToGlobal(ROOT::Math::XYZVector(0., 0., 0.), true).Phi();
        m_dataCache.phi = hit->getPosition().Phi();
        m_dataCache.theta = hit->getPosition().Theta();
        m_dataCache.uTime = hit->TimeU();
        m_dataCache.vTime = hit->TimeV();
      }
 
      const SpacePoint* getHit() const { return m_hit; }
 
      TrackFindingCDC::AutomatonCell& getAutomatonCell() { return m_automatonCell; }
 
 
      struct DataCache {
        float x = std::numeric_limits<float>::quiet_NaN();
        float y = std::numeric_limits<float>::quiet_NaN();
        float z = std::numeric_limits<float>::quiet_NaN();
        float xConformal = std::numeric_limits<float>::quiet_NaN();
        float yConformal = std::numeric_limits<float>::quiet_NaN();
        float localNormalizedu = std::numeric_limits<float>::quiet_NaN();
        float localNormalizedv = std::numeric_limits<float>::quiet_NaN();
        float sensorCenterPhi = std::numeric_limits<float>::quiet_NaN();
        float phi = std::numeric_limits<float>::quiet_NaN();
        float theta = std::numeric_limits<float>::quiet_NaN();
        float uTime = std::numeric_limits<float>::quiet_NaN();
        float vTime = std::numeric_limits<float>::quiet_NaN();
        unsigned short layer = std::numeric_limits<unsigned short>::quiet_NaN();
        unsigned short ladder = std::numeric_limits<unsigned short>::quiet_NaN();
        VxdID sensorID = 0;
      };
      const DataCache getDataCache() const { return m_dataCache; }
 
    protected:
      DataCache m_dataCache;
 
    private:
      const SpacePoint* m_hit = nullptr;
      TrackFindingCDC::AutomatonCell m_automatonCell;
    };
 
  }
}