EFexEMEnergyWeightedClusterTool.h

Go to the documentation of this file.

// Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration

 
#ifndef LVL1__EFEXEMENERGYWEIGHTEDCLUSTERTOOL
#define LVL1__EFEXEMENERGYWEIGHTEDCLUSTERTOOL
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "CaloEvent/CaloConstCellContainer.h"
#include "xAODTrigL1Calo/TriggerTowerContainer.h"
 
#include <string>
#include <vector>
 
 
namespace LVL1
{
   class EFexEMEnergyWeightedClusterTool : public AthAlgTool
   {
   public:
      EFexEMEnergyWeightedClusterTool(const std::string& type, const std::string& name, const IInterface* parent);
 
      // find the cells at EMB2 and EMEC2 above a given threshold in an event to be used as seeds to build electron/tau candidates
      void findCellsAbove_EMB2_EMEC2(const CaloConstCellContainer*, const float& thr, std::vector<const CaloCell*>& out) const;
 
      // finds cells around a seed cell. These cells will be part of the cluster
      void findCellsAround(const CaloConstCellContainer*, const CaloCell* cell, std::vector<const CaloCell*>& out, 
                           const float deta, const float dphi) const;
 
      void findCellsAround(const CaloConstCellContainer*, const float eta, const float phi, std::vector<const CaloCell*>& out,
                           const float deta, const float dphi) const ;
 

      void findTTsAround(const xAOD::TriggerTowerContainer*, const float eta, const float phi, std::vector<const xAOD::TriggerTower*>& out) const ;

      bool isCellEmMaximum(const std::vector<const CaloCell*>& scells, const CaloCell* cell) const ;
      float sumEmCells(const std::vector<const CaloCell*>& scells) const ;
      float sumEmCells2nd(const std::vector<const CaloCell*>& scells) const ;
      float sumHadCells(const std::vector<const CaloCell*>& scells) const ;
      float sumHadTTs(const std::vector<const xAOD::TriggerTower*>& scells) const ;

      void findCluster(const std::vector<const CaloCell*>& scells, float &etaCluster, float &phiCluster) const;
 
 
 
   private:
 
      // find the cells above a given threshold in an event to be used as seeds to build electron candidates
      void findCellsAbove(const CaloConstCellContainer*, const float& thr, std::vector<const CaloCell*>& out) const;
 

   const float m_detaTT{0.125};  
   const float m_dphiTT{0.15};  
   };
}
 
#endif