EfficiencyPtPlots.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#include <utility>
 
#include "EfficiencyPtPlots.h"
                       
namespace Tau{
 
  EfficiencyPtPlots::EfficiencyPtPlots(PlotBase* pParent, const std::string& sDir, std::string sTauJetContainerName):
    PlotBase(pParent, sDir),
    m_eff_pt_jetRNNloose(nullptr),
    m_eff_pt_jetRNNmed(nullptr),
    m_eff_pt_jetRNNtight(nullptr),
    m_eff_pt_jetRNNlooseHighPt(nullptr),
    m_eff_pt_jetRNNmedHighPt(nullptr),
    m_eff_pt_jetRNNtightHighPt(nullptr),  
    m_eff_jetRNNloose(nullptr),
    m_eff_jetRNNmed(nullptr),
    m_eff_jetRNNtight(nullptr),
    m_eff_pt_jetGNTauloose(nullptr),
    m_eff_pt_jetGNTaumed(nullptr),
    m_eff_pt_jetGNTautight(nullptr),
    m_eff_pt_jetGNTaulooseHighPt(nullptr),
    m_eff_pt_jetGNTaumedHighPt(nullptr),
    m_eff_pt_jetGNTautightHighPt(nullptr),
    m_eff_jetGNTauloose(nullptr),
    m_eff_jetGNTaumed(nullptr),
    m_eff_jetGNTautight(nullptr),
    m_sTauJetContainerName(std::move(sTauJetContainerName))
  { 
  }
 
  EfficiencyPtPlots::~EfficiencyPtPlots()
  {
  }
 
  void EfficiencyPtPlots::initializePlots()
  {
    m_eff_pt_jetRNNloose       = BookTProfile("Eff_Pt_jetRNNloose"," Matched Tau loose RNN eff in pt; pt; eff", 20, 0., 150.0);
    m_eff_pt_jetRNNmed         = BookTProfile("Eff_Pt_jetRNNmed","Matched Tau med RNN eff in pt; pt; eff", 20, 0.0, 150.0);
    m_eff_pt_jetRNNtight       = BookTProfile("Eff_Pt_jetRNNtight","Matched Tau tight RNN eff in pt; pt; eff", 20, 0.0, 150.0);
    m_eff_pt_jetRNNlooseHighPt = BookTProfile("Eff_Pt_jetRNNlooseHightPt"," Matched Tau loose RNN eff in pt; pt; eff", 20, 0.0, 1500.0);
    m_eff_pt_jetRNNmedHighPt   = BookTProfile("Eff_Pt_jetRNNmedHightPt","Matched Tau med RNN eff in pt; pt; eff", 20, 0.0, 1500.0);
    m_eff_pt_jetRNNtightHighPt = BookTProfile("Eff_Pt_jetRNNtightHightPt","Matched Tau tight RNN eff in pt; pt; eff", 20, 0.0, 1500.0);
    m_eff_jetRNNloose          = BookTProfile("Eff_jetRNNloose"," Matched Tau loose RNN eff total; bin; eff",3,-1.5,1.5);
    m_eff_jetRNNmed            = BookTProfile("Eff_jetRNNmed","Matched Tau med RNN eff total; bin; eff",3,-1.5,1.5);
    m_eff_jetRNNtight          = BookTProfile("Eff_jetRNNtight","Matched Tau tight RNN eff total; bin; eff",3,-1.5,1.5);
 
    m_eff_pt_jetGNTauloose       = BookTProfile("Eff_Pt_jetGNTauloose"," Matched Tau loose GNTau eff in pt; pt; eff", 20, 0., 150.0);
    m_eff_pt_jetGNTaumed         = BookTProfile("Eff_Pt_jetGNTaumed","Matched Tau med GNTau eff in pt; pt; eff", 20, 0.0, 150.0);
    m_eff_pt_jetGNTautight       = BookTProfile("Eff_Pt_jetGNTautight","Matched Tau tight GNTau eff in pt; pt; eff", 20, 0.0, 150.0);
    m_eff_pt_jetGNTaulooseHighPt = BookTProfile("Eff_Pt_jetGNTaulooseHightPt"," Matched Tau loose GNTau eff in pt; pt; eff", 20, 0.0, 1500.0);
    m_eff_pt_jetGNTaumedHighPt   = BookTProfile("Eff_Pt_jetGNTaumedHightPt","Matched Tau med GNTau eff in pt; pt; eff", 20, 0.0, 1500.0);
    m_eff_pt_jetGNTautightHighPt = BookTProfile("Eff_Pt_jetGNTautightHightPt","Matched Tau tight GNTau eff in pt; pt; eff", 20, 0.0, 1500.0);
    m_eff_jetGNTauloose          = BookTProfile("Eff_jetGNTauloose"," Matched Tau loose GNTau eff total; bin; eff",3,-1.5,1.5);
    m_eff_jetGNTaumed            = BookTProfile("Eff_jetGNTaumed","Matched Tau med GNTau eff total; bin; eff",3,-1.5,1.5);
    m_eff_jetGNTautight          = BookTProfile("Eff_jetGNTautight","Matched Tau tight GNTau eff total; bin; eff",3,-1.5,1.5);
 
 
  }
 
