ExtHit.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2010 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Leo Piilonen                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <framework/datastore/RelationsObject.h>
#include <framework/gearbox/Const.h>
#include <TVector3.h>
#include <TMatrixDSym.h>
#include <G4ThreeVector.hh>
#include <G4ErrorSymMatrix.hh>
 
namespace Belle2 {
  enum ExtHitStatus { EXT_FIRST = -1, EXT_ENTER, EXT_EXIT, EXT_STOP, EXT_ESCAPE,
                      EXT_ECLCROSS, EXT_ECLDL, EXT_ECLNEAR
                    };
 
  class ExtHit : public RelationsObject {
 
  public:
 
    ExtHit();
 
    ExtHit(int pdgCode, Const::EDetector detector, int element, ExtHitStatus status, bool backwardPropagation, double t,
           const TVector3& r, const TVector3& p, const TMatrixDSym& e);
 
    ExtHit(int pdgCode, Const::EDetector detector, int element, ExtHitStatus status, bool backwardPropagation, double t,
           const G4ThreeVector& r, const G4ThreeVector& p, const G4ErrorSymMatrix& e);
 
    ExtHit(const ExtHit&);
 
    ExtHit& operator=(const ExtHit&);
 
    virtual ~ExtHit() {}
 
    void setPDGCode(int pdgCode) { m_PdgCode = pdgCode; }
 
    void setDetectorID(Const::EDetector detectorID) { m_DetectorID = detectorID; }
 
    void setCopyID(int copyID) { m_CopyID = copyID; }
 
    void setBackwardPropagation(bool backwardPropagation) { m_BackwardPropagation = backwardPropagation; }
 
    void setStatus(ExtHitStatus extHitStatus) { m_Status = extHitStatus; }
 
    void setTOF(double tof) { m_TOF = tof; }
 
    void setPosition(const TVector3& position)
    {
      m_Position[0] = position.x();
      m_Position[1] = position.y();
      m_Position[2] = position.z();
    }
 
    void setMomentum(const TVector3& momentum)
    {
      m_Momentum[0] = momentum.x();
      m_Momentum[1] = momentum.y();
      m_Momentum[2] = momentum.z();
    }
 
    void setCovariance(double covArray[21])
    {
      for (int k = 0; k < 21; ++k) {
        m_Covariance[k] = covArray[k];
      }
    }
 
    int getPdgCode() const { return m_PdgCode; }
 
    Const::EDetector getDetectorID() const { return m_DetectorID; }
 
    int getCopyID() const { return m_CopyID; }
 
    ExtHitStatus getStatus() const { return m_Status; }
 
    bool isBackwardPropagated() const { return m_BackwardPropagation; }
 
    double getTOF() const { return m_TOF; }
 
    double getLength() const { return m_TOF; }
 
    TVector3 getPosition() const { return TVector3(m_Position[0], m_Position[1], m_Position[2]); }
 
    TVector3 getMomentum() const { return TVector3(m_Momentum[0], m_Momentum[1], m_Momentum[2]); }
 
    TMatrixDSym getCovariance() const;
 
    void update(ExtHitStatus status, double t,
                const G4ThreeVector& r, const G4ThreeVector& p, const G4ErrorSymMatrix& e);
 
    double getErrorTheta() const;
 
    double getErrorPhi() const;
 
  private:
 
    double getPolarCovariance(int i) const;
 
    int m_PdgCode;
 
    Const::EDetector m_DetectorID;
 
    int m_CopyID;
 
    ExtHitStatus m_Status;
 
    bool m_BackwardPropagation;
 
    float m_TOF;
 
    float m_Position[3];
 
    float m_Momentum[3];
 
    float m_Covariance[21];
 
    ClassDef(ExtHit, 8);
 
  };
}