/cvmfs/atlas.cern.ch/repo/sw/ASG/AnalysisBase/2.4.31/PhotonEfficiencyCorrection/PhotonEfficiencyCorrection/AsgPhotonEfficiencyCorrectionTool.h

// Dear emacs, this is -*-c++-*-

#ifndef __ASGPHOTONEFFICIENCYCORRECTIONTOOL__
#define __ASGPHOTONEFFICIENCYCORRECTIONTOOL__

// STL includes
#include <vector>
#include <string>
#include <fstream>      // std::ifstream

// Utility includes
#include "boost/unordered_map.hpp"
#include "boost/algorithm/string.hpp" // this one to replace std::string names

// Include the return object and the underlying ROOT tool
#include "PATCore/TResult.h"

//xAOD includes
#include "AsgTools/AsgTool.h"
#include "PATInterfaces/ISystematicsTool.h"
#include "PATInterfaces/SystematicRegistry.h"
#include "PATInterfaces/CorrectionCode.h"
#include "PhotonEfficiencyCorrection/TPhotonEfficiencyCorrectionTool.h"
#include "PhotonEfficiencyCorrection/IAsgPhotonEfficiencyCorrectionTool.h"

#include "xAODEgamma/Egamma.h"




class AsgPhotonEfficiencyCorrectionTool
  : virtual public IAsgPhotonEfficiencyCorrectionTool,
    virtual public CP::ISystematicsTool,
            public asg::AsgTool
{
  ASG_TOOL_CLASS2(AsgPhotonEfficiencyCorrectionTool, IAsgPhotonEfficiencyCorrectionTool, CP::ISystematicsTool )

public:
  AsgPhotonEfficiencyCorrectionTool ( const std::string myname );

  virtual ~AsgPhotonEfficiencyCorrectionTool();


public:
  virtual StatusCode initialize();

  virtual StatusCode finalize();


  // Main methods from IUserDataCalcTool
public:
  const Root::TResult& calculate( const xAOD::IParticle* part ) const;
  const Root::TResult& calculate( const xAOD::Egamma* egam ) const;
  const Root::TResult& calculate( const xAOD::Egamma& egam ) const{  
            return calculate(&egam);} // pass the Egamma obj by reference

  virtual CP::CorrectionCode getEfficiencyScaleFactor(const xAOD::Egamma& inputObject, double& efficiencyScaleFactor) const;
  virtual CP::CorrectionCode getEfficiencyScaleFactorError(const xAOD::Egamma& inputObject, double& efficiencyScaleFactorError) const;
  virtual CP::CorrectionCode applyEfficiencyScaleFactor(xAOD::Egamma& inputObject) const;


  virtual bool isAffectedBySystematic( const CP::SystematicVariation& systematic ) const;
  
  virtual CP::SystematicSet affectingSystematics() const ;
  
  virtual CP::SystematicSet recommendedSystematics() const ;
  
  const CP::SystematicSet& appliedSystematics() const {
    assert (m_appliedSystematics != nullptr);
    return *m_appliedSystematics;
  }
  
  virtual CP::SystematicCode applySystematicVariation ( const CP::SystematicSet& systConfig );

  CP::SystematicCode registerSystematics();

  // Private member variables
private:
  Root::TPhotonEfficiencyCorrectionTool* m_rootTool_unc;
  Root::TPhotonEfficiencyCorrectionTool* m_rootTool_con;
  Root::TPhotonEfficiencyCorrectionTool* m_rootTool_uncRadZ;
  Root::TPhotonEfficiencyCorrectionTool* m_rootTool_conRadZ;
  
  Root::TResult m_resultDummy;

  boost::unordered_map<CP::SystematicSet, CP::SystematicSet> m_systFilter;
  
  CP::SystematicSet* m_appliedSystematics = nullptr;
  
  // The prefix for the systematic name
  std::string m_sysSubstring;
  std::string m_sysSubstringRadZ;
  
  // Get the correction filename from the map
  std::string getFileName(std::string isoWP, bool isConv, std::string sufix);
  
  // Properties
  
  std::string m_corrFileNameConv;
  std::string m_corrFileNameUnconv;
 
  std::string m_resultPrefix;

  std::string m_resultName;

  int m_dataTypeOverwrite;
  
  std::string m_isoWP;
  
  std::string m_mapFile;  
  
  float m_Threshold_lowPT;
  float m_Threshold_highPT;
  bool  m_UseRadZ_mediumPT;



}; // End: class definition


#endif

---