#include <SiPixelPhase1Summary.h>

Inheritance diagram for SiPixelPhase1Summary:

Public Member Functions
	SiPixelPhase1Summary (const edm::ParameterSet &conf)

	~SiPixelPhase1Summary () override

Public Member Functions inherited from DQMEDHarvester
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

void	beginJob () override

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

void	beginRun (edm::Run const &run, edm::EventSetup const &) override

virtual void	dqmAnalyze (DQMStore::IBooker &, DQMStore::IGetter &, edm::Event const &, edm::EventSetup const &)

	DQMEDHarvester (edm::ParameterSet const &iConfig)

	DQMEDHarvester ()

virtual void	dqmEndRun (DQMStore::IBooker &, DQMStore::IGetter &, edm::Run const &, edm::EventSetup const &)

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

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &es) final

void	endProcessBlockProduce (edm::ProcessBlock &) final

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &es) final

	~DQMEDHarvester () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::EndProcessBlockProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) 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

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	beginRun (edm::Run const &run, edm::EventSetup const &eSetup) override

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

void	dqmEndLuminosityBlock (DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &c) override

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Protected Attributes
bool	runOnEndJob_

bool	runOnEndLumi_

std::string	topFolderName_

Protected Attributes inherited from DQMEDHarvester
DQMStore *	dqmstore_

edm::GetterOfProducts< DQMToken >	jobmegetter_

edm::EDPutTokenT< DQMToken >	jobToken_

edm::GetterOfProducts< DQMToken >	lumimegetter_

edm::EDPutTokenT< DQMToken >	lumiToken_

edm::GetterOfProducts< DQMToken >	runmegetter_

edm::EDPutTokenT< DQMToken >	runToken_

Private Types
enum	trendPlots { offline, fpix, bpix, layer1, layer2, layer3, layer4, ring1, ring2 }

Private Member Functions
void	bookSummaries (DQMStore::IBooker &iBooker)

void	bookTrendPlots (DQMStore::IBooker &iBooker)

void	fillSummaries (DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, int lumiSeg=0)

void	fillTrendPlots (DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, int lumiSeg=0)

Private Attributes
edm::ParameterSet	conf_

MonitorElement *	deadROCSummary

std::vector< double >	deadRocThresholds_

std::map< trendPlots, MonitorElement * >	deadROCTrends_

std::vector< double >	deadRocWarnThresholds_

bool	firstLumi

std::map< trendPlots, MonitorElement * >	ineffROCTrends_

MonitorElement *	reportSummary

edm::InputTag	src_

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

std::map< std::string, std::string >	summaryPlotName_

Additional Inherited Members
Public Types inherited from DQMEDHarvester
typedef dqm::harvesting::DQMStore	DQMStore

typedef dqm::harvesting::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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)

Detailed Description

Description: Summary map generation for the Phase 1 pixel

Usage: <usage>

Definition at line 35 of file SiPixelPhase1Summary.h.

Member Enumeration Documentation

◆ trendPlots

enum SiPixelPhase1Summary::trendPlots

private

Enumerator
offline
fpix
bpix
layer1
layer2
layer3
layer4
ring1
ring2

Definition at line 59 of file SiPixelPhase1Summary.h.

59 { offline, fpix, bpix, layer1, layer2, layer3, layer4, ring1, ring2 };

SiPixelPhase1Summary::bpix

Definition: SiPixelPhase1Summary.h:59

SiPixelPhase1Summary::ring1

Definition: SiPixelPhase1Summary.h:59

SiPixelPhase1Summary::layer1

Definition: SiPixelPhase1Summary.h:59

SiPixelPhase1Summary::layer2

Definition: SiPixelPhase1Summary.h:59

SiPixelPhase1Summary::fpix

Definition: SiPixelPhase1Summary.h:59

SiPixelPhase1Summary::offline

Definition: SiPixelPhase1Summary.h:59

SiPixelPhase1Summary::layer4

Definition: SiPixelPhase1Summary.h:59

SiPixelPhase1Summary::layer3

Definition: SiPixelPhase1Summary.h:59

SiPixelPhase1Summary::ring2

Definition: SiPixelPhase1Summary.h:59

Constructor & Destructor Documentation

◆ SiPixelPhase1Summary()

SiPixelPhase1Summary::SiPixelPhase1Summary ( const edm::ParameterSet & conf )

