MuonHitDiffSummaryPlots.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonHistUtils/MuonHitDiffSummaryPlots.h"
 
#include "AthContainers/ConstAccessor.h"
#include "xAODTracking/TrackingPrimitives.h"
 
namespace Muon {
 
MuonHitDiffSummaryPlots::MuonHitDiffSummaryPlots(PlotBase* pParent,
                                                 const std::string& sDir)
    : PlotBase(pParent, sDir),
      innerSmallHits(this, "diff_innersmallhits",
                     "diff(reco - truth) inner small hits", -10, 10),
      innerLargeHits(this, "diff_innerlargehits",
                     "diff(reco - truth) inner large hits", -10, 10),
      middleSmallHits(this, "diff_middlesmallhits",
                      "diff(reco - truth) middle small hits", -10, 10),
      middleLargeHits(this, "diff_middlelargehits",
                      "diff(reco - truth) middle large hits", -10, 10),
      outerSmallHits(this, "diff_outersmallhits",
                     "diff(reco - truth) outer small hits", -10, 10),
      outerLargeHits(this, "diff_outerlargehits",
                     "diff(reco - truth) outer large hits", -10, 10),
      extendedSmallHits(this, "diff_extendedsmallhits",
                        "diff(reco - truth) extended small hits", -10, 10),
      extendedLargeHits(this, "diff_extendedlargehits",
                        "diff(reco - truth) extended large hits", -10, 10),
 
      phiLayer1Hits(this, "diff_phiLayer1hits",
                    "diff(reco - truth) phi Layer 1 hits", -10, 10),
      phiLayer2Hits(this, "diff_phiLayer2hits",
                    "diff(reco - truth) phi Layer 2 hits", -10, 10),
      phiLayer3Hits(this, "diff_phiLayer3hits",
                    "diff(reco - truth) phi Layer 3 hits", -10, 10),
      phiLayer4Hits(this, "diff_phiLayer4hits",
                    "diff(reco - truth) phi Layer 4 hits", -10, 10),
 
      etaLayer1Hits(this, "diff_etaLayer1hits",
                    "diff(reco - truth) eta Layer 1 hits", -10, 10),
      etaLayer2Hits(this, "diff_etaLayer2hits",
                    "diff(reco - truth) eta Layer 2 hits", -10, 10),
      etaLayer3Hits(this, "diff_etaLayer3hits",
                    "diff(reco - truth) eta Layer 3 hits", -10, 10),
      etaLayer4Hits(this, "diff_etaLayer4hits",
                    "diff(reco - truth) eta Layer 4 hits", -10, 10) {}
 
void MuonHitDiffSummaryPlots::fill(const xAOD::Muon& muon,
                                   const xAOD::TruthParticle& truthprt,
                                   float weight) {
    fillPlot(innerSmallHits, xAOD::innerSmallHits, muon, "innerSmallHits",
             truthprt, weight);
    fillPlot(innerLargeHits, xAOD::innerLargeHits, muon, "innerLargeHits",
             truthprt, weight);
    fillPlot(middleSmallHits, xAOD::middleSmallHits, muon, "middleSmallHits",
             truthprt, weight);
    fillPlot(middleLargeHits, xAOD::middleLargeHits, muon, "middleLargeHits",
             truthprt, weight);
    fillPlot(outerSmallHits, xAOD::outerSmallHits, muon, "outerSmallHits",
             truthprt, weight);
    fillPlot(outerLargeHits, xAOD::outerLargeHits, muon, "outerLargeHits",
             truthprt, weight);
    fillPlot(extendedSmallHits, xAOD::extendedSmallHits, muon,
             "extendedSmallHits", truthprt, weight);
    fillPlot(extendedLargeHits, xAOD::extendedLargeHits, muon,
             "extendedLargeHits", truthprt, weight);
    fillPlot(phiLayer1Hits, xAOD::phiLayer1Hits, muon, "phiLayer1Hits",
             truthprt, weight);
    fillPlot(phiLayer2Hits, xAOD::phiLayer2Hits, muon, "phiLayer2Hits",
             truthprt, weight);
    fillPlot(phiLayer3Hits, xAOD::phiLayer3Hits, muon, "phiLayer3Hits",
             truthprt, weight);
    fillPlot(phiLayer4Hits, xAOD::phiLayer4Hits, muon, "phiLayer4Hits",
             truthprt, weight);
    fillPlot(etaLayer1Hits, xAOD::etaLayer1Hits, muon, "etaLayer1Hits",
             truthprt, weight);
    fillPlot(etaLayer2Hits, xAOD::etaLayer2Hits, muon, "etaLayer2Hits",
             truthprt, weight);
    fillPlot(etaLayer3Hits, xAOD::etaLayer3Hits, muon, "etaLayer3Hits",
             truthprt, weight);
    fillPlot(etaLayer4Hits, xAOD::etaLayer4Hits, muon, "etaLayer4Hits",
             truthprt, weight);
}
 
void MuonHitDiffSummaryPlots::fillPlot(Trk::HitTypePlots& hitPlots,
                                       xAOD::MuonSummaryType info,
                                       const xAOD::Muon& muon,
                                       const std::string& sInfo,
                                       const xAOD::TruthParticle& truthprt,
                                       float weight) {
    uint8_t hitval = 0;
    if (!muon.summaryValue(hitval, info)) {
        return;
    }
    static const SG::ConstAccessor<uint8_t> acc(sInfo);
    if (!acc.isAvailable(truthprt)) {
        return;
    }
    hitPlots.fill(hitval - acc(truthprt), muon.eta(), muon.phi(), weight);
}
 
}  // namespace Muon