#include <MultShiftMETcorrInputProducer.h>

Inheritance diagram for MultShiftMETcorrInputProducer:

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

	~MultShiftMETcorrInputProducer () 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_

std::vector< int >	counts_

std::vector< double >	etaMax_

std::vector< double >	etaMin_

std::vector< std::unique_ptr< TF1 > >	formula_x_

std::vector< std::unique_ptr< TF1 > >	formula_y_

std::string	moduleLabel_

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

std::vector< double >	sumPt_

std::vector< int >	type_

std::vector< int >	varType_

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

Definition at line 31 of file MultShiftMETcorrInputProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 34 of file MultShiftMETcorrInputProducer.cc.

References mps_setup::append, cfgCorrParameters_, counts_, etaMax_, etaMin_, formula_x_, formula_y_, edm::ParameterSet::getParameter(), mps_fire::i, moduleLabel_, pflow_, AlCaHLTBitMon_QueryRunRegistry::string, sumPt_, type_, findQualityFiles::v, varType_, and vertices_.

     : moduleLabel_(cfg.getParameter<std::string>("@module_label")) {
   pflow_ = consumes<edm::View<reco::Candidate> >(cfg.getParameter<edm::InputTag>("srcPFlow"));
   vertices_ = consumes<edm::View<reco::Vertex> >(cfg.getParameter<edm::InputTag>("vertexCollection"));
 
   cfgCorrParameters_ = cfg.getParameter<std::vector<edm::ParameterSet> >("parameters");
   etaMin_.clear();
   etaMax_.clear();
   type_.clear();
   varType_.clear();
 
   produces<CorrMETData>();
 
   for (std::vector<edm::ParameterSet>::const_iterator v = cfgCorrParameters_.begin(); v != cfgCorrParameters_.end();
        v++) {
     TString corrPxFormula = v->getParameter<std::string>("fx");
     TString corrPyFormula = v->getParameter<std::string>("fy");
     std::vector<double> corrPxParams = v->getParameter<std::vector<double> >("px");
     std::vector<double> corrPyParams = v->getParameter<std::vector<double> >("py");
 
     formula_x_.push_back(std::unique_ptr<TF1>(new TF1(
         std::string(moduleLabel_).append("_").append(v->getParameter<std::string>("name")).append("_corrPx").c_str(),
         v->getParameter<std::string>("fx").c_str())));
     formula_y_.push_back(std::unique_ptr<TF1>(new TF1(
         std::string(moduleLabel_).append("_").append(v->getParameter<std::string>("name")).append("_corrPy").c_str(),
         v->getParameter<std::string>("fy").c_str())));
 
     for (unsigned i = 0; i < corrPxParams.size(); i++)
       formula_x_.back()->SetParameter(i, corrPxParams[i]);
     for (unsigned i = 0; i < corrPyParams.size(); i++)
       formula_y_.back()->SetParameter(i, corrPyParams[i]);
 
     counts_.push_back(0);
     sumPt_.push_back(0.);
     etaMin_.push_back(v->getParameter<double>("etaMin"));
     etaMax_.push_back(v->getParameter<double>("etaMax"));
     type_.push_back(v->getParameter<int>("type"));
     varType_.push_back(v->getParameter<int>("varType"));
   }
 }

MultShiftMETcorrInputProducer::~MultShiftMETcorrInputProducer ( )

override

Definition at line 75 of file MultShiftMETcorrInputProducer.cc.

75 {}

Member Function Documentation

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

overrideprivate

Definition at line 77 of file MultShiftMETcorrInputProducer.cc.

References funct::abs(), HltBtagPostValidation_cff::c, cfgCorrParameters_, counts_, DEFINE_FWK_MODULE, reco::Candidate::eta(), etaMax_, etaMin_, formula_x_, formula_y_, edm::Event::getByToken(), metFilters_cff::goodVertices, mps_fire::i, edm::HandleBase::isValid(), dqmiolumiharvest::j, eostools::move(), pfLinker_cff::particleFlow, reco::Candidate::pdgId(), pflow_, position, reco::Candidate::pt(), edm::Event::put(), rho, sumPt_, translateTypeToAbsPdgId(), type_, findQualityFiles::v, heppy_batch::val, varType_, vertices_, and z.

                                                                                   {
   //get primary vertices
   edm::Handle<edm::View<reco::Vertex> > hpv;
   evt.getByToken(vertices_, hpv);
   if (!hpv.isValid()) {
     edm::LogError("MultShiftMETcorrInputProducer::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();
 
   for (unsigned i = 0; i < counts_.size(); i++)
     counts_[i] = 0;
   for (unsigned i = 0; i < sumPt_.size(); i++)
     sumPt_[i] = 0.;
 
   edm::Handle<edm::View<reco::Candidate> > particleFlow;
   evt.getByToken(pflow_, particleFlow);
   for (unsigned i = 0; i < particleFlow->size(); ++i) {
     const reco::Candidate& c = particleFlow->at(i);
     for (unsigned j = 0; j < type_.size(); j++) {
       if (abs(c.pdgId()) == translateTypeToAbsPdgId(reco::PFCandidate::ParticleType(type_[j]))) {
         if ((c.eta() > etaMin_[j]) and (c.eta() < etaMax_[j])) {
           counts_[j] += 1;
           sumPt_[j] += c.pt();
           continue;
         }
       }
     }
   }
 
   //MM: loop over all constituent types and sum each correction
   std::unique_ptr<CorrMETData> metCorr(new CorrMETData());
 
   double corx = 0.;
   double cory = 0.;
 
   for (std::vector<edm::ParameterSet>::const_iterator v = cfgCorrParameters_.begin(); v != cfgCorrParameters_.end();
        v++) {
     unsigned j = v - cfgCorrParameters_.begin();
 
     double val(0.);
     if (varType_[j] == 0) {
       val = counts_[j];
     }
     if (varType_[j] == 1) {
       val = ngoodVertices;
     }
     if (varType_[j] == 2) {
       val = sumPt_[j];
     }
 
     corx -= formula_x_[j]->Eval(val);
     cory -= formula_y_[j]->Eval(val);
 
   }  //end loop over corrections
 
   metCorr->mex = corx;
   metCorr->mey = cory;
   evt.put(std::move(metCorr), "");
 }

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

staticprivate

Definition at line 12 of file MultShiftMETcorrInputProducer.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;
   }
 }