explicit

Definition at line 49 of file SiPixelPhase1Summary.cc.

References conf_, deadRocThresholds_, deadRocWarnThresholds_, edm::ParameterSet::getParameter(), runOnEndJob_, runOnEndLumi_, AlCaHLTBitMon_QueryRunRegistry::string, summaryPlotName_, and topFolderName_.

     : DQMEDHarvester(iConfig), conf_(iConfig), firstLumi(true) {
   LogInfo("PixelDQM") << "SiPixelPhase1Summary::SiPixelPhase1Summary: Got DQM BackEnd interface" << endl;
   topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
   runOnEndLumi_ = conf_.getParameter<bool>("RunOnEndLumi");
   runOnEndJob_ = conf_.getParameter<bool>("RunOnEndJob");
 
   std::vector<edm::ParameterSet> mapPSets = conf_.getParameter<std::vector<edm::ParameterSet> >("SummaryMaps");
 
   //Go through the configuration file and add in
   for (auto const& mapPSet : mapPSets) {
     summaryPlotName_[mapPSet.getParameter<std::string>("MapName")] = mapPSet.getParameter<std::string>("MapHist");
   }
   deadRocThresholds_ = conf_.getParameter<std::vector<double> >("DeadROCErrorThreshold");
   deadRocWarnThresholds_ = conf_.getParameter<std::vector<double> >("DeadROCWarningThreshold");
 }

◆ ~SiPixelPhase1Summary()

SiPixelPhase1Summary::~SiPixelPhase1Summary ( )

override

Definition at line 66 of file SiPixelPhase1Summary.cc.

                                             {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   LogInfo("PixelDQM") << "SiPixelPhase1Summary::~SiPixelPhase1Summary: Destructor" << endl;
 }

Member Function Documentation

◆ beginRun()

void SiPixelPhase1Summary::beginRun	(	edm::Run const &	run,
		edm::EventSetup const &	eSetup
	)

overrideprotected

Definition at line 72 of file SiPixelPhase1Summary.cc.

72 {}

◆ bookSummaries()

void SiPixelPhase1Summary::bookSummaries ( DQMStore::IBooker & iBooker )

private

Definition at line 114 of file SiPixelPhase1Summary.cc.

References dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookFloat(), dqm::implementation::NavigatorBase::cd(), deadROCSummary, dqm::impl::MonitorElement::Fill(), mps_fire::i, dqmiolumiharvest::j, Skims_PA_cff::name, reportSummary, runOnEndLumi_, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), summaryMap_, and summaryPlotName_.

