/cvmfs/atlas.cern.ch/repo/sw/ASG/AnalysisBase/2.4.31/FsrUtils/FsrUtils/FsrPhotonTool.h

#ifndef FSRUTILS_FsrPhotonTool_H
#define FSRUTILS_FsrPhotonTool_H

// Framework include(s):
#include "AsgTools/AsgTool.h"
#include "AsgTools/ToolHandle.h"

// Local include(s):
#include "FsrUtils/IFsrPhotonTool.h"

// Forward declarations
namespace CP
{
    class IIsolationSelectionTool;
}

namespace FSR {


   class FsrPhotonTool : public virtual IFsrPhotonTool,
                         public asg::AsgTool {

      ASG_TOOL_CLASS( FsrPhotonTool, FSR::IFsrPhotonTool )

   public:
      FsrPhotonTool( const std::string& name );


      virtual StatusCode initialize();



      virtual CP::CorrectionCode getFsrPhoton(const xAOD::IParticle* part, FsrCandidate& candidate,
                                              const xAOD::PhotonContainer* photons = 0,
                                              const xAOD::ElectronContainer* electrons = 0);

      virtual std::vector<FsrCandidate>* getFsrCandidateList(const xAOD::IParticle* part,
                                                             const xAOD::PhotonContainer* photons = 0,
                                                             const xAOD::ElectronContainer* electrons = 0);

      virtual std::vector<FsrCandidate>* getFarFsrCandidateList(const xAOD::IParticle* part,
                                                                const xAOD::PhotonContainer* photons_cont);

      virtual std::vector<FsrCandidate>* getNearFsrCandidateList(const xAOD::Muon* part,
                                                                const xAOD::PhotonContainer* photons_cont,
                                                                const xAOD::ElectronContainer* electrons_cont);

   private:
      std::vector<FsrCandidate>* sortFsrCandidates( std::vector< std::pair <const xAOD::IParticle*, double> > FsrCandList,
                                                     std::string option="ET");
      bool isOverlap(const xAOD::Electron_v1* electron, std::vector< std::pair <const xAOD::IParticle*, double> > phfsr,
              unsigned int nofPhFsr);
      double deltaR(float muonEta, float muonPhi, float phEta, float phPhi) const;
      double deltaPhi(float phi1, float phi2) const;
      static bool compareEt(FsrCandidate c1, FsrCandidate c2) { return (c1.Et > c2.Et); }

      //const xAOD::PhotonContainer* m_photons;
      //const xAOD::ElectronContainer* m_electrons;

      double m_high_et_min;
      double m_overlap_el_ph;
      double m_overlap_el_mu;

      double m_far_fsr_drcut;
      double m_far_fsr_etcut;
      std::string m_far_fsr_isoWorkingPoint;

      double m_drcut;
      double m_etcut;
      double m_f1cut;
      double m_topo_drcut;
      double m_topo_f1cut;


      std::vector<FsrCandidate> m_fsrPhotons;
      FsrCandidate::FsrType     m_fsr_type;

      ToolHandle<CP::IIsolationSelectionTool> m_isoSelTool;

   }; // class FsrPhotonTool

} // namespace FSR

#endif // FSRUTILS_FsrPhotonTool_H

---