---

Cluster.h

Go to the documentation of this file.

// ================================================
// Cluster class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
// Namespace Atlfast::
//
// class: Cluster
//
// Authors: P.Clarke, H.Phillips, E.Richter-Was, P.Sherwood, R.Steward
//
// Description: 
//
//   Class to represent Clusters constructed from calorimeter Cells
//
//   It implements an "IKinematic" interface which is
//   designed to standardise kinematic attribute query method names.
//
//   Ongoing points during development:
//
//   1 It was developed during the first phase of ATLFAST++ re-engineering
//     and inclusion in the ATLAS (Gaudi/Athena) analysis framework. 
//     It is therefore one of the first implementations of such analysis classes,
//     and as such should be seen as a prototype which may be superceeded
//     in concept.
// 
//   3 At the present time we have not solved the association/relationship
//     to other entities. This is a deep problem which is wider than
//     Atlfast which is under study.
//  
//
//..................................................


#ifndef ATLFAST_CLUSTER_H
#define ATLFAST_CLUSTER_H

// Kinematic classes.
#include "CLHEP/Vector/LorentzVector.h"
#include "AtlfastEvent/IKinematic.h"

// Management of Associations
#include "AtlfastEvent/AssociationManager.h"

// Gaudi includes
#include "GaudiKernel/ContainedObject.h"
//#include "GaudiKernel/ObjectVector.h"


//Class ID for Cluster class
static const CLID CLID_ATLFAST_CLUSTER=2302 ;


//******************************************
//              Cluster class
//
// See description at head of this file 
//******************************************


namespace Atlfast {
  using ::HepLorentzVector;
  class Cluster : 
    public IKinematic, 
    public ContainedObject, 
    public AssociationManager{
      
    private:
    
    // Kinematic information assocated with this cluster  
    HepLorentzVector m_momentum ;
    
    public:
      //.............
      // Constructors
      
      template <class CellIterator>
        Cluster( const HepLorentzVector& vec, 
                 CellIterator first, 
                 CellIterator last  ) : 
        IKinematic(),
        ContainedObject(),
        AssociationManager(), 
        m_momentum(vec)
        {
          // Add associations to Cells from which the cluster was made
          // Only needs to be one direction Clusters->Cells
          for ( CellIterator i=first; i != last; ++i ) 
            {
              this->associate(*i);
            }
        }
      
      Cluster( const Cluster& src ) : 
        IKinematic(),
        ContainedObject(), 
        AssociationManager( src ),
        m_momentum(src.m_momentum)
        {}
      
      
      Cluster() : 
        IKinematic(),
        ContainedObject(), 
        AssociationManager(),
        m_momentum(HepLorentzVector(0.0,0.0,0.0,0.0))
        {}
      
      virtual ~Cluster() {}
      
      
      //..............................................
      // Methods required by IKinematic interface to kinematic quantitites
      // See interface documentation for description.
      virtual HepLorentzVector momentum() const {return m_momentum; } 
      virtual double  eta() const { return m_momentum.pseudoRapidity() ;}
      virtual double  phi() const { return m_momentum.phi() ;  } 
      virtual double  pT()  const { return m_momentum.perp() ;}
      //  FIXME: not sure if this is correct
      virtual double  eT()  const { 
        return  m_momentum.e()*m_momentum.perp()/m_momentum.rho() ;}  
      virtual double  mT()  const { return m_momentum.mt(); }   
      
      
      //......................
      // Athena mandatory methods
      static  const CLID& classID()       { return CLID_ATLFAST_CLUSTER; }
      virtual const CLID& clID()  const   {  return CLID_ATLFAST_CLUSTER; }
      
  };
  
  
} // end of namespace 
#endif

---