Referenced by dqmEndJob(), and dqmEndLuminosityBlock().

                                                                  {
   iBooker.cd();
 
   std::vector<std::string> xAxisLabels_ = {"BMO",
                                            "BMI",
                                            "BPO ",
                                            "BPI",
                                            "HCMO_1",
                                            "HCMO_2",
                                            "HCMI_1",
                                            "HCMI_2",
                                            "HCPO_1",
                                            "HCPO_2",
                                            "HCPI_1",
                                            "HCPI_2"};  // why not having a global variable !?!?!?!
   std::vector<std::string> yAxisLabels_ = {
       "1",
       "2",
       "3",
       "4"};  // why not having a global variable ?!?!?!!? - I originally did, but was told not to by David Lange!
 
   iBooker.setCurrentFolder("PixelPhase1/Summary");
   //Book the summary plots for the variables as described in the config file
   for (const auto& mapInfo : summaryPlotName_) {
     auto name = mapInfo.first;
     summaryMap_[name] = iBooker.book2D("pixel" + name + "Summary", "Pixel " + name + " Summary", 12, 0, 12, 4, 0, 4);
   }
   //Make the new 6 bin ROC summary
   deadROCSummary =
       iBooker.book2D("deadROCSummary", "Percentage of filled ROCs (with digis) per layer/ring", 2, 0, 2, 4, 0, 4);
   std::vector<std::string> xAxisLabelsReduced_ = {"Barrel", "Forward"};
   deadROCSummary->setAxisTitle("Subdetector", 1);
   for (unsigned int i = 0; i < xAxisLabelsReduced_.size(); i++) {
     deadROCSummary->setBinLabel(i + 1, xAxisLabelsReduced_[i]);
   }
 
   //Book the summary plot
   iBooker.setCurrentFolder("PixelPhase1/EventInfo");
 
   if (runOnEndLumi_) {
     //New less granular summary plot - this is currently only done online
     summaryMap_["Grand"] = iBooker.book2D("reportSummaryMap", "Pixel Summary Map", 2, 0, 2, 4, 0, 4);
     summaryMap_["Grand"]->setAxisTitle("Subdetector", 1);
     for (unsigned int i = 0; i < xAxisLabelsReduced_.size(); i++) {
       summaryMap_["Grand"]->setBinLabel(i + 1, xAxisLabelsReduced_[i]);
       for (unsigned int j = 0; j < 4; j++) {
         summaryMap_["Grand"]->Fill(i, j, -1);
       }
     }
   } else {
     //Book the original summary plot, for now juts doing this one offline.
     summaryMap_["Grand"] = iBooker.book2D("reportSummaryMap", "Pixel Summary Map", 12, 0, 12, 4, 0, 4);
   }
 
   reportSummary = iBooker.bookFloat("reportSummary");
 
   //Now set up axis and bin labels
   for (const auto& summaryMapEntry : summaryMap_) {
     if (summaryMapEntry.first == "Grand")
       continue;
     auto summaryMap = summaryMapEntry.second;
     for (unsigned int i = 0; i < xAxisLabels_.size(); i++) {
       summaryMap->setBinLabel(i + 1, xAxisLabels_[i], 1);
     }
     for (unsigned int i = 0; i < yAxisLabels_.size(); i++) {
       summaryMap->setBinLabel(i + 1, yAxisLabels_[i], 2);
     }
     summaryMap->setAxisTitle("Subdetector", 1);
     summaryMap->setAxisTitle("Layer/disk", 2);
     for (int i = 0; i < summaryMap->getTH1()->GetXaxis()->GetNbins(); i++) {    // !??!?!? xAxisLabels_.size() ?!?!
       for (int j = 0; j < summaryMap->getTH1()->GetYaxis()->GetNbins(); j++) {  // !??!?!? yAxisLabels_.size() ?!?!?!
         summaryMap->Fill(i, j, -1.);
       }
     }
   }
   reportSummary->Fill(-1.);
 
   //Reset the iBooker
   iBooker.setCurrentFolder("PixelPhase1/");
 }

◆ bookTrendPlots()

void SiPixelPhase1Summary::bookTrendPlots ( DQMStore::IBooker & iBooker )

private

Definition at line 198 of file SiPixelPhase1Summary.cc.

References dqm::implementation::IBooker::bookProfile(), deadROCTrends_, mps_fire::i, ineffROCTrends_, layer1, layer2, layer3, layer4, offline, ring1, ring2, runOnEndLumi_, dqm::implementation::NavigatorBase::setCurrentFolder(), and multiplicitycorr_cfi::yMax.

Referenced by dqmEndJob(), and dqmEndLuminosityBlock().

                                                                   {
   //We need different plots depending on if we're online (runOnEndLumi) or offline (!runOnEndLumi)
   iBooker.setCurrentFolder("PixelPhase1/");
   std::vector<string> binAxisLabels = {"Layer 1", "Layer 2", "Layer 3", "Layer 4", "Ring 1", "Ring 2"};
   if (runOnEndLumi_) {
     std::vector<trendPlots> histoOrder = {layer1, layer2, layer3, layer4, ring1, ring2};
     std::vector<string> varName = {"Layer_1", "Layer_2", "Layer_3", "Layer_4", "Ring_1", "Ring_2"};
     std::vector<int> yMax = {1536, 3584, 5632, 8192, 4224, 6528};
     for (unsigned int i = 0; i < histoOrder.size(); i++) {
       string varNameStr = "deadRocTrend" + varName[i];
       string varTitle = binAxisLabels[i] + " dead ROC trend";
       deadROCTrends_[histoOrder[i]] = iBooker.bookProfile(varNameStr, varTitle, 500, 0., 5000, 0., yMax[i], "");
       varNameStr = "ineffRocTrend" + varName[i];
       varTitle = binAxisLabels[i] + " inefficient ROC trend";
       ineffROCTrends_[histoOrder[i]] = iBooker.bookProfile(varNameStr, varTitle, 500, 0., 5000, 0., yMax[i], "");
       deadROCTrends_[histoOrder[i]]->setAxisTitle("Lumisection", 1);
       ineffROCTrends_[histoOrder[i]]->setAxisTitle("Lumisection", 1);
     }
   } else {
     deadROCTrends_[offline] = iBooker.bookProfile("deadRocTotal", "N dead ROCs", 6, 0., 6, 0., 8192, "");
     ineffROCTrends_[offline] = iBooker.bookProfile("ineffRocTotal", "N inefficient ROCs", 6, 0., 6, 0., 8192, "");
     deadROCTrends_[offline]->setAxisTitle("Subdetector", 1);
     ineffROCTrends_[offline]->setAxisTitle("Subdetector", 1);
     for (unsigned int i = 1; i <= binAxisLabels.size(); i++) {
       deadROCTrends_[offline]->setBinLabel(i, binAxisLabels[i - 1]);
       ineffROCTrends_[offline]->setBinLabel(i, binAxisLabels[i - 1]);
     }
   }
 }

