MuonAlignmentUncertTool.h

Go to the documentation of this file.

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

//   MuonAlignmentUncertTool
//      Helper classes to read the scattering and alignment uncertainties
//      of the muons from the histograms
 
#ifndef MUONCOMBINEDBASETOOLS_MUONALIGNMENTUNCERTTOOL_H
#define MUONCOMBINEDBASETOOLS_MUONALIGNMENTUNCERTTOOL_H
 
#include <AthenaBaseComps/AthAlgTool.h>
#include <GaudiKernel/ServiceHandle.h>
#include <MuidInterfaces/IMuonAlignmentUncertTool.h>
#include <TrkDetDescrInterfaces/ITrackingVolumesSvc.h>
#include <TrkGeometry/TrackingVolume.h>
#include <TrkTrack/Track.h>
#include <TrkTrack/TrackStateOnSurface.h>
 
#include "GaudiKernel/ToolHandle.h"
#include <memory>
#include <TH1.h>
 
#include <functional>
#include <map>
#include <string>
#include <vector>
 
namespace Muon {
    class MuonAlignmentUncertTool : public AthAlgTool, virtual public IMuonAlignmentUncertTool {
    public:
        MuonAlignmentUncertTool(const std::string& type, const std::string& name, const IInterface* parent);
        ~MuonAlignmentUncertTool() = default;
 
        StatusCode initialize() override;
 
        float get_uncertainty(const Trk::Track* trk) const override;

        void get_track_state_measures(const Trk::Track* track, const Trk::TrackStateOnSurface*& id_exit,
                                      const Trk::TrackStateOnSurface*& calo_entrance, const Trk::TrackStateOnSurface*& calo_exit,
                                      const Trk::TrackStateOnSurface*& ms_entrance) const override;
 
    private:
        std::function<int(const Trk::Track* trk)> make_reader(const TAxis* axis) const;

        ServiceHandle<Trk::ITrackingVolumesSvc> m_trackingVolumesSvc{this, "TrackingVolumesSvc", "Trk::TrackingVolumesSvc/TrackingVolumesSvc"};
 
        Gaudi::Property<std::string> m_in_file{this, "InFile",
                                               "MuonCombinedBaseTools/AlignmentUncertainties/201029_initial/ID_MS_Uncertainties.root",
                                               "Location of the calibration file"};
        Gaudi::Property<std::string> m_histo_name{this, "HistoName", "", "Name of the histogram in the file"};

        std::unique_ptr<TH1> m_histo;
        unsigned int m_dim;
        std::function<int(const Trk::Track* trk)> m_x_reader;
        std::function<int(const Trk::Track* trk)> m_y_reader;
        std::function<int(const Trk::Track* trk)> m_z_reader;
 
        std::unique_ptr<const Trk::Volume> m_calorimeterVolume;
        std::unique_ptr<const Trk::Volume> m_indetVolume;
    };
 
}  // namespace Muon
 
#endif