#include <MixCollectionValidation.h>

Inheritance diagram for MixCollectionValidation:

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

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

	MixCollectionValidation (const edm::ParameterSet &)

	~MixCollectionValidation () 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 () 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 Member Functions
template<class T1 , class T2 >
void	fillCaloHitTime (T1 &theItr_, T2 &theColl_, MonitorElement *theProfile_)

template<class T1 , class T2 >
void	fillGenParticleMulti (T1 &theItr_, T2 &theColl_, MonitorElement *theProfile_)

template<class T1 , class T2 >
void	fillMultiplicity (T1 &theItr_, T2 &theColl_, MonitorElement *theProfile_)

template<class T1 , class T2 >
void	fillSimHitTime (T1 &theItr_, T2 &theColl_, MonitorElement *theProfile_)

Private Attributes
std::map< std::string, MonitorElement * >	CaloHitNrmap_

std::map< std::string, MonitorElement * >	CaloHitTimemap_

edm::EDGetTokenT< CrossingFrame< edm::HepMCProduct > >	crossingFrame_Hep_Token_

std::vector< edm::EDGetTokenT< CrossingFrame< PCaloHit > > >	crossingFrame_PCaloHit_Tokens_

std::vector< edm::EDGetTokenT< CrossingFrame< PSimHit > > >	crossingFrame_PSimHit_Tokens_

edm::EDGetTokenT< CrossingFrame< SimTrack > >	crossingFrame_SimTr_Token_

edm::EDGetTokenT< CrossingFrame< SimVertex > >	crossingFrame_SimVtx_Token_

std::vector< edm::InputTag >	HepMCProductTags_

int	maxbunch_

int	minbunch_

edm::ParameterSet	mixObjextsSet_

int	nbin_

MonitorElement *	nrHepMCProductH_

MonitorElement *	nrSimTrackH_

MonitorElement *	nrSimVertexH_

std::string	outputFile_

std::vector< edm::InputTag >	PCaloHitTags_

std::vector< edm::InputTag >	PSimHitTags_

std::map< std::string, MonitorElement * >	SimHitNrmap_

std::map< std::string, MonitorElement * >	SimHitTimemap_

std::vector< edm::InputTag >	SimTrackTags_

std::vector< edm::InputTag >	SimVertexTags_

bool	verbose_

Additional Inherited Members
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

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 34 of file MixCollectionValidation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 22 of file MixCollectionValidation.cc.

References edm::ParameterSet::getParameter(), and mixObjextsSet_.

                                                                               :
   minbunch_(iConfig.getParameter<int>("minBunch")),
   maxbunch_(iConfig.getParameter<int>("maxBunch")),
   verbose_(iConfig.getUntrackedParameter<bool>("verbose",false)),
   nbin_(maxbunch_-minbunch_+1)
 {
   // Histograms will be defined according to the configuration
   ParameterSet mixObjextsSet_ = iConfig.getParameter<ParameterSet>("mixObjects");
 }

MixCollectionValidation::~MixCollectionValidation ( )

override

Definition at line 32 of file MixCollectionValidation.cc.

 {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

void MixCollectionValidation::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iConfig
	)

override

Definition at line 128 of file MixCollectionValidation.cc.