◆ dqmEndJob()

void SiPixelPhase1Summary::dqmEndJob	(	DQMStore::IBooker &	iBooker,
		DQMStore::IGetter &	iGetter
	)

overrideprotectedvirtual

Implements DQMEDHarvester.

Definition at line 96 of file SiPixelPhase1Summary.cc.

References bookSummaries(), bookTrendPlots(), fillSummaries(), fillTrendPlots(), firstLumi, runOnEndJob_, and runOnEndLumi_.

                                                                                        {
   if (firstLumi) {  //Book the plots in the (maybe possible?) case that they aren't booked in the dqmEndLuminosityBlock method
     bookSummaries(iBooker);
     bookTrendPlots(iBooker);
     firstLumi = false;
   }
   if (runOnEndJob_) {
     if (!runOnEndLumi_)
       fillTrendPlots(
           iBooker,
           iGetter);  //If we're filling these plots at the end lumi step, it doesn't really make sense to also do them at the end job
     fillSummaries(iBooker, iGetter);
   }
 }

◆ dqmEndLuminosityBlock()

void SiPixelPhase1Summary::dqmEndLuminosityBlock	(	DQMStore::IBooker &	iBooker,
		DQMStore::IGetter &	iGetter,
		edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	c
	)

overrideprotectedvirtual

Reimplemented from DQMEDHarvester.

Definition at line 74 of file SiPixelPhase1Summary.cc.

References bookSummaries(), bookTrendPlots(), fillSummaries(), fillTrendPlots(), firstLumi, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), and runOnEndLumi_.

                                                                          {
   if (firstLumi) {
     bookSummaries(iBooker);
     bookTrendPlots(iBooker);
     firstLumi = false;
   }
 
   if (runOnEndLumi_) {
     int lumiSec = lumiSeg.id().luminosityBlock();
     fillTrendPlots(iBooker, iGetter, lumiSec);
     fillSummaries(iBooker, iGetter, lumiSec);
   }
 
   //  iBooker.cd();
 }

◆ fillSummaries()

void SiPixelPhase1Summary::fillSummaries	(	DQMStore::IBooker &	iBooker,
		DQMStore::IGetter &	iGetter,
		int	lumiSeg = `0`
	)

private

Definition at line 231 of file SiPixelPhase1Summary.cc.

References deadROCSummary, deadROCTrends_, relativeConstraints::empty, dqm::impl::MonitorElement::Fill(), dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getBinContent(), mps_fire::i, dqmiolumiharvest::j, layer1, layer2, layer3, layer4, hlt_dqm_clientPB-live_cfg::me, reco::HaloData::minus, dqmiodumpmetadata::n, Skims_PA_cff::name, offline, SurfaceOrientation::outer, reportSummary, ring1, ring2, runOnEndLumi_, dqm::impl::MonitorElement::setBinContent(), jetUpdater_cfi::sort, AlCaHLTBitMon_QueryRunRegistry::string, summaryMap_, summaryPlotName_, topFolderName_, photonAnalyzer_cfi::xBin, and photonAnalyzer_cfi::yBin.

