#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::EDAnalyzer
typedef EDAnalyzer	ModuleType

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

void	beginJob ()

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

	DuplicationChecker (const edm::ParameterSet &)

void	endJob ()

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

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

virtual	~DuplicationChecker ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

Additional Inherited Members
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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 38 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 dbe, generatedCollection_, generatedCollectionToken_, edm::ParameterSet::getParameter(), lheEventProduct_, lheEventProductToken_, cppFunctionSkipper::operator, 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");
   }
   dbe = 0;
   dbe = edm::Service<DQMStore>().operator->();
 
   xBjorkenHistory.clear();
 
   if (searchForLHE_) lheEventProductToken_=consumes<LHEEventProduct>(lheEventProduct_);
   else generatedCollectionToken_=consumes<HepMCProduct>(generatedCollection_);
 
 }

DuplicationChecker::~DuplicationChecker ( )

virtual

Definition at line 30 of file DuplicationChecker.cc.

References xBjorkenHistory.

 {
   xBjorkenHistory.clear();
 }

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 46 of file DuplicationChecker.cc.

References MonitorElement::Fill(), generatedCollectionToken_, edm::Event::getByToken(), hepeup_, edm::EventBase::id(), lheEventProductToken_, lhef::HEPEUP::PUP, searchForLHE_, WeightManager::weight(), histoStyle::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::beginJob ( void )

virtual

Setting the DQM top directories

Booking the ME's

Reimplemented from edm::EDAnalyzer.

Definition at line 35 of file DuplicationChecker.cc.

References DQMStore::book1D(), dbe, DQMStore::setCurrentFolder(), and xBjorkenME.

 {
   if(dbe){
     dbe->setCurrentFolder("Generator/DuplicationCheck");
     
     xBjorkenME = dbe->book1D("xBjorkenME", "x Bjorken ratio", 1000000, 0., 1.);
   }
 }

virtual void DuplicationChecker::beginRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

inlinevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 51 of file DuplicationChecker.h.

51 { return;}

void DuplicationChecker::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 106 of file DuplicationChecker.cc.

References findValuesAssociatedWithKey(), i, edm::second(), and xBjorkenHistory.

 {
 
   itemList theObjects;
   theObjects.reserve(10);
 
   for (associationMap::iterator it = xBjorkenHistory.begin(); it != xBjorkenHistory.end(); it++) {
     double theKey = (*it).first;
 
     findValuesAssociatedWithKey(xBjorkenHistory, theKey, theObjects);
 
     if (theObjects.size() > 1) {
       edm::LogWarning("DuplicatedEventFound") << "Duplicated events found with xBjorken = " << std::fixed << std::setw(16) << std::setprecision(14) << theKey; 
       for (unsigned int i = 0; i < theObjects.size(); i++) {
         edm::LogPrint("DuplicatedEventList") << "Event = " << (*theObjects[i]).second;
       }
     }
 
     theObjects.clear();
  
   }
 
 }

virtual void DuplicationChecker::endRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

inlinevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 52 of file DuplicationChecker.h.

52 { return;}

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

Definition at line 86 of file DuplicationChecker.cc.

Referenced by endJob().

 {
   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);
 }  

Member Data Documentation

DQMStore* DuplicationChecker::dbe

private

Definition at line 65 of file DuplicationChecker.h.

Referenced by beginJob(), and DuplicationChecker().

edm::InputTag DuplicationChecker::generatedCollection_

private

Definition at line 58 of file DuplicationChecker.h.

Referenced by DuplicationChecker().

edm::EDGetTokenT<edm::HepMCProduct> DuplicationChecker::generatedCollectionToken_

private

Definition at line 70 of file DuplicationChecker.h.

Referenced by analyze(), and DuplicationChecker().

edm::InputTag DuplicationChecker::lheEventProduct_

private

Definition at line 59 of file DuplicationChecker.h.

Referenced by DuplicationChecker().

edm::EDGetTokenT<LHEEventProduct> DuplicationChecker::lheEventProductToken_

private

Definition at line 69 of file DuplicationChecker.h.

Referenced by analyze(), and DuplicationChecker().

bool DuplicationChecker::searchForLHE_

private

Definition at line 61 of file DuplicationChecker.h.

Referenced by analyze(), and DuplicationChecker().

WeightManager DuplicationChecker::wmanager_

private

Definition at line 56 of file DuplicationChecker.h.

Referenced by analyze().

associationMap DuplicationChecker::xBjorkenHistory

private

Definition at line 63 of file DuplicationChecker.h.

Referenced by analyze(), DuplicationChecker(), endJob(), and ~DuplicationChecker().

MonitorElement* DuplicationChecker::xBjorkenME

private

Definition at line 67 of file DuplicationChecker.h.

Referenced by analyze(), and beginJob().

Private Attributes
DQMStore *	dbe

edm::InputTag	generatedCollection_

edm::EDGetTokenT < edm::HepMCProduct >	generatedCollectionToken_

edm::InputTag	lheEventProduct_

edm::EDGetTokenT< LHEEventProduct >	lheEventProductToken_

bool	searchForLHE_

WeightManager	wmanager_

associationMap	xBjorkenHistory

MonitorElement *	xBjorkenME

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation