PFEnergyCalibrationHF.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_PFClusterTools_PFEnergyCalibrationHF_h
#define RecoParticleFlow_PFClusterTools_PFEnergyCalibrationHF_h 

#include "DataFormats/ParticleFlowReco/interface/PFCluster.h"

// -*- C++ -*-
//
// Package:    PFClusterTools
// Class:      PFEnergyCalibrationHF
// 
#include <iostream>
#include "FWCore/ParameterSet/interface/ParameterSet.h"

class PFEnergyCalibrationHF 
{
 public:
  PFEnergyCalibrationHF(); // default constructor;
                         // needed by PFRootEvent
 explicit PFEnergyCalibrationHF(bool calibHF_use, 
                         const std::vector<double>& calibHF_eta_step,
                         const std::vector<double>& calibHF_a_EMonly,
                         const std::vector<double>& calibHF_b_HADonly,
                         const std::vector<double>& calibHF_a_EMHAD,
                         const std::vector<double>& calibHF_b_EMHAD);
  
  ~PFEnergyCalibrationHF();


  double energyEm(double uncalibratedEnergyECAL, 
          double eta, double phi);

  // HCAL only calibration
  double energyHad(double uncalibratedEnergyHCAL, 
           double eta, double phi) ;
  
  
  // ECAL+HCAL (abc) calibration
  double energyEmHad(double uncalibratedEnergyECAL, 
             double uncalibratedEnergyHCAL, 
             double eta, double phi) ;


  
  friend std::ostream& operator<<(std::ostream& out, 
                  const PFEnergyCalibrationHF& calib);

  const bool& getcalibHF_use() const
    { return calibHF_use_;}
  const std::vector<double>& getcalibHF_eta_step() const
    {return calibHF_eta_step_;}
  const std::vector<double>& getcalibHF_a_EMonly() const
    {return calibHF_a_EMonly_;}
  const std::vector<double>& getcalibHF_b_HADonly() const
    {return calibHF_b_HADonly_;}
  const std::vector<double>& getcalibHF_a_EMHAD() const
    {return calibHF_a_EMHAD_;}
  const std::vector<double>& getcalibHF_b_EMHAD() const
    {return calibHF_b_EMHAD_;}


 protected:

 private:
  bool calibHF_use_;
  std::vector<double>  calibHF_eta_step_;
  std::vector<double>  calibHF_a_EMonly_;
  std::vector<double>  calibHF_b_HADonly_;
  std::vector<double>  calibHF_a_EMHAD_;
  std::vector<double>  calibHF_b_EMHAD_ ;
};

#endif