HiggsTo2GammaSkim.cc

Go to the documentation of this file.

 
/* \class HiggsTo2GammaSkim
 *
 * Consult header file for description
 *
 * author:  Kati Lassila-Perini Helsinki Institute of Physics
 *
 */
 
// system include files
#include <HiggsAnalysis/Skimming/interface/HiggsTo2GammaSkim.h>
 
// User include files
#include <FWCore/ParameterSet/interface/ParameterSet.h>
 
// Message logger
#include <FWCore/MessageLogger/interface/MessageLogger.h>
 
// Photons:
#include <DataFormats/EgammaCandidates/interface/Photon.h>
 
// C++
#include <iostream>
#include <vector>
 
using namespace std;
using namespace edm;
using namespace reco;
 
// Constructor
HiggsTo2GammaSkim::HiggsTo2GammaSkim(const edm::ParameterSet& pset) {
  // Local Debug flag
  debug = pset.getParameter<bool>("DebugHiggsTo2GammaSkim");
 
  // Reconstructed objects
  thePhotonToken = consumes<reco::PhotonCollection>(pset.getParameter<edm::InputTag>("PhotonCollectionLabel"));
 
  // Minimum Pt for photons for skimming
  photon1MinPt = pset.getParameter<double>("photon1MinimumPt");
  nPhotonMin = pset.getParameter<int>("nPhotonMinimum");
 
  nEvents = 0;
  nSelectedEvents = 0;
}
 
// Destructor
HiggsTo2GammaSkim::~HiggsTo2GammaSkim() {
  edm::LogVerbatim("HiggsTo2GammaSkim") << " Number_events_read " << nEvents << " Number_events_kept "
                                        << nSelectedEvents << " Efficiency         "
                                        << ((double)nSelectedEvents) / ((double)nEvents + 0.01) << std::endl;
}
 
// Filter event
bool HiggsTo2GammaSkim::filter(edm::Event& event, const edm::EventSetup& setup) {
  nEvents++;
 
  using reco::PhotonCollection;
 
  bool keepEvent = false;
  int nPhotons = 0;
 
  // Look at photons:
 
  // Get the photon collection from the event
  edm::Handle<reco::PhotonCollection> photonHandle;
 
  event.getByToken(thePhotonToken, photonHandle);
 
  if (photonHandle.isValid()) {
    const reco::PhotonCollection* phoCollection = photonHandle.product();
 
    reco::PhotonCollection::const_iterator photons;
 
    // Loop over photon collections and count how many photons there are,
    // and how many are above the thresholds
 
    // Question: do we need to take the reconstructed primary vertex at this point?
    // Here, I assume that the et is taken with respect to the nominal vertex (0,0,0).
    for (photons = phoCollection->begin(); photons != phoCollection->end(); ++photons) {
      float et_p = photons->et();
      if (et_p > photon1MinPt)
        nPhotons++;
    }
  }
 
  // Make decision:
  if (nPhotons >= nPhotonMin)
    keepEvent = true;
 
  if (keepEvent)
    nSelectedEvents++;
 
  return keepEvent;
}