gFEXaltMetAlgo.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
//***************************************************************************
//    gFEXaltMetAlgo - Noise cut and Rho+RMS algorithm for gFEX MET
//                              -------------------
//     begin                : 31 03 2022
//     email                : cecilia.tosciri@cern.ch
//***************************************************************************
 
#ifndef gFEXaltMetAlgo_H
#define gFEXaltMetAlgo_H
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "L1CaloFEXToolInterfaces/IgFEXaltMetAlgo.h"
#include "AthenaKernel/CLASS_DEF.h"
#include "L1CaloFEXSim/gTowerContainer.h"
#include "L1CaloFEXSim/FEXAlgoSpaceDefs.h"
 
 
 
namespace LVL1 {
 
  class gFEXaltMetAlgo : public extends<AthAlgTool, IgFEXaltMetAlgo> {
 
  public:
    gFEXaltMetAlgo(const std::string& type, const std::string& name, const IInterface* parent);

    virtual StatusCode initialize() override;
 
                                 
    virtual void setAlgoConstant(std::vector<int>&& A_thr,
                             std::vector<int>&& B_thr,
                             std::vector<int>&& C_thr,
                             const int rhoPlusThr) override;
 
    virtual void altMetAlgo(const gTowersCentral &Atwr,
                        const gTowersCentral &Btwr,
                        const gTowersCentral &Ctwr,
                        std::array<uint32_t, 4> & outTOB) const override;
 
  private:
 
    std::array<std::vector<int>, 3> m_etaThr;   // A, B, C
    int m_rhoPlusThr{};
    
    
    void metFPGA(const gTowersCentral &twrs, int & MET_x, int & MET_y, const unsigned short FPGA_NO) const;
 
    void metTotal(const int A_MET_x, const int A_MET_y,
              const int B_MET_x, const int B_MET_y,
              const int C_MET_x, const int C_MET_y,
              int & MET_x, int & MET_y, int & MET) const;
 
    int get_rho(const gTowersCentral &twrs) const;
 
    int get_sigma(const gTowersCentral &twrs) const;
 
    void rho_MET(const gTowersCentral &twrs, int & MET_x, int & MET_y, const int rho, const int sigma) const;
 
    int sumEtFPGAnc(const gTowersCentral &twrs, const unsigned short FPGA_NO) const;
 
    int sumEtFPGArms(const gTowersCentral &twrs, const int sigma) const;
 
    int sumEt(const int A_sumEt, const int B_sumEt, const int C_sumEt) const;
 
    float sinLUT(const unsigned int phiIDX, const unsigned int aw) const;
 
    float cosLUT(const unsigned int phiIDX, const unsigned int aw) const;
 
  };
 
} // end of namespace
 
 
#endif