Referenced by dqmEndJob(), and dqmEndLuminosityBlock().

                                                                                                         {
   //Firstly, we will fill the regular summary maps.
   for (const auto& mapInfo : summaryPlotName_) {
     auto name = mapInfo.first;
     std::ostringstream histNameStream;
     std::string histName;
 
     for (int i = 0; i < 12; i++) {   // !??!?!? xAxisLabels_.size() ?!?!
       for (int j = 0; j < 4; j++) {  // !??!?!? yAxisLabels_.size() ?!?!?!
         if (i > 3 && j == 3)
           continue;
         bool minus = i < 2 || (i > 3 && i < 8);  // bleah !
         int iOver2 = floor(i / 2.);
         bool outer = (i > 3) ? iOver2 % 2 == 0 : i % 2 == 0;
         //Complicated expression that creates the name of the histogram we are interested in.
         histNameStream.str("");
         histNameStream << topFolderName_.c_str() << "PX" << ((i > 3) ? "Forward" : "Barrel") << "/"
                        << ((i > 3) ? "HalfCylinder" : "Shell") << "_" << (minus ? "m" : "p") << ((outer) ? "O" : "I")
                        << "/" << ((i > 3) ? ((i % 2 == 0) ? "PXRing_1/" : "PXRing_2/") : "")
                        << summaryPlotName_[name].c_str() << "_PX" << ((i > 3) ? "Disk" : "Layer") << "_"
                        << ((i > 3) ? ((minus) ? "-" : "+") : "") << (j + 1);
         histName = histNameStream.str();
         MonitorElement* me = iGetter.get(histName);
 
         if (!me) {
           edm::LogWarning("SiPixelPhase1Summary") << "ME " << histName << " is not available !!";
           continue;  // Ignore non-existing MEs, as this can cause the whole thing to crash
         }
 
         if (summaryMap_[name] == nullptr) {
           edm::LogWarning("SiPixelPhase1Summary") << "Summary map " << name << " is not available !!";
           continue;  // Based on reported errors it seems possible that we're trying to access a non-existant summary map, so if the map doesn't exist but we're trying to access it here we'll skip it instead.
         }
         if (!(me->getQReports()).empty())
           summaryMap_[name]->setBinContent(i + 1, j + 1, (me->getQReports())[0]->getQTresult());
         else
           summaryMap_[name]->setBinContent(i + 1, j + 1, -1);
       }
     }
   }
 
   //Fill the dead ROC summary
   std::vector<trendPlots> trendOrder = {layer1, layer2, layer3, layer4, ring1, ring2};
   std::vector<int> nRocsPerTrend = {1536, 3584, 5632, 8192, 4224, 6528};
   for (unsigned int i = 0; i < trendOrder.size(); i++) {
     int xBin = i < 4 ? 1 : 2;
     int yBin = i % 4 + 1;
     float numDeadROCs = 0.;
     float numFilledROCs = 0.;
     float fracFilledROCs = 1.;
     unsigned int lastTrendBin = 0;
     int numEmptyBins = 0;
     float sumOf5Bins = 0.;
     std::vector<float> last5TrendBinsValue = {0., 0., 0., 0., 0.};
     if (runOnEndLumi_) {  //Online case
       TH1* tempProfile = deadROCTrends_[trendOrder[i]]->getTH1();
 
       // Old strategy: simply pick the last bin value (deadROCTrend: 1 bin every 10 LS)
       //nROCs = tempProfile->GetBinContent(tempProfile->FindLastBinAbove());
 
       // New strategy, check last 5 bins and decide
       lastTrendBin = lumiSec / 10;
       if (lastTrendBin >= 5) {
         // Fill vector with value of last 5 bins of deadROCTrend
         for (unsigned int n = 0; n < 5; n++)
           last5TrendBinsValue[n] = tempProfile->GetBinContent(lastTrendBin - 5 + n + 1);
       } else {
         // Fill vector with bins from 1 to last available bin
         for (unsigned int n = 0; n < lastTrendBin; n++)
           last5TrendBinsValue[n] = tempProfile->GetBinContent(n + 1);
       }
 
       for (unsigned int n = 0; n < 5; n++)
         if (last5TrendBinsValue[n] == 0)
           numEmptyBins++;
       if (numEmptyBins == 0) {
         // If the last 5 bins are filled: take the median value (to suppress effect of peaks)
         std::sort(last5TrendBinsValue.begin(), last5TrendBinsValue.end());
         numDeadROCs = last5TrendBinsValue[2];
       } else if (numEmptyBins != 5) {
         //If there are some empty bins: take the average of the non-empty bins
         for (unsigned int n = 0; n < 5; n++)
           sumOf5Bins += last5TrendBinsValue[n];
         numDeadROCs = sumOf5Bins / (5 - numEmptyBins);
       }
 
     } else {  //Offline case
       TH1* tempProfile = deadROCTrends_[offline]->getTH1();
       numDeadROCs = tempProfile->GetBinContent(i + 1);
     }
     // Filled ROCs = Total number - dead ROCs
     numFilledROCs = nRocsPerTrend[i] - numDeadROCs;
     //Fill with fraction of filled ROCs (with digis)
     fracFilledROCs = numFilledROCs / nRocsPerTrend[i];
     deadROCSummary->setBinContent(xBin, yBin, fracFilledROCs);
     deadROCSummary->setBinContent(2, 3, -1);
     deadROCSummary->setBinContent(2, 4, -1);
   }
 
   //Sum of non-negative bins for the reportSummary
   float sumOfNonNegBins = 0.;
   //Now we will use the other summary maps to create the overall map.
   //For now we only want to do this offline
   if (!runOnEndLumi_) {
     for (int i = 0; i < 12; i++) {  // !??!?!? xAxisLabels_.size() ?!?!
       if (summaryMap_["Grand"] == nullptr) {
         edm::LogWarning("SiPixelPhase1Summary") << "Grand summary does not exist!";
         break;
       }
       for (int j = 0; j < 4; j++) {  // !??!?!? yAxisLabels_.size() ?!?!?!
         summaryMap_["Grand"]->setBinContent(
             i + 1,
             j + 1,
             1);  // This resets the map to be good. We only then set it to 0 if there has been a problem in one of the other summaries.
         for (auto const& mapInfo : summaryPlotName_) {  //Check summary maps
           auto name = mapInfo.first;
           if (summaryMap_[name] == nullptr) {
             edm::LogWarning("SiPixelPhase1Summary") << "Summary " << name << " does not exist!";
             continue;
           }
           if (summaryMap_["Grand"]->getBinContent(i + 1, j + 1) > summaryMap_[name]->getBinContent(i + 1, j + 1))
             summaryMap_["Grand"]->setBinContent(i + 1, j + 1, summaryMap_[name]->getBinContent(i + 1, j + 1));
         }
         if (summaryMap_["Grand"]->getBinContent(i + 1, j + 1) > -0.1)
           sumOfNonNegBins += summaryMap_["Grand"]->getBinContent(i + 1, j + 1);
       }
     }
     reportSummary->Fill(sumOfNonNegBins / 40.);  // The average of the 40 useful bins in the summary map.
   }
 
   //Fill the new overall map
   //  if (!runOnEndLumi_) return;
   else {  //Do this for online only
     for (int i = 0; i < 2; i++) {
       if (summaryMap_["Grand"] == nullptr) {
         edm::LogWarning("SiPixelPhase1Summary") << "Grand summary does not exist!";
         break;
       }
       for (int j = 0; j < 4; j++) {  // !??!?!? yAxisLabels_.size() ?!?!?!
         //Ignore the bins without detectors in them
         if (i == 1 && j > 1) {
           summaryMap_["Grand"]->setBinContent(i + 1, j + 1, -1);
         } else {
           //Copy dead ROC Summary into Online Pixel Summary
           summaryMap_["Grand"]->setBinContent(i + 1, j + 1, deadROCSummary->getBinContent(i + 1, j + 1));
 
           sumOfNonNegBins += summaryMap_["Grand"]->getBinContent(i + 1, j + 1);
         }
       }
     }
     reportSummary->Fill(sumOfNonNegBins / 6.);  // The average of the 6 useful bins in the summary map.
   }
 }

