Had3ProngPlots.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#include "Had3ProngPlots.h"
#include "AthContainers/ConstAccessor.h"
 
namespace Tau{
 
Had3ProngPlots::Had3ProngPlots(PlotBase* pParent, const std::string& sDir, const std::string& sTauJetContainerName):
   PlotBase(pParent, sDir),
   m_oGeneralTauPlots(this, "", sTauJetContainerName),
   m_tauCoreFrac(nullptr),
   m_tauEoverPTrk(nullptr),
   m_tauTrkAvgDist(nullptr),
   m_tauDRMax(nullptr),
   m_tauSflight(nullptr),
   m_tauMtrks(nullptr),
   m_tauIpSig(nullptr),
   m_innerTrkAvgDist(nullptr),
   m_ptRatioEflowApprox(nullptr),
   m_mEflowApprox(nullptr),
   m_ChPiEMEOverCaloEME(nullptr),
   m_EMPOverTrkSysP(nullptr),
   m_sTauJetContainerName(sTauJetContainerName)
{
}
    
Had3ProngPlots::~Had3ProngPlots()
{
}
 
void Had3ProngPlots::initializePlots(){
 
  //  m_oGeneralTauPlots.initialize();
  m_tauCoreFrac   = Book1D("CentFrac",m_sTauJetContainerName + " Tau CentFrac; CentFrac; # Taus",40,0.,2.);
  m_tauEoverPTrk  = Book1D("EoverPTrk",m_sTauJetContainerName + " Tau E over pLeadTrk; EoverP; # Taus",50,0.,4.);
  m_tauTrkAvgDist = Book1D("TrkAvgDist",m_sTauJetContainerName + " Tau track avg dist; TrkAvgDist; # Taus",10,-0.2,.8);
  m_tauDRMax      = Book1D("DRMax",m_sTauJetContainerName + " Tau DR Max track-seed; DRMax; # Taus",20,0.,0.5);       
  m_tauSflight    = Book1D("Sflight",m_sTauJetContainerName + " Tau flight sign. ; Sflight; # Taus",100,-10.,20.);
  m_tauMtrks      = Book1D("Mtracks",m_sTauJetContainerName + " Tau tracks mass ; Mtrks; # Taus",100,-10.,20.);
  m_tauIpSig      = Book1D("IpSig",m_sTauJetContainerName + " Tau lead track IP signif. ; IpSig; # Taus",30,-15.,15.);
  m_innerTrkAvgDist    = Book1D("innerTrkAvgDist",m_sTauJetContainerName + "Tau innerTrkAvgDist; innerTrkAvgDist; # of Taus", 10,-0.2,.8);
  m_ptRatioEflowApprox = Book1D("ptRatioEflowApprox",m_sTauJetContainerName + "Tau ptRatioEflowApprox; ptRatioEflowApprox; # of Taus", 10, 0, 10);
  m_mEflowApprox       = Book1D("mEflowApprox",m_sTauJetContainerName + "Tau mEflowApprox; mEflowApprox; #Taus",30,0.0,3000.0);
  m_ChPiEMEOverCaloEME = Book1D("ChPiEMEOverCaloEME",m_sTauJetContainerName + "Tau ChPiEMEOverCaloEME; ChPiEMEOverCaloEME; # of Taus", 30, -15, 15);
  m_EMPOverTrkSysP     = Book1D("EMPOverTrkSysP",m_sTauJetContainerName + "Tau EMPOverTrkSysP; EMPOverTrkSysP; # of Taus", 81, -1, 80);
 
}
 
void Had3ProngPlots::fill(const xAOD::TauJet& tau, float weight) {
  m_oGeneralTauPlots.fill(tau, weight);
 
  float avariable = 0.;
 
  bool test = tau.detail(xAOD::TauJetParameters::centFrac, avariable);
  if(test) m_tauCoreFrac->Fill(avariable, weight);
 
  test = tau.detail(xAOD::TauJetParameters::etOverPtLeadTrk, avariable);
  if(test) m_tauEoverPTrk->Fill(avariable, weight);
 
  test = tau.detail(xAOD::TauJetParameters::trkAvgDist, avariable);
  if(test) m_tauTrkAvgDist->Fill(avariable, weight);
 
  test = tau.detail(xAOD::TauJetParameters::dRmax, avariable);
  if(test) m_tauDRMax->Fill(avariable, weight);
 
  test = tau.detail(xAOD::TauJetParameters::trFlightPathSig, avariable);
  if(test) m_tauSflight->Fill(avariable, weight);
 
  test = tau.detail(xAOD::TauJetParameters::massTrkSys, avariable);
  if(test) m_tauMtrks->Fill(avariable/1000., weight);
 
  static const SG::ConstAccessor<float> d0SigTJVAAcc("d0SigTJVA");
  if(tau.nTracks()>0 && d0SigTJVAAcc.isAvailable(*tau.track(0))) {
      m_tauIpSig->Fill(tau.track(0)->d0SigTJVA(), weight);
  }
 
  test = tau.detail(xAOD::TauJetParameters::innerTrkAvgDist, avariable);
  if(test) m_innerTrkAvgDist->Fill(avariable, weight);
 
  test = tau.detail(xAOD::TauJetParameters::ptRatioEflowApprox, avariable);
  if(test) m_ptRatioEflowApprox->Fill(avariable, weight);
 
  test = tau.detail(xAOD::TauJetParameters::mEflowApprox, avariable);
  if(test) m_mEflowApprox->Fill(avariable, weight);
 
  test = tau.detail(xAOD::TauJetParameters::ChPiEMEOverCaloEME, avariable);
  if(test) m_ChPiEMEOverCaloEME->Fill(avariable, weight);
 
  test = tau.detail(xAOD::TauJetParameters::EMPOverTrkSysP, avariable);
  if(test) m_EMPOverTrkSysP->Fill(avariable, weight);
}
 
}