P4EEtaPhiM.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef FOURMOM_P4EETAPHIM_H
#define FOURMOM_P4EETAPHIM_H
 
#include "FourMom/P4EEtaPhiMBase.h"
//#include "SealBase/DebugAids.h"
#include <cmath>
 
// CLHEP includes
#include "CLHEP/Units/SystemOfUnits.h"

 
class P4EEtaPhiM : public virtual P4EEtaPhiMBase
{
 public:

  P4EEtaPhiM(const double e, const double eta, const double phi, const double m);

  P4EEtaPhiM();

  P4EEtaPhiM(const CLHEP::HepLorentzVector & theHlv);

  P4EEtaPhiM(const I4Momentum & theI4M );

  P4EEtaPhiM(const I4Momentum  * const pI4M );
  
  
  virtual double e() const;
  virtual double eta() const;
  virtual double phi() const;
  virtual double m() const;

  virtual void setE (double theE);
  virtual void setEta (double theEta);
  virtual void setPhi (double thePhi);
  virtual void setM (double theM);
  
  virtual void set4Mom (const I4Momentum & theI4Mom );  
  virtual void set4Mom (const I4Momentum * const theI4Mom );  
  virtual void set4Mom (const CLHEP::HepLorentzVector & theHlv );  
  
  virtual const I4MomentumError* errors() const;
 
 protected:
 
 
 
  double m_e;
  double m_eta;
  double m_phi;
  double m_m;
  
};
 
inline P4EEtaPhiM::P4EEtaPhiM() :
  P4EEtaPhiMBase(),
  m_e   ( 0.*CLHEP::GeV ),
  m_eta ( 0. ),
  m_phi ( 0. ),
  m_m   ( 0.*CLHEP::GeV )
{}
 
inline P4EEtaPhiM::P4EEtaPhiM(const CLHEP::HepLorentzVector & theHlv)
{
  this->set4Mom(theHlv);
}
 
inline P4EEtaPhiM::P4EEtaPhiM(const I4Momentum & theI4Mom ) 
{this->set4Mom(theI4Mom); }
 
inline P4EEtaPhiM::P4EEtaPhiM(const I4Momentum  * const  theI4Mom )
{this->set4Mom(*theI4Mom); }
 
 
inline double P4EEtaPhiM::e() const 
{ return m_e;}
 
inline double P4EEtaPhiM::eta() const
{ return m_eta;}
 
inline double P4EEtaPhiM::phi() const
{ return m_phi;}
 
inline double P4EEtaPhiM::m() const
{ return m_m;}
 
inline void P4EEtaPhiM::setE(double theE) 
{ m_e = theE;}
 
inline void P4EEtaPhiM::setEta(double theEta)
{ m_eta = theEta;}
 
inline void P4EEtaPhiM::setPhi(double thePhi)
{ m_phi = thePhi;}
 
inline void P4EEtaPhiM::setM(double theM)
{ m_m = theM;}
 
inline void P4EEtaPhiM::set4Mom(const I4Momentum & theI4Mom )
{ 
  m_e=theI4Mom.e();
  m_eta=theI4Mom.eta();
  m_phi=theI4Mom.phi();
  m_m=theI4Mom.m();
}
 
inline void P4EEtaPhiM::set4Mom(const I4Momentum * const theI4Mom )
{ 
  this->set4Mom(*theI4Mom);
}
 
inline void P4EEtaPhiM::set4Mom(const CLHEP::HepLorentzVector & theHlv)
{
  m_e=theHlv.e();
  m_eta=theHlv.eta();
  //FIXME protect against negative energy
  //assert(m_e >= 0,"P4EEtaPhiM::set4Mom cannot have negative energy");
  //assert(m_e >= 0);
  // FIXME of the FIXME in fact it is not necessary to prtoect against negative energy
  // and besides Seal assert does not work
  //  ASSERT( m_e >= 0 );
  m_phi=theHlv.phi();
  m_m=theHlv.m();
 
 
}
 
 
#endif