References CaloHitNrmap_, CaloHitTimemap_, crossingFrame_Hep_Token_, crossingFrame_PCaloHit_Tokens_, crossingFrame_PSimHit_Tokens_, crossingFrame_SimTr_Token_, crossingFrame_SimVtx_Token_, fillCaloHitTime(), fillGenParticleMulti(), fillMultiplicity(), fillSimHitTime(), edm::Event::getByToken(), HepMCProductTags_, mps_fire::i, instance, createfilelist::int, nrHepMCProductH_, nrSimTrackH_, nrSimVertexH_, PCaloHitTags_, edm::Handle< T >::product(), PSimHitTags_, SimHitNrmap_, trackerHits::simHits, SimHitTimemap_, tkConvValidator_cfi::simTracks, SimTrackTags_, SimVertexTags_, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   using namespace edm;
 
   if ( !HepMCProductTags_.empty() ) {
     bool gotHepMCProduct;
     edm::Handle<CrossingFrame<HepMCProduct> > crossingFrame;
     gotHepMCProduct = iEvent.getByToken(crossingFrame_Hep_Token_, crossingFrame);
 
     if (gotHepMCProduct){
       std::auto_ptr<MixCollection<HepMCProduct> >
           hepMCProduct (new MixCollection<HepMCProduct>(crossingFrame.product ()));
       MixCollection<HepMCProduct>::MixItr hitItr;
 
       fillGenParticleMulti(hitItr, hepMCProduct, nrHepMCProductH_);
     }
   }
 
   if ( !SimTrackTags_.empty() ) {
     bool gotSimTrack;
     edm::Handle<CrossingFrame<SimTrack> > crossingFrame;
     gotSimTrack = iEvent.getByToken(crossingFrame_SimTr_Token_,crossingFrame);
 
     if (gotSimTrack){
       std::auto_ptr<MixCollection<SimTrack> >
           simTracks (new MixCollection<SimTrack>(crossingFrame.product ()));
       MixCollection<SimTrack>::MixItr hitItr;
 
       fillMultiplicity(hitItr, simTracks, nrSimTrackH_);
     }
   }
 
   if ( !SimVertexTags_.empty() ) {
     bool gotSimVertex;
     edm::Handle<CrossingFrame<SimVertex> > crossingFrame;
     std::string SimVertexLabel = SimVertexTags_[0].label();
     gotSimVertex = iEvent.getByToken(crossingFrame_SimVtx_Token_, crossingFrame);
 
     if (gotSimVertex){
       std::auto_ptr<MixCollection<SimVertex> >
           simVerteces (new MixCollection<SimVertex>(crossingFrame.product ()));
       MixCollection<SimVertex>::MixItr hitItr;
 
       fillMultiplicity(hitItr, simVerteces, nrSimVertexH_);
     }
   }
 
   if ( !PSimHitTags_.empty() ) {
 
     edm::Handle<CrossingFrame<PSimHit> > crossingFrame;
 
     for ( int i = 0; i < (int)PSimHitTags_.size(); i++ ) {
       bool gotPSimHit;
       gotPSimHit = iEvent.getByToken(crossingFrame_PSimHit_Tokens_[i], crossingFrame);
 
       if (gotPSimHit){
         std::auto_ptr<MixCollection<PSimHit> >
             simHits (new MixCollection<PSimHit>(crossingFrame.product ()));
 
         MixCollection<PSimHit>::MixItr hitItr;
 
         fillMultiplicity(hitItr, simHits, SimHitNrmap_[PSimHitTags_[i].instance()]);
 
         fillSimHitTime(hitItr, simHits, SimHitTimemap_[PSimHitTags_[i].instance()]);
       }
     }
   }
 
   if ( !PCaloHitTags_.empty() ) {
 
     edm::Handle<CrossingFrame<PCaloHit> > crossingFrame;
 
     for ( int i = 0; i < (int)PCaloHitTags_.size(); i++ ) {
       bool gotPCaloHit;
       std::string PCaloHitLabel = PCaloHitTags_[i].label()+PCaloHitTags_[i].instance();
       gotPCaloHit = iEvent.getByToken(crossingFrame_PCaloHit_Tokens_[i], crossingFrame);
 
       if (gotPCaloHit){
         std::auto_ptr<MixCollection<PCaloHit> >
             caloHits (new MixCollection<PCaloHit>(crossingFrame.product ()));
 
         MixCollection<PCaloHit>::MixItr hitItr;
 
         fillMultiplicity(hitItr, caloHits, CaloHitNrmap_[PCaloHitTags_[i].instance()]);
 
         fillCaloHitTime(hitItr, caloHits, CaloHitTimemap_[PCaloHitTags_[i].instance()]);
       }
     }
   }
 
 }

void MixCollectionValidation::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	iRun,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 38 of file MixCollectionValidation.cc.

