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 ()
ModuleDescription const &	moduleDescription () const
virtual	~EDFilter ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Attributes
bool	debug
int	nEvents
int	nPhotonMin
int	nSelectedEvents
float	photon1MinPt
float	photon2MinPt
edm::EDGetTokenT < reco::PhotonCollection >	thePhotonToken

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 27 of file HiggsTo2GammaSkim.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 33 of file HiggsTo2GammaSkim.cc.

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

                                                                {

  // 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;

}

HiggsTo2GammaSkim::~HiggsTo2GammaSkim

(

)

Definition at line 53 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 64 of file HiggsTo2GammaSkim.cc.

References edm::HandleBase::isValid(), nEvents, interactiveExample::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.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;
}

Member Data Documentation

bool HiggsTo2GammaSkim::debug

private

Definition at line 44 of file HiggsTo2GammaSkim.h.

Referenced by rrapi.RRApi::dprint(), and rrapi.RRApi::get().

int HiggsTo2GammaSkim::nEvents

private

Definition at line 41 of file HiggsTo2GammaSkim.h.

Referenced by looper.Looper::loop().

int HiggsTo2GammaSkim::nPhotonMin

private

Definition at line 47 of file HiggsTo2GammaSkim.h.

int HiggsTo2GammaSkim::nSelectedEvents

private

Definition at line 41 of file HiggsTo2GammaSkim.h.

float HiggsTo2GammaSkim::photon1MinPt

private

Definition at line 45 of file HiggsTo2GammaSkim.h.

float HiggsTo2GammaSkim::photon2MinPt

private

Definition at line 46 of file HiggsTo2GammaSkim.h.

edm::EDGetTokenT<reco::PhotonCollection> HiggsTo2GammaSkim::thePhotonToken

private

Definition at line 50 of file HiggsTo2GammaSkim.h.