#include <METProducer.h>

Inheritance diagram for cms::METProducer:

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

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

virtual	~METProducer ()

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

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_CaloMET (edm::Event &event)

void	produce_else (edm::Event &event)

void	produce_GenMET (edm::Event &event)

void	produce_PFClusterMET (edm::Event &event)

void	produce_PFMET (edm::Event &event)

void	produce_TCMET (edm::Event &event, const edm::EventSetup &setup)

Private Attributes
std::string	alias

bool	applyFiducialThresholdForFractions

bool	calculateSignificance_

double	globalThreshold

edm::InputTag	inputLabel

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

std::string	inputType

edm::InputTag	jetsLabel_

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

std::string	METtype

bool	noHF

bool	onlyFiducial

metsig::SignAlgoResolutions *	resolutions_

TCMETAlgo	tcMetAlgo_

bool	usePt

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

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Definition at line 46 of file METProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 55 of file METProducer.cc.

References alias, applyFiducialThresholdForFractions, calculateSignificance_, TCMETAlgo::configure(), edm::EDConsumerBase::consumesCollector(), edm::ParameterSet::getParameter(), jetsLabel_, jetToken_, METtype, noHF, onlyFiducial, resolutions_, tcMetAlgo_, and usePt.

     : inputLabel(iConfig.getParameter<edm::InputTag>("src"))
     , inputType(iConfig.getParameter<std::string>("InputType"))
     , METtype(iConfig.getParameter<std::string>("METType"))
     , alias(iConfig.getParameter<std::string>("alias"))
     , inputToken_(consumes<edm::View<reco::Candidate> >(inputLabel))
     , calculateSignificance_(false)
     , resolutions_(0)
     , globalThreshold(iConfig.getParameter<double>("globalThreshold"))
   {
 
 
     if( METtype == "CaloMET" ) 
       {
     noHF = iConfig.getParameter<bool>("noHF");
     produces<reco::CaloMETCollection>().setBranchAlias(alias.c_str()); 
     calculateSignificance_ = iConfig.getParameter<bool>("calculateSignificance");
       }
     else if( METtype == "GenMET" )  
       {
     onlyFiducial = iConfig.getParameter<bool>("onlyFiducialParticles");
         usePt = iConfig.getParameter<bool>("usePt");
         applyFiducialThresholdForFractions = iConfig.getParameter<bool>("applyFiducialThresholdForFractions");
     produces<reco::GenMETCollection>().setBranchAlias(alias.c_str());
       }
     else if( METtype == "PFMET" )
       {
     produces<reco::PFMETCollection>().setBranchAlias(alias.c_str()); 
 
     calculateSignificance_ = iConfig.getParameter<bool>("calculateSignificance");
 
     if(calculateSignificance_)
       {
         jetsLabel_ = iConfig.getParameter<edm::InputTag>("jets");
         jetToken_ = consumes<edm::View<reco::PFJet> >(iConfig.getParameter<edm::InputTag>("jets"));
       }
 
       }
     else if( METtype == "PFClusterMET" )
       {
     produces<reco::PFClusterMETCollection>().setBranchAlias(alias.c_str()); 
       }
     else if (METtype == "TCMET" )
       {
     produces<reco::METCollection>().setBranchAlias(alias.c_str());
     tcMetAlgo_.configure(iConfig, consumesCollector());
       }
     else                            
       produces<reco::METCollection>().setBranchAlias(alias.c_str()); 
 
     if (calculateSignificance_ && ( METtype == "CaloMET" || METtype == "PFMET"))
       {
     resolutions_ = new metsig::SignAlgoResolutions(iConfig);
       }
   }

virtual cms::METProducer::~METProducer ( )

inlinevirtual

Definition at line 50 of file METProducer.h.

50 { }

Member Function Documentation

void METProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	setup
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 112 of file METProducer.cc.

References METtype, produce_CaloMET(), produce_else(), produce_GenMET(), produce_PFClusterMET(), produce_PFMET(), and produce_TCMET().

   {
     if( METtype == "CaloMET" ) 
       {
     produce_CaloMET(event);
     return;
       }
 
     if( METtype == "TCMET" )
       {
     produce_TCMET(event, setup);
     return;
       }
 
     if( METtype == "PFMET" )
       {
     produce_PFMET(event);
     return;
       }
 
     if( METtype == "PFClusterMET" )
       {
     produce_PFClusterMET(event);
     return;
       }
 
     if( METtype == "GenMET" ) 
       {
     produce_GenMET(event);
     return;
       }
 
     produce_else(event);
   }

void METProducer::produce_CaloMET ( edm::Event & event )

private

Definition at line 147 of file METProducer.cc.

References CaloSpecificAlgo::addInfo(), SignCaloSpecificAlgo::calculateBaseCaloMET(), calculateSignificance_, SignCaloSpecificAlgo::getSignificance(), SignCaloSpecificAlgo::getSignificanceMatrix(), globalThreshold, input, inputToken_, noHF, edm::Handle< T >::product(), resolutions_, METAlgo::run(), reco::CaloMET::SetMetSignificance(), and reco::MET::setSignificanceMatrix().

Referenced by produce().

   {
     edm::Handle<edm::View<reco::Candidate> > input;
     event.getByToken(inputToken_, input);
 
     METAlgo algo;
     CommonMETData commonMETdata = algo.run(*input.product(), globalThreshold);
 
     CaloSpecificAlgo calospecalgo;
     reco::CaloMET calomet = calospecalgo.addInfo(input, commonMETdata, noHF, globalThreshold);
 
     if( calculateSignificance_ ) 
       {
     SignCaloSpecificAlgo signcalospecalgo;
     signcalospecalgo.calculateBaseCaloMET(input, commonMETdata, *resolutions_, noHF, globalThreshold);
     calomet.SetMetSignificance(signcalospecalgo.getSignificance() );
     calomet.setSignificanceMatrix(signcalospecalgo.getSignificanceMatrix());
       }
 
     std::auto_ptr<reco::CaloMETCollection> calometcoll;
     calometcoll.reset(new reco::CaloMETCollection);
     calometcoll->push_back( calomet ) ;
     event.put( calometcoll );  
   }

