TauJetAuxContainer_v2.h

Go to the documentation of this file.

// Dear emacs, this is -*- c++ -*-
 
/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
// $Id: TauJetAuxContainer_v2.h 750177 2016-05-27 08:12:45Z krasznaa $
#ifndef XAODTAU_VERSIONS_TAUJETAUXCONTAINER_V2_H
#define XAODTAU_VERSIONS_TAUJETAUXCONTAINER_V2_H
 
// System include(s):
#include <vector>
extern "C" {
#   include <stdint.h>
}
 
//local includes
#include "xAODTau/TauDefs.h"
 
// EDM include(s):
#include "xAODCore/AuxContainerBase.h"
#include "xAODTracking/TrackParticleContainer.h"
#include "xAODTracking/VertexContainer.h"
#include "xAODJet/JetContainer.h"
#include "xAODPFlow/PFOContainer.h"
 
namespace xAOD {
 

  class TauJetAuxContainer_v2 : public AuxContainerBase {
 
  public:
    TauJetAuxContainer_v2();
 
  private:
    std::vector< float > pt;
    std::vector< float > eta;
    std::vector< float > phi;
    std::vector< float > m;
 
    std::vector< float > ptJetSeed;
    std::vector< float > etaJetSeed;
    std::vector< float > phiJetSeed;
    std::vector< float > mJetSeed;
 
    std::vector< float > ptDetectorAxis;
    std::vector< float > etaDetectorAxis;
    std::vector< float > phiDetectorAxis;
    std::vector< float > mDetectorAxis;
 
    std::vector< float > ptIntermediateAxis;
    std::vector< float > etaIntermediateAxis;
    std::vector< float > phiIntermediateAxis;
    std::vector< float > mIntermediateAxis;
 
    std::vector< float > ptTauEnergyScale;
    std::vector< float > etaTauEnergyScale;
    std::vector< float > phiTauEnergyScale;
    std::vector< float > mTauEnergyScale;
 
    std::vector< float > ptTauEtaCalib;
    std::vector< float > etaTauEtaCalib;
    std::vector< float > phiTauEtaCalib;
    std::vector< float > mTauEtaCalib;
   
    std::vector< float > ptPanTauCellBasedProto;
    std::vector< float > etaPanTauCellBasedProto;
    std::vector< float > phiPanTauCellBasedProto;
    std::vector< float > mPanTauCellBasedProto;
   
    std::vector< float > ptPanTauCellBased;
    std::vector< float > etaPanTauCellBased;
    std::vector< float > phiPanTauCellBased;
    std::vector< float > mPanTauCellBased;
 
    std::vector< float > ptTrigCaloOnly;
    std::vector< float > etaTrigCaloOnly;
    std::vector< float > phiTrigCaloOnly;
    std::vector< float > mTrigCaloOnly;
   
        
             
 
 
    std::vector< unsigned int > ROIWord;
    std::vector< float > charge;
 
    std::vector< uint32_t > vetoFlags;
    std::vector< uint32_t > isTauFlags;
    
    std::vector< float > Likelihood;
    std::vector< float > SafeLikelihood;
    std::vector< float > BDTJetScore;
    std::vector< float > BDTEleScore;
 
    //additional discriminant output
    std::vector< float > BDTJetScoreSigTrans;
    std::vector< float > BDTJetScoreBkgTrans;
 
    typedef std::vector< ElementLink< TrackParticleContainer > > TrackLink_t;
    std::vector< TrackLink_t > trackLinks;
    std::vector< TrackLink_t > wideTrackLinks;
    std::vector< TrackLink_t > otherTrackLinks;
 
 
    std::vector< std::vector < uint8_t > > trackFlags;
    std::vector< int > trackFilterProngs;
    std::vector< int > trackFilterQuality;             
    std::vector< float > pi0ConeDR;             
 
