PatElectronHitFitTranslator.cc

Go to the documentation of this file.

 //
// $Id: PatElectronHitFitTranslator.cc,v 1.1 2011/05/26 09:47:00 mseidel Exp $
//

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

namespace hitfit {

template<>
LeptonTranslatorBase<pat::Electron>::LeptonTranslatorBase()
{

    std::string CMSSW_BASE(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(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