SimpleElectron.cc

Go to the documentation of this file.

#ifndef SimpleElectron_STANDALONE
#include "EgammaAnalysis/ElectronTools/interface/SimpleElectron.h"
SimpleElectron::SimpleElectron(const reco::GsfElectron &in, unsigned int runNumber, bool isMC) :
        run_(runNumber),
        eClass_(int(in.classification())), 
        r9_(in.full5x5_r9()),
        scEnergy_(in.superCluster()->rawEnergy() + in.isEB() ? 0 : in.superCluster()->preshowerEnergy()), 
        scEnergyError_(-999.),  // FIXME???
        trackMomentum_(in.trackMomentumAtVtx().R()), 
        trackMomentumError_(in.trackMomentumError()), 
        regEnergy_(in.correctedEcalEnergy()), 
        regEnergyError_(in.correctedEcalEnergyError()), 
        eta_(in.superCluster()->eta()), 
        isEB_(in.isEB()), 
        isMC_(isMC), 
        isEcalDriven_(in.ecalDriven()), 
        isTrackerDriven_(in.trackerDrivenSeed()), 
        newEnergy_(regEnergy_), 
        newEnergyError_(regEnergyError_),
        combinedMomentum_(in.p4(reco::GsfElectron::P4_COMBINATION).P()), 
        combinedMomentumError_(in.p4Error(reco::GsfElectron::P4_COMBINATION)),
        scale_(1.0), smearing_(0.0)
{
}


void SimpleElectron::writeTo(reco::GsfElectron & out) const 
{
    math::XYZTLorentzVector oldMomentum = out.p4();
    math::XYZTLorentzVector newMomentum = math::XYZTLorentzVector(oldMomentum.x()*getCombinedMomentum()/oldMomentum.t(),
                                                                  oldMomentum.y()*getCombinedMomentum()/oldMomentum.t(),
                                                                  oldMomentum.z()*getCombinedMomentum()/oldMomentum.t(),
                                                                  getCombinedMomentum());
    out.setCorrectedEcalEnergy(getNewEnergy());
    out.setCorrectedEcalEnergyError(getNewEnergyError());
    out.correctMomentum(newMomentum, getTrackerMomentumError(), getCombinedMomentumError());
}
#endif