ClusEcalEtSumVisitor.cxx

Go to the documentation of this file.

#include "AtlfastUtils/ClusEcalEtSumVisitor.h"

#include "AtlfastEvent/ReconstructedParticle.h" 
#include "AtlfastEvent/Cell.h"
#include "AtlfastEvent/TwoCptCell.h"
#include "AtlfastEvent/Cluster.h"
#include "AtlfastEvent/KtCluster.h"
#include "AtlfastEvent/Jet.h"
#include "AtlfastEvent/SimpleKinematic.h"
#include "AtlfastEvent/Track.h"
#include "AtlfastEvent/IKinematic.h"


namespace Atlfast{
  void ClusEcalEtSumVisitor::operator()(const IAOO* iaoo){iaoo->accept(this);}

  void ClusEcalEtSumVisitor::process(const Cell* t){
    const IKinematic* ik = t;
    m_cellEt += ik->eT();
  }
  void ClusEcalEtSumVisitor::process(const Cluster* t){
    const IKinematic* ik = t;
    m_clusterEt += ik->eT();
  }
  void ClusEcalEtSumVisitor::process(const KtCluster* t){
    const IKinematic* ik = t;
    m_ktClusterEt += ik->eT();
  }
  void ClusEcalEtSumVisitor::process(const Jet* t){
    const IKinematic* ik = t;
    m_jetEt += ik->eT();
  }
  void ClusEcalEtSumVisitor::process(const ReconstructedParticle* t){
    const IKinematic* ik = t;
    m_reconstructedParticleEt += ik->eT();
  }
  void ClusEcalEtSumVisitor::process(const SimpleKinematic* t){
    const IKinematic* ik = t;
    m_simpleKinematicEt += ik->eT();
  }
  void ClusEcalEtSumVisitor::process(const Track* t){
    const IKinematic* ik = t;
    m_trackEt += ik->eT();
  }
  void ClusEcalEtSumVisitor::process(const TwoCptCell* t){
    m_coneEt += t->ecalEt();
    const ICell* ic = t;
    double dR = std::sqrt(std::pow(ic->eta()-m_eta,2)+std::pow(ic->phi()-m_phi,2));
     if (dR < m_coreRadius) m_coreEt += t->ecalEt();
 }

  double ClusEcalEtSumVisitor::getSum( const Cell& )                   const {
    return m_cellEt;
  } 
  double ClusEcalEtSumVisitor::getSum( const TwoCptCell& )             const {
    return m_coneEt;
  }
  double ClusEcalEtSumVisitor::getSum( const Cluster& )                const {
    return m_clusterEt;
  }
  double ClusEcalEtSumVisitor::getSum( const KtCluster& )              const {
    return m_ktClusterEt;
  }
  double ClusEcalEtSumVisitor::getSum( const Jet& )                    const {
    return m_jetEt;
  }
  double ClusEcalEtSumVisitor::getSum( const SimpleKinematic& )        const {
    return m_simpleKinematicEt;
  }
  double ClusEcalEtSumVisitor::getSum( const Track& )                  const {
    return m_trackEt;
  }
  double  ClusEcalEtSumVisitor::getSum( const ReconstructedParticle& ) const {
    return m_reconstructedParticleEt;
  }


  double ClusEcalEtSumVisitor::getCore(const TwoCptCell&)              const {
    return m_coreEt; 
  }
  double ClusEcalEtSumVisitor::getIsolRing(const TwoCptCell&)          const {
    return  m_coneEt - m_coreEt;
  }


}//namespace;

---