#include <ShiftedMETcorrInputProducer.h>

Inheritance diagram for ShiftedMETcorrInputProducer:

Classes
struct	binningEntryType

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

	~ShiftedMETcorrInputProducer ()

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 Types
typedef std::vector < edm::InputTag >	vInputTag

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
std::vector< binningEntryType * >	binning_

std::string	moduleLabel_

double	shiftBy_

vInputTag	src_

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

Vary px, py and sumEt of "unclustered energy" (PFJets of Pt < 10 GeV plus PFCandidates not within jets) by +/- 1 standard deviation, in order to estimate resulting uncertainty on MET

Author: Christian Veelken, LLR

Version

Revision:: 1.2

Id:: ShiftedMETcorrInputProducer.h,v 1.2 2011/10/14 12:00:17 veelken Exp

Definition at line 25 of file ShiftedMETcorrInputProducer.h.

Member Typedef Documentation

typedef std::vector<edm::InputTag> ShiftedMETcorrInputProducer::vInputTag

private

Definition at line 38 of file ShiftedMETcorrInputProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 7 of file ShiftedMETcorrInputProducer.cc.

References binning_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), shiftBy_, and src_.

   : moduleLabel_(cfg.getParameter<std::string>("@module_label"))
 {
   src_ = cfg.getParameter<vInputTag>("src");
 
 //--- check that all InputTags refer to the same module label
 //   (i.e. differ by instance label only)
   for ( vInputTag::const_iterator src_ref = src_.begin();
     src_ref != src_.end(); ++src_ref ) {
     for ( vInputTag::const_iterator src_test = src_ref;
       src_test != src_.end(); ++src_test ) {
       if ( src_test->label() != src_ref->label() )
     throw cms::Exception("ShiftedMETcorrInputProducer") 
       << "InputTags specified by 'src' Configuration parameter must not refer to different module labels !!\n";
     }
   }
 
   shiftBy_ = cfg.getParameter<double>("shiftBy");
 
   if ( cfg.exists("binning") ) {
     typedef std::vector<edm::ParameterSet> vParameterSet;
     vParameterSet cfgBinning = cfg.getParameter<vParameterSet>("binning");
     for ( vParameterSet::const_iterator cfgBinningEntry = cfgBinning.begin();
       cfgBinningEntry != cfgBinning.end(); ++cfgBinningEntry ) {
       binning_.push_back(new binningEntryType(*cfgBinningEntry));
     }
   } else {
     double uncertainty = cfg.getParameter<double>("uncertainty");
     binning_.push_back(new binningEntryType(uncertainty));
   }
   
   for ( vInputTag::const_iterator src_i = src_.begin();
     src_i != src_.end(); ++src_i ) {
     for ( std::vector<binningEntryType*>::const_iterator binningEntry = binning_.begin();
       binningEntry != binning_.end(); ++binningEntry ) {
       produces<CorrMETData>((*binningEntry)->getInstanceLabel_full(src_i->instance()));
     }
   }
 }

ShiftedMETcorrInputProducer::~ShiftedMETcorrInputProducer ( )

Definition at line 47 of file ShiftedMETcorrInputProducer.cc.

References binning_.

 {
   for ( std::vector<binningEntryType*>::const_iterator it = binning_.begin();
     it != binning_.end(); ++it ) {
     delete (*it);
   }
 }

Member Function Documentation

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

privatevirtual

Implements edm::EDProducer.

Definition at line 55 of file ShiftedMETcorrInputProducer.cc.

References binning_, edm::Event::getByLabel(), edm::Event::put(), edm::shift, shiftBy_, and src_.

 {
   for ( vInputTag::const_iterator src_i = src_.begin();
     src_i != src_.end(); ++src_i ) {
     for ( std::vector<binningEntryType*>::iterator binningEntry = binning_.begin();
       binningEntry != binning_.end(); ++binningEntry ) {
       edm::Handle<CorrMETData> originalObject;
       evt.getByLabel(edm::InputTag(src_i->label(), (*binningEntry)->getInstanceLabel_full(src_i->instance())), originalObject);
   
       double shift = shiftBy_*(*binningEntry)->binUncertainty_;
       
       std::auto_ptr<CorrMETData> shiftedObject(new CorrMETData(*originalObject));
 //--- MET balances momentum of reconstructed particles,
 //    hence variations of "unclustered energy" and MET are opposite in sign
       shiftedObject->mex   = -shift*originalObject->mex;
       shiftedObject->mey   = -shift*originalObject->mey;
       shiftedObject->sumet = shift*originalObject->sumet;
       
       evt.put(shiftedObject, (*binningEntry)->getInstanceLabel_full(src_i->instance()));
     }
   }
 }