Slimmer of PAT METs. More...

Inheritance diagram for pat::PATMETSlimmer:

Classes
class	OneMETShift

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

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

	~PATMETSlimmer () override

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

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

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	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	maybeReadShifts (const edm::ParameterSet &basePSet, const std::string &name, pat::MET::METCorrectionType level, bool readFromMiniAOD=false)

Private Attributes
const bool	onMiniAOD_

std::vector< OneMETShift >	shifts_

const edm::EDGetTokenT< pat::METCollection >	src_

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::global::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

Slimmer of PAT METs.

Definition at line 19 of file PATMETSlimmer.cc.

Constructor & Destructor Documentation

pat::PATMETSlimmer::PATMETSlimmer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 53 of file PATMETSlimmer.cc.

References pat::MET::Calo, pat::MET::Chs, maybeReadShifts(), pat::MET::None, onMiniAOD_, pat::MET::Smear, pat::MET::T0, pat::MET::T1, pat::MET::Trk, pat::MET::TXY, pat::MET::TXYForRaw, pat::MET::TXYForT01, pat::MET::TXYForT01Smear, and pat::MET::TXYForT1Smear.

                                                                :
   src_(consumes<pat::METCollection>(iConfig.getParameter<edm::InputTag>("src"))),
   onMiniAOD_(iConfig.existsAs<bool>("runningOnMiniAOD") ?  iConfig.getParameter<bool>("runningOnMiniAOD") : false)
 {
   maybeReadShifts( iConfig, "rawVariation", pat::MET::None );
   maybeReadShifts( iConfig, "t1Uncertainties", pat::MET::T1 );
   maybeReadShifts( iConfig, "t01Variation", pat::MET::T0, onMiniAOD_ );
   maybeReadShifts( iConfig, "t1SmearedVarsAndUncs", pat::MET::Smear );
   
   maybeReadShifts( iConfig, "tXYUncForRaw", pat::MET::TXYForRaw );
   maybeReadShifts( iConfig, "tXYUncForT1", pat::MET::TXY );
   maybeReadShifts( iConfig, "tXYUncForT01", pat::MET::TXYForT01 ); 
   maybeReadShifts( iConfig, "tXYUncForT1Smear", pat::MET::TXYForT1Smear );
   maybeReadShifts( iConfig, "tXYUncForT01Smear", pat::MET::TXYForT01Smear );
   maybeReadShifts( iConfig, "caloMET", pat::MET::Calo );
   maybeReadShifts( iConfig, "chsMET", pat::MET::Chs );
   maybeReadShifts( iConfig, "trkMET", pat::MET::Trk );
 
   produces<std::vector<pat::MET> >();
 }

pat::PATMETSlimmer::~PATMETSlimmer ( )

inlineoverride

Definition at line 22 of file PATMETSlimmer.cc.

References iEvent, and produce().

22 { }

Member Function Documentation

void pat::PATMETSlimmer::maybeReadShifts	(	const edm::ParameterSet &	basePSet,
		const std::string &	name,
		pat::MET::METCorrectionType	level,
		bool	readFromMiniAOD = `false`
	)

private

Definition at line 74 of file PATMETSlimmer.cc.

References edm::EDConsumerBase::consumesCollector(), pat::MET::ElectronEnDown, pat::MET::ElectronEnUp, Exception, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), pat::MET::JetEnDown, pat::MET::JetEnUp, pat::MET::JetResDown, pat::MET::JetResUp, pat::MET::MuonEnDown, pat::MET::MuonEnUp, dataset::name, pat::MET::NoShift, pat::MET::PhotonEnDown, pat::MET::PhotonEnUp, shifts_, pat::MET::Smear, pat::MET::T1, pat::MET::TauEnDown, pat::MET::TauEnUp, pat::MET::UnclusteredEnDown, and pat::MET::UnclusteredEnUp.

Referenced by PATMETSlimmer().

                                                                                                                                                   {
  
     if (basePSet.existsAs<edm::ParameterSet>(name)) {
         throw cms::Exception("Unsupported", "Reading PSets not supported, for now just use input tag");
     } else if (basePSet.existsAs<edm::InputTag>(name) ) {
         const edm::InputTag & baseTag = basePSet.getParameter<edm::InputTag>(name);
 
     if(level==pat::MET::T1) {
       shifts_.push_back(OneMETShift(pat::MET::NoShift, level, baseTag, consumesCollector(), readFromMiniAOD, true, false, false));
       shifts_.push_back(OneMETShift(pat::MET::NoShift, level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::JetResUp,   level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::JetResDown, level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::JetEnUp,   level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::JetEnDown, level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::MuonEnUp,   level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::MuonEnDown, level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::ElectronEnUp,   level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::ElectronEnDown, level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::PhotonEnUp,   level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::PhotonEnDown, level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::TauEnUp,   level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::TauEnDown, level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::UnclusteredEnUp,   level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
       shifts_.push_back(OneMETShift(pat::MET::UnclusteredEnDown, level, baseTag, consumesCollector(), readFromMiniAOD, false, true));
     }
     else if(level==pat::MET::Smear) {
       shifts_.push_back(OneMETShift(pat::MET::NoShift, level, baseTag, consumesCollector(), readFromMiniAOD, true, false, true));
       shifts_.push_back(OneMETShift(pat::MET::JetResUp,   level, baseTag, consumesCollector(), readFromMiniAOD, false, true, true));
       shifts_.push_back(OneMETShift(pat::MET::JetResDown, level, baseTag, consumesCollector(), readFromMiniAOD, false, true, true));
     }
     else {
       shifts_.push_back(OneMETShift(pat::MET::NoShift, level, baseTag, consumesCollector(), readFromMiniAOD, true, false));
     }
     }
 
 }

void pat::PATMETSlimmer::produce	(	edm::StreamID	,
		edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		const

override

Definition at line 141 of file PATMETSlimmer.cc.

References edm::Event::getByToken(), RazorAnalyzer::met, eostools::move(), MillePedeFileConverter_cfg::out, edm::Event::put(), pat::PATMETSlimmer::OneMETShift::shift, shifts_, TrackRefitter_38T_cff::src, and src_.

Referenced by ~PATMETSlimmer().

                                                                                             {
     using namespace edm;
     using namespace std;
 
     Handle<pat::METCollection>  src;
     iEvent.getByToken(src_, src);
     if (src->size() != 1) throw cms::Exception("CorruptData", "More than one MET in the collection");
 
     auto out = std::make_unique<std::vector<pat::MET>>(1, src->front());
     pat::MET & met = out->back();
 
     for (const OneMETShift &shift : shifts_) {
         shift.readAndSet(iEvent, met);
     }
 
     iEvent.put(std::move(out));
 
 }