#include <EmbeddingKineReweightProducer.h>

Inheritance diagram for EmbeddingKineReweightProducer:

Classes
struct	lutEntryType

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

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

	~EmbeddingKineReweightProducer ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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 Attributes
std::vector< lutEntryType * >	lutEntries_

double	maxWeight_

double	minWeight_

edm::InputTag	srcGenDiTaus_

int	verbosity_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	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::EDProducer
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

Compute reweighting factors compensating (small) biases of the embedding procedure

Authors: Christian Veelken

Version

Revision:: 1.3

Id:: EmbeddingKineReweightProducer.h,v 1.3 2013/04/17 07:11:37 veelken Exp

Definition at line 35 of file EmbeddingKineReweightProducer.h.

Constructor & Destructor Documentation

EmbeddingKineReweightProducer::EmbeddingKineReweightProducer ( const edm::ParameterSet & cfg )

explicit

Definition at line 15 of file EmbeddingKineReweightProducer.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNamesForType(), analyzePatCleaning_cfg::inputFile, writeAntiElectronDiscrMVA_cfi::inputFileName, lutEntries_, maxWeight_, minWeight_, srcGenDiTaus_, AlCaHLTBitMon_QueryRunRegistry::string, edm::FileInPath::Unknown, makeLayoutFileForGui::variables, and verbosity_.

 {
   srcGenDiTaus_ = cfg.getParameter<edm::InputTag>("srcGenDiTaus"); 
 
   edm::FileInPath inputFileName = cfg.getParameter<edm::FileInPath>("inputFileName");
   if ( inputFileName.location() == edm::FileInPath::Unknown) 
     throw cms::Exception("EmbeddingReweightProducer") 
       << " Failed to find File = " << inputFileName << " !!\n";
   std::auto_ptr<TFile> inputFile(new TFile(inputFileName.fullPath().data()));
 
   edm::ParameterSet cfgLUTs = cfg.getParameter<edm::ParameterSet>("lutNames"); 
   std::vector<std::string> variables = cfgLUTs.getParameterNamesForType<std::string>();
   for ( std::vector<std::string>::const_iterator variable = variables.begin(); 
         variable != variables.end(); ++variable ) {
     std::string lutName = cfgLUTs.getParameter<std::string>(*variable);
     lutEntryType* lutEntry = 
       new lutEntryType(*inputFile, *variable, lutName);
     lutEntries_.push_back(lutEntry);
   }
   
   minWeight_ = cfg.getParameter<double>("minWeight"); 
   maxWeight_ = cfg.getParameter<double>("maxWeight");
 
   verbosity_ = ( cfg.exists("verbosity") ) ?
     cfg.getParameter<int>("verbosity") : 0;
 
   for ( std::vector<lutEntryType*>::iterator lutEntry = lutEntries_.begin();
     lutEntry != lutEntries_.end(); ++lutEntry ) {
     produces<double>((*lutEntry)->variableName_);
   }
 }

EmbeddingKineReweightProducer::~EmbeddingKineReweightProducer ( )

Definition at line 47 of file EmbeddingKineReweightProducer.cc.

References lutEntries_.

 {
   for ( std::vector<lutEntryType*>::iterator it = lutEntries_.begin();
     it != lutEntries_.end(); ++it ) {
     delete (*it);
   }
 }

Member Function Documentation

void EmbeddingKineReweightProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

virtual

Implements edm::EDProducer.

Definition at line 55 of file EmbeddingKineReweightProducer.cc.

References gather_cfg::cout, reco::Candidate::eta(), edm::Event::getByLabel(), lutEntries_, reco::Candidate::mass(), maxWeight_, minWeight_, reco::Candidate::phi(), reco::Candidate::pt(), edm::Event::put(), srcGenDiTaus_, verbosity_, and histoStyle::weight.

 {
   if ( verbosity_ ) {
     std::cout << "<EmbeddingKineReweightProducer::produce>:" << std::endl;
   }
 
   typedef edm::View<reco::Candidate> CandidateView;
   edm::Handle<CandidateView> genDiTaus;
   evt.getByLabel(srcGenDiTaus_, genDiTaus);
   if ( genDiTaus->size() != 1 )
     throw cms::Exception("EmbeddingKineReweightProducer") 
       << "Failed to find unique genDiTau object !!\n";  
   const reco::Candidate& genDiTau = genDiTaus->front();
   if ( verbosity_ ) {
     std::cout << "diTau: Pt = " << genDiTau.pt() << ", eta = " << genDiTau.eta() << ", phi = " << genDiTau.phi() << ", mass = " << genDiTau.mass() << std::endl;
   }
 
   for ( std::vector<lutEntryType*>::const_iterator lutEntry = lutEntries_.begin();
     lutEntry != lutEntries_.end(); ++lutEntry ) {
     double weight = (**lutEntry)(genDiTau);
     if ( weight < minWeight_ ) weight = minWeight_;
     if ( weight > maxWeight_ ) weight = maxWeight_;
     if ( verbosity_ ) {
       std::cout << " " << (*lutEntry)->variableName_ << " = " << weight << std::endl;
     }
     std::auto_ptr<double> weightPtr(new double(weight));
     evt.put(weightPtr, (*lutEntry)->variableName_);
   }
 }