  void EfficiencyPtPlots::fill(const xAOD::TauJet& tau, float weight)
  {
    if ( tau.isTau(xAOD::TauJetParameters::JetRNNSigLoose) ) {
      m_eff_pt_jetRNNloose      ->Fill(tau.pt()/1000., 1., weight);
      m_eff_pt_jetRNNlooseHighPt->Fill(tau.pt()/1000., 1., weight);
      m_eff_jetRNNloose         ->Fill(0., 1., weight);
    }
    else {
      m_eff_pt_jetRNNloose      ->Fill(tau.pt()/1000., 0., weight);
      m_eff_pt_jetRNNlooseHighPt->Fill(tau.pt()/1000., 0., weight);
      m_eff_jetRNNloose         ->Fill(0., 0., weight);
    }
   
    if ( tau.isTau(xAOD::TauJetParameters::JetRNNSigMedium) ) {
      m_eff_pt_jetRNNmed      ->Fill(tau.pt()/1000., 1., weight);
      m_eff_pt_jetRNNmedHighPt->Fill(tau.pt()/1000., 1., weight);
      m_eff_jetRNNmed         ->Fill(0., 1., weight);
    }
    else {
      m_eff_pt_jetRNNmed      ->Fill(tau.pt()/1000., 0., weight);
      m_eff_pt_jetRNNmedHighPt->Fill(tau.pt()/1000., 0., weight);
      m_eff_jetRNNmed         ->Fill(0., 0., weight);
    }
 
    if ( tau.isTau(xAOD::TauJetParameters::JetRNNSigTight) ) {
      m_eff_pt_jetRNNtight      ->Fill(tau.pt()/1000., 1., weight);
      m_eff_pt_jetRNNtightHighPt->Fill(tau.pt()/1000., 1., weight);
      m_eff_jetRNNtight         ->Fill(0., 1., weight);
    }
    else {
      m_eff_pt_jetRNNtight      ->Fill(tau.pt()/1000., 0., weight);
      m_eff_pt_jetRNNtightHighPt->Fill(tau.pt()/1000., 0., weight);
      m_eff_jetRNNtight         ->Fill(0., 0., weight);
    }
 
 
    static const SG::ConstAccessor<char> acc_GNTauL("GNTauL_v0prune");
    double pass_loose = acc_GNTauL.withDefault(tau,false);
    m_eff_pt_jetGNTauloose      ->Fill(tau.pt()/1000., pass_loose, weight);
 
 
    static const SG::ConstAccessor<char> acc_GNTauM("GNTauM_v0prune");
    double pass_medium = acc_GNTauM.withDefault(tau,false);
    m_eff_pt_jetGNTaumed      ->Fill(tau.pt()/1000., pass_medium, weight);
 
 
    static const SG::ConstAccessor<char> acc_GNTauT("GNTauT_v0prune");
    double pass_tight = acc_GNTauT.withDefault(tau,false);
    m_eff_pt_jetGNTautight      ->Fill(tau.pt()/1000., pass_tight, weight);
 
 
  }
  
 
}