---

Cell.cxx

Go to the documentation of this file.

// Cell implementation
//

#include "AtlfastEvent/Cell.h"
#include "HepMC/GenParticle.h"
#include "AtlfastEvent/EPileupDeposit.h"
static ContainedObjectFactory < Atlfast::Cell > s_AtlfastCellFactory;
const IFactory& AtlfastCellFactory = s_AtlfastCellFactory;


namespace Atlfast {

  Cell::Cell(): 
    ContainedObject(), 
    AssociationManager(),
    m_id(),
    m_ptSum(0.),
    m_momentum(0.),
    m_particles(){
  }
  Cell::Cell(const CellDescriptor& cellID) :
    ContainedObject(),
    AssociationManager(),
    m_id(cellID),
    m_ptSum(0.),
    m_momentum(0.),
    m_particles(){
  }
  
  Cell& Cell::operator=(const Cell& otherCell){

    if (this!=&otherCell){
      this->m_id=otherCell.m_id;
      this->m_ptSum=otherCell.m_ptSum;
      this->m_particles=otherCell.m_particles;//m_particles is a 
      //ContainedObject therefore do not deep copy
    }
    
    return *this;
  }
  
  Cell::Cell(const Cell& otherCell) : 
    ContainedObject(), 
    IKinematic(),
    AssociationManager(){
    *this=otherCell; //no deep copy of vector of pointers as is a 
    //ContainedObject
  }
  
  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;
  }
  
  
  void Cell::newHit(const HepMC::GenParticle* part){
    m_ptSum += part->momentum().perp();
    m_particles.push_back(part);
  }
  void Cell::newHit(const EPileupDeposit* part){
    m_ptSum += part->energyDeposit();
  }
  
  void Cell::reset(){
    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 {
    std::vector<const GenParticle*> particles;
    //std::copy(m_particles.begin(),m_particles.end(),particles.begin());
    //for(DataVector<const GenParticle>::const_iterator itr = m_particles.begin();
    //  itr != m_particles.end();++itr) particles.push_back(*itr);
    return m_particles;
    }

} // end of namespace bracket

---