ChargeDepParamPlots.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonHistUtils/ChargeDepParamPlots.h"
 
#include <utility>
 
namespace Muon {
 
ChargeDepParamPlots::ChargeDepParamPlots(PlotBase* pParent,
                                         const std::string& sDir,
                                         std::string sType)
    : PlotBase(pParent, sDir), m_sType(std::move(sType)) {}
 
void ChargeDepParamPlots::initializePlots() {
    q_times_eta =
        Book1D("qTimesEta_" + m_sType,
               "qTimesEta_" + m_sType + ";charge*Pseudo-Rapidity;Entries", 128,
               -3.2, 3.2);
    q_over_pt =
        Book1D("qOverPt_" + m_sType,
               "qOverPt_" + m_sType +
                   ";charge over Transverse Momentum [GeV^{-1}];Entries",
               200, -0.5, 0.5);
    q_over_p =
        Book1D("qOverP_" + m_sType,
               "qOverP_" + m_sType + ";charge over Momentum [GeV^{-1}];Entries",
               200, -0.5, 0.5);
 
    // charge comparison of ID, MS, ME tracks
    qMS_times_qCB =
        Book1D("qMS_times_qCB_" + m_sType,
               "qMS_times_qCB_" + m_sType + ";qMS*qCB;Entries", 3, -1.5, 1.5);
    qID_times_qCB =
        Book1D("qID_times_qCB_" + m_sType,
               "qID_times_qCB_" + m_sType + ";qID*qCB;Entries", 3, -1.5, 1.5);
}
 
void ChargeDepParamPlots::fill(const xAOD::Muon& mu, float weight) {
    const xAOD::TrackParticle* trkprt = mu.trackParticle(xAOD::Muon::Primary);
    if (!trkprt) {
        return;
    }
 
    q_times_eta->Fill(trkprt->charge() * trkprt->eta(), weight);
    q_over_pt->Fill(trkprt->charge() / trkprt->pt() / 0.001, weight);
    q_over_p->Fill(trkprt->qOverP() / 0.001, weight);
 
    const xAOD::TrackParticle* trkprtID =
        mu.trackParticle(xAOD::Muon::InnerDetectorTrackParticle);
    const xAOD::TrackParticle* trkprtMS =
        mu.trackParticle(xAOD::Muon::MuonSpectrometerTrackParticle);
 
    if (trkprtID) {
        qID_times_qCB->Fill(trkprtID->charge() * trkprt->charge(), weight);
    }
 
    if (trkprtMS) {
        qMS_times_qCB->Fill(trkprtMS->charge() * trkprt->charge(), weight);
    }
}
 
}  // namespace Muon