    std::vector< std::vector < float > > trackEtaStrip;
    std::vector< std::vector < float > > trackPhiStrip;
 
 
    typedef ElementLink< xAOD::JetContainer > JetLink_t;
    std::vector< JetLink_t > jetLink;
 
    typedef ElementLink< xAOD::VertexContainer > VertexLink_t;
    std::vector< VertexLink_t > vertexLink;
    std::vector< VertexLink_t > secondaryVertexLink;
 
 
    // //pfo links for substructure
    typedef std::vector< ElementLink< PFOContainer > > PFOLink_t;
 
    std::vector< PFOLink_t > chargedPFOLinks;// (Links to default charged PFOs)
    std::vector< PFOLink_t > neutralPFOLinks;// (Links to default all neutral PFOs)
    std::vector< PFOLink_t > pi0PFOLinks;// (Links to default pi0 neutral PFOs)
    std::vector< PFOLink_t > protoChargedPFOLinks;// (Links to cell-based charged PFOs)
    std::vector< PFOLink_t > protoNeutralPFOLinks;// (Links to cell-based all neutral PFOs)
    std::vector< PFOLink_t > protoPi0PFOLinks;// (Links to cell-based pi0 neutral PFOs)
    std::vector< PFOLink_t > shotPFOLinks;// (Links to shots)
    std::vector< PFOLink_t > hadronicPFOLinks;// (Links to hadronic pfos)
 
 
    //set of details variables
    std::vector< float > ipZ0SinThetaSigLeadTrk;
    std::vector< float > etOverPtLeadTrk;
    std::vector< float > leadTrkPt;
    std::vector< float > ipSigLeadTrk;
    std::vector< float > massTrkSys;
    std::vector< float > trkWidth2;
    std::vector< float > trFlightPathSig;
    std::vector< float > ele_E237E277;
    std::vector< float > ele_PresamplerFraction;
    std::vector< float > ele_ECALFirstFraction;
    //for topocluster ID variables
    std::vector< int > numCells;
    std::vector< int > numTopoClusters;
    std::vector< float > numEffTopoClusters;
    std::vector< float > topoInvMass;
    std::vector< float > effTopoInvMass;
    std::vector< float > topoMeanDeltaR;
    std::vector< float > effTopoMeanDeltaR;
    std::vector< float > EMRadius;
    std::vector< float > hadRadius;
    std::vector< float > etEMAtEMScale;
    std::vector< float > etHadAtEMScale;
    std::vector< float > isolFrac;
    std::vector< float > centFrac;
    std::vector< float > stripWidth2;
    std::vector< int > nStrip;
    std::vector< float > etEMCalib;
    std::vector< float > etHadCalib;
    std::vector< float > seedCalo_eta;
    std::vector< float > seedCalo_phi;
    std::vector< float > trkAvgDist;
    std::vector< float > trkRmsDist;
    std::vector< float > lead2ClusterEOverAllClusterE;
    std::vector< float > lead3ClusterEOverAllClusterE;
    std::vector< float > caloIso;
    std::vector< float > caloIsoCorrected;
    std::vector< float > dRmax;
    
    std::vector< float > secMaxStripEt;
    std::vector< float > sumEMCellEtOverLeadTrkPt;
    std::vector< float > hadLeakEt;

    std::vector< float > EM_TES_scale;
    std::vector< float > LC_TES_precalib;
 
    std::vector< float > TESOffset;
    std::vector< float > TESCalibConstant;

    //
    // Cell-based energies in rings getters/setters
    //
    std::vector< float > cellBasedEnergyRing1;
    std::vector< float > cellBasedEnergyRing2;
    std::vector< float > cellBasedEnergyRing3;
    std::vector< float > cellBasedEnergyRing4;
    std::vector< float > cellBasedEnergyRing5;
    std::vector< float > cellBasedEnergyRing6;
    std::vector< float > cellBasedEnergyRing7;
    std::vector< float > TRT_NHT_OVER_NLT;
    std::vector< float > TauJetVtxFraction;
 
 
    //generic substructure details
    std::vector< float > etEflow;
    std::vector< float > mEflow;
    std::vector< float > ptRatioEflow;
    std::vector< int > nPi0;
    std::vector< int > nCharged;
 