◆ fillTrendPlots()

void SiPixelPhase1Summary::fillTrendPlots	(	DQMStore::IBooker &	iBooker,
		DQMStore::IGetter &	iGetter,
		int	lumiSeg = `0`
	)

private

Definition at line 387 of file SiPixelPhase1Summary.cc.

References dqm::implementation::NavigatorBase::cd(), deadROCTrends_, dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getBinContent(), dqm::legacy::MonitorElement::getTH1(), mps_fire::i, ineffROCTrends_, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmiolumiharvest::j, layer1, layer2, layer3, layer4, offline, dqm::impl::MonitorElement::Reset(), ring1, ring2, runOnEndLumi_, AlCaHLTBitMon_QueryRunRegistry::string, and to_string().

Referenced by dqmEndJob(), and dqmEndLuminosityBlock().

                                                                                                          {
   // If we're running in online mode and the lumi section is not modulo 10, return. Offline running always uses lumiSec=0, so it will pass this test.
   if (lumiSec % 10 != 0)
     return;
 
   if (runOnEndLumi_) {
     MonitorElement* nClustersAll = iGetter.get("PixelPhase1/Phase1_MechanicalView/num_clusters_per_Lumisection_PXAll");
     if (nClustersAll == nullptr) {
       edm::LogWarning("SiPixelPhase1Summary") << "All pixel cluster trend plot not available!!";
       return;
     }
     if (nClustersAll->getTH1()->GetBinContent(lumiSec) < 100)
       return;
   }
 
   std::string histName;
 
   //Find the total number of filled bins and hi efficiency bins
   std::vector<trendPlots> trendOrder = {layer1, layer2, layer3, layer4, ring1, ring2};
   std::vector<int> nFilledROCs(trendOrder.size(), 0);
   std::vector<int> hiEffROCs(trendOrder.size(), 0);
   std::vector<int> nRocsPerTrend = {1536, 3584, 5632, 8192, 4224, 6528};
   std::vector<string> trendNames = {};
 
   for (auto it : {1, 2, 3, 4}) {
     histName = "PXBarrel/digi_occupancy_per_SignedModuleCoord_per_SignedLadderCoord_PXLayer_" + std::to_string(it);
     trendNames.push_back(histName);
   }
   for (auto it : {1, 2}) {
     histName = "PXForward/digi_occupancy_per_SignedDiskCoord_per_SignedBladePanelCoord_PXRing_" + std::to_string(it);
     trendNames.push_back(histName);
   }
   //Loop over layers. This will also do the rings, but we'll skip the ring calculation for
   for (unsigned int trendIt = 0; trendIt < trendOrder.size(); trendIt++) {
     iGetter.cd();
     histName = "PixelPhase1/Phase1_MechanicalView/" + trendNames[trendIt];
     MonitorElement* tempLayerME = iGetter.get(histName);
     if (tempLayerME == nullptr)
       continue;
     float lowEffValue = 0.25 * (tempLayerME->getTH1()->Integral() / nRocsPerTrend[trendIt]);
     for (int i = 1; i <= tempLayerME->getTH1()->GetXaxis()->GetNbins(); i++) {
       for (int j = 1; j <= tempLayerME->getTH1()->GetYaxis()->GetNbins(); j++) {
         if (tempLayerME->getBinContent(i, j) > 0.)
           nFilledROCs[trendIt]++;
         if (tempLayerME->getBinContent(i, j) > lowEffValue)
           hiEffROCs[trendIt]++;
       }
     }
     if (runOnEndLumi_) {
       tempLayerME->Reset();  //If we're doing online monitoring, reset the digi maps.
     }
   }  // Close layers/ring loop
 
   if (!runOnEndLumi_) {  //offline
     for (unsigned int i = 0; i < trendOrder.size(); i++) {
       deadROCTrends_[offline]->Fill(i, nRocsPerTrend[i] - nFilledROCs[i]);
       ineffROCTrends_[offline]->Fill(i, nFilledROCs[i] - hiEffROCs[i]);
     }
   } else {  //online
     for (unsigned int i = 0; i < trendOrder.size(); i++) {
       deadROCTrends_[trendOrder[i]]->Fill(lumiSec - 1, nRocsPerTrend[i] - nFilledROCs[i]);
       ineffROCTrends_[trendOrder[i]]->Fill(lumiSec - 1, nFilledROCs[i] - hiEffROCs[i]);
     }
   }
 
   if (!runOnEndLumi_)
     return;  // The following only occurs in the online
   //Reset some MEs every 10LS here
   for (auto it : {1, 2, 3, 4}) {  //PXBarrel
     histName = "PixelPhase1/Phase1_MechanicalView/PXBarrel/clusterposition_zphi_PXLayer_" + std::to_string(it);
     MonitorElement* toReset = iGetter.get(histName);
     if (toReset != nullptr) {
       toReset->Reset();
     }
   }
   for (auto it : {"-3", "-2", "-1", "+1", "+2", "+3"}) {  //PXForward
     histName = "PixelPhase1/Phase1_MechanicalView/PXForward/clusterposition_xy_PXDisk_" + std::string(it);
     MonitorElement* toReset = iGetter.get(histName);
     if (toReset != nullptr) {
       toReset->Reset();
     }
   }
 }

