CKFState.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Nils Braun, Christian Wessel                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include <tracking/trackFindingCDC/ca/AutomatonCell.h>
 
#include <genfit/MeasuredStateOnPlane.h>
#include <framework/logging/Logger.h>
 
#include <vxd/dataobjects/VxdID.h>
 
namespace Belle2 {
  template<class ASeed, class AHit>
  class CKFState {
  public:
    using Seed = ASeed;
    using Hit = AHit;
 
    CKFState() = default;
    CKFState(CKFState&&) = default;
    CKFState(const CKFState&) = delete;
    CKFState& operator=(const CKFState&) = delete;
    CKFState& operator=(CKFState&&) = default;
    ~CKFState() = default;
 
    explicit CKFState(const Seed* seed) : m_seed(seed)
    {
    }
 
    explicit CKFState(const Hit* hit) : m_hit(hit)
    {
    }
 
    const Seed* getSeed() const
    {
      return m_seed;
    }
 
    const Hit* getHit() const
    {
      return m_hit;
    }
 
    double getChi2() const
    {
      B2ASSERT("Chi2 is not set", isFitted());
      return m_chi2;
    }
 
    void setChi2(double chi2)
    {
      m_chi2 = chi2;
      m_isFitted = true;
    }
 
    void setMeasuredStateOnPlane(const genfit::MeasuredStateOnPlane& mSoP)
    {
      m_measuredStateOnPlane = mSoP;
      m_hasMSoP = true;
    }
 
    const genfit::MeasuredStateOnPlane& getMeasuredStateOnPlane() const
    {
      B2ASSERT("You are asking for an invalid variable!", mSoPSet());
      return m_measuredStateOnPlane;
    }
 
    bool isFitted() const
    {
      return m_isFitted;
    }
 
    bool mSoPSet() const
    {
      return m_hasMSoP;
    }
 
    void reset()
    {
      m_isFitted = false;
      m_hasMSoP = false;
    }
 
    TrackFindingCDC::AutomatonCell& getAutomatonCell()
    {
      return m_automatonCell;
    }
 
    struct stateCache {
      bool isHitState = false;
      VxdID sensorID = VxdID(0);
      unsigned short geoLayer = std::numeric_limits<unsigned short>::quiet_NaN();
      unsigned short ladder = std::numeric_limits<unsigned short>::quiet_NaN();
      double sensorCenterPhi = std::numeric_limits<double>::quiet_NaN();
      double phi = std::numeric_limits<double>::quiet_NaN();
      double theta = std::numeric_limits<double>::quiet_NaN();
      double localNormalizedu = std::numeric_limits<double>::quiet_NaN();
      double localNormalizedv = std::numeric_limits<double>::quiet_NaN();
    };
 
  protected:
    stateCache m_stateCache;
 
 
  private:
    const Seed* m_seed = nullptr;
    const Hit* m_hit = nullptr;
    double m_chi2 = 0;
    bool m_isFitted = false;
    bool m_hasMSoP = false;
    genfit::MeasuredStateOnPlane m_measuredStateOnPlane;
    TrackFindingCDC::AutomatonCell m_automatonCell;
  };
}