EgammaHLTRemoveDuplicatedSC.cc

Go to the documentation of this file.

#include "RecoEgamma/EgammaHLTProducers/interface/EgammaHLTRemoveDuplicatedSC.h"
#include "DataFormats/EgammaReco/interface/BasicCluster.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"

#include <string>

EgammaHLTRemoveDuplicatedSC::EgammaHLTRemoveDuplicatedSC(const edm::ParameterSet& ps)
{
 
 
  // the input producers
  sCInputProducer_ = consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("L1NonIsoUskimmedSC"));
  alreadyExistingSC_= consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("L1IsoSC"));

  // set the producer parameters
  outputCollection_ = ps.getParameter<std::string>("L1NonIsoSkimmedCollection");
  produces<reco::SuperClusterCollection>(outputCollection_);
 
}

EgammaHLTRemoveDuplicatedSC::~EgammaHLTRemoveDuplicatedSC()
{}

void EgammaHLTRemoveDuplicatedSC::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {

  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>(("L1NonIsoUskimmedSC"), edm::InputTag("hltCorrectedMulti5x5EndcapSuperClustersWithPreshowerL1NonIsolatedTemp"));
  desc.add<edm::InputTag>(("L1IsoSC"), edm::InputTag("hltCorrectedMulti5x5EndcapSuperClustersWithPreshowerL1Isolated"));
  desc.add<std::string>(("L1NonIsoSkimmedCollection"), ""); 
  descriptions.add(("hltEgammaHLTRemoveDuplicatedSC"), desc);  
}


void
EgammaHLTRemoveDuplicatedSC::produce(edm::Event& evt, const edm::EventSetup& es) {
 
  using namespace edm;
  
  // Get raw SuperClusters from the event    
  Handle<reco::SuperClusterCollection> UnskimmedSuperCluster;
  evt.getByToken(sCInputProducer_, UnskimmedSuperCluster);
  Handle<reco::SuperClusterCollection> L1IsoSuperCluster;
  evt.getByToken(alreadyExistingSC_, L1IsoSuperCluster);

  // Create a pointer to the existing SuperClusters
    const reco::SuperClusterCollection *UnskimmedL1NonIsoSC = UnskimmedSuperCluster.product();
    const reco::SuperClusterCollection *L1IsoSC = L1IsoSuperCluster.product();


  // Define a collection of corrected SuperClusters to put back into the event
  std::auto_ptr<reco::SuperClusterCollection> corrClusters(new reco::SuperClusterCollection);
  
  //  Loop over raw clusters and make corrected ones
  reco::SuperClusterCollection::const_iterator aClus;
  reco::SuperClusterCollection::const_iterator cit;
  for(aClus = UnskimmedL1NonIsoSC->begin(); aClus != UnskimmedL1NonIsoSC->end(); aClus++)
    {
      bool AlreadyThere = false;
      //reco::SuperCluster newClus;
      for(cit = L1IsoSC->begin(); cit != L1IsoSC->end(); cit++){
    if( fabs(aClus->energy()- cit->energy()) < 0.5 && fabs(aClus->eta()- cit->eta()) < 0.0175 ){
      float deltaphi=fabs( aClus->phi() - cit->phi() );
      if(deltaphi>6.283185308) deltaphi -= 6.283185308;
      if(deltaphi>3.141592654) deltaphi = 6.283185308-deltaphi;

      if( deltaphi < 0.035 ){AlreadyThere = true; break;}
    }
      }
      // if(AlreadyThere){std::cout<<"AAAA: SC removed: "<<aClus->energy()<<" "<<aClus->eta()<<" "<<aClus->phi()<<std::endl;}
      if(!AlreadyThere){corrClusters->push_back(*aClus);}
    }

  // Put collection of corrected SuperClusters into the event
  evt.put(corrClusters, outputCollection_);   
  
}