{

public:

/// \brief Default constructor

Hit() = default;

/// \brief Hit constructor

///

/// Fully defining information of the PHOS Hit (position,

/// momentum, energy, track, ...)

///

/// \param trackID Index of the track entered PHOS

/// \param detID ID of the detector segment

/// \param pos Position vector of the Hit

/// \param mom Momentum vector for the particle at the Hit

/// \param initialEnergy Energy of the primary particle enering the EMCAL

/// \param tof Time of the hit

/// \param length Length of the segment

Hit(Int_t trackID, Int_t detID, const math_utils::Point3D<float>& pos, const math_utils::Vector3D<float>& mom, Double_t totE, Double_t tof,

Double_t eLoss)

: o2::BasicXYZEHit<float>(pos.X(), pos.Y(), pos.Z(), tof, eLoss, trackID, detID),

mPvector(mom),

mInitialEnergy(totE)

{

}

Hit& operator=(const Hit& hit) = default;

/// \brief Check whether the points are from the same SuperParent and in the same detector volume

/// \return True if points are the same (origin and detector), false otherwise

Bool_t operator==(const Hit& rhs) const;

/// \brief Sorting points according to parent particle and detector volume

/// \return True if this Hit is smaller, false otherwise

Bool_t operator<(const Hit& rhs) const;

/// \brief Adds energy loss from the other Hit to this Hit

/// \param rhs phos::Hit to add to this Hit

/// \return This Hit with the summed energy loss

Hit& operator+=(const Hit& rhs);

/// \brief Creates a new Hit based on this Hit but adding the energy loss of the right hand side

/// \param

/// \return New Hit based on this Hit

Hit operator+(const Hit& rhs) const;

/// \brief Destructor

~Hit() = default;

/// \brief Get the initial energy of the primary particle entering PHOS

/// \return Energy of the primary particle entering PHOS

Double_t getInitialEnergy() const { return mInitialEnergy; }

void addEnergyLoss(Double_t eloss) { SetEnergyLoss(GetEnergyLoss() + eloss); }

/// \brief Writing Hit information to an output stream;

/// \param stream target output stream

void PrintStream(std::ostream& stream) const;

private:

math_utils::Vector3D<float> mPvector; // Momentum Vector

Double32_t mInitialEnergy; // Energy of the parent particle that entered the PHOS front surface

ClassDefNV(Hit, 1);