MuidCaloEnergyMeas.h

Go to the documentation of this file.

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


 
#ifndef MUIDCALOENERGYTOOLS_MUIDCALOENERGYMEAS_H
#define MUIDCALOENERGYTOOLS_MUIDCALOENERGYMEAS_H
 
//<<<<<< INCLUDES                                                       >>>>>>
 
#include <atomic>
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "CaloConditions/CaloNoise.h"
#include "CaloEvent/CaloCellContainer.h"
#include "GaudiKernel/ToolHandle.h"
#include "MuidInterfaces/IMuidCaloEnergyMeas.h"
#include "StoreGate/ReadHandleKey.h"
#include "StoreGate/ReadCondHandleKey.h"
#include "CaloDetDescr/CaloDetDescrManager.h"
 
//<<<<<< CLASS DECLARATIONS                                             >>>>>>
 
class CaloCellContainer;
class Identifier;
class LArEM_ID;
class LArHEC_ID;
class TileID;
 
namespace Rec {
 
    class caloMeas;
    class IMuidCaloEnergyParam;
 
    class MuidCaloEnergyMeas : public AthAlgTool, virtual public IMuidCaloEnergyMeas {
    public:
        MuidCaloEnergyMeas(const std::string& type, const std::string& name, const IInterface* parent);
        virtual ~MuidCaloEnergyMeas();  // destructor
 
        StatusCode initialize() override;
        StatusCode finalize() override;

        std::unique_ptr<CaloMeas> energyMeasurement(const EventContext& ctx, double etaEM, double phiEM, double etaHad,
                                                    double phiHad) const override;
 
    private:
        enum SubCaloId {
            TILE = 0,
            LARHEC = 1,
            LAREM = 2,
        };
        // private methods
        void energyInCalo(CaloMeas& caloMeas, const CaloCellContainer* cellContainer, const CaloDetDescrManager* detMgr, 
                          const CaloNoise* noiseCDO, double eta, double phi, int iSubCalo) const;
 
        void isolationEnergy(CaloMeas& caloMeas, const CaloCellContainer* cellContainer, const CaloDetDescrManager* detMgr, 
                             const CaloNoise* noiseCDO, double eta, double phi, int iSubCalo) const;
 
        double energyInTile(const CaloCellContainer* cellContainer, const CaloDetDescrManager* detMgr, const CaloNoise* noiseCDO, 
                            double eta, double phi, int, int) const;
 
        double energyInLArHEC(const CaloCellContainer* cellContainer, const CaloDetDescrManager* detMgr, const CaloNoise* noiseCDO, 
                              double eta, double phi, int, int) const;
 
        double energyInLArEM(const CaloCellContainer* cellContainer, const CaloDetDescrManager* detMgr, const CaloNoise* noiseCDO, 
                             double eta, double phi, int, int) const;
 
        int cellCounting(const CaloCellContainer* cellContainer, const CaloDetDescrManager* detMgr, const CaloNoise* noiseCDO, 
                             double eta, double phi) const;
 
        //
        int samplingID(const CaloCell* cell, int iSubCalo) const;
 
        // helpers, managers, tools
        SG::ReadCondHandleKey<CaloNoise> m_noiseCDOKey{
            this,
            "CaloNoiseKey",
            "totalNoise",
            "SG Key of CaloNoise data object",
        };
        ToolHandle<IMuidCaloEnergyParam> m_caloParamTool{
            this,
            "CaloParamTool",
            "",
        };
 
        const TileID* m_tileID;
        const LArEM_ID* m_emID;
        const LArHEC_ID* m_hecID;
        SG::ReadHandleKey<CaloCellContainer> m_cellContainerLocation{
            this,
            "CellContainerLocation",
            "AllCalo",
            "calo cell container location",
        };
        
        SG::ReadCondHandleKey<CaloDetDescrManager> m_caloMgrKey{this,"CaloDetDescrManager", "CaloDetDescrManager"};
        
        double m_measurementConeTile;
        double m_measurementConeLArHEC;
        double m_measurementConeLArEM;
        double m_isolationConeTile;
        double m_isolationConeLArHEC;
        double m_isolationConeLArEM;
 
        double m_sigmasAboveNoise;      // The minimum sigmas above the noise tool rms
        double m_sigmasAboveNoiseCore;  // The minimum sigmas above the noise tool rms
 
        mutable std::atomic_int m_totalCoreCellsEM{0};
        mutable std::atomic_int m_totalCoreCellsHEC{0};
        mutable std::atomic_int m_totalCoreCellsTile{0};
        mutable std::atomic_int m_totalSelectedEM{0};
        mutable std::atomic_int m_totalSelectedHEC{0};
        mutable std::atomic_int m_totalSelectedTile{0};
    };
 
}  // namespace Rec
 
#endif  // MUIDCALOENERGYTOOLS_MUIDCALOENERGYMEAS_H