jTower.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
//***************************************************************************
//                           jTower.h  -  description
//                              -------------------
//     begin                : 19 02 2019
//     email                : Alan.Watson@cern.ch, jacob.julian.kempster@cern.ch
//  ***************************************************************************/
 
 
#ifndef jTower_H
#define jTower_H
 
#include <vector>
#include "AthenaKernel/CLASS_DEF.h"
#include "CaloEvent/CaloCellContainer.h"
#include "CaloIdentifier/CaloIdManager.h"
#include "CaloIdentifier/CaloCell_SuperCell_ID.h"
 
#include "xAODBase/IParticle.h"
#include "xAODCore/AuxStoreAccessorMacros.h"
 
namespace LVL1 {
  
  //Doxygen class description below:
  class jTower {
    
  public:
    
    jTower();
    jTower(int ieta, int iphi, int towerid, int posneg, float centre_eta = -1.0, float centre_phi = -1.0, int fcal_layer = -1);
    
    virtual ~jTower() = default;
 
    void clearET();
    
    void clear_EM_scIDs();
    void clear_HAD_scIDs();
    
    void setEMSat(){m_EM_sat = true;};
    bool getEMSat() const {return m_EM_sat; };
    
    void setHADSat(){m_HAD_sat = true;};
    bool getHADSat() const {return m_HAD_sat;};
    
    bool getTowerSat() const { return (m_EM_sat || m_HAD_sat);};
 
    void addET(float et, int cell);
 
    void recordMD_ET(float et, int cell);
 
    void setCentreEta(float ieta);
    void setCentrePhi(float iphi);
    void setiEta(int ieta);
    void setiPhi(int iphi);
    
    int iEta() const;
    int iPhi() const; 
       
    float centreEta()      const {return m_centre_eta;}
    float centrePhi()      const {return m_centre_phi;}
    float centrephi_toPI() const {return m_centre_phi_toPI;}    
 
    void setTTowerArea(int area,int layer);
    int getTTowerArea(int layer) const;
    
    void setTTowerAreaInv(int area,int layer);
    int getTTowerAreaInv(int layer) const;
 
    void setMinEtforPileup(int etval){m_minEt_pileup_thr=etval;};
    int getMinEtforPileup() const {return m_minEt_pileup_thr;};
    void setMaxEtforPileup(int etval){m_maxEt_pileup_thr=etval;};
    int getMaxEtforPileup() const {return m_maxEt_pileup_thr;};
 
 
    
 
    int fcalLayer() const {return m_fcal_layer;}
    
    // jTower ID Online and Offline
    void setOnlineID(int tower_id_online);
    
    int OnlineID() {return m_tower_id_online;};
    int OnlineID() const {return m_tower_id_online;}
    int OfflineID() {return m_tower_id;};
    int OfflineID() const {return m_tower_id;}
 
    float constid() const {return m_tower_id;};
 
    int getET(unsigned int layer, int cell = 0) const;
    
    int getTotalET() const;
 
    int getLayerTotalET(unsigned int layer) const;
 
    std::vector<int> getLayerETvec(unsigned int layer) const;
    
    std::vector<int> getETs() const {return m_et;};
    
    void set_Et(int layer, int et);
    
    int getET_EM() const {return m_et[0];};
    
    int getET_HAD() const {return m_et[1];};
 
    float getET_float(unsigned int layer, int cell = 0) const;
 
    float getTotalET_float() const;
 
    float getLayerTotalET_float(unsigned int layer) const;
 
    std::vector<float> getLayerETvec_float(unsigned int layer) const;
 
    std::vector<float> getETs_float() const {return m_et_float_raw;};
 
    void set_TileCal_Et(int layer, int et);
 
    void set_LAr_Et(Identifier ID, int cell, float et, int layer);
    
    void Do_LAr_encoding();
    
    void  setNoiseForMet(int noiseVal,int layer);
    int getNoiseForMet(int layer)const;
    void  setNoiseForJet(int noiseVal,int layer);
    int getNoiseForJet(int layer)const;
 
    std::vector<Identifier> getEMSCIDs() const { return m_EM_scID; }
    std::vector<Identifier> getHADSCIDs() const { return m_HAD_scID; }
 
    Identifier getEMSCID(int cell) const { return m_EM_scID[cell]; }
    Identifier getHADSCID(int cell) const { return m_HAD_scID[cell]; }
 
    std::vector<Identifier> getLayerSCIDs(unsigned int layer) const;
 
    void setPosNeg(int posneg);
 
    inline int getPosNeg() const {return m_posneg;}
 
  private:
        
    int m_iEta=0;
    int m_iPhi=0;
    int m_tower_id = 0;
    int m_tower_id_online = 0;
    int m_posneg = 0;
    float m_centre_eta =0;
    float m_centre_phi =0;
    float m_centre_phi_toPI=0;
    
    std::vector<Identifier> m_EM_scID;
    std::vector<Identifier> m_HAD_scID;
    std::vector<int> m_et; // Real energy from TILE and the decoded energy from LATOME
    std::vector<float> m_et_float_raw;  // Raw imput energy from LATOME (not encoded/decoded) and TILE.
    std::vector<int> m_TTowerArea{ 1, 1};
    std::vector<int> m_TTowerAreaInv{ 1, 1};
    int m_fcal_layer = -1;
    int m_NoiseForMet[2] = {0};
    int m_NoiseForJet[2] = {0};
    int m_minEt_pileup_thr = -999;
    int m_maxEt_pileup_thr = -999;
    bool m_EM_sat = false;
    bool m_HAD_sat = false;
    
    
  };
  
} // end of namespace
 
CLASS_DEF( LVL1::jTower , 41848655 , 1 )
 
#endif