SVDRecoHit.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 SVDRECOHIT_H_
#define SVDRECOHIT_H_
 
#include <vxd/dataobjects/VxdID.h>
#include <svd/dataobjects/SVDTrueHit.h>
#include <svd/dataobjects/SVDCluster.h>
 
// ROOT includes
#include <TMatrixD.h>
 
// GenFit includes
#include <genfit/PlanarMeasurement.h>
#include <genfit/HMatrixU.h>
#include <genfit/HMatrixV.h>
#include <genfit/TrackCandHit.h>
 
namespace Belle2 {
  class SVDRecoHit: public genfit::PlanarMeasurement {
  public:
    SVDRecoHit();
 
    SVDRecoHit(const SVDTrueHit* hit, bool uDirection, float sigma = -1);
 
    explicit SVDRecoHit(const SVDCluster* hit, const genfit::TrackCandHit* trackCandHit = nullptr);
 
    virtual ~SVDRecoHit() {}
 
    genfit::AbsMeasurement* clone() const override;
 
    VxdID getSensorID() const { return m_sensorID; }
 
    const SVDTrueHit* getTrueHit() const { return m_trueHit; }
 
    const SVDCluster* getCluster() const { return m_cluster; }
 
    bool isU() const { return m_isU; }
 
    float getPosition() const { return rawHitCoords_(0); }
 
    float getPositionVariance() const { return rawHitCov_(0, 0); }
 
    float getEnergyDep() const { return m_energyDep; }
 
    float getRotation() const { return m_rotationPhi; }
 
    //virtual genfit::SharedPlanePtr constructPlane(const genfit::StateOnPlane&) const;
    virtual std::vector<genfit::MeasurementOnPlane*> constructMeasurementsOnPlane(const genfit::StateOnPlane& state) const override;
 
    // TODO: use HMatrixPhi for wedge sensors instead of rotating the plane!
    virtual const genfit::AbsHMatrix* constructHMatrix(const genfit::AbsTrackRep*) const override { if (m_isU) return new genfit::HMatrixU(); else return new genfit::HMatrixV(); }
 
  private:
 
    enum { HIT_DIMENSIONS = 1  };
 
    unsigned short m_sensorID; 
    const SVDTrueHit* m_trueHit; 
 
    const SVDCluster* m_cluster; 
    bool m_isU; 
    float m_energyDep; 
    //float m_energyDepError; /**< error in dep. energy.*/
    float m_rotationPhi; 
    void setDetectorPlane();
 
    TVectorD applyPlanarDeformation(TVectorD rawHit, std::vector<double> planarParameters, const genfit::StateOnPlane& state) const;
 
    ClassDefOverride(SVDRecoHit, 7)
  };
 
} // namespace Belle2
 
#endif /* SVDRECOHIT_H_ */