/cvmfs/atlas.cern.ch/repo/sw/ASG/AnalysisBase/2.4.31/xAODTau/xAODTau/versions/TauJetAuxContainer_v1.h

// Dear emacs, this is -*- c++ -*-
// $Id: TauJetAuxContainer_v1.h 629816 2014-11-20 14:55:09Z janus $
#ifndef XAODTAU_VERSIONS_TAUJETAUXCONTAINER_V1_H
#define XAODTAU_VERSIONS_TAUJETAUXCONTAINER_V1_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_v1 : public AuxContainerBase {

  public:
    TauJetAuxContainer_v1();

  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 > ptPanTauEFlowRecProto;
    // std::vector< float > etaPanTauEFlowRecProto;
    // std::vector< float > phiPanTauEFlowRecProto;
    // std::vector< float > mPanTauEFlowRecProto;
   
    // std::vector< float > ptPanTauEFlowRec;
    // std::vector< float > etaPanTauEFlowRec;
    // std::vector< float > phiPanTauEFlowRec;
    // std::vector< float > mPanTauEFlowRec;
   
    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< 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;
    std::vector< float > PanTauScore;


    //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 > conversionTrackLinks;
    std::vector< TrackLink_t > wideTrackLinks;
    std::vector< TrackLink_t > otherTrackLinks;


    std::vector< std::vector < int > > trackFilterPass;
    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 > charged_PFOLinks;// (Links to default charged PFOs)
    std::vector< PFOLink_t > neutral_PFOLinks;// (Links to default all neutral PFOs)
    std::vector< PFOLink_t > pi0_PFOLinks;// (Links to default pi0 neutral PFOs)
    std::vector< PFOLink_t > cellBased_Charged_PFOLinks;// (Links to cell-based charged PFOs)
    std::vector< PFOLink_t > cellBased_Neutral_PFOLinks;// (Links to cell-based all neutral PFOs)
    std::vector< PFOLink_t > cellBased_Pi0_PFOLinks;// (Links to cell-based pi0 neutral PFOs)
    std::vector< PFOLink_t > eflowRec_Charged_PFOLinks;// (Links to eflowRec charged PFOs)
    std::vector< PFOLink_t > eflowRec_Neutral_PFOLinks;// (Links to eflowRec all neutral PFOs)
    std::vector< PFOLink_t > eflowRec_Pi0_PFOLinks;// (Links to eflowRec pi0 neutral PFOs)
    std::vector< PFOLink_t > shot_PFOLinks;// (Links to shots)
    std::vector< PFOLink_t > hadronic_PFOLinks;// (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;


    //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;
    
    // //!PanTau variables when using eflowRec pfos
    // //Flag whether this seed has pantau info
    // std::vector< int > pantau_eflowRecInput_isPanTauCandidate;
    // //decay modes (input mode [proto] and output mode )
    // std::vector< int > pantau_eflowRecInput_DecayModeProto;
    // std::vector< int > pantau_eflowRecInput_DecayMode;
    // //BDT output distributions
    // std::vector< float > pantau_eflowRecInput_BDTValue_1p0n_vs_1p1n;
    // std::vector< float > pantau_eflowRecInput_BDTValue_1p1n_vs_1pXn;
    // std::vector< float > pantau_eflowRecInput_BDTValue_3p0n_vs_3pXn;
    // //Variables used in BDTs
    // //NOTE: They are different from the CellBased ones in general!
    // std::vector< int > pantau_eflowRecInput_BDTVar_Basic_NPi0NeutConsts;
    // std::vector< int > pantau_eflowRecInput_BDTVar_Basic_NNeutralConsts;
    // std::vector< float > pantau_eflowRecInput_BDTVar_Charged_HLV_SumPt;
    // std::vector< float > pantau_eflowRecInput_BDTVar_Charged_Ratio_EtOverEtAllConsts;
    // std::vector< float > pantau_eflowRecInput_BDTVar_Neutral_HLV_SumM;
    // std::vector< float > pantau_eflowRecInput_BDTVar_Neutral_PID_BDTValues_EtSort_1;
    // std::vector< float > pantau_eflowRecInput_BDTVar_Neutral_PID_BDTValues_BDTSort_2;
    // std::vector< float > pantau_eflowRecInput_BDTVar_Neutral_Ratio_EtOverEtAllConsts;
    // std::vector< float > pantau_eflowRecInput_BDTVar_Neutral_Mean_DRToLeading_WrtEtAllConsts;
    // std::vector< float > pantau_eflowRecInput_BDTVar_Combined_DeltaR1stNeutralTo1stCharged;

    
    
    
  }; // class TauJetAuxContainer_v1

} // namespace xAOD

// Set up a CLID and StoreGate inheritance for the class:
#ifndef XAOD_STANDALONE
#include "SGTools/BaseInfo.h"
SG_BASE( xAOD::TauJetAuxContainer_v1, xAOD::AuxContainerBase );
#endif // not XAOD_STANDALONE


#endif // XAODTAU_VERSIONS_TAUJETAUXCONTAINER_V1_H

---