References DQMStore::IBooker::bookProfile(), CaloHitNrmap_, CaloHitTimemap_, edm::EDConsumerBase::consumes(), crossingFrame_Hep_Token_, crossingFrame_PCaloHit_Tokens_, crossingFrame_PSimHit_Tokens_, crossingFrame_SimTr_Token_, crossingFrame_SimVtx_Token_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNames(), HepMCProductTags_, cuy::ii, diffTwoXMLs::label, maxbunch_, minbunch_, mixObjextsSet_, dataset::name, cscdqm::h::names, nbin_, nrHepMCProductH_, nrSimTrackH_, nrSimVertexH_, resolutioncreator_cfi::object, PCaloHitTags_, muonDTDigis_cfi::pset, PSimHitTags_, DQMStore::IBooker::setCurrentFolder(), SimHitNrmap_, SimHitTimemap_, SimTrackTags_, SimVertexTags_, AlCaHLTBitMon_QueryRunRegistry::string, fireworks::subdets, and fftjetcommon_cfi::title.

 {
   iBooker.setCurrentFolder("MixingV/Mixing");
 
   std::vector<std::string> names = mixObjextsSet_.getParameterNames();
 
   for (std::vector<std::string>::iterator it = names.begin();it!= names.end();++it)
   {
     ParameterSet pset = mixObjextsSet_.getParameter<ParameterSet>((*it));
     if (!pset.exists("type"))  continue; //to allow replacement by empty pset
     std::string object = pset.getParameter<std::string>("type");
     std::vector<InputTag> tags = pset.getParameter<std::vector<InputTag> >("input");
 
     if ( object == "HepMCProduct" ) {
 
       std::string title = "Log10 Number of GenParticle in " + object;
       std::string name = "NumberOf" + object;
       nrHepMCProductH_ = iBooker.bookProfile(name,title,nbin_,minbunch_,maxbunch_+1,40,0.,40.);
 
       HepMCProductTags_ = tags;
       if (!HepMCProductTags_.empty()) {
         crossingFrame_Hep_Token_ = consumes<CrossingFrame<HepMCProduct> >(
             edm::InputTag("mix", HepMCProductTags_[0].label()));
       }
     }
     else if ( object == "SimTrack" ) {
 
       std::string title = "Log10 Number of " + object;
       std::string name = "NumberOf" + object;
       nrSimTrackH_ = iBooker.bookProfile(name,title,nbin_,minbunch_,maxbunch_+1,40,0.,40.);
 
       SimTrackTags_ = tags;
       if (!SimTrackTags_.empty()) {
         crossingFrame_SimTr_Token_ = consumes<CrossingFrame<SimTrack> >(
             edm::InputTag("mix", SimTrackTags_[0].label()));
       }
     }
     else if ( object == "SimVertex" )  {
 
       std::string title = "Log10 Number of " + object;
       std::string name = "NumberOf" + object;
       nrSimVertexH_ = iBooker.bookProfile(name,title,nbin_,minbunch_,maxbunch_+1,40,0.,40.);
 
       SimVertexTags_ = tags;
       if (!SimVertexTags_.empty()) {
         crossingFrame_SimVtx_Token_ = consumes<CrossingFrame<SimVertex> >(
             edm::InputTag("mix", SimVertexTags_[0].label()));
       }
     }
     else if ( object == "PSimHit" ) {
       std::vector<std::string> subdets=pset.getParameter<std::vector<std::string> >("subdets");
       for (unsigned int ii=0;ii<subdets.size();ii++) {
 
         std::string title = "Log10 Number of " + subdets[ii];
         std::string name = "NumberOf" + subdets[ii];
         SimHitNrmap_[subdets[ii]] = iBooker.bookProfile(name,title,nbin_,minbunch_,maxbunch_+1,40,0.,40.);
 
         title = "Time of " + subdets[ii];
         name = "TimeOf" + subdets[ii];
         SimHitTimemap_[subdets[ii]] = iBooker.bookProfile(name,title,nbin_,minbunch_,maxbunch_+1,40,-125.,375.);
 
       }
 
       PSimHitTags_ = tags;
       for (auto const & it : PSimHitTags_)
         crossingFrame_PSimHit_Tokens_.push_back(consumes<CrossingFrame<PSimHit> >(
             edm::InputTag("mix", it.label() + it.instance())));
     }
     else if ( object == "PCaloHit" ) {
       std::vector<std::string> subdets=pset.getParameter<std::vector<std::string> >("subdets");
       for (unsigned int ii=0;ii<subdets.size();ii++) {
 
         std::string title = "Log10 Number of " + subdets[ii];
         std::string name = "NumberOf" + subdets[ii];
         CaloHitNrmap_[subdets[ii]] = iBooker.bookProfile(name,title,nbin_,minbunch_,maxbunch_+1,40,0.,40.);
 
         title = "Time of " + subdets[ii];
         name = "TimeOf" + subdets[ii];
         CaloHitTimemap_[subdets[ii]] = iBooker.bookProfile(name,title,nbin_,minbunch_,maxbunch_+1,40,-125.,375.);
 
       }
 
       PCaloHitTags_ = tags;
       for (auto const & it : PCaloHitTags_)
         crossingFrame_PCaloHit_Tokens_.push_back(consumes<CrossingFrame<PCaloHit> >(
             edm::InputTag("mix", it.label() + it.instance())));
     }
   }
 }

