ReconstructedParticle.cxx

Go to the documentation of this file.

// ReconstructedParticle implementation
//

#include "AtlfastEvent/ReconstructedParticle.h"
#include "AtlfastEvent/IAOOvisitor.h"

namespace Atlfast {
  ReconstructedParticle::ReconstructedParticle(const int pdg_id, 
                                               const HepLorentzVector& vec,
                                               const HepMC::GenParticle* 
                                               progenitor ): 
    IKinematic(),
    AODNavigationImp(),m_pdg_id(pdg_id), 
    m_momentum(vec), m_truth(progenitor),
    m_halo(0.),
    m_halovector(HepLorentzVector(0.,0.,0.,0.)) {
  }
  ReconstructedParticle::ReconstructedParticle(const ReconstructedParticle& 
                                               src 
                                               ): 
    Atlfast::IAOO(src),
    INavigable(src),
    I4Momentum(src),
    INavigable4Momentum(src),
    Atlfast::IKinematic(src),
    Atlfast::IAODNavigation(src),
    AODNavigationImp( src ){
    m_pdg_id     = src.m_pdg_id;
    m_momentum   = src.m_momentum;
    m_truth      = src.m_truth;
    m_halo       = src.m_halo;
    m_halovector = src.m_halovector;
  }
  ReconstructedParticle::ReconstructedParticle(
                                               const ReconstructedParticle* 
                                               src 
                                               ):AODNavigationImp( *src ){
    m_pdg_id     = src->m_pdg_id; 
    m_momentum   = src->m_momentum; 
    m_truth      = src->m_truth;
    m_halo       = src->m_halo;
    m_halovector = src->m_halovector;
  }
  ReconstructedParticle& 
  ReconstructedParticle::operator=(const ReconstructedParticle& src ){

    if(this == &src) return *this;
    
    AODNavigationImp::operator=( src );
    m_pdg_id     = src.m_pdg_id; 
    m_momentum   = src.m_momentum; 
    m_truth      = src.m_truth;
    m_halo       = src.m_halo;
    m_halovector = src.m_halovector;
    return *this;
  }
  
  
  ReconstructedParticle::ReconstructedParticle():
    AODNavigationImp(),
    m_pdg_id(0), 
    m_momentum(HepLorentzVector(0.0,0.0,0.0,0.0)), 
    m_truth(0),
    m_halo(0.),
    m_halovector(HepLorentzVector(0.,0.,0.,0.)){
  }
  HepLorentzVector ReconstructedParticle::momentum() const {
    return m_momentum;
  } 
  double  ReconstructedParticle::eta() const { 
    return m_momentum.pseudoRapidity() ;
  }
  double  ReconstructedParticle::phi() const { 
    return m_momentum.phi() ;  
  } 
  double  ReconstructedParticle::pT()  const {
    return m_momentum.perp() ;
  }
  double  ReconstructedParticle::pt()  const {
    return m_momentum.perp() ;
  }
  double  ReconstructedParticle::eT()  const {
      return m_momentum.e()*m_momentum.perp()/m_momentum.rho() ;
  }
  double  ReconstructedParticle::mT()  const {
    return m_momentum.mt() ;
  }
  IKinematic*  ReconstructedParticle::clone()  const {
    IKinematic* ik =  new ReconstructedParticle(*this) ;
    return ik;
  }
  void ReconstructedParticle::accept(IAOOvisitor* iaPtr) const{
    iaPtr->process(this);
  }
  int ReconstructedParticle::pdg_id() const {
    return m_pdg_id ; 
  }
  const HepMC::GenParticle* ReconstructedParticle::truth() const {
    return m_truth ; 
  }
  void ReconstructedParticle::print(const std::string& coor, 
                                    std::string t)  const{
    AODNavigationImp::print(coor, t);
    std::cout<<t<<std::string(73,' ')<<"ReconstructedParticle"<<std::endl;
  }
  
} // end of namespace bracket

---