void METProducer::produce_else ( edm::Event & event )

private

Definition at line 244 of file METProducer.cc.

References globalThreshold, input, inputToken_, CaloMET_cfi::met, CommonMETData::met, CommonMETData::mex, CommonMETData::mey, p4, edm::Handle< T >::product(), METAlgo::run(), and CommonMETData::sumet.

Referenced by produce().

   {
     edm::Handle<edm::View<reco::Candidate> > input;
     event.getByToken(inputToken_, input);
 
     METAlgo algo;
     CommonMETData commonMETdata = algo.run(*input.product(), globalThreshold);
 
     math::XYZTLorentzVector p4( commonMETdata.mex, commonMETdata.mey, 0.0, commonMETdata.met);
     math::XYZPoint vtx(0,0,0);
     reco::MET met( commonMETdata.sumet, p4, vtx );
     std::auto_ptr<reco::METCollection> metcoll;
     metcoll.reset(new reco::METCollection);
     metcoll->push_back( met );
     event.put( metcoll );
   }

void METProducer::produce_GenMET ( edm::Event & event )

private

Definition at line 230 of file METProducer.cc.

References GenSpecificAlgo::addInfo(), applyFiducialThresholdForFractions, relval_steps::gen(), globalThreshold, input, inputToken_, onlyFiducial, and usePt.

Referenced by produce().

   {
     edm::Handle<edm::View<reco::Candidate> > input;
     event.getByToken(inputToken_, input);
 
     CommonMETData commonMETdata;
 
     GenSpecificAlgo gen;
     std::auto_ptr<reco::GenMETCollection> genmetcoll;
     genmetcoll.reset (new reco::GenMETCollection);
     genmetcoll->push_back( gen.addInfo(input, &commonMETdata, globalThreshold, onlyFiducial,applyFiducialThresholdForFractions, usePt) );
     event.put( genmetcoll );
   }

void METProducer::produce_PFClusterMET ( edm::Event & event )

private

Definition at line 214 of file METProducer.cc.

References PFClusterSpecificAlgo::addInfo(), globalThreshold, input, inputToken_, edm::Handle< T >::product(), and METAlgo::run().

Referenced by produce().

   {
     edm::Handle<edm::View<reco::Candidate> > input;
     event.getByToken(inputToken_, input);
 
     METAlgo algo;
     CommonMETData commonMETdata = algo.run(*input.product(), globalThreshold);
 
     PFClusterSpecificAlgo pfcluster;
     std::auto_ptr<reco::PFClusterMETCollection> pfclustermetcoll;
     pfclustermetcoll.reset (new reco::PFClusterMETCollection);
     
     pfclustermetcoll->push_back( pfcluster.addInfo(input, commonMETdata) );
     event.put( pfclustermetcoll );
   }

void METProducer::produce_PFMET ( edm::Event & event )

private

Definition at line 180 of file METProducer.cc.

References metsig::SignPFSpecificAlgo::addPFJets(), calculateSignificance_, globalThreshold, input, inputToken_, fwrapper::jets, jetToken_, CommonMETData::met, CommonMETData::mex, CommonMETData::mey, metsig::SignPFSpecificAlgo::mkSignifMatrix(), p4, edm::Handle< T >::product(), resolutions_, PFSpecificAlgo::run(), METAlgo::run(), metsig::SignPFSpecificAlgo::setResolutions(), reco::MET::setSignificanceMatrix(), timingPdfMaker::specific, and CommonMETData::sumet.

Referenced by produce().

   {
     edm::Handle<edm::View<reco::Candidate> > input;
     event.getByToken(inputToken_, input);
 
     METAlgo algo;
     CommonMETData commonMETdata = algo.run(*input.product(), globalThreshold);
 
     const math::XYZTLorentzVector p4(commonMETdata.mex, commonMETdata.mey, 0.0, commonMETdata.met);
     const math::XYZPoint vtx(0.0, 0.0, 0.0);
 
     PFSpecificAlgo pf;
     SpecificPFMETData specific = pf.run(*input.product());
 
     reco::PFMET pfmet(specific, commonMETdata.sumet, p4, vtx);
 
     if(calculateSignificance_)
       {
     metsig::SignPFSpecificAlgo pfsignalgo;
     pfsignalgo.setResolutions(resolutions_);
 
     edm::Handle<edm::View<reco::PFJet> > jets;
     event.getByToken(jetToken_, jets);
     pfsignalgo.addPFJets(jets.product());
     pfmet.setSignificanceMatrix(pfsignalgo.mkSignifMatrix(input));
       }
 
     std::auto_ptr<reco::PFMETCollection> pfmetcoll;
     pfmetcoll.reset(new reco::PFMETCollection);
 
     pfmetcoll->push_back(pfmet);
     event.put(pfmetcoll);
   }

void METProducer::produce_TCMET	(	edm::Event &	event,
		const edm::EventSetup &	setup
	)

private

Definition at line 172 of file METProducer.cc.

References TCMETAlgo::CalculateTCMET(), and tcMetAlgo_.

Referenced by produce().

   {
     std::auto_ptr<reco::METCollection> tcmetcoll;
     tcmetcoll.reset(new reco::METCollection);
     tcmetcoll->push_back( tcMetAlgo_.CalculateTCMET(event, setup ) ) ;
     event.put( tcmetcoll );
   }