MuidCaloEnergyTool.h

Go to the documentation of this file.

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


 
#ifndef MUIDCALOENERGYTOOLS_MUIDCALOENERGYTOOL_H
#define MUIDCALOENERGYTOOLS_MUIDCALOENERGYTOOL_H
 
//<<<<<< INCLUDES                                                       >>>>>>
 
#include <atomic>
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "GaudiKernel/ToolHandle.h"
#include "MuidInterfaces/IMuidCaloEnergy.h"
#include "MuidInterfaces/IMuidCaloEnergyParam.h"
#include "MuidInterfaces/IMuidTrackIsolation.h"
 
//<<<<<< CLASS DECLARATIONS                                             >>>>>>
 
namespace Rec {
 
    class CaloMeas;
    class IMuidCaloEnergyMeas;
    class IMuidCaloEnergyParam;
    class IMuidTrackIsolation;
 
    class MuidCaloEnergyTool : public AthAlgTool, virtual public IMuidCaloEnergy {
    public:
        MuidCaloEnergyTool(const std::string& type, const std::string& name, const IInterface* parent);
        ~MuidCaloEnergyTool(void);  // destructor
 
        StatusCode initialize() override;
        StatusCode finalize() override;

        std::unique_ptr<CaloEnergy> energyLoss(const EventContext& ctx, double trackMomentum, double eta, double phi) const override;

        std::unique_ptr<Trk::TrackStateOnSurface> trackStateOnSurface(const EventContext& ctx, const Trk::TrackParameters& middleParameters,
                                                                      const Trk::TrackParameters* innerParameters,
                                                                      const Trk::TrackParameters* outerParameters) const override;
 
    private:
        // private methods
        std::unique_ptr<CaloEnergy> measurement(const EventContext& ctx, double trackMomentum, double eta, double phi,
                                                const CaloMeas& caloMeas) const;
        static double muSpecResolParam(double trackMomentum, double eta) ;
        static double paramCorrection(double trackMomentum, double eta, double MopLoss, double MopSigma) ;
        static double landau(double x, double mpv, double sigma, bool norm) ;
 
        // helpers, managers, tools
        ToolHandle<IMuidCaloEnergyMeas> m_caloMeasTool{
            this,
            "CaloMeasTool",
            "Rec::MuidCaloEnergyMeas/MuidCaloEnergyMeas",
        };
        ToolHandle<IMuidCaloEnergyParam> m_caloParamTool{
            this,
            "CaloParamTool",
            "Rec::MuidCaloEnergyParam/MuidCaloEnergyParam",
        };
        ToolHandle<IMuidTrackIsolation> m_trackIsolationTool{
            this,
            "TrackIsolationTool",
            "Rec::MuidTrackIsolation/MuidTrackIsolation",
        };
 
        // configurable options
        bool m_cosmics;
        bool m_energyLossMeasurement;
        bool m_forceIsolationFailure;
        bool m_FSRtreatment;
        bool m_MOPparametrization;
        bool m_trackIsolation;
 
        // thresholds for use of energy measurement
        double m_emEtCut;         // minimum Et in em for FSR treatment
        double m_emF1Cut;         // minimum F1 in em for FSR treatment
        double m_emMinEnergy;     // mininum energy in the EM
        double m_hecMinEnergy;    // minimum energy in the HEC
        int m_maxNTracksIso;      // max #tracks in the isolation cone
        double m_minFinalEnergy;  // minimum measured final energy
        double m_minMuonPt;       // minimum pt of the muon
 
        // counters (for finalize)
        mutable std::atomic_int m_countMean{0};         // number of tracks using mean
        mutable std::atomic_int m_countMeasurement{0};  // number of tracks using measurement
        mutable std::atomic_int m_countMop{0};          // number of tracks using mop
    };
 
}  // namespace Rec
 
#endif  // MUIDCALOENERGYTOOLS_MUIDCALOENERGYTOOL_H