InDetGlobalPrimaryVertexMonTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef InDetGlobalPrimaryVertexMonTool_H
#define InDetGlobalPrimaryVertexMonTool_H
 
#include <string>
#include "GaudiKernel/StatusCode.h"
#include "StoreGate/ReadHandleKey.h"
#include "AthenaMonitoring/ManagedMonitorToolBase.h"
#include "xAODTracking/VertexContainer.h"
#include "xAODTracking/Vertex.h"
 
class TH1F_LW;
class TH2F_LW;
class TProfile_LW;
 
namespace Trk {
  class VxCandidate;
}
 
class InDetGlobalPrimaryVertexMonTool : public ManagedMonitorToolBase {
 
public:
    InDetGlobalPrimaryVertexMonTool( const std::string & type, const std::string & name, const IInterface* parent );
 
    virtual ~InDetGlobalPrimaryVertexMonTool();
 
        virtual StatusCode initialize();
        virtual StatusCode bookHistogramsRecurrent();
    virtual StatusCode fillHistograms();
 
protected:
        // basic monitoring
        TH1F_LW* m_hPvN;
        TH1F_LW* m_hPvNPriVtx;
        TH1F_LW* m_hPvNPileupVtx;
 
        TProfile_LW* m_hPvN_LB;
    
        TH2F_LW* m_hPvNaveMu;
    
        TH1F_LW* m_hPvX;
        TH1F_LW* m_hPvY;
        TH1F_LW* m_hPvZ;
        TH1F_LW* m_hPvErrX;
        TH1F_LW* m_hPvErrY;
        TH1F_LW* m_hPvErrZ;
        TH1F_LW* m_hPvChiSqDoF;
 
        TH1F_LW* m_hPvNTracks;
        TH1F_LW* m_hPvTrackPt;
        TH1F_LW* m_hPvTrackEta;
 
        // enhnaced monitoring
        TH2F_LW* m_hVrt_XpullVsNtrkAverage_split;
        TH2F_LW* m_hVrt_YpullVsNtrkAverage_split;
        TH2F_LW* m_hVrt_ZpullVsNtrkAverage_split;
 
        TH2F_LW* m_hVrt_XpullVsPt2Average_split;
        TH2F_LW* m_hVrt_YpullVsPt2Average_split;
        TH2F_LW* m_hVrt_ZpullVsPt2Average_split;
 
        TH2F_LW* m_hVrt_Xerr_vs_ntrk;
        TH2F_LW* m_hVrt_Yerr_vs_ntrk;
        TH2F_LW* m_hVrt_Zerr_vs_ntrk;
 
        TH2F_LW* m_hVrt_Xerr_vs_pt2;
        TH2F_LW* m_hVrt_Yerr_vs_pt2;
        TH2F_LW* m_hVrt_Zerr_vs_pt2;
 
    // reconstruction efficiency plots
    TH1F_LW* m_hVrt_split_tag_ntrk; 
    TH1F_LW* m_hVrt_split_probe_ntrk; 
    TH1F_LW* m_hVrt_split_matched_tag_ntrk; 
    TH1F_LW* m_hVrt_split_matched_probe_ntrk; 
    TH1F_LW* m_hVrt_split_dist_tag; 
    TH1F_LW* m_hVrt_split_dist_probe; 
 
private:
        TH1F_LW* makeAndRegisterTH1F(MonGroup& mon, const char* name, const std::string& title, int nBins, float minX, float maxX);
        TH2F_LW* makeAndRegisterTH2F(MonGroup& mon, const char* hName, const std::string& hTitle,
                                     int nBinsX, float minX, float maxX,
                                     int nBinsY, float minY, float maxY);
        TH2F_LW* makeAndRegisterTH2FVariableXBins(MonGroup& mon, const char* hName, const std::string& hTitle,
                                     int nBinsX, double* xRange,
                                     int nBinsY, float minY, float maxY);
        // athena monitoring needs dummy registration of histos nad graphs which are created during post processing by a macro
        void makeAndRegisterDummyTGraph(MonGroup& mon, const char* hName);
 
    // Properties
        std::string m_histFolder;
        SG::ReadHandleKey<xAOD::VertexContainer> m_vxContainerName{this,"vxContainerName","VxPrimaryCandidate","Primary Vertices for Global Monitoring"};
        SG::ReadHandleKey<xAOD::VertexContainer> m_vxContainerNameWithoutBeamConstraint{this,"vxContainerNameWithOutBeamConstraint","VxCandidatesWithoutBeamConstraint","Vertices without beam constraint for Global Monitoring"};
        SG::ReadHandleKey<xAOD::VertexContainer> m_vxContainerNameSplit{this,"vxContainerNameSplit","SplitVxCandidates","Split Vertices for Global Monitoring"};
    int         m_splitVertexTrkInvFraction; 
    float       m_distanceSplitVxMatch; 
 
    int         m_splitMatchingMetric; 
        bool m_doEnhancedMonitoring; // trigger the enhanced monitoring (to be activated by InDetFlags.doMonitoringPrimaryVertexingEnhanced because split vertexing and vertexing with no beam constraint need to be run as well)
    // Helper functions
    double GetSplitMatchDistance(const xAOD::Vertex* splitVx, const xAOD::Vertex* originalVx);
    
};
 
#endif