FakeRatePlotTool.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file FakeRatePlotTool.hpp

// This file is part of the Acts project.
//
// Copyright (C) 2019 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.

#pragma once

#include "Acts/EventData/TrackParameters.hpp"
#include "Acts/Utilities/Logger.hpp"
#include "ActsExamples/Utilities/Helpers.hpp"
#include "ActsFatras/EventData/Particle.hpp"

#include <map>
#include <memory>
#include <string>

namespace ActsExamples {

// Tools to make fake rate plots to show tracking fake rate.
//
// The fake rate is investigated for all reco tracks. A track is 'fake' if it's
// not matched with truth.
class FakeRatePlotTool {
 public:
  struct Config {
    std::map<std::string, PlotHelpers::Binning> varBinning = {
        {"Eta", PlotHelpers::Binning("#eta", 40, -4, 4)},
        {"Phi", PlotHelpers::Binning("#phi", 100, -3.15, 3.15)},
        {"Pt", PlotHelpers::Binning("pT [GeV/c]", 40, 0, 100)},
        {"Num", PlotHelpers::Binning("N", 30, -0.5, 29.5)}};
  };

  struct FakeRatePlotCache {
    TH2F* nReco_vs_pT;          
    TH2F* nTruthMatched_vs_pT;  

    TH2F* nFake_vs_pT;   

    TH2F* nReco_vs_eta;  
    TH2F* nTruthMatched_vs_eta;  

    TH2F* nFake_vs_eta;  

    TEfficiency* fakeRate_vs_pT;   
    TEfficiency* fakeRate_vs_eta;  
    TEfficiency* fakeRate_vs_phi;  
  };

  FakeRatePlotTool(const Config& cfg, Acts::Logging::Level lvl);

  void book(FakeRatePlotCache& fakeRatePlotCache) const;

  void fill(FakeRatePlotCache& fakeRatePlotCache,
            const Acts::BoundTrackParameters& fittedParameters,
            bool status) const;

  void fill(FakeRatePlotCache& fakeRatePlotCache,
            const ActsFatras::Particle& truthParticle,
            size_t nTruthMatchedTracks, size_t nFakeTracks) const;

  void write(const FakeRatePlotCache& fakeRatePlotCache) const;

  void clear(FakeRatePlotCache& fakeRatePlotCache) const;

 private:
  Config m_cfg;                                  
  std::unique_ptr<const Acts::Logger> m_logger;  

  const Acts::Logger& logger() const { return *m_logger; }
};

}  // namespace ActsExamples