RecoMuonTrackPlots.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonTrackMonitoring/RecoMuonTrackPlots.h"
 
RecoMuonTrackPlots::RecoMuonTrackPlots(PlotBase* pParent, const std::string& sDir):PlotBase(pParent, sDir),
m_oAllPlots(this, "/", "Reco Muon"),
m_oImpactPlots(this, "/"),
m_oTrkRecoInfoPlots(this, "/"),
m_pt_broad(nullptr),
m_eta_phi_broad(nullptr)
{}
 
void RecoMuonTrackPlots::initializePlots(){
  //booking histograms
  m_pt_broad = Book1D("_pt_broad", "High p_{T} Distribution", 70, 100, 1500);
  m_eta_phi_broad = Book2D("_eta_phi_broad", "High p_{T} Muon #eta #phi Distribution;#eta;#phi", 52, -2.6, 2.6, 64, -3.2, 3.2);
}
 
//when the plot function called with a Muon Container
//loop through each muon, get the corresponding link and fill it
 
//when the plot function called with a Muon
//get's the corresponding link and fill it
void RecoMuonTrackPlots::fill(const xAOD::Muon& mu, int component){
  if (component > 1) {
      return;
  }
  xAOD::Muon::TrackParticleType tType = component == 0 ? xAOD::Muon::TrackParticleType::MuonSpectrometerTrackParticle
                                                       : xAOD::Muon::TrackParticleType::ExtrapolatedMuonSpectrometerTrackParticle;
  const xAOD::TrackParticle* trk = mu.trackParticle(tType);
  if (trk) {
      fill(*trk);
  } 
}
 
 
void RecoMuonTrackPlots::fill(const xAOD::TrackParticle& muTP){
  m_oAllPlots.fill(muTP);
  m_oImpactPlots.fill(muTP);
  m_oTrkRecoInfoPlots.fill(muTP);
  
  if (muTP.pt()/1000.0 > 100) {//ony for high pt muons
    m_pt_broad->Fill(muTP.pt()/1000.0);
    m_eta_phi_broad->Fill(muTP.eta(), muTP.phi());
  }
  
}