HIJetClusterSubtractorTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
// HIJetClusterSubtractorTool.h
 
#ifndef HIJETREC_HIJETCLUSTERSUBTRACTORTOOL_H
#define HIJETREC_HIJETCLUSTERSUBTRACTORTOOL_H
 
 
#include "HIJetRec/HIJetSubtractorToolBase.h"
#include <map>
#include "AsgTools/AsgTool.h"
 
 
class TH3F;
 
class HIJetClusterSubtractorTool : public HIJetSubtractorToolBase
{
  ASG_TOOL_CLASS(HIJetClusterSubtractorTool, IHISubtractorTool);
 
  // FIX
  //needs lookup for static geometric factors!
 
public:
  HIJetClusterSubtractorTool(const std::string& myname);
  virtual void subtract(xAOD::IParticle::FourMom_t&, const xAOD::IParticle*, const xAOD::HIEventShapeContainer*, const HIEventShapeIndex*, const ToolHandle<IHIUEModulatorTool>&, const xAOD::HIEventShape* eshape) const override;
  virtual void subtractWithMoments(xAOD::CaloCluster*, const xAOD::HIEventShapeContainer*, const HIEventShapeIndex* index, const ToolHandle<IHIUEModulatorTool>&, const xAOD::HIEventShape* eshape ) const override;
  virtual void updateUsingCluster(xAOD::HIEventShapeContainer* shape, const HIEventShapeIndex* index, const xAOD::CaloCluster* cl) const override;
 
  virtual StatusCode initialize() override;
  virtual StatusCode initializeTool();
 
  private:
 
  bool m_init;
 
  float getWeight(float eta, float phi, int sample) const;
  float getWeightEta(float eta, float phi, int sample) const;
  float getWeightPhi(float eta, float phi, int sample) const;
  void updateSlice(xAOD::HIEventShape* slice, float ET, float phi0, float area_cluster) const;
 
  Gaudi::Property< std::string > m_inputFile { this, "InputFile", "cluster.geo.HIJING_2018.root", "File containing cluster geometric moments." };
 
  Gaudi::Property< std::string > m_configDir { this, "ConfigDir", "HIJetCorrection/", "Directory containing configuration file." };
 
  TH3F* m_h3W;
  TH3F* m_h3Eta;
  TH3F* m_h3Phi;
 
  Gaudi::Property< bool > m_useSamplings { this, "UseSamplings", true, "Boolean for samplings use" };
 
};
 
#endif