Tau1P3PDetails.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/

 
#include "tauEvent/Tau1P3PDetails.h"
#include "AthenaKernel/ClassName.h"
 
namespace Analysis
{
    Tau1P3PDetails :: Tau1P3PDetails():
        TauDetails(),
        m_numStripCells( TauDetails::DEFAULT_INT ),
        m_stripWidth2( TauDetails::DEFAULT ),
        m_emRadius( TauDetails::DEFAULT ),
        m_ET12Frac( TauDetails::DEFAULT ),
        m_etIsolHAD( TauDetails::DEFAULT ),
        m_etIsolEM( TauDetails::DEFAULT ),
        m_etChrgHAD( TauDetails::DEFAULT ),
        m_nAssocTracksCore( TauDetails::DEFAULT_INT ),
        m_nAssocTracksIsol( TauDetails::DEFAULT_INT ),
        m_signD0Trk3P( TauDetails::DEFAULT ),
        m_massTrk3P( TauDetails::DEFAULT ),
        m_rWidth2Trk3P( TauDetails::DEFAULT ),
        m_z0SinThetaSig( TauDetails::DEFAULT ),
        m_etChrgHADoverPttot( TauDetails::DEFAULT ),
        m_etIsolFrac( TauDetails::DEFAULT ),
        m_etHadAtEMScale( TauDetails::DEFAULT ),
        m_etEMAtEMScale( TauDetails::DEFAULT ),
        m_etEMCL( TauDetails::DEFAULT ),
        m_etChrgEM( TauDetails::DEFAULT ),
        m_etNeuEM( TauDetails::DEFAULT ),
        m_etResNeuEM( TauDetails::DEFAULT ),
        m_trFlightPathSig(TauDetails::DEFAULT),
        m_secVertex(nullptr),
        m_etEflow(TauDetails::DEFAULT),
        m_sumEtCellsLArOverLeadTrackPt(TauDetails::DEFAULT),
        m_hadronicLeak(TauDetails::DEFAULT),
        m_secondaryMax(TauDetails::DEFAULT)
    {
    }
 
    Tau1P3PDetails::~Tau1P3PDetails() = default;
 
    Tau1P3PDetails::Tau1P3PDetails(const Tau1P3PDetails& rhs):
      TauDetails(),
      m_numStripCells(rhs.m_numStripCells),
      m_stripWidth2(rhs.m_stripWidth2),
      m_emRadius(rhs.m_emRadius),
      m_ET12Frac(rhs.m_ET12Frac),
      m_etIsolHAD(rhs.m_etIsolHAD),
      m_etIsolEM(rhs.m_etIsolEM),
      m_etChrgHAD(rhs.m_etChrgHAD),
      m_nAssocTracksCore(rhs.m_nAssocTracksCore),
      m_nAssocTracksIsol(rhs.m_nAssocTracksIsol),
      m_signD0Trk3P(rhs.m_signD0Trk3P),
      m_massTrk3P(rhs.m_massTrk3P),
      m_rWidth2Trk3P(rhs.m_rWidth2Trk3P),
      m_z0SinThetaSig(rhs.m_z0SinThetaSig),
      m_etChrgHADoverPttot(rhs.m_etChrgHADoverPttot),
      m_etIsolFrac(rhs.m_etIsolFrac),
      m_cellEM012Cluster(rhs.m_cellEM012Cluster),
      m_etHadAtEMScale(rhs.m_etHadAtEMScale),
      m_etEMAtEMScale(rhs.m_etEMAtEMScale),
      m_etEMCL(rhs.m_etEMCL),
      m_etChrgEM(rhs.m_etChrgEM),
      m_etNeuEM(rhs.m_etNeuEM),
      m_etResNeuEM(rhs.m_etResNeuEM),
      m_trFlightPathSig(rhs.m_trFlightPathSig),
      m_etChrgEM01Trk(rhs.m_etChrgEM01Trk),
      m_etResChrgEMTrk(rhs.m_etResChrgEMTrk),
      m_sumEM(rhs.m_sumEM),
      m_secVertex(rhs.m_secVertex ? std::make_unique<Trk::RecVertex>(*rhs.m_secVertex) : nullptr),
      m_pi0(rhs.m_pi0),
      m_etEflow(rhs.m_etEflow),
      m_sumEtCellsLArOverLeadTrackPt(rhs.m_sumEtCellsLArOverLeadTrackPt),
      m_hadronicLeak(rhs.m_hadronicLeak),
      m_secondaryMax(rhs.m_secondaryMax){
    }
 
    Tau1P3PDetails&
    Tau1P3PDetails::operator=(const Tau1P3PDetails& rhs){
      if (this != &rhs) {
        Tau1P3PDetails tmp(rhs);
        std::swap(*this, tmp);
      }
      return *this;
    }
 
    double Tau1P3PDetails :: parameter( TauJetParameters :: Detail detail ) const
    {
        switch( detail )
        {
            case TauJetParameters :: emRadius:
                return emRadius();
            case TauJetParameters :: isolationFraction:
                return isolationFraction();
            case TauJetParameters :: stripWidth2:
                return stripWidth2();
            case TauJetParameters :: annularIsolationFraction:
                return isolationFraction();
            case TauJetParameters :: etCaloAtEMScale:
                return etEMAtEMScale() + etHadAtEMScale();
            case TauJetParameters :: etChargedHadCells:
                return etChrgHAD();
            case TauJetParameters :: etOtherEMCells:
                return etIsolEM();
            case TauJetParameters :: etOtherHadCells:
                return etIsolHAD();
                //            case TauJetParameters :: sumPTtracks:
                //                return sumPTTracks();
            case TauJetParameters :: nAssocTracksCore:
                return (double)nAssocTracksCore();
            case TauJetParameters :: nAssocTracksIsol:
                return (double)nAssocTracksIsol();
            case TauJetParameters :: rWidth2Trk3P:
                return rWidth2Trk3P();
            case TauJetParameters :: massTrk3P:
                return massTrk3P();
            case TauJetParameters :: signDOTrk3P:
                return signD0Trk3P();
            default:
                return TauDetails::DEFAULT;
        }
    }
 
    bool Tau1P3PDetails :: hasParameter( TauJetParameters :: Detail detail ) const
    {
        switch( detail )
        {
            case TauJetParameters :: emRadius:
            case TauJetParameters :: isolationFraction:
            case TauJetParameters :: stripWidth2:
            case TauJetParameters :: annularIsolationFraction:
            case TauJetParameters :: etCaloAtEMScale:
            case TauJetParameters :: etChargedHadCells:
            case TauJetParameters :: etOtherEMCells:
            case TauJetParameters :: etOtherHadCells:
            case TauJetParameters :: nAssocTracksCore:
            case TauJetParameters :: nAssocTracksIsol:
            case TauJetParameters :: rWidth2Trk3P:
            case TauJetParameters :: massTrk3P:
            case TauJetParameters :: signDOTrk3P:
                return true;
            default:
                return false;
        }
    }
 
    const std::string Tau1P3PDetails::s_className =
        ClassName<Tau1P3PDetails>::name();
    const std::string& Tau1P3PDetails::className() const
    {
        return s_className;
    }
}