#include "PhysicsTools/PatAlgos/interface/PATMETProducer.h"

Inheritance diagram for pat::PATMETProducer:

Public Member Functions
	PATMETProducer (const edm::ParameterSet &iConfig)

virtual void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override

	~PATMETProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	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 (const std::string &iProcessName, std::vector< const char * > &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 ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
const reco::METCovMatrix	getMETCovMatrix (const edm::Event &event) const

Private Attributes
bool	addEfficiencies_

bool	addGenMET_

bool	addMuonCorr_

bool	addResolutions_

bool	calculateMETSignificance_

pat::helper::EfficiencyLoader	efficiencyLoader_

GreaterByEt< MET >	eTComparator_

edm::EDGetTokenT< edm::View < reco::GenMET > >	genMETToken_

edm::EDGetTokenT< edm::View < reco::Jet > >	jetToken_

std::vector< edm::EDGetTokenT < edm::View< reco::Candidate > > >	lepTokens_

metsig::METSignificance *	metSigAlgo_

edm::InputTag	metSrc_

edm::EDGetTokenT< edm::View < reco::MET > >	metToken_

edm::InputTag	muonSrc_

edm::EDGetTokenT< edm::View < reco::Candidate > >	pfCandToken_

pat::helper::KinResolutionsLoader	resolutionLoader_

pat::PATUserDataHelper< pat::MET >	userDataHelper_

bool	useUserData_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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

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

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

Produces the pat::MET.

The PATMETProducer produces the analysis-level pat::MET starting from a collection of objects of METType.

Author: Steven Lowette

Version

Id:: PATMETProducer.h,v 1.10 2009/06/25 23:49:35 gpetrucc Exp

Definition at line 38 of file PATMETProducer.h.

Constructor & Destructor Documentation

PATMETProducer::PATMETProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 18 of file PATMETProducer.cc.

References addEfficiencies_, addGenMET_, addResolutions_, calculateMETSignificance_, edm::EDConsumerBase::consumesCollector(), efficiencyLoader_, genMETToken_, edm::ParameterSet::getParameter(), jetToken_, lepTokens_, edm::EDConsumerBase::mayConsume(), metSigAlgo_, METSignificance_cfi::METSignificance, metSrc_, metToken_, pfCandToken_, resolutionLoader_, userDataHelper_, and useUserData_.

                                                              :
   useUserData_(iConfig.exists("userData"))
 {
   // initialize the configurables
   metSrc_         = iConfig.getParameter<edm::InputTag>("metSource");
   metToken_       = consumes<edm::View<reco::MET> >(metSrc_);
   addGenMET_      = iConfig.getParameter<bool>         ("addGenMET");
   genMETToken_    = mayConsume<edm::View<reco::GenMET> >(iConfig.getParameter<edm::InputTag>("genMETSource"));
   addResolutions_ = iConfig.getParameter<bool>         ("addResolutions");
 
   // Efficiency configurables
   addEfficiencies_ = iConfig.getParameter<bool>("addEfficiencies");
   if (addEfficiencies_) {
      efficiencyLoader_ = pat::helper::EfficiencyLoader(iConfig.getParameter<edm::ParameterSet>("efficiencies"), consumesCollector());
   }
 
   // Resolution configurables
   addResolutions_ = iConfig.getParameter<bool>("addResolutions");
   if (addResolutions_) {
      resolutionLoader_ = pat::helper::KinResolutionsLoader(iConfig.getParameter<edm::ParameterSet>("resolutions"));
   }
 
   // Check to see if the user wants to add user data
   if ( useUserData_ ) {
     userDataHelper_ = PATUserDataHelper<MET>(iConfig.getParameter<edm::ParameterSet>("userData"), consumesCollector());
   }
 
   // MET Significance
   calculateMETSignificance_ = iConfig.getParameter<bool>("computeMETSignificance");
   if(calculateMETSignificance_)
     {
       metSigAlgo_ = new metsig::METSignificance(iConfig);
       jetToken_ = mayConsume<edm::View<reco::Jet> >(iConfig.getParameter<edm::InputTag>("srcJets"));
       pfCandToken_ = mayConsume<edm::View<reco::Candidate> >(iConfig.getParameter<edm::InputTag>("srcPFCands"));
       std::vector<edm::InputTag> srcLeptonsTags = iConfig.getParameter< std::vector<edm::InputTag> >("srcLeptons");
       for(std::vector<edm::InputTag>::const_iterator it=srcLeptonsTags.begin();it!=srcLeptonsTags.end();it++) {
     lepTokens_.push_back( mayConsume<edm::View<reco::Candidate> >( *it ) );
       }
     }  
 
   // produces vector of mets
   produces<std::vector<MET> >();
 }

PATMETProducer::~PATMETProducer ( )

Definition at line 63 of file PATMETProducer.cc.

63 {

64 }

Member Function Documentation

void PATMETProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 131 of file PATMETProducer.cc.

References edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addOptional(), pat::helper::KinResolutionsLoader::fillDescription(), pat::PATUserDataHelper< ObjectType >::fillDescription(), HLT_25ns14e33_v1_cff::InputTag, edm::ParameterSetDescription::setAllowAnything(), and edm::ParameterSetDescription::setComment().

 {
   edm::ParameterSetDescription iDesc;
   iDesc.setComment("PAT MET producer module");
 
   // input source
   iDesc.add<edm::InputTag>("metSource", edm::InputTag("no default"))->setComment("input collection");
 
   // MC configurations
   iDesc.add<bool>("addGenMET", false);
   iDesc.add<edm::InputTag>("genMETSource", edm::InputTag("genMetCalo"));
 
   pat::helper::KinResolutionsLoader::fillDescription(iDesc);
 
   // Efficiency configurables
   edm::ParameterSetDescription efficienciesPSet;
   efficienciesPSet.setAllowAnything(); // TODO: the pat helper needs to implement a description.
   iDesc.add("efficiencies", efficienciesPSet);
   iDesc.add<bool>("addEfficiencies", false);
 
   // Check to see if the user wants to add user data
   edm::ParameterSetDescription userDataPSet;
   PATUserDataHelper<MET>::fillDescription(userDataPSet);
   iDesc.addOptional("userData", userDataPSet);
 
   // muon correction
   iDesc.add<bool>("addMuonCorrections", false);
   iDesc.add<edm::InputTag>("muonSource", edm::InputTag("muons"));
 
 }

const reco::METCovMatrix PATMETProducer::getMETCovMatrix ( const edm::Event & event ) const

private

Definition at line 163 of file PATMETProducer.cc.

References metsig::METSignificance::getCovariance(), HLT_25ns14e33_v1_cff::inputJets, jetToken_, lepTokens_, HLT_25ns14e33_v1_cff::leptons, metSigAlgo_, and pfCandToken_.

Referenced by produce().

                                                            {
   std::vector< edm::Handle<reco::CandidateView> > leptons;
   for ( std::vector<edm::EDGetTokenT<edm::View<reco::Candidate> > >::const_iterator srcLeptons_i = lepTokens_.begin();
     srcLeptons_i != lepTokens_.end(); ++srcLeptons_i ) {
     edm::Handle<reco::CandidateView> leptons_i;
     event.getByToken(*srcLeptons_i, leptons_i);
     leptons.push_back( leptons_i );
   }
   // jets
   edm::Handle<edm::View<reco::Jet> > inputJets;
   event.getByToken( jetToken_, inputJets );
 
   //candidates
   edm::Handle<edm::View<reco::Candidate> > inputCands;
   event.getByToken( pfCandToken_, inputCands );
 
   //Compute the covariance matrix and fill it
   reco::METCovMatrix cov = metSigAlgo_->getCovariance( *inputJets, leptons, *inputCands);
   return cov;
 }

void PATMETProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 67 of file PATMETProducer.cc.

References addGenMET_, calculateMETSignificance_, efficiencyLoader_, pat::helper::EfficiencyLoader::enabled(), pat::helper::KinResolutionsLoader::enabled(), edm::InputTag::encode(), genMETToken_, edm::Event::getByToken(), getMETCovMatrix(), metsig::METSignificance::getSignificance(), customizeTrackingMonitorSeedNumber::idx, metSigAlgo_, metSrc_, metToken_, pat::helper::EfficiencyLoader::newEvent(), pat::helper::KinResolutionsLoader::newEvent(), metProducer_cfi::patMETs, edm::Event::put(), resolutionLoader_, pat::helper::EfficiencyLoader::setEfficiencies(), pat::MET::setGenMET(), pat::MET::setMETSignificance(), pat::helper::KinResolutionsLoader::setResolutions(), reco::MET::setSignificanceMatrix(), userDataHelper_, and useUserData_.

                                                                             {
 
   // Get the vector of MET's from the event
   edm::Handle<edm::View<reco::MET> > mets;
   iEvent.getByToken(metToken_, mets);
 
   if (mets->size() != 1) throw cms::Exception("Corrupt Data") << "The input MET collection " << metSrc_.encode() << " has size " << mets->size() << " instead of 1 as it should.\n";
   if (efficiencyLoader_.enabled()) efficiencyLoader_.newEvent(iEvent);
   if (resolutionLoader_.enabled()) resolutionLoader_.newEvent(iEvent, iSetup);
 
   // Get the vector of generated met from the event if needed
   edm::Handle<edm::View<reco::GenMET> > genMETs;
   if (addGenMET_) {
     iEvent.getByToken(genMETToken_, genMETs);
   }
 
   // loop over mets
   std::vector<MET> * patMETs = new std::vector<MET>();
   for (edm::View<reco::MET>::const_iterator itMET = mets->begin(); itMET != mets->end(); itMET++) {
     // construct the MET from the ref -> save ref to original object
     unsigned int idx = itMET - mets->begin();
     edm::RefToBase<reco::MET> metsRef = mets->refAt(idx);
     edm::Ptr<reco::MET> metsPtr = mets->ptrAt(idx);
     MET amet(metsRef);
     // add the generated MET
     if (addGenMET_) amet.setGenMET((*genMETs)[idx]);
 
     //add the MET significance
     if(calculateMETSignificance_) {
       const reco::METCovMatrix& sigcov = getMETCovMatrix(iEvent);
       amet.setSignificanceMatrix(sigcov);
       double metSig=metSigAlgo_->getSignificance(sigcov, amet);
       amet.setMETSignificance(metSig);
     }
 
     if (efficiencyLoader_.enabled()) {
         efficiencyLoader_.setEfficiencies( amet, metsRef );
     }
 
     if (resolutionLoader_.enabled()) {
         resolutionLoader_.setResolutions(amet);
     }
 
 
     if ( useUserData_ ) {
       userDataHelper_.add( amet, iEvent, iSetup );
     }
 
 
     // correct for muons if demanded... never more: it's now done by JetMETCorrections
     // add the MET to the vector of METs
     patMETs->push_back(amet);
   }
 
   // sort MET in ET .. don't mess with this
   //  std::sort(patMETs->begin(), patMETs->end(), eTComparator_);
 
   // put genEvt object in Event
   std::auto_ptr<std::vector<MET> > myMETs(patMETs);
   iEvent.put(myMETs);
 
 }