#include <BDHadronTrackMonitoringHarvester.h>

Inheritance diagram for BDHadronTrackMonitoringHarvester:

Public Member Functions
	BDHadronTrackMonitoringHarvester (const edm::ParameterSet &pSet)

	~BDHadronTrackMonitoringHarvester ()

Public Member Functions inherited from DQMEDHarvester
virtual void	analyze (edm::Event const &, edm::EventSetup const &) final

virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

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

	DQMEDHarvester (void)

virtual void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	endJob () final

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

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

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginJob ()

void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override

Private Attributes
MonitorElement *	nTrk_absolute_bjet

MonitorElement *	nTrk_absolute_cjet

MonitorElement *	nTrk_absolute_dusgjet

MonitorElement *	nTrk_relative_bjet

MonitorElement *	nTrk_relative_cjet

MonitorElement *	nTrk_relative_dusgjet

MonitorElement *	nTrk_std_bjet

MonitorElement *	nTrk_std_cjet

MonitorElement *	nTrk_std_dusgjet

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Detailed Description

Top level steering routine for b tag performance harvesting.

Definition at line 24 of file BDHadronTrackMonitoringHarvester.h.

Constructor & Destructor Documentation

BDHadronTrackMonitoringHarvester::BDHadronTrackMonitoringHarvester ( const edm::ParameterSet & pSet )

explicit

Definition at line 9 of file BDHadronTrackMonitoringHarvester.cc.

10 {

11 }

BDHadronTrackMonitoringHarvester::~BDHadronTrackMonitoringHarvester ( )

Definition at line 14 of file BDHadronTrackMonitoringHarvester.cc.

15 {

16 }

Member Function Documentation

void BDHadronTrackMonitoringHarvester::beginJob ( void )

privatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 18 of file BDHadronTrackMonitoringHarvester.cc.

19 {

20 }

void BDHadronTrackMonitoringHarvester::dqmEndJob	(	DQMStore::IBooker &	ibook,
		DQMStore::IGetter &	iget
	)

overrideprivatevirtual

Implements DQMEDHarvester.

Definition at line 22 of file BDHadronTrackMonitoringHarvester.cc.

