PatMuonHitFitTranslator.cc

Go to the documentation of this file.

//
//
#include "TopQuarkAnalysis/TopHitFit/interface/LeptonTranslatorBase.h"
#include "DataFormats/PatCandidates/interface/Muon.h"

namespace hitfit {

  template <>
  LeptonTranslatorBase<pat::Muon>::LeptonTranslatorBase() {
    std::string CMSSW_BASE(std::getenv("CMSSW_BASE"));
    std::string resolution_filename =
        CMSSW_BASE + std::string("/src/TopQuarkAnalysis/PatHitFit/data/exampleMuonResolution.txt");
    resolution_ = EtaDepResolution(resolution_filename);

  }  // LeptonTranslatorBase<pat::Muon>::LeptonTranslatorBase()

  template <>
  LeptonTranslatorBase<pat::Muon>::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/exampleMuonResolution.txt");
    } else {
      resolution_filename = ifile;
    }

    resolution_ = EtaDepResolution(resolution_filename);

  }  // LeptonTranslatorBase<pat::Muon>::LeptonTranslatorBase(const std::string& s)

  template <>
  LeptonTranslatorBase<pat::Muon>::~LeptonTranslatorBase() {
  }  // LeptonTranslatorBase<pat::Muon>::~LeptonTranslatorBase()

  template <>
  Lepjets_Event_Lep LeptonTranslatorBase<pat::Muon>::operator()(const pat::Muon& lepton,
                                                                int type /*= hitfit::lepton_label */,
                                                                bool useObjEmbRes /* = false */) {
    Fourvec p(lepton.px(), lepton.py(), lepton.pz(), lepton.energy());

    double muon_eta = lepton.eta();
    Vector_Resolution muon_resolution = resolution_.GetResolution(muon_eta);

    Lepjets_Event_Lep muon(p, muon_label, muon_resolution);
    return muon;

  }  // Lepjets_Event_Lep LeptonTranslatorBase<pat::Muon>::operator()

  template <>
  const EtaDepResolution& LeptonTranslatorBase<pat::Muon>::resolution() const {
    return resolution_;
  }

  template <>
  bool LeptonTranslatorBase<pat::Muon>::CheckEta(const pat::Muon& lepton) const {
    return resolution_.CheckEta(lepton.eta());
  }

}  // namespace hitfit