#include <MultShiftMETcorrDBInputProducer.h>

Inheritance diagram for MultShiftMETcorrDBInputProducer:

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

	~MultShiftMETcorrDBInputProducer () 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

Static Private Member Functions
static int	translateTypeToAbsPdgId (reco::PFCandidate::ParticleType type)

Private Attributes
std::vector< edm::ParameterSet >	cfgCorrParameters_

int	counts_

std::vector< double >	etaMax_

std::vector< double >	etaMin_

std::unique_ptr< TF1 >	formula_x_

std::unique_ptr< TF1 >	formula_y_

bool	mIsData

edm::ESGetToken< MEtXYcorrectParametersCollection, MEtXYcorrectRecord >	mMEtXYcorParaColl_

std::string	moduleLabel_

std::string	mPayloadName

std::string	mSampleType

edm::EDGetTokenT< edm::View< reco::Candidate > >	pflow_

double	sumPt_

edm::EDGetTokenT< edm::View< reco::Vertex > >	vertices_

edm::EDGetTokenT< edm::ValueMap< float > >	weightsToken_

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: SangEun Lee,; Robert Schoefbeck, Vienna

Definition at line 36 of file MultShiftMETcorrDBInputProducer.h.

Constructor & Destructor Documentation

◆ MultShiftMETcorrDBInputProducer()

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

explicit

Definition at line 38 of file MultShiftMETcorrDBInputProducer.cc.

     : moduleLabel_(cfg.getParameter<std::string>("@module_label")) {
   mPayloadName = cfg.getParameter<std::string>("payloadName");
   mSampleType = (cfg.exists("sampleType")) ? cfg.getParameter<std::string>("sampleType") : "MC";
   mIsData = cfg.getParameter<bool>("isData");
  
   pflow_ = consumes<edm::View<reco::Candidate>>(cfg.getParameter<edm::InputTag>("srcPFlow"));
   vertices_ = consumes<edm::View<reco::Vertex>>(cfg.getParameter<edm::InputTag>("vertexCollection"));
   mMEtXYcorParaColl_ =
       esConsumes<MEtXYcorrectParametersCollection, MEtXYcorrectRecord>(edm::ESInputTag("", mPayloadName));
  
   edm::InputTag srcWeights = cfg.getParameter<edm::InputTag>("srcWeights");
   if (!srcWeights.label().empty())
     weightsToken_ = consumes<edm::ValueMap<float>>(srcWeights);
  
   etaMin_.clear();
   etaMax_.clear();
  
   produces<CorrMETData>();
 }

References looper::cfg, etaMax_, etaMin_, mIsData, mMEtXYcorParaColl_, mPayloadName, mSampleType, pflow_, HLT_FULL_cff::srcWeights, AlCaHLTBitMon_QueryRunRegistry::string, vertices_, and weightsToken_.

◆ ~MultShiftMETcorrDBInputProducer()

MultShiftMETcorrDBInputProducer::~MultShiftMETcorrDBInputProducer ( )

override

Definition at line 59 of file MultShiftMETcorrDBInputProducer.cc.

59 {}

Member Function Documentation

◆ produce()

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

overrideprivate

