/cvmfs/atlas.cern.ch/repo/sw/ASG/AnalysisBase/2.4.31/tauRecTools/tauRecTools/CaloClusterVariables.h

#ifndef CALOCLUSTERVARIABLES_H
#define CALOCLUSTERVARIABLES_H

#include <vector>
#ifndef XAOD_ANALYSIS
#include "CaloUtils/CaloVertexedCluster.h"
typedef xAOD::CaloVertexedCluster CaloVertexedClusterType;
#else
#include "xAODCaloEvent/CaloVertexedTopoCluster.h"
typedef xAOD::CaloVertexedTopoCluster CaloVertexedClusterType;
#endif
//#include "CaloEvent/CaloVertexedCluster.h"
#include "CxxUtils/fpcompare.h"
#include "xAODTau/TauJet.h"

#include "TLorentzVector.h"
#include "xAODCaloEvent/CaloVertexedTopoCluster.h"


class CaloClusterVariables {
public:

    static const double DEFAULT;

    CaloClusterVariables();

    ~CaloClusterVariables() {
    }

    bool update(const xAOD::TauJet& pTau, bool inAODmode=false); 

    void setVertexCorrection(bool flag) {m_doVertexCorrection=flag;}

    // ID Variables
    unsigned int numConstituents() { return (unsigned int) m_numConstit; }

    double totalMass()     { return m_totMass; }
    double effectiveMass() { return m_effMass; }
    
    double effectiveNumConstituents()  { return m_effNumConstit; }
    int effectiveNumConstituents_int() { return m_effNumConstit_int; }

    double averageEffectiveRadius() { return m_aveEffRadius; }
    double averageRadius()          { return m_aveRadius; }

    // Energy Variables
    double totalEnergy()     { return m_totEnergy; }
    double effectiveEnergy() { return m_effEnergy; }

    //cells
    unsigned int numCells() { return m_numCells; }

private:
    int m_numConstit;
    int m_effNumConstit_int;
    double m_effNumConstit;
    double m_aveRadius;
    double m_aveEffRadius;
    double m_totMass;
    double m_effMass;
    double m_totEnergy;
    double m_effEnergy;
    unsigned int m_numCells;

    TLorentzVector calculateTauCentroid(int nConst, const std::vector<CaloVertexedClusterType>& constituents);
    
    bool m_doVertexCorrection;
};

//-------------------------------------------------------------------------
//-------------------------------------------------------------------------
struct CaloClusterCompare { 
  bool operator()(const CaloVertexedClusterType& left, const CaloVertexedClusterType& right) {
    //volatile double leftE = left.e();
    //volatile double rightE = right.e();
    return CxxUtils::fpcompare::greater (left.e(),right.e());
  }
};

struct DefCaloClusterCompare { 
  bool operator()(const xAOD::CaloCluster* left, const xAOD::CaloCluster* right) {
    //volatile double leftE = left.e();
    //volatile double rightE = right.e();
    return CxxUtils::fpcompare::greater (left->e(),right->e());
  }
};

#endif

---