ClusEcalEtSumVisitor.h

Go to the documentation of this file.

#ifndef ATLFAST_CLUSECALETSUMVISITOR_H
#define ATLFAST_CLUSECALETSUMVISITOR_H


#ifndef ATLFAST_IAOO_H
#include "AtlfastEvent/IAOO.h"
#endif

#ifndef ATLFAST_IAOOVISITOR_H
#include "AtlfastEvent/IAOOvisitor.h"
#endif

#ifndef ATLFAST_CLUSSTER_H
#include "AtlfastEvent/Cluster.h"
#endif

#ifndef STD_VECTOR_H
#include <vector>
#define STD_VECTOR_H
#endif

namespace Atlfast{
  class ReconstructedParticle;
  class Cell;
  class IKinematic;
  class TwoCptCell;
  class Cluster;
  class KtCluster;
  class Jet;
  class SimpleKinematic;
  class Track;

  class IAOO;

  class ClusEcalEtSumVisitor: public IAOOvisitor{
  public:
    ClusEcalEtSumVisitor(const Cluster* cluster, const double coreRadius=0.2):
      m_phi(cluster->phi()),
      m_eta(cluster->eta()),
      m_cellEt(0), 
      m_coneEt(0), 
      m_clusterEt(0),
      m_ktClusterEt(0),
      m_reconstructedParticleEt(0),
      m_simpleKinematicEt(0), 
      m_trackEt(0),
      m_coreRadius(coreRadius), m_coreEt(0){}

    ~ClusEcalEtSumVisitor(){}
    void reset();
    void operator()(const IAOO* ia);

    virtual void process(const Cell*);
    virtual void process(const Cluster*);
    virtual void process(const KtCluster*);
    virtual void process(const Jet*);
    virtual void process(const ReconstructedParticle*);
    virtual void process(const SimpleKinematic*);
    virtual void process(const Track*);
    virtual void process(const TwoCptCell*);

    //
    
    double getSum(const Cell&)                   const;
    double getSum(const TwoCptCell&)             const;
    double getSum(const Cluster&)                const;
    double getSum(const KtCluster&)              const;
    double getSum(const Jet&)                    const;
    double getSum(const SimpleKinematic&)        const;
    double getSum(const Track&)                  const;
    double getSum(const ReconstructedParticle&)  const;


    // m_coreRadius < m_coneRadius <= cluster cone radius
    // Sum Et for  dR < m_coneRadius
    double getCore(const TwoCptCell&)            const;
    // Sum Et for  m_coreRadius < dR < m_coneRadius
    double getIsolRing(const TwoCptCell&)       const;

  private:

    double m_phi;
    double m_eta;
    double m_cellEt;
    double m_coneEt;
    double m_clusterEt;
    double m_ktClusterEt;
    double m_jetEt;
    double m_reconstructedParticleEt;
    double m_simpleKinematicEt;
    double m_trackEt;
    double m_coreRadius;
    double m_coreEt;
  };
}//namespace
#endif 
  
 

---