UniqueElectrons.h

Go to the documentation of this file.

#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
 
#include <vector>
 
std::vector<reco::GsfElectronRef> uniqueElectronFinder(edm::Handle<reco::GsfElectronCollection>& pElectrons) {
  const reco::GsfElectronCollection* electrons = pElectrons.product();
  //Remove duplicate electrons which share a supercluster
  std::vector<reco::GsfElectronRef> UniqueElectrons;
  int index = 0;
  for (reco::GsfElectronCollection::const_iterator elec = electrons->begin(); elec != electrons->end(); ++elec) {
    const reco::GsfElectronRef electronRef(pElectrons, index);
    reco::GsfElectronCollection::const_iterator BestDuplicate = elec;
    for (reco::GsfElectronCollection::const_iterator elec2 = electrons->begin(); elec2 != electrons->end(); ++elec2) {
      if (elec != elec2) {
        if (elec->superCluster() == elec2->superCluster()) {
          edm::LogDebug_("", "MySelection.cc", 122) << "e/p Best duplicate = " << BestDuplicate->eSuperClusterOverP()
                                                    << "\telec2 = " << elec2->eSuperClusterOverP();
          if (fabs(BestDuplicate->eSuperClusterOverP() - 1.) >= fabs(elec2->eSuperClusterOverP() - 1.)) {
            BestDuplicate = elec2;
            edm::LogDebug_("", "MySelection.cc", 122) << "elec2 is now best duplicate";
          } else
            edm::LogDebug_("", "MySelection.cc", 122) << "BestDuplicate remains best duplicate";
        }
      }
    }
    if (BestDuplicate == elec)
      UniqueElectrons.push_back(electronRef);
    ++index;
  }
  return UniqueElectrons;
}