Member Data Documentation

◆ conf_

edm::ParameterSet SiPixelPhase1Summary::conf_

private

Definition at line 60 of file SiPixelPhase1Summary.h.

Referenced by SiPixelPhase1Summary().

◆ deadROCSummary

MonitorElement* SiPixelPhase1Summary::deadROCSummary

private

Definition at line 65 of file SiPixelPhase1Summary.h.

Referenced by bookSummaries(), and fillSummaries().

◆ deadRocThresholds_

std::vector<double> SiPixelPhase1Summary::deadRocThresholds_

private

Definition at line 75 of file SiPixelPhase1Summary.h.

Referenced by SiPixelPhase1Summary().

◆ deadROCTrends_

std::map<trendPlots, MonitorElement*> SiPixelPhase1Summary::deadROCTrends_

private

Definition at line 71 of file SiPixelPhase1Summary.h.

Referenced by bookTrendPlots(), fillSummaries(), and fillTrendPlots().

◆ deadRocWarnThresholds_

std::vector<double> SiPixelPhase1Summary::deadRocWarnThresholds_

private

Definition at line 76 of file SiPixelPhase1Summary.h.

Referenced by SiPixelPhase1Summary().

◆ firstLumi

bool SiPixelPhase1Summary::firstLumi

private

Definition at line 62 of file SiPixelPhase1Summary.h.

