#include <MultShiftMETcorrDBInputProducer.h>

Inheritance diagram for MultShiftMETcorrDBInputProducer:

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

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

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

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

override

Definition at line 50 of file MultShiftMETcorrDBInputProducer.cc.

51 {

52 }

Member Function Documentation

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

overrideprivate

Definition at line 54 of file MultShiftMETcorrDBInputProducer.cc.

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