PhysValMET.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
// PhysValMET.h 
// Header file for class PhysValMET
// Author: Daniel Buescher <daniel.buescher@cern.ch>, Philipp Mogg <philipp.mogg@cern.ch>
#ifndef MISSINGETDQA_PHYSVALMET_H
#define MISSINGETDQA_PHYSVALMET_H 1
 
// STL includes
#include <string>
 
// FrameWork includes
#include "GaudiKernel/ServiceHandle.h"
 
// Local includes
#include "AthenaMonitoring/ManagedMonitorToolBase.h"
 
#include "MuonAnalysisInterfaces/IMuonSelectionTool.h"
#include "EgammaAnalysisInterfaces/IAsgElectronLikelihoodTool.h"
#include "EgammaAnalysisInterfaces/IAsgPhotonIsEMSelector.h"
#include "JetInterface/IJetUpdateJvt.h"
#include "METInterface/IMETMaker.h"
#include "TauAnalysisTools/ITauSelectionTool.h"
 
#include "AsgTools/ToolHandle.h"
 
// Root includes
#include "TH1.h"
 
// Forward declaration
class IMETMaker;
class IAsgElectronLikelihoodTool;
class IAsgPhotonIsEMSelector;
namespace TauAnalysisTools {
  class ITauSelectionTool;
}
 
 
namespace MissingEtDQA {
 
class PhysValMET
  : public ManagedMonitorToolBase
{ 
  // Public methods: 
 public: 
 
  // Copy constructor: 

  PhysValMET( const std::string& type,
          const std::string& name, 
          const IInterface* parent);

  virtual ~PhysValMET(); 
 
  // Athena algtool's Hooks
  virtual StatusCode initialize();
  virtual StatusCode bookHistograms();
  virtual StatusCode fillHistograms();
  virtual StatusCode procHistograms();
 

  // Const methods: 

  // Non-const methods: 
 

  // Private data: 
 private:
  bool m_doTruth;
  bool m_inputIsDAOD;
  bool m_doMETRefPlots;

  PhysValMET();
 
  // Containers
  std::string m_eleColl;
  std::string m_gammaColl;
  std::string m_tauColl;
  std::string m_muonColl;
  std::string m_mapname;
  std::string m_corename;
 
  // Methods
  bool Accept(const xAOD::Electron* el);
  bool Accept(const xAOD::Photon* ph);
  bool Accept(const xAOD::TauJet* tau);
  bool Accept(const xAOD::Muon* muon);
  bool Accept(const xAOD::Jet* jet, double JvtCut, ToolHandle<IJetUpdateJvt>* jvtTool);
 
  // vector of collections
  std::vector <std::string> m_types;
 
  // vector of terms
  std::vector <std::string> m_terms;
 
  // Map for names
  std::map <std::string,std::string> m_names;
 
  // Hists
  TH1D *m_MET_Track = nullptr, *m_MET_Track_x = nullptr, *m_MET_Track_y = nullptr, *m_MET_Track_phi = nullptr, *m_MET_Track_sum = nullptr;
  TH1D *m_MET_PVTrack_Nominal = nullptr, *m_MET_PVTrack_Nominal_x = nullptr, *m_MET_PVTrack_Nominal_y = nullptr, *m_MET_PVTrack_Nominal_phi = nullptr, *m_MET_PVTrack_Nominal_sum = nullptr;
  TH1D *m_MET_PVTrack_Pileup = nullptr, *m_MET_PVTrack_Pileup_x = nullptr, *m_MET_PVTrack_Pileup_y = nullptr, *m_MET_PVTrack_Pileup_phi = nullptr, *m_MET_PVTrack_Pileup_sum = nullptr;
  TH1D *m_MET_Calo = nullptr, *m_MET_Calo_x = nullptr, *m_MET_Calo_y = nullptr, *m_MET_Calo_phi = nullptr, *m_MET_Calo_sum = nullptr;
 
  //Maps
  std::map<std::string,std::vector<TH1D*> > m_MET_Ref;
  std::map<std::string,std::vector<TH1D*> > m_MET_Ref_x;
  std::map<std::string,std::vector<TH1D*> > m_MET_Ref_y;
  std::map<std::string,std::vector<TH1D*> > m_MET_Ref_phi;
  std::map<std::string,std::vector<TH1D*> > m_MET_Ref_sum;
  std::map<std::string,std::vector<TH1D*> > m_MET_Diff_Ref;
  std::map<std::string,std::vector<TH1D*> > m_MET_Diff_Ref_x;
  std::map<std::string,std::vector<TH1D*> > m_MET_Diff_Ref_y;
  std::map<std::string,std::vector<TH1D*> > m_MET_Diff_Ref_phi;
  std::map<std::string,std::vector<TH1D*> > m_MET_Diff_Ref_sum;
  std::map<std::string,std::vector<TH1D*> > m_MET_Cumu_Ref;
  std::map<std::string,std::vector<TH1D*> > m_MET_Resolution_Ref;
  std::map<std::string,std::vector<TH1D*> > m_MET_Significance_Ref;
  std::map<std::string,std::vector<TH1D*> > m_MET_dPhi_Ref;
  std::map<std::string,std::vector<TH2D*> > m_MET_CorrFinalTrk_Ref;
  std::map<std::string,std::vector<TH2D*> > m_MET_CorrFinalClus_Ref;
  std::map<std::string,std::vector<TH1D*> > m_MET_Reb;
  std::map<std::string,std::vector<TH1D*> > m_MET_Reb_x;
  std::map<std::string,std::vector<TH1D*> > m_MET_Reb_y;
  std::map<std::string,std::vector<TH1D*> > m_MET_Reb_phi;
  std::map<std::string,std::vector<TH1D*> > m_MET_Reb_sum;
  std::map<std::string,std::vector<TH1D*> > m_MET_Diff_Reb;
  std::map<std::string,std::vector<TH1D*> > m_MET_Diff_Reb_x;
  std::map<std::string,std::vector<TH1D*> > m_MET_Diff_Reb_y;
  std::map<std::string,std::vector<TH1D*> > m_MET_Diff_Reb_phi;
  std::map<std::string,std::vector<TH1D*> > m_MET_Diff_Reb_sum;
  std::map<std::string,std::vector<TH1D*> > m_MET_Cumu_Reb;
  std::map<std::string,std::vector<TH1D*> > m_MET_Resolution_Reb;
  std::map<std::string,std::vector<TH1D*> > m_MET_Significance_Reb;
  std::map<std::string,std::vector<TH1D*> > m_MET_dPhi_Reb;
  std::map<std::string,std::vector<TH2D*> > m_MET_CorrFinalTrk_Reb;
  std::map<std::string,std::vector<TH2D*> > m_MET_CorrFinalClus_Reb;
 
  std::vector<std::string> m_dir_met;
 
  ToolHandle<CP::IMuonSelectionTool> m_muonSelTool            {this, "MuonSelectionTool",        "", "Muon selection tool"};
  ToolHandle<IAsgElectronLikelihoodTool> m_elecSelLHTool      {this, "ElectronLHSelectionTool",  "", "Electron likelihood selection tool"};
  ToolHandle<IAsgPhotonIsEMSelector> m_photonSelIsEMTool      {this, "PhotonIsEMSelectionTool" , "", "Photon selection tool"};
  ToolHandle<IJetUpdateJvt> m_jvtToolEM                       {this, "JVTToolEMTopo",            "", "JVT tool for EMTopo jets"};
  ToolHandle<IJetUpdateJvt> m_jvtToolPFlow                    {this, "JVTToolEMPFlow",           "", "JVT tool forEMPFlow jets"};
  ToolHandle<IMETMaker> m_metmakerTopo                        {this, "METMakerTopo",             "", "METMaker for EMTopo jets"};
  ToolHandle<IMETMaker> m_metmakerPFlow                       {this, "METMakerPFlow",            "", "METMaker for EMPFlow jets"};
  ToolHandle<TauAnalysisTools::ITauSelectionTool> m_tauSelTool{this, "TauSelectionTool",         "", "Tau selection tool"};
  ToolHandle<IMETMaker>* m_metmaker;
}; 
 
}
 
#endif //> !MISSINGETDQA_PHYSVALMET_H