Referenced by dqmEndJob(), and dqmEndLuminosityBlock().

◆ ineffROCTrends_

std::map<trendPlots, MonitorElement*> SiPixelPhase1Summary::ineffROCTrends_

private

Definition at line 72 of file SiPixelPhase1Summary.h.

Referenced by bookTrendPlots(), and fillTrendPlots().

◆ reportSummary

MonitorElement* SiPixelPhase1Summary::reportSummary

private

Definition at line 66 of file SiPixelPhase1Summary.h.

Referenced by bookSummaries(), and fillSummaries().

◆ runOnEndJob_

bool SiPixelPhase1Summary::runOnEndJob_

protected

Definition at line 56 of file SiPixelPhase1Summary.h.

Referenced by dqmEndJob(), and SiPixelPhase1Summary().

◆ runOnEndLumi_

bool SiPixelPhase1Summary::runOnEndLumi_

protected

Definition at line 55 of file SiPixelPhase1Summary.h.

Referenced by bookSummaries(), bookTrendPlots(), dqmEndJob(), dqmEndLuminosityBlock(), fillSummaries(), fillTrendPlots(), and SiPixelPhase1Summary().

◆ src_

edm::InputTag SiPixelPhase1Summary::src_

private

Definition at line 61 of file SiPixelPhase1Summary.h.

◆ summaryMap_

std::map<std::string, MonitorElement*> SiPixelPhase1Summary::summaryMap_

private

Definition at line 64 of file SiPixelPhase1Summary.h.

Referenced by bookSummaries(), and fillSummaries().

◆ summaryPlotName_

std::map<std::string, std::string> SiPixelPhase1Summary::summaryPlotName_

private

Definition at line 68 of file SiPixelPhase1Summary.h.

Referenced by bookSummaries(), fillSummaries(), and SiPixelPhase1Summary().

◆ topFolderName_

std::string SiPixelPhase1Summary::topFolderName_

protected

Definition at line 54 of file SiPixelPhase1Summary.h.

Referenced by fillSummaries(), and SiPixelPhase1Summary().

Public Member Functions

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ trendPlots

Constructor & Destructor Documentation

◆ SiPixelPhase1Summary()

◆ ~SiPixelPhase1Summary()

Member Function Documentation

◆ beginRun()

◆ bookSummaries()

◆ bookTrendPlots()

◆ dqmEndJob()

◆ dqmEndLuminosityBlock()

◆ fillSummaries()

◆ fillTrendPlots()

Member Data Documentation

◆ conf_

◆ deadROCSummary

◆ deadRocThresholds_

◆ deadROCTrends_

◆ deadRocWarnThresholds_

◆ firstLumi

◆ ineffROCTrends_

◆ reportSummary

◆ runOnEndJob_

◆ runOnEndLumi_

◆ src_

◆ summaryMap_

◆ summaryPlotName_

◆ topFolderName_