References DQMStore::IBooker::book1D(), DEFINE_FWK_MODULE, DQMStore::IGetter::get(), MonitorElement::getMean(), MonitorElement::getRMS(), i, nTrk_absolute_bjet, nTrk_absolute_cjet, nTrk_absolute_dusgjet, nTrk_relative_bjet, nTrk_relative_cjet, nTrk_relative_dusgjet, nTrk_std_bjet, nTrk_std_cjet, nTrk_std_dusgjet, DQMStore::IGetter::removeElement(), MonitorElement::setAxisRange(), MonitorElement::setAxisTitle(), MonitorElement::setBinContent(), MonitorElement::setBinLabel(), RecoBTag::setTDRStyle(), and BDHadronTrackMonitoringAnalyzer::TrkHistCat.

 {
     // ***********************
     //
     // Book all new histograms.
     //
     // ***********************
     RecoBTag::setTDRStyle();
 
     // b jets
     // absolute average number of tracks
     nTrk_absolute_bjet = ibook.book1D("nTrk_absolute_bjet","absolute average number of tracks in b jets",6,-0.5,5.5);
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_absolute_bjet->setBinLabel(i+1,BDHadronTrackMonitoringAnalyzer::TrkHistCat[i],1);
     }
     nTrk_absolute_bjet->setAxisRange(0, 5, 2);
     nTrk_absolute_bjet->setAxisTitle("average number of tracks",2);
 
     // relative (in percent) average number of tracks
     nTrk_relative_bjet = ibook.book1D("nTrk_relative_bjet","relative average number of tracks in b jets",6,-0.5,5.5);
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_relative_bjet->setBinLabel(i+1,BDHadronTrackMonitoringAnalyzer::TrkHistCat[i],1);
     }
     nTrk_relative_bjet->setAxisRange(0, 1, 2);
     nTrk_relative_bjet->setAxisTitle("average fraction of tracks",2);
 
     // standard deviation of number of tracks
     nTrk_std_bjet = ibook.book1D("nTrk_std_bjet","RMS of number of tracks in b jets",6,-0.5,5.5);
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_std_bjet->setBinLabel(i+1,BDHadronTrackMonitoringAnalyzer::TrkHistCat[i],1);
     }
     nTrk_std_bjet->setAxisRange(0, 3, 2);
     nTrk_std_bjet->setAxisTitle("RMS of number of tracks",2);
 
 
 
 
     // c jets
     nTrk_absolute_cjet = ibook.book1D("nTrk_absolute_cjet","absolute average number of tracks in c jets",6,-0.5,5.5);
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_absolute_cjet->setBinLabel(i+1,BDHadronTrackMonitoringAnalyzer::TrkHistCat[i],1);
     }
     nTrk_absolute_cjet->setAxisRange(0, 5, 2);
     nTrk_absolute_cjet->setAxisTitle("average number of tracks",2);
 
     nTrk_relative_cjet = ibook.book1D("nTrk_relative_cjet","relative average number of tracks in c jets",6,-0.5,5.5);
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_relative_cjet->setBinLabel(i+1,BDHadronTrackMonitoringAnalyzer::TrkHistCat[i],1);
     }
     nTrk_relative_cjet->setAxisRange(0, 1, 2);
     nTrk_relative_cjet->setAxisTitle("average fraction of tracks",2);
 
     nTrk_std_cjet = ibook.book1D("nTrk_std_cjet","RMS of number of tracks in c jets",6,-0.5,5.5);
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_std_cjet->setBinLabel(i+1,BDHadronTrackMonitoringAnalyzer::TrkHistCat[i],1);
     }
     nTrk_std_cjet->setAxisRange(0, 3, 2);
     nTrk_std_cjet->setAxisTitle("RMS of number of tracks",2);
 
 
 
     // udsg jets
     nTrk_absolute_dusgjet = ibook.book1D("nTrk_absolute_dusgjet","absolute average number of tracks in dusg jets",6,-0.5,5.5);
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_absolute_dusgjet->setBinLabel(i+1,BDHadronTrackMonitoringAnalyzer::TrkHistCat[i],1);
     }
     nTrk_absolute_dusgjet->setAxisRange(0, 5, 2);
     nTrk_absolute_dusgjet->setAxisTitle("average number of tracks",2);
 
     nTrk_relative_dusgjet = ibook.book1D("nTrk_relative_dusgjet","relative average number of tracks in dusg jets",6,-0.5,5.5);
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_relative_dusgjet->setBinLabel(i+1,BDHadronTrackMonitoringAnalyzer::TrkHistCat[i],1);
     }
     nTrk_relative_dusgjet->setAxisRange(0, 1, 2);
     nTrk_relative_dusgjet->setAxisTitle("average fraction of tracks",2);
 
     nTrk_std_dusgjet = ibook.book1D("nTrk_std_dusgjet","RMS of number of tracks in dusg jets",6,-0.5,5.5);
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_std_dusgjet->setBinLabel(i+1,BDHadronTrackMonitoringAnalyzer::TrkHistCat[i],1);
     }
     nTrk_std_dusgjet->setAxisRange(0, 3, 2);
     nTrk_std_dusgjet->setAxisTitle("RMS of number of tracks",2);
 
 
     // ***********************
     //
     // get all the old histograms
     //
     // ***********************
 
     // b jets
     MonitorElement *nTrk_bjet[6];
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_bjet[i] = iget.get("nTrk_bjet_"+BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
     }
     MonitorElement *nTrkAll_bjet = iget.get("nTrkAll_bjet");
 
     // c jets
     MonitorElement *nTrk_cjet[6];
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_cjet[i] = iget.get("nTrk_cjet_"+BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
     }
     MonitorElement *nTrkAll_cjet = iget.get("nTrkAll_cjet");
 
     // dusg jets
     MonitorElement *nTrk_dusgjet[6];
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         nTrk_dusgjet[i] = iget.get("nTrk_dusgjet_"+BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
     }
     MonitorElement *nTrkAll_dusgjet = iget.get("nTrkAll_dusgjet");
 
 
 
     // ***********************
     //
     // Calculate contents of new histograms
     //
     // ***********************
 
     // b jets
     float mean_bjets[6];
     float std_bjets[6];
     float meanAll_bjets;
     meanAll_bjets = nTrkAll_bjet->getMean(1);
     for (unsigned int i = 0; i< BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++){
         mean_bjets[i] = nTrk_bjet[i]->getMean(1); // mean number of tracks per category
         std_bjets[i] = nTrk_bjet[i]->getRMS(1);
         nTrk_absolute_bjet->setBinContent(i+1,mean_bjets[i]);
         nTrk_relative_bjet->setBinContent(i+1,mean_bjets[i]/meanAll_bjets);
         nTrk_std_bjet->setBinContent(i+1,std_bjets[i]);
     }
 
     // c jets
     float mean_cjets[6];
     float std_cjets[6];
     float meanAll_cjets;
     meanAll_cjets = nTrkAll_cjet->getMean(1);
     for (unsigned int i = 0; i< BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++){
         mean_cjets[i] = nTrk_cjet[i]->getMean(1); // mean number of tracks per category
         std_cjets[i] = nTrk_cjet[i]->getRMS(1);
         nTrk_absolute_cjet->setBinContent(i+1,mean_cjets[i]);
         nTrk_relative_cjet->setBinContent(i+1,mean_cjets[i]/meanAll_cjets);
         nTrk_std_cjet->setBinContent(i+1,std_cjets[i]);
     }
 
 
     // dusg jets
     float mean_dusgjets[6];
     float std_dusgjets[6];
     float meanAll_dusgjets;
     meanAll_dusgjets = nTrkAll_dusgjet->getMean(1);
     for (unsigned int i = 0; i< BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++){
         mean_dusgjets[i] = nTrk_dusgjet[i]->getMean(1); // mean number of tracks per category
         std_dusgjets[i] = nTrk_dusgjet[i]->getRMS(1);
         nTrk_absolute_dusgjet->setBinContent(i+1,mean_dusgjets[i]);
         nTrk_relative_dusgjet->setBinContent(i+1,mean_dusgjets[i]/meanAll_dusgjets);
         nTrk_std_dusgjet->setBinContent(i+1,std_dusgjets[i]);
     }
 
 
     // ***********************
     //
     // Remove histograms with nTracks that we do not want to show
     //
     // ***********************
 
     // b jets
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         iget.removeElement("nTrk_bjet_"+BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
     }
     iget.removeElement("nTrkAll_bjet");
 
     // c jets
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         iget.removeElement("nTrk_cjet_"+BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
     }
     iget.removeElement("nTrkAll_cjet");
 
     // dusg jets
     for(unsigned int i = 0; i < BDHadronTrackMonitoringAnalyzer::TrkHistCat.size(); i++) {
         iget.removeElement("nTrk_dusgjet_"+BDHadronTrackMonitoringAnalyzer::TrkHistCat[i]);
     }
     iget.removeElement("nTrkAll_dusgjet");
 
 
 
 }