#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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

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_

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

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	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 32 of file MultShiftMETcorrDBInputProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 38 of file MultShiftMETcorrDBInputProducer.cc.

References etaMax_, etaMin_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), mIsData, mPayloadName, mSampleType, pflow_, AlCaHLTBitMon_QueryRunRegistry::string, and vertices_.

     : 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"));
 
   etaMin_.clear();
   etaMax_.clear();
 
   produces<CorrMETData>();
 }

MultShiftMETcorrDBInputProducer::~MultShiftMETcorrDBInputProducer ( )

override

Definition at line 53 of file MultShiftMETcorrDBInputProducer.cc.

53 {}

Member Function Documentation

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

overrideprivate

Definition at line 55 of file MultShiftMETcorrDBInputProducer.cc.

                                                                                     {
   // Get para.s from DB
   edm::ESHandle<MEtXYcorrectParametersCollection> MEtXYcorParaColl;
   es.get<MEtXYcorrectRecord>().get(mPayloadName, MEtXYcorParaColl);
 
   // 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);
 
   //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))) {
           counts_ += 1;
           sumPt_ += c.pt();
           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_.reset(new TF1("corrPx", MEtXYcorParams.definitions().formula().c_str()));
     formula_y_.reset(new 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), "");
 }

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

staticprivate

Definition at line 16 of file MultShiftMETcorrDBInputProducer.cc.

References MillePedeFileConverter_cfg::e, CustomPhysics_cfi::gamma, h, amptDefaultParameters_cff::mu, and X.

Referenced by produce().

                                                                                              {
   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;
   }
 }