#include <RemovePileUpDominatedEventsGen/RemovePileUpDominatedEventsGen/plugins/RemovePileUpDominatedEventsGen.cc>

Inheritance diagram for RemovePileUpDominatedEventsGen:

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

	~RemovePileUpDominatedEventsGen () override

Public Member Functions inherited from edm::stream::EDFilter<>
	EDFilter ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Public Member Functions inherited from edm::stream::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

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

	~ProducerBase () noexcept(false) override

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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
bool	filter (edm::Event &, const edm::EventSetup &) override

Private Attributes
unsigned int	bunchCrossing

const edm::EDGetTokenT< GenEventInfoProduct >	generatorInfo_

const edm::EDGetTokenT< std::vector< PileupSummaryInfo > >	pileupSummaryInfos_

Additional Inherited Members
Public Types inherited from edm::stream::EDFilter<>
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::EDFilterBase
typedef EDFilterAdaptorBase	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

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

Description: [one line class summary]

This EDFilter select events having the generator pt-hat greater than the pt-had of the pile-up collisions. This code is used by STEAM-TSG in order to estimate the HLT rate using MC sample, especially to avoid a rate double counting due to pile-up. For more information on the theory of this method see Appendix B of https://www.dropbox.com/home/TesiPhD_Silvio?preview=PhD_thesis.pdf .

Implementation: [Notes on implementation]

Definition at line 34 of file RemovePileUpDominatedEventsGen.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 48 of file RemovePileUpDominatedEventsGen.cc.

References bunchCrossing, and ~RemovePileUpDominatedEventsGen().

     : pileupSummaryInfos_(
           consumes<std::vector<PileupSummaryInfo> >(iConfig.getParameter<edm::InputTag>("pileupSummaryInfos"))),
       generatorInfo_(consumes<GenEventInfoProduct>(iConfig.getParameter<edm::InputTag>("generatorInfo"))) {
   bunchCrossing = 0;
   produces<float>();
 }

RemovePileUpDominatedEventsGen::~RemovePileUpDominatedEventsGen ( )

overridedefault

Referenced by RemovePileUpDominatedEventsGen().

Member Function Documentation

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

static

Definition at line 110 of file RemovePileUpDominatedEventsGen.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addDefault(), DEFINE_FWK_MODULE, and HLT_2018_cff::InputTag.

                                                                                                 {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("pileupSummaryInfos", edm::InputTag("addPileupInfo"));
   desc.add<edm::InputTag>("generatorInfo", edm::InputTag("generator"));
   descriptions.addDefault(desc);
 }

bool RemovePileUpDominatedEventsGen::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 58 of file RemovePileUpDominatedEventsGen.cc.

References bunchCrossing, newFWLiteAna::found, generatorInfo_, edm::Event::getByToken(), PileupSummaryInfo::getPU_pT_hats(), eostools::move(), pileupSummaryInfos_, edm::Handle< T >::product(), edm::Event::put(), and GenEventInfoProduct::qScale().

                                                                                          {
   using namespace edm;
   using namespace std;
 
   edm::Handle<GenEventInfoProduct> generatorInfo;
   iEvent.getByToken(generatorInfo_, generatorInfo);
 
   edm::Handle<std::vector<PileupSummaryInfo> > pileupSummaryInfos;
   iEvent.getByToken(pileupSummaryInfos_, pileupSummaryInfos);
 
   //find in-time pile-up
   if (bunchCrossing >= pileupSummaryInfos.product()->size() ||
       pileupSummaryInfos.product()->at(bunchCrossing).getBunchCrossing() != 0) {
     bool found = false;
     for (bunchCrossing = 0; bunchCrossing < pileupSummaryInfos.product()->size() && !found; ++bunchCrossing) {
       if (pileupSummaryInfos.product()->at(bunchCrossing).getBunchCrossing() == 0) {
         found = true;
         bunchCrossing--;
       }
     }
     if (!found) {
       edm::LogInfo("RemovePileUpDominatedEventsGen") << "In-time pile-up not found!" << endl;
       return true;
     }
   }
 
   //   cout << "Using "<<bunchCrossing<<endl;
   //   cout << "pileupSummaryInfos.product()->at(bunchCrossing).getBunchCrossing() "<<pileupSummaryInfos.product()->at(bunchCrossing).getBunchCrossing()<<endl;
 
   //get the PU pt-hat max
   float signal_pT_hat = -1;
   float pu_pT_hat_max = -1;
 
   PileupSummaryInfo puSummary_onTime = pileupSummaryInfos.product()->at(bunchCrossing);
   for (const auto& pu_pT_hat : puSummary_onTime.getPU_pT_hats())
     if (pu_pT_hat > pu_pT_hat_max)
       pu_pT_hat_max = pu_pT_hat;
 
   //get the signal pt-hat
   signal_pT_hat = generatorInfo->qScale();
 
   //save PU - signal pt-hat
   std::unique_ptr<float> pOut(new float());
   *pOut = signal_pT_hat - pu_pT_hat_max;
   iEvent.put(std::move(pOut));
 
   //filter the event
   if (signal_pT_hat > pu_pT_hat_max)
     return true;
   return false;
 }