PatElectronHitFitTranslator.cc

Go to the documentation of this file.

//
//
 
#include "TopQuarkAnalysis/TopHitFit/interface/LeptonTranslatorBase.h"
#include "DataFormats/PatCandidates/interface/Electron.h"
 
namespace hitfit {
 
  template <>
  LeptonTranslatorBase<pat::Electron>::LeptonTranslatorBase() {
    std::string CMSSW_BASE(std::getenv("CMSSW_BASE"));
    std::string resolution_filename =
        CMSSW_BASE + std::string("/src/TopQuarkAnalysis/PatHitFit/data/exampleElectronResolution.txt");
    resolution_ = EtaDepResolution(resolution_filename);
 
  }  // LeptonTranslatorBase<pat::Electron>::LeptonTranslatorBase()
 
  template <>
  LeptonTranslatorBase<pat::Electron>::LeptonTranslatorBase(const std::string& ifile) {
    std::string CMSSW_BASE(std::getenv("CMSSW_BASE"));
    std::string resolution_filename;
 
    if (ifile.empty()) {
      resolution_filename =
          CMSSW_BASE + std::string("/src/TopQuarkAnalysis/PatHitFit/data/exampleElectronResolution.txt");
    } else {
      resolution_filename = ifile;
    }
 
    resolution_ = EtaDepResolution(resolution_filename);
 
  }  // LeptonTranslatorBase<pat::Electron>::LeptonTranslatorBase(const std::string& ifile)
 
  template <>
  LeptonTranslatorBase<pat::Electron>::~LeptonTranslatorBase() {}
 
  template <>
  Lepjets_Event_Lep LeptonTranslatorBase<pat::Electron>::operator()(const pat::Electron& lepton,
                                                                    int type /* = hitfit::lepton_label */,
                                                                    bool useObjEmbRes /* = false */) {
    Fourvec p(lepton.px(), lepton.py(), lepton.pz(), lepton.energy());
 
    double electron_eta = lepton.superCluster()->eta();
    Vector_Resolution electron_resolution = resolution_.GetResolution(electron_eta);
 
    Lepjets_Event_Lep electron(p, electron_label, electron_resolution);
    return electron;
 
  }  // Lepjets_Event_Lep LeptonTranslatorBase<pat::Electron>::operator()
 
  template <>
  const EtaDepResolution& LeptonTranslatorBase<pat::Electron>::resolution() const {
    return resolution_;
  }
 
  template <>
  bool LeptonTranslatorBase<pat::Electron>::CheckEta(const pat::Electron& lepton) const {
    double electron_eta = lepton.superCluster()->eta();
    return resolution_.CheckEta(electron_eta);
  }
 
}  // namespace hitfit