#include <MultShiftMETcorrInputProducer.h>

Inheritance diagram for MultShiftMETcorrInputProducer:

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

	~MultShiftMETcorrInputProducer ()

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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 bitset_utilities::append(), cfgCorrParameters_, counts_, etaMax_, etaMin_, formula_x_, formula_y_, edm::ParameterSet::getParameter(), 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 ( )

Definition at line 63 of file MultShiftMETcorrInputProducer.cc.

64 {

65 }

Member Function Documentation

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

overrideprivatevirtual

Implements edm::stream::EDProducerBase.

Definition at line 67 of file MultShiftMETcorrInputProducer.cc.

References funct::abs(), EnergyCorrector::c, cfgCorrParameters_, counts_, reco::Candidate::eta(), etaMax_, etaMin_, formula_x_, formula_y_, edm::Event::getByToken(), metFilters_cff::goodVertices, i, edm::HandleBase::isValid(), j, pfLinker_cff::particleFlow, reco::Candidate::pdgId(), pflow_, position, reco::Candidate::pt(), edm::Event::put(), rho, sumPt_, translateTypeToAbsPdgId(), type_, findQualityFiles::v, 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::auto_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(metCorr, "");
   
 }

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

staticprivate

Definition at line 12 of file MultShiftMETcorrInputProducer.cc.

References alignCSCRings::e, 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;
   }
 }