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 == 0 ) {
    const ElementLink<xAOD::TrackParticleContainer>& Mu_MStrack = mu.muonSpectrometerTrackParticleLink();
    if(Mu_MStrack.isValid()){
      const xAOD::TrackParticle* trk = *Mu_MStrack;
      fill(*trk);
    }
  }
 
  if (component == 1 ){
 
    const ElementLink<xAOD::TrackParticleContainer>& Mu_metrack = mu.trackParticleLink(xAOD::Muon::TrackParticleType::ExtrapolatedMuonSpectrometerTrackParticle);
    if(Mu_metrack.isValid()){
      const xAOD::TrackParticle* trk = *Mu_metrack;
      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());
  }
  
}