#include <SysShiftMETcorrInputProducer.h>

Inheritance diagram for SysShiftMETcorrInputProducer:

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

	~SysShiftMETcorrInputProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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

Private Attributes
TFormula *	corrPx_

TFormula *	corrPy_

std::string	moduleLabel_

edm::EDGetTokenT< edm::View< reco::MET > >	token_

bool	useNvtx

edm::EDGetTokenT< reco::VertexCollection >	verticesToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Compute MET correction to compensate systematic shift of MET in x/y-direction (cf. https://indico.cern.ch/getFile.py/access?contribId=1&resId=0&materialId=slides&confId=174318 )

Authors: Christian Veelken, LLR

Definition at line 28 of file SysShiftMETcorrInputProducer.h.

Constructor & Destructor Documentation

◆ SysShiftMETcorrInputProducer()

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

explicit

Definition at line 11 of file SysShiftMETcorrInputProducer.cc.

     : moduleLabel_(cfg.getParameter<std::string>("@module_label")), useNvtx(false), corrPx_(nullptr), corrPy_(nullptr) {
   token_ = consumes<edm::View<reco::MET> >(cfg.getParameter<edm::InputTag>("src"));
  
   edm::ParameterSet cfgCorrParameter = cfg.getParameter<edm::ParameterSet>("parameter");
   TString corrPxFormula = cfgCorrParameter.getParameter<std::string>("px");
   TString corrPyFormula = cfgCorrParameter.getParameter<std::string>("py").data();
   if (corrPxFormula.Contains("Nvtx") || corrPyFormula.Contains("Nvtx")) {
     useNvtx = true, verticesToken_ = consumes<reco::VertexCollection>(cfg.getParameter<edm::InputTag>("srcVertices"));
   }
  
   corrPxFormula.ReplaceAll("sumEt", "x");
   corrPxFormula.ReplaceAll("Nvtx", "y");
   std::string corrPxName = std::string(moduleLabel_).append("_corrPx");
   corrPx_ = new TFormula(corrPxName.data(), corrPxFormula.Data());
  
   corrPyFormula.ReplaceAll("sumEt", "x");
   corrPyFormula.ReplaceAll("Nvtx", "y");
   std::string corrPyName = std::string(moduleLabel_).append("_corrPy");
   corrPy_ = new TFormula(corrPyName.data(), corrPyFormula.Data());
  
   produces<CorrMETData>();
 }

References looper::cfg, corrPx_, corrPy_, data, edm::ParameterSet::getParameter(), moduleLabel_, AlCaHLTBitMon_QueryRunRegistry::string, token_, useNvtx, and verticesToken_.

◆ ~SysShiftMETcorrInputProducer()

SysShiftMETcorrInputProducer::~SysShiftMETcorrInputProducer ( )

override

Definition at line 35 of file SysShiftMETcorrInputProducer.cc.

                                                             {
   // nothing to be done yet...
 }

Member Function Documentation

◆ produce()

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

overrideprivate

Definition at line 39 of file SysShiftMETcorrInputProducer.cc.

                                                                                  {
   //std::cout << "<SysShiftMETcorrInputProducer::produce>:" << std::endl;
  
   typedef edm::View<reco::MET> METView;
   edm::Handle<METView> met;
   evt.getByToken(token_, met);
   if (met->size() != 1)
     throw cms::Exception("SysShiftMETcorrInputProducer::produce") << "Failed to find unique MET object !!\n";
  
   double sumEt = met->front().sumEt();
   //std::cout << " sumEt = " << sumEt << std::endl;
  
   size_t Nvtx = 0;
   if (useNvtx) {
     edm::Handle<reco::VertexCollection> vertices;
     evt.getByToken(verticesToken_, vertices);
     Nvtx = vertices->size();
   }
   //std::cout << " Nvtx = " << Nvtx << std::endl;
  
   std::unique_ptr<CorrMETData> metCorr(new CorrMETData());
   metCorr->mex = -corrPx_->Eval(sumEt, Nvtx);
   metCorr->mey = -corrPy_->Eval(sumEt, Nvtx);
   //std::cout << "--> metCorr: Px = " << metCorr->mex << ", Py = " << metCorr->mey << std::endl;
  
   evt.put(std::move(metCorr));
 }