#include <DuplicationChecker.h>

Inheritance diagram for DuplicationChecker:

Public Types
typedef std::multimap< double, edm::EventID >	associationMap

typedef std::vector< associationMap::iterator >	itemList

Public Types inherited from edm::one::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

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override

	DuplicationChecker (const edm::ParameterSet &)

void	findValuesAssociatedWithKey (associationMap &mMap, double &key, itemList &theObjects)

	~DuplicationChecker () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

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

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

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

	~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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::InputTag	generatedCollection_

edm::EDGetTokenT< edm::HepMCProduct >	generatedCollectionToken_

edm::InputTag	lheEventProduct_

edm::EDGetTokenT< LHEEventProduct >	lheEventProductToken_

bool	searchForLHE_

WeightManager	wmanager_

associationMap	xBjorkenHistory

MonitorElement *	xBjorkenME

Additional Inherited Members
Static Public Member Functions inherited from edm::one::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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 39 of file DuplicationChecker.h.

Member Typedef Documentation

typedef std::multimap<double,edm::EventID> DuplicationChecker::associationMap

Definition at line 43 of file DuplicationChecker.h.

typedef std::vector<associationMap::iterator> DuplicationChecker::itemList

Definition at line 44 of file DuplicationChecker.h.

Constructor & Destructor Documentation

DuplicationChecker::DuplicationChecker ( const edm::ParameterSet & iPSet )

explicit

Definition at line 12 of file DuplicationChecker.cc.

References generatedCollection_, generatedCollectionToken_, edm::ParameterSet::getParameter(), lheEventProduct_, lheEventProductToken_, searchForLHE_, and xBjorkenHistory.

                                                                   :
   wmanager_(iPSet,consumesCollector()),
   generatedCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection")),
   searchForLHE_(iPSet.getParameter<bool>("searchForLHE"))
 { 
   if (searchForLHE_) {
     lheEventProduct_ = iPSet.getParameter<edm::InputTag>("lheEventProduct");
   }
   xBjorkenHistory.clear();
 
   if (searchForLHE_) lheEventProductToken_=consumes<LHEEventProduct>(lheEventProduct_);
   else generatedCollectionToken_=consumes<HepMCProduct>(generatedCollection_);
 
 }

DuplicationChecker::~DuplicationChecker ( )

override

Definition at line 27 of file DuplicationChecker.cc.

References xBjorkenHistory.

 {
   xBjorkenHistory.clear();
 }

Member Function Documentation

void DuplicationChecker::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 40 of file DuplicationChecker.cc.

References MonitorElement::Fill(), generatedCollectionToken_, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), LHEEventProduct::hepeup(), hepeup_, edm::EventBase::id(), lheEventProductToken_, lhef::HEPEUP::PUP, searchForLHE_, WeightManager::weight(), wmanager_, xBjorkenHistory, and xBjorkenME.

 {
     
   double bjorken = 0;
  
   double weight = 1.;
 
   if (searchForLHE_) {
 
     Handle<LHEEventProduct> evt;
     iEvent.getByToken(lheEventProductToken_, evt);
 
     const lhef::HEPEUP hepeup_ = evt->hepeup();
 
     const std::vector<lhef::HEPEUP::FiveVector> pup_ = hepeup_.PUP;
 
     double pz1=(pup_[0])[3];
     double pz2=(pup_[1])[3];
     bjorken+=(pz1/(pz1+pz2));
   }
   else {
     //change teh weight in this case
     weight = wmanager_.weight(iEvent);
 
     edm::Handle<HepMCProduct> evt;
     iEvent.getByToken(generatedCollectionToken_, evt);
 
     const HepMC::PdfInfo *pdf = evt->GetEvent()->pdf_info();    
     if(pdf){
       bjorken = ((pdf->x1())/((pdf->x1())+(pdf->x2())));
     }
 
   }
 
   xBjorkenHistory.insert(std::pair<double,edm::EventID>(bjorken,iEvent.id()));
 
   xBjorkenME->Fill(bjorken,weight);
 
 }//analyze

void DuplicationChecker::bookHistograms	(	DQMStore::IBooker &	i,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Setting the DQM top directories

Booking the ME's

Implements DQMEDAnalyzer.

Definition at line 32 of file DuplicationChecker.cc.

References DQMHelper::book1dHisto(), DQMStore::IBooker::setCurrentFolder(), and xBjorkenME.

                                                                                                 {
   DQMHelper dqm(&i); i.setCurrentFolder("Generator/DuplicationCheck");
   
   xBjorkenME = dqm.book1dHisto("xBjorkenME", "x Bjorken ratio", 1000000, 0., 1.);
 }

void DuplicationChecker::findValuesAssociatedWithKey	(	associationMap &	mMap,
		double &	key,
		itemList &	theObjects
	)

Definition at line 80 of file DuplicationChecker.cc.

 {
   associationMap::iterator itr;
   associationMap::iterator lastElement;
         
   theObjects.clear();
 
   // locate an iterator to the first pair object associated with key
   itr = mMap.find(key);
   if (itr == mMap.end())
     return; // no elements associated with key, so return immediately
 
   // get an iterator to the element that is one past the last element associated with key
   lastElement = mMap.upper_bound(key);
 
   // for each element in the sequence [itr, lastElement)
   for ( ; itr != lastElement; ++itr)
     theObjects.push_back(itr);
 }  