P4IPtCotThPhiM.h

Go to the documentation of this file.

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

 
 
class P4IPtCotThPhiM : public virtual P4IPtCotThPhiMBase
 
{
 public:
  P4IPtCotThPhiM(const double iPt, const double cotTh, const double phi, const double m);

  P4IPtCotThPhiM();

  P4IPtCotThPhiM(const CLHEP::HepLorentzVector & theHlv);
  
  P4IPtCotThPhiM(const I4Momentum & theI4M );

  P4IPtCotThPhiM(const I4Momentum  * const pI4M );

  virtual double iPt() const;
  virtual double cotTh() const;
  virtual double phi() const;
  virtual double m() const;

  virtual void setIPt(double theIPt) ;
  virtual void setCotTh(double theCotTh) ;
  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_iPt;
  double m_cotTh;
  double m_phi;
  double m_m;
  
};
 
inline P4IPtCotThPhiM::P4IPtCotThPhiM() :
  P4IPtCotThPhiMBase( ),
  m_iPt   ( 0./CLHEP::GeV ),
  m_cotTh ( 0. ),
  m_phi   ( 0. ),
  m_m     ( 0.*CLHEP::GeV )
{}
 
inline P4IPtCotThPhiM::P4IPtCotThPhiM(const CLHEP::HepLorentzVector & theHlv)
{
  this->set4Mom(theHlv);
}
 
inline P4IPtCotThPhiM::P4IPtCotThPhiM(const I4Momentum & theI4Mom ) 
{this->set4Mom(theI4Mom); }
 
inline P4IPtCotThPhiM::P4IPtCotThPhiM(const I4Momentum  * const  theI4Mom )
{this->set4Mom(*theI4Mom); }
 
 
inline double P4IPtCotThPhiM::iPt() const 
{ return m_iPt;}
 
inline double P4IPtCotThPhiM::cotTh() const
{ return m_cotTh;}
 
inline double P4IPtCotThPhiM::phi() const
{ return m_phi;}
 
inline double P4IPtCotThPhiM::m() const
{ return m_m;}
 
inline void P4IPtCotThPhiM::setIPt(double theIPt) 
{ m_iPt=theIPt;}
 
inline void P4IPtCotThPhiM::setCotTh(double theCotTh) 
{ m_cotTh=theCotTh;}
 
inline void P4IPtCotThPhiM::setPhi(double thePhi) 
{ m_phi=thePhi;}
 
inline void P4IPtCotThPhiM::setM(double theM) 
{ m_m=theM;}
 
 
inline void P4IPtCotThPhiM::set4Mom(const I4Momentum & theI4Mom )
{
  m_iPt=theI4Mom.iPt();
  m_cotTh=theI4Mom.cotTh();
  m_phi=theI4Mom.phi();
  m_m=theI4Mom.m();
}
 
inline void P4IPtCotThPhiM::set4Mom(const I4Momentum * const theI4Mom )
{
  this->set4Mom(*theI4Mom);
}
 
 
inline void P4IPtCotThPhiM::set4Mom(const CLHEP::HepLorentzVector & theHlv)
{
  m_iPt=1./theHlv.perp();
  m_cotTh=theHlv.pz()*m_iPt;
  m_phi=theHlv.phi();
  m_m=theHlv.m();
}
 
 
#endif