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