#include <GenFilterEfficiencyProducer.h>

Inheritance diagram for GenFilterEfficiencyProducer:

Public Member Functions
	GenFilterEfficiencyProducer (const edm::ParameterSet &)

	~GenFilterEfficiencyProducer ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks >
	EDProducer ()=default

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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 (const std::string &iProcessName, std::vector< const char * > &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 Member Functions
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, const edm::EventSetup &) override

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, const edm::EventSetup &) override

virtual void	endLuminosityBlockProduce (edm::LuminosityBlock &, const edm::EventSetup &) override

virtual void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::string	filterPath

edm::EDGetTokenT < GenEventInfoProduct >	genEventInfoToken_

unsigned int	numEventsPassNeg_

unsigned int	numEventsPassPos_

unsigned int	numEventsTotalNeg_

unsigned int	numEventsTotalPos_

unsigned int	pathIndex

double	sumpass_w2_

double	sumpass_w_

double	sumtotal_w2_

double	sumtotal_w_

std::string	thisProcess

edm::service::TriggerNamesService *	tns_

edm::EDGetTokenT < edm::TriggerResults >	triggerResultsToken_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	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::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 36 of file GenFilterEfficiencyProducer.h.

Constructor & Destructor Documentation

GenFilterEfficiencyProducer::GenFilterEfficiencyProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 4 of file GenFilterEfficiencyProducer.cc.

References filterPath, genEventInfoToken_, edm::service::TriggerNamesService::getProcessName(), edm::service::TriggerNamesService::getTrigPaths(), i, HLT_25ns14e33_v1_cff::InputTag, cppFunctionSkipper::operator, pathIndex, thisProcess, tns_, and triggerResultsToken_.

                                                                                        :
   filterPath(iConfig.getParameter<std::string>("filterPath")),
   tns_(),
   thisProcess(),pathIndex(100000),
   numEventsPassPos_(0),
   numEventsPassNeg_(0),
   numEventsTotalPos_(0),
   numEventsTotalNeg_(0),
   sumpass_w_(0.),
   sumpass_w2_(0.),
   sumtotal_w_(0.),
   sumtotal_w2_(0.)
 {
    //now do what ever initialization is needed
   if (edm::Service<edm::service::TriggerNamesService>().isAvailable()) {
     // get tns pointer
     tns_ = edm::Service<edm::service::TriggerNamesService>().operator->();
     if (tns_!=0) {
       thisProcess = tns_->getProcessName();
       std::vector<std::string> theNames = tns_->getTrigPaths();
       for ( unsigned int i = 0; i < theNames.size(); i++ ) {
         if ( theNames[i] == filterPath ) { pathIndex = i; continue; }
       }
     }
     else
       edm::LogError("ServiceNotAvailable") << "TriggerNamesServive not available, no filter information stored";
   }
 
   triggerResultsToken_ = consumes<edm::TriggerResults>(edm::InputTag("TriggerResults","",thisProcess));
   genEventInfoToken_ = consumes<GenEventInfoProduct>(edm::InputTag("generator",""));
   produces<GenFilterInfo, edm::InLumi>(); 
 
   
 }

GenFilterEfficiencyProducer::~GenFilterEfficiencyProducer ( )

Definition at line 40 of file GenFilterEfficiencyProducer.cc.

 {
 
 } 

Member Function Documentation

void GenFilterEfficiencyProducer::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		const edm::EventSetup &
	)

overrideprivatevirtual

Definition at line 103 of file GenFilterEfficiencyProducer.cc.

References numEventsPassNeg_, numEventsPassPos_, numEventsTotalNeg_, numEventsTotalPos_, sumpass_w2_, sumpass_w_, sumtotal_w2_, and sumtotal_w_.

                                                                                                  {
 
   numEventsPassPos_=0;
   numEventsPassNeg_=0;
   numEventsTotalPos_=0;
   numEventsTotalNeg_=0;
   sumpass_w_=0;
   sumpass_w2_=0;
   sumtotal_w_=0;
   sumtotal_w2_=0;
  
 }

void GenFilterEfficiencyProducer::endLuminosityBlock	(	edm::LuminosityBlock const &	iLumi,
		const edm::EventSetup &
	)

overrideprivatevirtual

Definition at line 116 of file GenFilterEfficiencyProducer.cc.

116 {

117 }

void GenFilterEfficiencyProducer::endLuminosityBlockProduce	(	edm::LuminosityBlock &	iLumi,
		const edm::EventSetup &
	)

overrideprivatevirtual

Definition at line 120 of file GenFilterEfficiencyProducer.cc.

References numEventsPassNeg_, numEventsPassPos_, numEventsTotalNeg_, numEventsTotalPos_, edm::LuminosityBlock::put(), sumpass_w2_, sumpass_w_, sumtotal_w2_, and sumtotal_w_.

                                                                                                        {
 
   std::auto_ptr<GenFilterInfo> thisProduct(new GenFilterInfo(
                                  numEventsPassPos_,
                                  numEventsPassNeg_,
                                  numEventsTotalPos_,
                                  numEventsTotalNeg_,
                                  sumpass_w_,
                                  sumpass_w2_,
                                  sumtotal_w_,
                                  sumtotal_w2_
                                  ));
   iLumi.put(thisProduct);
 }

void GenFilterEfficiencyProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDProducerBase.

Definition at line 52 of file GenFilterEfficiencyProducer.cc.

References genEventInfoToken_, edm::Event::getByToken(), edm::HandleBase::isValid(), numEventsPassNeg_, numEventsPassPos_, numEventsTotalNeg_, numEventsTotalPos_, pathIndex, sumpass_w2_, sumpass_w_, sumtotal_w2_, sumtotal_w_, triggerResultsToken_, and puppiForMET_cff::weight.

 {
 
   edm::Handle<edm::TriggerResults> trigR;
   iEvent.getByToken(triggerResultsToken_,trigR); 
   edm::Handle<GenEventInfoProduct>    genEventScale;
   iEvent.getByToken(genEventInfoToken_,genEventScale);
   if (!genEventScale.isValid()) return;
   double weight = genEventScale->weight();
   
   unsigned int nSize = (*trigR).size();
   // std::cout << "Number of paths in TriggerResults = " << nSize  << std::endl;
   if ( nSize >= pathIndex ) {
 
     if (!trigR->wasrun(pathIndex))return;
     if ( trigR->accept(pathIndex) ) { 
       sumpass_w_ += weight;
       sumpass_w2_+= weight*weight;
       
       sumtotal_w_ += weight;
       sumtotal_w2_+= weight*weight;
 
       if(weight > 0)
     {
       numEventsPassPos_++;
       numEventsTotalPos_++;
     }
       else
     {
       numEventsPassNeg_++;
       numEventsTotalNeg_++;
     }
 
     }
     else // if fail the filter
       {
     sumtotal_w_ += weight;
     sumtotal_w2_+= weight*weight;
 
     if(weight > 0)
       numEventsTotalPos_++;
     else
       numEventsTotalNeg_++;
       }
     //    std::cout << "Total events = " << numEventsTotal << " passed = " << numEventsPassed << std::endl;
 
   }
   
 }