#include <DuplicationChecker.h>

Inheritance diagram for DuplicationChecker:

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

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

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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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

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

override

Setting the DQM top directories

Booking the ME's

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