#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	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 27 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 63 of file MultShiftMETcorrInputProducer.cc.

64 {

65 }

Member Function Documentation

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

overrideprivate

Definition at line 67 of file MultShiftMETcorrInputProducer.cc.

References funct::abs(), EnergyCorrector::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(), 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, RPCpg::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;
   }
 }