L1AnalysisRecoElectron.cc

Go to the documentation of this file.

#include "L1Trigger/L1TNtuples/interface/L1AnalysisRecoElectron.h"

using namespace std;

L1Analysis::L1AnalysisRecoElectron::L1AnalysisRecoElectron() {}

L1Analysis::L1AnalysisRecoElectron::~L1AnalysisRecoElectron() {}

void L1Analysis::L1AnalysisRecoElectron::SetElectron(
    const edm::Event& event,
    const edm::EventSetup& setup,
    //const edm::Handle<edm::View<reco::GsfElectron>>& electrons,
    edm::Handle<reco::GsfElectronCollection> electrons,
    std::vector<edm::Handle<edm::ValueMap<bool> > > eleVIDDecisionHandles,
    const unsigned& maxElectron)

{
  recoElectron_.nElectrons = 0;

  for (reco::GsfElectronCollection::const_iterator el = electrons->begin();
       el != electrons->end() && recoElectron_.nElectrons < maxElectron;
       ++el) {
    recoElectron_.e.push_back(el->energy());
    recoElectron_.pt.push_back(el->pt());
    recoElectron_.et.push_back(el->et());
    recoElectron_.eta.push_back(el->eta());
    recoElectron_.phi.push_back(el->phi());
    recoElectron_.eta_SC.push_back((el->superClusterPosition()).eta());
    recoElectron_.phi_SC.push_back((el->superClusterPosition()).phi());
    recoElectron_.e_ECAL.push_back(el->ecalEnergy());
    recoElectron_.e_SC.push_back(el->superCluster()->energy());
    recoElectron_.charge.push_back(el->charge());

    edm::Ref<reco::GsfElectronCollection> electronEdmRef(electrons, recoElectron_.nElectrons);

    recoElectron_.isVetoElectron.push_back((*(eleVIDDecisionHandles[0]))[electronEdmRef]);
    recoElectron_.isLooseElectron.push_back((*(eleVIDDecisionHandles[1]))[electronEdmRef]);
    recoElectron_.isMediumElectron.push_back((*(eleVIDDecisionHandles[2]))[electronEdmRef]);
    recoElectron_.isTightElectron.push_back((*(eleVIDDecisionHandles[3]))[electronEdmRef]);

    double iso = (el->pfIsolationVariables().sumChargedHadronPt +
                  max(el->pfIsolationVariables().sumNeutralHadronEt + el->pfIsolationVariables().sumPhotonEt -
                          0.5 * el->pfIsolationVariables().sumPUPt,
                      0.0)) /
                 el->pt();

    recoElectron_.iso.push_back(iso);
    // cout<<ConversionTools::hasMatchedConversion(*el,conversions,theBeamSpot->position())<<endl;

    recoElectron_.nElectrons++;
  }
  // for(reco::GsfElectronCollection::const_iterator it=electrons->begin();
  //     it!=electrons->end() && recoElectron_.nElectrons < maxElectron;
  //     ++it) {

  //   recoElectron_.e.push_back(it->energy());
  //   recoElectron_.pt.push_back(it->pt());
  //   recoElectron_.et.push_back(it->et());
  //   recoElectron_.eta.push_back(it->eta());
  //   recoElectron_.phi.push_back(it->phi());

  //   // cout<<it->superCluster().position().eta()<<endl;
  //   isVetoElectronCustom(*it, vertices, conversions, theBeamSpot, Rho);

  //   // cout<<ConversionTools::hasMatchedConversion(*it,conversions,theBeamSpot->position())<<endl;

  //   recoElectron_.nElectrons++;

  // }
}