SVDRecoHit2D.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2010 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Peter Kvasnicka, Martin Ritter, Moritz Nadler            *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef SVDRECOHIT2D_H_
#define SVDRECOHIT2D_H_
 
#include <vxd/dataobjects/VxdID.h>
#include <svd/dataobjects/SVDTrueHit.h>
#include <svd/dataobjects/SVDCluster.h>
#include <svd/reconstruction/SVDRecoHit.h>
 
// ROOT includes
#include <TMatrixD.h>
 
// GenFit includes
#include <genfit/PlanarMeasurement.h>
#include <genfit/TrackCandHit.h>
 
namespace Belle2 {
  class SVDRecoHit2D: public genfit::PlanarMeasurement {
  public:
    SVDRecoHit2D();
 
    explicit SVDRecoHit2D(const SVDTrueHit* hit, const genfit::TrackCandHit* trackCandHit = nullptr, float sigmaU = -1,
                          float sigmaV = -1);
 
 
    SVDRecoHit2D(VxdID::baseType vxdid, const double u, const double v, double sigmaU = -1, double sigmaV = -1);
 
    SVDRecoHit2D(const SVDCluster& uHit, const SVDCluster& vHit);
 
    SVDRecoHit2D(const SVDRecoHit& uRecoHit, const SVDRecoHit& vRecoHit);
 
    virtual ~SVDRecoHit2D() {}
 
    genfit::AbsMeasurement* clone() const override;
 
    VxdID getSensorID() const { return m_sensorID; }
 
    const SVDTrueHit* getTrueHit() const { return m_trueHit; }
    const SVDCluster* getUCluster() const {return m_uCluster; }
    const SVDCluster* getVCluster() const {return m_vCluster; }
 
    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; }
 
    //virtual genfit::SharedPlanePtr constructPlane(const genfit::StateOnPlane&) const;
    virtual std::vector<genfit::MeasurementOnPlane*> constructMeasurementsOnPlane(const genfit::StateOnPlane& state) const override;
 
  private:
 
    enum { HIT_DIMENSIONS = 2  };
 
    unsigned short m_sensorID; 
    const SVDTrueHit* m_trueHit; 
    const SVDCluster* m_uCluster; 
    const SVDCluster* m_vCluster; 
    float m_energyDep; 
    //float m_energyDepError; /**< error in dep. energy.*/
 
    void setDetectorPlane();
 
    TVectorD applyPlanarDeformation(TVectorD rawHit, std::vector<double> planarParameters, const genfit::StateOnPlane& state) const;
 
    ClassDefOverride(SVDRecoHit2D, 7)
  };
 
} // namespace Belle2
 
#endif /* SVDRECOHIT_H_ */