Definition at line 61 of file MultShiftMETcorrDBInputProducer.cc.

                                                                                     {
   // Get para.s from DB
   edm::ESHandle<MEtXYcorrectParametersCollection> MEtXYcorParaColl = es.getHandle(mMEtXYcorParaColl_);
  
   // get the sections from Collection (pair of section and METCorr.Par class)
   std::vector<MEtXYcorrectParametersCollection::key_type> keys;
   // save level to keys for each METParameter in METParameter collection
   MEtXYcorParaColl->validKeys(keys);
  
   //get primary vertices
   edm::Handle<edm::View<reco::Vertex>> hpv;
   evt.getByToken(vertices_, hpv);
   if (!hpv.isValid()) {
     edm::LogError("MultShiftMETcorrDBInputProducer::produce") << "could not find vertex collection ";
   }
   std::vector<reco::Vertex> goodVertices;
   for (unsigned i = 0; i < hpv->size(); i++) {
     if ((*hpv)[i].ndof() > 4 && (fabs((*hpv)[i].z()) <= 24.) && (fabs((*hpv)[i].position().rho()) <= 2.0))
       goodVertices.push_back((*hpv)[i]);
   }
   int ngoodVertices = goodVertices.size();
  
   edm::Handle<edm::View<reco::Candidate>> particleFlow;
   evt.getByToken(pflow_, particleFlow);
  
   edm::Handle<edm::ValueMap<float>> weights;
   if (!weightsToken_.isUninitialized())
     evt.getByToken(weightsToken_, weights);
  
   //loop over all constituent types and sum each correction
   //std::unique_ptr<CorrMETData> metCorr(new CorrMETData());
   std::unique_ptr<CorrMETData> metCorr(new CorrMETData());
  
   double corx = 0.;
   double cory = 0.;
  
   // check DB
   for (std::vector<MEtXYcorrectParametersCollection::key_type>::const_iterator ikey = keys.begin(); ikey != keys.end();
        ++ikey) {
     if (mIsData) {
       if (!MEtXYcorParaColl->isShiftData(*ikey))
         throw cms::Exception("MultShiftMETcorrDBInputProducer::produce")
             << "DB is not for Data. Set proper option: \"corrPfMetXYMultDB.isData\" !!\n";
     } else {
       if (MEtXYcorParaColl->isShiftData(*ikey))
         throw cms::Exception("MultShiftMETcorrDBInputProducer::produce")
             << "DB is for Data. Set proper option: \"corrPfMetXYMultDB.isData\" !!\n";
     }
   }
  
   for (std::vector<MEtXYcorrectParametersCollection::key_type>::const_iterator ikey = keys.begin(); ikey != keys.end();
        ++ikey) {
     if (!mIsData) {
       if (mSampleType == "MC") {
         if (!MEtXYcorParaColl->isShiftMC(*ikey))
           continue;
       } else if (mSampleType == "DY") {
         if (!MEtXYcorParaColl->isShiftDY(*ikey))
           continue;
       } else if (mSampleType == "TTJets") {
         if (!MEtXYcorParaColl->isShiftTTJets(*ikey))
           continue;
       } else if (mSampleType == "WJets") {
         if (!MEtXYcorParaColl->isShiftWJets(*ikey))
           continue;
       } else
         throw cms::Exception("MultShiftMETcorrDBInputProducer::produce")
             << "SampleType: " << mSampleType << " is not reserved !!!\n";
     }
  
     std::string sectionName = MEtXYcorParaColl->findLabel(*ikey);
     MEtXYcorrectParameters const& MEtXYcorParams = (*MEtXYcorParaColl)[*ikey];
  
     counts_ = 0;
     sumPt_ = 0;
  
     for (unsigned i = 0; i < particleFlow->size(); ++i) {
       const reco::Candidate& c = particleFlow->at(i);
       if (abs(c.pdgId()) ==
           translateTypeToAbsPdgId(reco::PFCandidate::ParticleType(MEtXYcorParams.definitions().PtclType()))) {
         if ((c.eta() > MEtXYcorParams.record(0).xMin(0)) and (c.eta() < MEtXYcorParams.record(0).xMax(0))) {
           float weight = (!weightsToken_.isUninitialized()) ? (*weights)[particleFlow->ptrAt(i)] : 1.0;
           counts_ += (weight > 0);
           sumPt_ += c.pt() * weight;
           continue;
         }
       }
     }
     double val(0.);
     unsigned parVar = MEtXYcorParams.definitions().parVar(0);
  
     if (parVar == 0) {
       val = counts_;
  
     } else if (parVar == 1) {
       val = ngoodVertices;
  
     } else if (parVar == 2) {
       val = sumPt_;
  
     } else {
       throw cms::Exception("MultShiftMETcorrDBInputProducer::produce")
           << "parVar: " << parVar << " is not reserved !!!\n";
     }
  
     formula_x_ = std::make_unique<TF1>("corrPx", MEtXYcorParams.definitions().formula().c_str());
     formula_y_ = std::make_unique<TF1>("corrPy", MEtXYcorParams.definitions().formula().c_str());
  
     for (unsigned i(0); i < MEtXYcorParams.record(0).nParameters(); i++) {
       formula_x_->SetParameter(i, MEtXYcorParams.record(0).parameter(i));
     }
     for (unsigned i(0); i < MEtXYcorParams.record(1).nParameters(); i++) {
       formula_y_->SetParameter(i, MEtXYcorParams.record(1).parameter(i));
     }
  
     corx -= formula_x_->Eval(val);
     cory -= formula_y_->Eval(val);
  
   }  //end loop over corrections
  
   metCorr->mex = corx;
   metCorr->mey = cory;
   evt.put(std::move(metCorr), "");
 }

◆ translateTypeToAbsPdgId()

int MultShiftMETcorrDBInputProducer::translateTypeToAbsPdgId ( reco::PFCandidate::ParticleType type )

staticprivate

Definition at line 16 of file MultShiftMETcorrDBInputProducer.cc.

                                                                                              {
   switch (type) {
     case reco::PFCandidate::ParticleType::h:
       return 211;  // pi+
     case reco::PFCandidate::ParticleType::e:
       return 11;
     case reco::PFCandidate::ParticleType::mu:
       return 13;
     case reco::PFCandidate::ParticleType::gamma:
       return 22;
     case reco::PFCandidate::ParticleType::h0:
       return 130;  // K_L0
     case reco::PFCandidate::ParticleType::h_HF:
       return 1;  // dummy pdg code
     case reco::PFCandidate::ParticleType::egamma_HF:
       return 2;  // dummy pdg code
     case reco::PFCandidate::ParticleType::X:
     default:
       return 0;
   }
 }