Cell.cxx

Go to the documentation of this file.

// Cell implementation
//

#include "AtlfastEvent/Cell.h"
#include "AtlfastEvent/IAOOvisitor.h"
#include <iomanip>

namespace Atlfast {
  Cell::Cell(): 
    AODNavigationImp(),
    m_id(),
    m_ptSum(0.),
    m_momentum(0.),
    m_particles(){
  }
  Cell::Cell(const CellDescriptor& cellID) :
    AODNavigationImp(),
    m_id(cellID),
    m_ptSum(0.),
    m_momentum(0.),
    m_particles(){

  }
  
  Cell& Cell::operator=(const Cell& otherCell){

    if (this!=&otherCell){
      AODNavigationImp::operator=(otherCell);
      this->m_id=otherCell.m_id;
      this->m_ptSum=otherCell.m_ptSum;
      this->m_particles=otherCell.m_particles;//shallow copy
    }
    
    return *this;
  }
  
  Cell::Cell(const Cell& otherCell): 
    Atlfast::IAOO(otherCell),
    INavigable(otherCell),
    I4Momentum(otherCell),
    INavigable4Momentum(otherCell),
    Atlfast::IKinematic(otherCell),
    Atlfast::IAODNavigation(otherCell),
    Atlfast::ICell(otherCell),
    AODNavigationImp(otherCell){
    *this=otherCell; //shallow copy
  }
  
  HepLorentzVector Cell::momentum() const{
    //sum momentum calculated as if all particles hit cell center 
    double px,py,pz,e;
    px=pT()*cos(m_id.phi());
    py=pT()*sin(m_id.phi());
    pz=pT()*sinh(m_id.eta());
    e=pT()*cosh(m_id.eta());
    HepLorentzVector vec(px,py,pz,e);
    return vec;
  }
  //IAOO interface
  void Cell::accept(IAOOvisitor* iaPtr) const{
    iaPtr->process(this);
  }
  void Cell::newHit(const HepMC::GenParticle* part){
    m_ptSum += part->momentum().perp();
    this->addParticle(part);
  }    
    void Cell::newHit(const HepLorentzVector momentum){
    m_ptSum += momentum.perp();
  }    
  void Cell::newHit(double pT){
    //    cerr<<"Cell::newHit(double pT) "<<pT<<endl;
    m_ptSum += pT;
  }    

  void Cell::newHit(const EPileupDeposit* part){
    m_ptSum += part->energyDeposit();
  }
  void Cell::resetCell(){
    m_ptSum=0.0;
    m_particles.erase(m_particles.begin(), m_particles.end());
  }
  void Cell::setPt(HepLorentzVector& vec){
    m_ptSum=vec.perp(); 
  } 

  std::vector<const GenParticle*> Cell::particles() const {
    return m_particles;
  }
  double Cell::eta() const {return m_id.eta();}
  double Cell::phi() const {return m_id.phi();} 
  double Cell::pT()  const {return m_ptSum;}
  double Cell::eT() const {return m_ptSum;} 
  double Cell::mT() const {return 0;}
  IKinematic* Cell::clone() const{
    IKinematic* ik = new Cell(*this);
    return ik;
  }
  void Cell::addParticle(const HepMC::GenParticle* part){
    if(find(m_particles.begin(), m_particles.end(), part)==m_particles.end()){
      m_particles.push_back(part);
    }    
  }
  void Cell::print(const std::string& coor, std::string t) const{
    AODNavigationImp::print(coor, t);
    std::cout<<t<<"Number of particles: "
             <<std::setw(6)<<m_particles.size()
             <<std::string(45,' ')
             <<" Cell"<<std::endl;
  }
} // end of namespace bracket

---