TauJetAuxContainer_v3.h

Go to the documentation of this file.

// Dear emacs, this is -*- c++ -*-
 
/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef XAODTAU_VERSIONS_TAUJETAUXCONTAINER_V3_H
#define XAODTAU_VERSIONS_TAUJETAUXCONTAINER_V3_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 "xAODTau/TauTrackContainer.h"
#include "xAODTracking/TrackParticleContainer.h"
#include "xAODTracking/VertexContainer.h"
#include "xAODJet/JetContainer.h"
#include "xAODPFlow/PFOContainer.h"
 
namespace xAOD {
 
 
  class TauJetAuxContainer_v3 : public AuxContainerBase {
 
  public:
    TauJetAuxContainer_v3();
 
  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; // same as IntermediateAxis
    std::vector< float > phiTauEnergyScale; // as IntermediateAxis
    std::vector< float > mTauEnergyScale; // is 0
 
    std::vector< float > ptTauEtaCalib;
    std::vector< float > etaTauEtaCalib;
    std::vector< float > phiTauEtaCalib; // same as IntermediateAxis
    std::vector< float > mTauEtaCalib; // is 0
   
    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< float > ptFinalCalib;
    std::vector< float > etaFinalCalib;
    std::vector< float > phiFinalCalib;
    std::vector< float > mFinalCalib;
 
 
    std::vector< unsigned int > ROIWord;
    std::vector< float > charge;
 
    //    std::vector< uint32_t > vetoFlags;//r21 cleanup
    std::vector< uint32_t > isTauFlags;
    
    // std::vector< float > Likelihood;
    // std::vector< float > SafeLikelihood;
    std::vector< float > BDTJetScore;
    std::vector< float > BDTEleScore;
    std::vector< float > EleMatchLikelihoodScore;//from TauEleOLRDecorator
 
    //additional discriminant output
    std::vector< float > BDTJetScoreSigTrans;
    // std::vector< float > BDTJetScoreBkgTrans;
 
    typedef std::vector< ElementLink< TauTrackContainer > > TauTrackLink_t;
    std::vector< TauTrackLink_t > tauTrackLinks;
 
    //simply nTracks(TauJetParameters::TauTrackFlag::classifiedCharged)
    //after the first pass of TauTrackFinder
    std::vector< int > nChargedTracks;
    std::vector< int > nIsolatedTracks;
    std::vector< int > nAllTracks;
    std::vector< int > nLargeRadiusTracks;
 
    typedef std::vector< ElementLink< IParticleContainer > > IParticleLink_t;
    std::vector< IParticleLink_t > clusterLinks;//actually xAOD::CaloClusters
    std::vector< IParticleLink_t > pi0Links;//acutally xAOD::Particles
 
 
    std::vector< int > trackFilterProngs;
    std::vector< int > trackFilterQuality;             
    std::vector< float > pi0ConeDR;             
 
    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; //r21 cleanup
    // std::vector< float > ele_PresamplerFraction; //r21 cleanup
    // std::vector< float > ele_ECALFirstFraction; //r21 cleanup
    //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;//r21 cleanup
    //    std::vector< float > etHadCalib;//r21 cleanup
    //std::vector< float > seedCalo_eta; //r21 cleanup DetectorAxis eta
    //    std::vector< float > seedCalo_phi;//r21 cleanup DetectorAxis 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;//r21 cleanup
    //std::vector< float > LC_TES_precalib; //r21 cleanup
 
    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;//r21 cleanup
    //    std::vector< float > mEflow;//r21 cleanup
    //    std::vector< float > ptRatioEflow;//r21 cleanup
    //    std::vector< int > nPi0;//r21 cleanup
    std::vector< int > nCharged; //used in PanTau/PanTauAlgs/src/Tool_DetailsArranger.cxx
 
    //    std::vector< float > etEflowTopo;//r21 cleanup
    //    std::vector< float > mEflowTopo;//r21 cleanup
    //    std::vector< float > ptRatioEflowTopo;//r21 cleanup
    //    std::vector< int > nPi0Topo;//r21 cleanup
    //    std::vector< int > nChargedTopo;//r21 cleanup
 
    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_isPanTauCandidate;
    //decay modes (input mode [proto] and output mode )
    std::vector< int > PanTau_DecayModeProto;
    std::vector< int > PanTau_DecayMode;
    //BDT output distributions
    std::vector< float > PanTau_BDTValue_1p0n_vs_1p1n;
    std::vector< float > PanTau_BDTValue_1p1n_vs_1pXn;
    std::vector< float > PanTau_BDTValue_3p0n_vs_3pXn;
    //Variables used in BDTs
    std::vector< int > PanTau_BDTVar_Basic_NNeutralConsts;
    std::vector< float > PanTau_BDTVar_Charged_JetMoment_EtDRxTotalEt;
    std::vector< float > PanTau_BDTVar_Charged_StdDev_Et_WrtEtAllConsts;
    std::vector< float > PanTau_BDTVar_Neutral_HLV_SumM;
    std::vector< float > PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_1;
    std::vector< float > PanTau_BDTVar_Neutral_PID_BDTValues_BDTSort_2;
    std::vector< float > PanTau_BDTVar_Neutral_Ratio_1stBDTEtOverEtAllConsts;
    std::vector< float > PanTau_BDTVar_Neutral_Ratio_EtOverEtAllConsts;
    std::vector< float > PanTau_BDTVar_Neutral_Shots_NPhotonsInSeed;
    std::vector< float > PanTau_BDTVar_Combined_DeltaR1stNeutralTo1stCharged;
    std::vector< float > PanTau_BDTVar_Charged_HLV_SumM;
 
    // RNN tau ID
    std::vector< float > RNNJetScore;
    std::vector< float > RNNJetScoreSigTrans;
 
    // RNN tau eVeto
    std::vector< float > RNNEleScore;
    std::vector< float > RNNEleScoreSigTrans;
    
  }; // class TauJetAuxContainer_v3
 
} // namespace xAOD
 
// Set up the StoreGate inheritance for the class:
#include "xAODCore/BaseInfo.h"
SG_BASE( xAOD::TauJetAuxContainer_v3, xAOD::AuxContainerBase );
 
#endif // XAODTAU_VERSIONS_TAUJETAUXCONTAINER_V3_H