HiggsTo2GammaSkim Class Reference

#include <HiggsTo2GammaSkim.h>

Inheritance diagram for HiggsTo2GammaSkim:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)
	Get event properties to send to builder to fill seed collection. More...
	HiggsTo2GammaSkim (const edm::ParameterSet &)
	~HiggsTo2GammaSkim ()
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
virtual	~EDFilter ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

Private Attributes
bool	debug
int	nEvents
int	nPhotonMin
int	nSelectedEvents
float	photon1MinPt
float	photon2MinPt
edm::InputTag	thePhotonLabel

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType
typedef WorkerT< EDFilter >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 26 of file HiggsTo2GammaSkim.h.

Constructor & Destructor Documentation

HiggsTo2GammaSkim::HiggsTo2GammaSkim ( const edm::ParameterSet &  pset )

explicit

Definition at line 34 of file HiggsTo2GammaSkim.cc.

References debug, edm::ParameterSet::getParameter(), and nEvents.

                                                                {

  // Local Debug flag
  debug              = pset.getParameter<bool>("DebugHiggsTo2GammaSkim");

  // Reconstructed objects
  thePhotonLabel     = 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;

}

HiggsTo2GammaSkim::~HiggsTo2GammaSkim

(

)

Definition at line 54 of file HiggsTo2GammaSkim.cc.

References nEvents.

                                      {

  edm::LogVerbatim("HiggsTo2GammaSkim") 
  << " Number_events_read " << nEvents          
  << " Number_events_kept " << nSelectedEvents 
  << " Efficiency         " << ((double)nSelectedEvents)/((double) nEvents + 0.01) << std::endl;
}

Member Function Documentation

bool HiggsTo2GammaSkim::filter ( edm::Event &  event, const edm::EventSetup &  setup  )

virtual

Get event properties to send to builder to fill seed collection.

Implements edm::EDFilter.

Definition at line 65 of file HiggsTo2GammaSkim.cc.

References edm::HandleBase::isValid(), nEvents, ExpressReco_HICollisions_FallBack::photons, and edm::Handle< T >::product().

                                                                           {

  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.getByLabel(thePhotonLabel.label(),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;
}

Member Data Documentation

bool HiggsTo2GammaSkim::debug

private

Definition at line 43 of file HiggsTo2GammaSkim.h.

int HiggsTo2GammaSkim::nEvents

private

Definition at line 40 of file HiggsTo2GammaSkim.h.

int HiggsTo2GammaSkim::nPhotonMin

private

Definition at line 46 of file HiggsTo2GammaSkim.h.

int HiggsTo2GammaSkim::nSelectedEvents

private

Definition at line 40 of file HiggsTo2GammaSkim.h.

float HiggsTo2GammaSkim::photon1MinPt

private

Definition at line 44 of file HiggsTo2GammaSkim.h.

float HiggsTo2GammaSkim::photon2MinPt

private

Definition at line 45 of file HiggsTo2GammaSkim.h.

edm::InputTag HiggsTo2GammaSkim::thePhotonLabel

private

Definition at line 49 of file HiggsTo2GammaSkim.h.