#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

	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_

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_

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

Definition at line 32 of file MultShiftMETcorrInputProducer.h.

Constructor & Destructor Documentation

◆ MultShiftMETcorrInputProducer()

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

explicit

Definition at line 37 of file MultShiftMETcorrInputProducer.cc.

References mps_setup::append, looper::cfg, cfgCorrParameters_, counts_, etaMax_, etaMin_, formula_x_, formula_y_, mps_fire::i, moduleLabel_, pflow_, HLT_2023v12_cff::srcWeights, AlCaHLTBitMon_QueryRunRegistry::string, sumPt_, type_, findQualityFiles::v, varType_, vertices_, and weightsToken_.

     : 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"));
 
   edm::InputTag srcWeights = cfg.getParameter<edm::InputTag>("srcWeights");
   if (!srcWeights.label().empty())
     weightsToken_ = consumes<edm::ValueMap<float>>(srcWeights);
 
   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::make_unique<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::make_unique<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::~MultShiftMETcorrInputProducer ( )

override

Definition at line 82 of file MultShiftMETcorrInputProducer.cc.

82 {}

Member Function Documentation

◆ produce()

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

overrideprivate

Definition at line 84 of file MultShiftMETcorrInputProducer.cc.

References funct::abs(), HltBtagPostValidation_cff::c, cfgCorrParameters_, counts_, etaMax_, etaMin_, spr::find(), formula_x_, formula_y_, edm::Event::getByToken(), mps_fire::i, edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), dqmiolumiharvest::j, eostools::move(), TrackCollections2monitor_cff::nVertices, pfLinker_cff::particleFlow, pflow_, position, edm::Event::put(), rho, sumPt_, translateTypeToAbsPdgId(), type_, parallelization::uint, findQualityFiles::v, heppy_batch::val, varType_, vertices_, mps_merge::weight, hltDeepSecondaryVertexTagInfosPFPuppi_cfi::weights, weightsToken_, 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 ";
   }
   uint ngoodVertices = 0;
   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))
       ngoodVertices += 1;
   }
   uint nVertices = hpv->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);
 
   edm::Handle<edm::ValueMap<float>> weights;
   if (!weightsToken_.isUninitialized())
     evt.getByToken(weightsToken_, weights);
   if (std::find(varType_.begin(), varType_.end(), 0) != varType_.end() ||
       std::find(varType_.begin(), varType_.end(), 2) != varType_.end()) {
     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])) {
             float weight = (!weightsToken_.isUninitialized()) ? (*weights)[particleFlow->ptrAt(i)] : 1.0;
             counts_[j] += (weight > 0);
             sumPt_[j] += c.pt() * weight;
             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];
     } else if (varType_[j] == 1) {
       val = ngoodVertices;
     } else if (varType_[j] == 2) {
       val = sumPt_[j];
     } else if (varType_[j] == 3) {
       val = nVertices;
     }
 
     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), "");
 }

◆ translateTypeToAbsPdgId()

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

staticprivate

Definition at line 15 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;
   }
 }