PXDRecoHit.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 <vxd/dataobjects/VxdID.h>
#include <pxd/dataobjects/PXDTrueHit.h>
#include <pxd/dataobjects/PXDCluster.h>
 
// ROOT includes
#include <TMatrixD.h>
 
// GenFit includes
#include <genfit/PlanarMeasurement.h>
#include <genfit/HMatrixUV.h>
#include <genfit/TrackCandHit.h>
 
namespace Belle2 {
  class PXDRecoHit: public genfit::PlanarMeasurement {
  public:

    PXDRecoHit();

    explicit PXDRecoHit(const PXDTrueHit* hit, const genfit::TrackCandHit* trackCandHit = NULL, float sigmaU = -1, float sigmaV = -1);

    explicit PXDRecoHit(const PXDCluster* hit, float sigmaU, float sigmaV, float covUV);

    explicit PXDRecoHit(const PXDCluster* hit, const genfit::TrackCandHit* trackCandHit = NULL);

    genfit::AbsMeasurement* clone() const override;

    /* This method allows to provide hit position dependent on track direction. */
    virtual std::vector<genfit::MeasurementOnPlane*> constructMeasurementsOnPlane(const genfit::StateOnPlane& state) const override;

    VxdID getSensorID() const { return m_sensorID; }

    const PXDTrueHit* getTrueHit() const { return m_trueHit; }
    const PXDCluster* getCluster() const { return m_cluster; }

    float getU() const { return rawHitCoords_(0); }
    float getV() const { return rawHitCoords_(1); }

    float getUVariance() const { return rawHitCov_(0, 0); }
    float getVVariance() const { return rawHitCov_(1, 1); }
    float getUVCov() const { return rawHitCov_(0, 1); }

    float getEnergyDep() const { return m_energyDep; }

    //float getEnergyDepError() const { return m_energyDepError; }

    float getShapeLikelyhood(const genfit::StateOnPlane& state) const;

    virtual const genfit::AbsHMatrix* constructHMatrix(const genfit::AbsTrackRep*) const override { return new genfit::HMatrixUV(); };
 
  private:
 
    enum { HIT_DIMENSIONS = 2  };

    const PXDTrueHit* m_trueHit; 
    const PXDCluster* m_cluster; 
    float m_energyDep; 
    //float m_energyDepError; /**< error in dep. energy.*/
    unsigned short m_sensorID;

    void setDetectorPlane();

    TVectorD applyPlanarDeformation(TVectorD hitCoords, std::vector<double> planarParameters, const genfit::StateOnPlane& state) const;
 
    ClassDefOverride(PXDRecoHit, 9)
  };

} // namespace Belle2