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  };
 
    unsigned short m_sensorID; 
    const PXDTrueHit* m_trueHit; 
    const PXDCluster* m_cluster; 
    float m_energyDep; 
    //float m_energyDepError; /**< error in dep. energy.*/
 
    void setDetectorPlane();
 
    TVectorD applyPlanarDeformation(TVectorD hitCoords, std::vector<double> planarParameters, const genfit::StateOnPlane& state) const;
 
    ClassDefOverride(PXDRecoHit, 8)
  };
 
} // namespace Belle2