/cvmfs/atlas.cern.ch/repo/sw/ASG/AnalysisBase/2.4.31/JetMomentTools/JetMomentTools/JetBadChanCorrTool.h

// this file is -*- C++ -*-
#ifndef XAOD_ANALYSIS

#ifndef JETMOMENTTOOLS_JETBADCHANCORRTOOL_H
#define JETMOMENTTOOLS_JETBADCHANCORRTOOL_H

#include "JetRec/JetModifierBase.h"

#include "JetRecCalo/MissingCellListTool.h"

#include "CaloIdentifier/CaloCell_ID.h"

#include "TH1D.h"


class Identifier;
class CaloDetDescrManager;

class ITileBadChanTool;
class ILArBadChanTool;
class ITHistSvc;

class JetBadChanCorrTool: public JetModifierBase
{
  ASG_TOOL_CLASS0(JetBadChanCorrTool)
public:
  JetBadChanCorrTool(const std::string& name);
  
  virtual ~JetBadChanCorrTool();
  
  virtual StatusCode initialize();
 
  virtual int modifyJet( xAOD::Jet& pJet) const ;
  
  virtual StatusCode setupEvent();


protected:
  int correctionFromCellsInJet( xAOD::Jet* j, const jet::CaloCellFastMap * badCellMap) const;
  int correctionFromCellsInCone( xAOD::Jet* j, const jet::CaloCellFastMap * badCellMap) const;

  int correctionFromClustersBadCells( xAOD::Jet* j) const;

 private:

  ServiceHandle<ITHistSvc> m_thistSvc;

  // const CaloDetDescrManager* m_caloDDM;
  // const CaloCell_ID*       m_calo_id;

  // limit to calculate moments
  int m_nBadCellLimit;

  

  

  // for jet-level correction
  std::string m_streamName;
  std::string m_profileName;
  std::string m_profileTag;

  bool m_useCone;
  //  double m_coneDr;
  
  // for cell-level correction
  bool m_useCalibScale;
  //ToolHandle<IHadronicCalibrationTool> m_calibTool;

  //ToolHandle<IMissingCellListTool> m_missingCellToolHandle;
  //MissingCellListTool* m_missingCellTool;
  std::string m_missingCellMapName;

  bool m_forceMissingCellCheck;

  bool m_useClusters;

  // jet profiles
  class ProfileData {
  public:
    ProfileData(TH1D* th, int sample,
        double ptMin=0, double ptMax=9999,
        double etaMin=0, double etaMax=5.0,
        double phiMin=-M_PI, double phiMax=M_PI):
      m_th(th), m_sample(sample), 
      m_ptMin(ptMin), m_ptMax(ptMax), 
      m_etaMin(etaMin), m_etaMax(etaMax),
      m_phiMin(phiMin), m_phiMax(phiMax) {}

    virtual ~ProfileData(){}

    bool match(double pt, int sample, double eta, double phi) const {
      return ( pt>=m_ptMin && pt<m_ptMax 
           && sample==m_sample
           && fabs(eta)>=m_etaMin && fabs(eta)<m_etaMax 
           && phi>=m_phiMin && phi<m_phiMax);
    }
    
    double frac(double dr) const {
      int idr = m_th->FindBin(dr);
      return m_th->GetBinContent(idr);
    }
  private:
    TH1D* m_th;
    int m_sample;
    double m_ptMin;
    double m_ptMax;
    double m_etaMin;
    double m_etaMax;
    double m_phiMin;
    double m_phiMax;
  };
  std::vector<ProfileData> m_profileDatas[CaloCell_ID::Unknown];//24
  
  double getProfile(double pt, double dr, int sample, double eta, double phi) const;

};
#endif 

#endif 

// DoxygenDocumentation

---