template<class T1 , class T2 >

void MixCollectionValidation::fillCaloHitTime	(	T1 &	theItr_,
		T2 &	theColl_,
		MonitorElement *	theProfile_
	)

private

Definition at line 271 of file MixCollectionValidation.cc.

References MonitorElement::Fill(), minbunch_, nbin_, and ntuplemaker::time.

Referenced by analyze().

                                                                                                                                     {
 
   for ( theItr_ = theColl_->begin() ; theItr_ != theColl_->end() ; ++theItr_) {
 
     int bunch = (*theItr_).eventId().bunchCrossing();
     float time = (*theItr_).time();
     int index = bunch - minbunch_;
     if ( index >= 0 && index < nbin_ ) { theProfile_->Fill(float(bunch+0.5),time); }
     else { edm::LogWarning("MixCollectionValidation") << "fillCaloHitTime: bunch number " << bunch << " out of range"; }
 
   }
 
 }

template<class T1 , class T2 >

void MixCollectionValidation::fillGenParticleMulti	(	T1 &	theItr_,
		T2 &	theColl_,
		MonitorElement *	theProfile_
	)

private

Definition at line 239 of file MixCollectionValidation.cc.

References MonitorElement::Fill(), mps_fire::i, SiStripPI::max, minbunch_, and nbin_.

Referenced by analyze().

                                                                                                                                          {
 
   std::vector<int> theMult(nbin_);
 
   for ( theItr_ = theColl_->begin() ; theItr_ != theColl_->end() ; ++theItr_) {
 
     int bunch = theItr_.bunch();
     int index = bunch - minbunch_;
     if ( index >= 0 && index < nbin_ ) { theMult[index] += (*theItr_).GetEvent()->particles_size(); }
     else { edm::LogWarning("MixCollectionValidation") << "fillMultiplicity: bunch number " << bunch << " out of range"; }
 
   }
 
   for ( int i = 0; i < nbin_; i++ ) {
     theProfile_->Fill(float(i+minbunch_+0.5),std::log10(std::max(float(0.1),float(theMult[i]))));
   }
 }

template<class T1 , class T2 >

void MixCollectionValidation::fillMultiplicity	(	T1 &	theItr_,
		T2 &	theColl_,
		MonitorElement *	theProfile_
	)

private

Definition at line 220 of file MixCollectionValidation.cc.

References MonitorElement::Fill(), mps_fire::i, SiStripPI::max, minbunch_, and nbin_.

Referenced by analyze().

                                                                                                                                      {
 
   std::vector<int> theMult(nbin_);
 
   for ( theItr_ = theColl_->begin() ; theItr_ != theColl_->end() ; ++theItr_) {
 
     int bunch = (*theItr_).eventId().bunchCrossing();
     int index = bunch - minbunch_;
     if ( index >= 0 && index < nbin_ ) { theMult[index] += 1; }
     else { edm::LogWarning("MixCollectionValidation") << "fillMultiplicity: bunch number " << bunch << " out of range"; }
 
   }
 
   for ( int i = 0; i < nbin_; i++ ) {
     theProfile_->Fill(float(i+minbunch_+0.5),std::log10(std::max(float(0.1),float(theMult[i]))));
   }
 }

template<class T1 , class T2 >

void MixCollectionValidation::fillSimHitTime	(	T1 &	theItr_,
		T2 &	theColl_,
		MonitorElement *	theProfile_
	)

private

Definition at line 257 of file MixCollectionValidation.cc.

References MonitorElement::Fill(), minbunch_, nbin_, and ntuplemaker::time.

Referenced by analyze().

                                                                                                                                    {
 
   for ( theItr_ = theColl_->begin() ; theItr_ != theColl_->end() ; ++theItr_) {
 
     int bunch = (*theItr_).eventId().bunchCrossing();
     float time = (*theItr_).timeOfFlight();
     int index = bunch - minbunch_;
     if ( index >= 0 && index < nbin_ ) { theProfile_->Fill(float(bunch+0.5),time); }
     else { edm::LogWarning("MixCollectionValidation") << "fillSimHitTime: bunch number " << bunch << " out of range"; }
 
   }
 
 }