MultShiftMETcorrDBInputProducer Class Reference

#include <MultShiftMETcorrDBInputProducer.h>

Inheritance diagram for MultShiftMETcorrDBInputProducer:

Public Member Functions
	MultShiftMETcorrDBInputProducer (const edm::ParameterSet &)
	~MultShiftMETcorrDBInputProducer ()
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)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
virtual	~ProducerBase () noexcept(false)
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

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
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

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 31 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

(

)

Definition at line 50 of file MultShiftMETcorrDBInputProducer.cc.

{
}

Member Function Documentation

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

overrideprivate

Definition at line 54 of file MultShiftMETcorrDBInputProducer.cc.

References funct::abs(), EnergyCorrector::c, counts_, DEFINE_FWK_MODULE, MEtXYcorrectParameters::definitions(), reco::Candidate::eta(), Exception, MEtXYcorrectParametersCollection::findLabel(), MEtXYcorrectParameters::Definitions::formula(), formula_x_, formula_y_, edm::EventSetup::get(), edm::Event::getByToken(), metFilters_cff::goodVertices, mps_fire::i, MEtXYcorrectParametersCollection::isShiftData(), MEtXYcorrectParametersCollection::isShiftDY(), MEtXYcorrectParametersCollection::isShiftMC(), MEtXYcorrectParametersCollection::isShiftTTJets(), MEtXYcorrectParametersCollection::isShiftWJets(), edm::HandleBase::isValid(), relativeConstraints::keys, mIsData, eostools::move(), mPayloadName, mSampleType, MEtXYcorrectParameters::Record::nParameters(), MEtXYcorrectParameters::Record::parameter(), pfLinker_cff::particleFlow, MEtXYcorrectParameters::Definitions::parVar(), reco::Candidate::pdgId(), pflow_, position, reco::Candidate::pt(), MEtXYcorrectParameters::Definitions::PtclType(), edm::Event::put(), MEtXYcorrectParameters::record(), rho, AlCaHLTBitMon_QueryRunRegistry::string, sumPt_, translateTypeToAbsPdgId(), heppy_batch::val, MEtXYcorrectParametersCollection::validKeys(), vertices_, MEtXYcorrectParameters::Record::xMax(), MEtXYcorrectParameters::Record::xMin(), and z.

{
  // 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::auto_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, 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;
  }
}

Member Data Documentation

std::vector<edm::ParameterSet> MultShiftMETcorrDBInputProducer::cfgCorrParameters_

private

Definition at line 52 of file MultShiftMETcorrDBInputProducer.h.

int MultShiftMETcorrDBInputProducer::counts_

private

Definition at line 55 of file MultShiftMETcorrDBInputProducer.h.

Referenced by produce().

std::vector<double> MultShiftMETcorrDBInputProducer::etaMax_

private

Definition at line 54 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer().

std::vector<double> MultShiftMETcorrDBInputProducer::etaMin_

private

Definition at line 54 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer().

std::unique_ptr< TF1 > MultShiftMETcorrDBInputProducer::formula_x_

private

Definition at line 57 of file MultShiftMETcorrDBInputProducer.h.

Referenced by produce().

std::unique_ptr< TF1 > MultShiftMETcorrDBInputProducer::formula_y_

private

Definition at line 58 of file MultShiftMETcorrDBInputProducer.h.

Referenced by produce().

bool MultShiftMETcorrDBInputProducer::mIsData

private

Definition at line 50 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().

std::string MultShiftMETcorrDBInputProducer::moduleLabel_

private

Definition at line 47 of file MultShiftMETcorrDBInputProducer.h.

Referenced by Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), and Modules.ESProducer::nameInProcessDesc_().

std::string MultShiftMETcorrDBInputProducer::mPayloadName

private

Definition at line 48 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().

std::string MultShiftMETcorrDBInputProducer::mSampleType

private

Definition at line 49 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().

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

private

Definition at line 45 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().

double MultShiftMETcorrDBInputProducer::sumPt_

private

Definition at line 56 of file MultShiftMETcorrDBInputProducer.h.

Referenced by produce().

edm::EDGetTokenT<edm::View<reco::Vertex> > MultShiftMETcorrDBInputProducer::vertices_

private

Definition at line 46 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().