#ifndef __JPOINT4E__
#define __JPOINT4E__

#include "JFit/JPoint4D.hh"


/**
 * \author vcarretero
 */

namespace JFIT {}
namespace JPP { using namespace JFIT; }

namespace JFIT {

  using JFIT::JPoint4D;
    

  /**
   * Data structure for vertex fit.
   */
  class JPoint4E :
    public JPoint4D
  {
  public:
    /**
     * Default constructor.
     */
    JPoint4E() :
      JPoint4D(),
      __E(0.0)
    {}


    /**
     * Constructor
     *
     * \param vertex      vertex position and time
     * \param E           issotropic energy
     */
    JPoint4E(const JPoint4D& vertex,
		      const double     E) :
      JPoint4D(vertex),
      __E(E)
    {}

    /**
     * Get energy.
     *
     * \return               energy
     */
    inline double getE() const
    {
      return __E;
    }
    
    /**
     * Set energy.
     *
     * \param  E             energy
     */
    inline void setE(const double E)
    {
      __E = E;
    }
   
    typedef double JPoint4E::*parameter_type;

    static parameter_type pX() { return &JPoint4E::__x; }
    static parameter_type pY() { return &JPoint4E::__y; }
    static parameter_type pZ() { return &JPoint4E::__z; }
    static parameter_type pT() { return &JPoint4E::__t; }
    static parameter_type pE() { return &JPoint4E::__E; }
  protected:
    double __E;
  };
}
#endif