#include <GenFilterEfficiencyProducer.h>

Inheritance diagram for GenFilterEfficiencyProducer:

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

	~GenFilterEfficiencyProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

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 Member Functions
virtual void	beginLuminosityBlock (edm::LuminosityBlock &, const edm::EventSetup &)

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

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

Private Attributes
std::string	filterPath

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_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
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 36 of file GenFilterEfficiencyProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 5 of file GenFilterEfficiencyProducer.cc.

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

                                                                                        :
   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";
   }
 
   produces<GenFilterInfo, edm::InLumi>(); 
 
 }

GenFilterEfficiencyProducer::~GenFilterEfficiencyProducer ( )

Definition at line 38 of file GenFilterEfficiencyProducer.cc.

 {
 
 } 

Member Function Documentation

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

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 101 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 &	iLumi,
		const edm::EventSetup &
	)

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 115 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
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 50 of file GenFilterEfficiencyProducer.cc.

References edm::Event::getByLabel(), numEventsPassNeg_, numEventsPassPos_, numEventsTotalNeg_, numEventsTotalPos_, pathIndex, sumpass_w2_, sumpass_w_, sumtotal_w2_, sumtotal_w_, thisProcess, and CommonMethods::weight().

 {
 
   edm::InputTag theTrig("TriggerResults","",thisProcess);
   edm::Handle<edm::TriggerResults> trigR;
   iEvent.getByLabel(theTrig,trigR); 
   edm::Handle<GenEventInfoProduct>    genEventScale;
   if (!iEvent.getByLabel("generator", genEventScale))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;
   }
 
 }