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

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

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

	~DQMEDAnalyzer () override=default

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_

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

override

Definition at line 38 of file MixCollectionValidation.cc.

References DQMStore::IBooker::bookProfile(), CaloHitNrmap_, CaloHitTimemap_, 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, 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"; }
 
   }
 
 }