    std::vector< float > etEflowTopo;
    std::vector< float > mEflowTopo;
    std::vector< float > ptRatioEflowTopo;
    std::vector< int > nPi0Topo;
    std::vector< int > nChargedTopo;
 
    std::vector< float > mEflowApprox;
    std::vector< float > ptRatioEflowApprox;
 
    std::vector< float > innerTrkAvgDist;
 
    std::vector< float > SumPtTrkFrac;

    std::vector< float > etOverPtLeadTrkCorrected;
    std::vector< float > ipSigLeadTrkCorrected;
    std::vector< float > trFlightPathSigCorrected;
    std::vector< float > massTrkSysCorrected;
    std::vector< float > dRmaxCorrected;
    std::vector< float > ChPiEMEOverCaloEMECorrected;
    std::vector< float > EMPOverTrkSysPCorrected;
    std::vector< float > ptRatioEflowApproxCorrected;
    std::vector< float > mEflowApproxCorrected;
    std::vector< float > centFracCorrected;
    std::vector< float > innerTrkAvgDistCorrected;
    std::vector< float > SumPtTrkFracCorrected;
 
 
 
    //trigger details
    std::vector< float > PSSFraction;
    std::vector< float > ChPiEMEOverCaloEME;
    std::vector< float > EMPOverTrkSysP;
 
    //panTau details
    //Flag whether this seed has pantau info
    std::vector< int > pantau_CellBasedInput_isPanTauCandidate;
    //decay modes (input mode [proto] and output mode )
    std::vector< int > pantau_CellBasedInput_DecayModeProto;
    std::vector< int > pantau_CellBasedInput_DecayMode;
    //BDT output distributions
    std::vector< float > pantau_CellBasedInput_BDTValue_1p0n_vs_1p1n;
    std::vector< float > pantau_CellBasedInput_BDTValue_1p1n_vs_1pXn;
    std::vector< float > pantau_CellBasedInput_BDTValue_3p0n_vs_3pXn;
    //Variables used in BDTs
    std::vector< int > pantau_CellBasedInput_BDTVar_Basic_NNeutralConsts;
    std::vector< float > pantau_CellBasedInput_BDTVar_Charged_JetMoment_EtDRxTotalEt;
    std::vector< float > pantau_CellBasedInput_BDTVar_Charged_StdDev_Et_WrtEtAllConsts;
    std::vector< float > pantau_CellBasedInput_BDTVar_Neutral_HLV_SumM;
    std::vector< float > pantau_CellBasedInput_BDTVar_Neutral_PID_BDTValues_BDTSort_1;
    std::vector< float > pantau_CellBasedInput_BDTVar_Neutral_PID_BDTValues_BDTSort_2;
    std::vector< float > pantau_CellBasedInput_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts;
    std::vector< float > pantau_CellBasedInput_BDTVar_Neutral_Ratio_EtOverEtAllConsts;
    std::vector< float > pantau_CellBasedInput_BDTVar_Neutral_Shots_NPhotonsInSeed;
    std::vector< float > pantau_CellBasedInput_BDTVar_Combined_DeltaR1stNeutralTo1stCharged;
    std::vector< float > pantau_CellBasedInput_BDTVar_Charged_HLV_SumM;
 
    
    
    
  }; // class TauJetAuxContainer_v2
 
} // namespace xAOD
 
// Set up the StoreGate inheritance for the class:
#include "xAODCore/BaseInfo.h"
SG_BASE( xAOD::TauJetAuxContainer_v2, xAOD::AuxContainerBase );
 
#endif // XAODTAU_VERSIONS_TAUJETAUXCONTAINER_V2_H