---

ClusterConeStrategy.h

Go to the documentation of this file.

// ================================================
// ClusterConeStrategy class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
//
// Namespace Atlfast::
//
// class: ClusterConeStrategy
//
// Authors: P.Clarke, H.Phillips, E.Richter-Was, P.Sherwood, R.Steward
//
// ................................................................
//


#ifndef __Atlfast_ClusterConeStrategy__
#define __Atlfast_ClusterConeStrategy__

// STL
#include <string>
#include <vector>
#include <list>

// Gaudi 
#include "GaudiKernel/MsgStream.h"
// Other
#include "CLHEP/Vector/LorentzVector.h"
#include "HepMC/GenParticle.h"

// Atlfast
#include "AtlfastCode/ReconstructedParticle.h"
#include "AtlfastCode/ReconstructedParticleCollection.h"
#include "AtlfastCode/Cell.h"
#include "AtlfastCode/CellCollection.h"
#include "AtlfastCode/Cluster.h"
#include "AtlfastCode/ClusterCollection.h"

#include "AtlfastCode/IKinematic.h"
#include "AtlfastCode/KinematicHelper.h"
#include "AtlfastCode/IClusterStrategy.h"


//*************************************
//Default parameters for this algorithm
// (descriptions in class declaration)
//*************************************

//***************************************************************
//                 ClusterConeStrategy class declaration
//****************************************************************

namespace Atlfast {
  typedef std::vector<Cell*>                       localCellCollection ;
  typedef localCellCollection::iterator            localCellIterator ;
  class ClusterConeStrategy:  public IClusterStrategy {
  public:
   
    
    ClusterConeStrategy(
                        double rConeBarrel,
                        double rConeForward,
                        double minInitiatorET, 
                        double minClusterET
                        ):
      m_rConeBarrel(rConeBarrel),
      m_rConeForward(rConeForward),
      m_minInitiatorET(minInitiatorET), 
      m_minClusterET(minClusterET){}
    
    virtual void makeClusters(
                              MsgStream& log, 
                              std::vector<Cell*> storedCells,
                              CellCollection* unusedCells,
                              ClusterCollection* clusters) const;
 
  private:
    
    //------------------------------------
    // Types used internally by this class
    //------------------------------------
    // For local mutable copies of collections needed for algorithm operation
    std::list<Cell*> m_cells;       
    
    
    
    //---------------------------------------------------
    // Parameters of this algorithm 
    //---------------------------------------------------
    
    double        m_rConeBarrel;   
    double        m_rConeForward;   
    
    double        m_minInitiatorET ;
    
    double        m_minClusterET;     
    
    //-------------------------------
    // Objects used by this class to do its job
    //-------------------------------
    
    KinematicHelper m_kinehelp ; 
    
    /*#  Cell lnkCell; */
    //-------------------------------
    // Private methods
    //-------------------------------
    
    double rCone() const;
    //-------------------------------
    // Private helper class
    // For accumulating weighted kinematic quantities
    //-------------------------------
    
    class PreCluster {
      
    private:
      
      double m_eT_total ;
      double m_eta_weighted ;
      HepLorentzVector m_momentum_sum ;
      
    public:
      
      // Construct
      PreCluster( localCellIterator start, localCellIterator end )  ;
      
      double eT()  ;
      double eta() ;
      double phi() ;
      
      operator HepLorentzVector () ;
    };
    
    
    Cluster * lnkCluster;
  };
  
} // end of namespace bracket


#endif

---