VXDTrueHit.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
#ifndef VXD_DATAOBJECTS_VXDTRUEHIT_H
#define VXD_DATAOBJECTS_VXDTRUEHIT_H
 
#include <vxd/dataobjects/VxdID.h>
#include <framework/datastore/RelationsObject.h>
#include <Math/Vector3D.h>
#include <algorithm>
 
namespace Belle2 {
  class VXDTrueHit : public RelationsObject {
 
  public:
    VXDTrueHit(): RelationsObject(), m_positionEntry(),
      m_positionMidPoint(), m_positionExit(), m_momentumEntry(),
      m_momentumMidPoint(), m_momentumExit(), m_energyDep(0), m_globalTime(0), m_sensorID(0)
    {}

    VXDTrueHit(
      VxdID sensorID, float* posEntry, float* posMidPoint, float* posExit,
      float* momEntry, float* momMidPoint, float* momExit, float energyDep, float globalTime):
      RelationsObject(), m_energyDep(energyDep), m_globalTime(globalTime), m_sensorID(sensorID)
    {
      std::copy_n(posEntry, 3, m_positionEntry);
      std::copy_n(posMidPoint, 3, m_positionMidPoint);
      std::copy_n(posExit, 3, m_positionExit);
      std::copy_n(momEntry, 3, m_momentumEntry);
      std::copy_n(momMidPoint, 3, m_momentumMidPoint);
      std::copy_n(momExit, 3, m_momentumExit);
    }

    VxdID getSensorID() const { return m_sensorID; }
    unsigned int getRawSensorID() const {return m_sensorID; }
    float getU() const { return m_positionMidPoint[0]; }
    float getV() const { return m_positionMidPoint[1]; }
    float getW() const { return m_positionMidPoint[2]; }
    float getEntryU() const { return m_positionEntry[0]; }
    float getEntryV() const { return m_positionEntry[1]; }
    float getEntryW() const { return m_positionEntry[2]; }
    float getExitU() const { return m_positionExit[0]; }
    float getExitV() const { return m_positionExit[1]; }
    float getExitW() const { return m_positionExit[2]; }
    float getEnergyDep() const { return m_energyDep; }
    float getGlobalTime() const { return m_globalTime; }
    ROOT::Math::XYZVector getMomentum() const { return ROOT::Math::XYZVector(m_momentumMidPoint[0], m_momentumMidPoint[1], m_momentumMidPoint[2]); }
    ROOT::Math::XYZVector getEntryMomentum() const { return ROOT::Math::XYZVector(m_momentumEntry[0], m_momentumEntry[1], m_momentumEntry[2]);}
    ROOT::Math::XYZVector getExitMomentum() const { return ROOT::Math::XYZVector(m_momentumExit[0], m_momentumExit[1], m_momentumExit[2]); }

    virtual void shiftInTime(float delta) { m_globalTime += delta; }
 
  protected:
    float m_positionEntry[3];
    float m_positionMidPoint[3];
    float m_positionExit[3];
    float m_momentumEntry[3];
    float m_momentumMidPoint[3];
    float m_momentumExit[3];
    float m_energyDep;
    float m_globalTime;
    int m_sensorID;
 
    ClassDef(VXDTrueHit, 7)
  };

} // end namespace Belle2
 
#endif