#include <DTnoiseDBValidation.h>

Inheritance diagram for DTnoiseDBValidation:

Public Types
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

void	beginRun (const edm::Run &run, const edm::EventSetup &setup) override
	Operations. More...

	DTnoiseDBValidation (const edm::ParameterSet &pset)
	Constructor. More...

void	endJob () override

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

	~DTnoiseDBValidation () override
	Destructor. More...

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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

ESProxyIndex const *	esGetTokenIndices (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::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)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	bookHisto (const DTChamberId &)

Private Attributes
DQMStore *	dbe_

MonitorElement *	diffHisto_

std::string	diffTestName_

edm::ESHandle< DTGeometry >	dtGeom_

std::string	labelDB_

std::string	labelDBRef_

MonitorElement *	layerHisto_

std::string	layerTestName_

std::map< DTChamberId, MonitorElement * >	noiseHistoMap_

const DTStatusFlag *	noiseMap_

const DTStatusFlag *	noiseRefMap_

int	noisyCellsRef_

int	noisyCellsValid_

std::string	outputFileName_

bool	outputMEsInRootFile_

MonitorElement *	sectorHisto_

std::string	sectorTestName_

MonitorElement *	stationHisto_

std::string	stationTestName_

MonitorElement *	wheelHisto_

std::string	wheelTestName_

Additional Inherited Members
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

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

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

Plot the noise from the DB comparaison

Author: G. Mila - INFN Torino

Definition at line 28 of file DTnoiseDBValidation.h.

Member Typedef Documentation

typedef dqm::legacy::DQMStore DTnoiseDBValidation::DQMStore

Definition at line 31 of file DTnoiseDBValidation.h.

typedef dqm::legacy::MonitorElement DTnoiseDBValidation::MonitorElement

Definition at line 30 of file DTnoiseDBValidation.h.

Constructor & Destructor Documentation

DTnoiseDBValidation::DTnoiseDBValidation ( const edm::ParameterSet & pset )

Constructor.

Definition at line 39 of file DTnoiseDBValidation.cc.

References dbe_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), Utilities::operator, dqm::dqmstoreimpl::DQMStore::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                  {
   LogVerbatim("NoiseDBValidation") << "[DTnoiseDBValidation] Constructor called!";
 
   // Get the DQM needed services
   dbe_ = edm::Service<DQMStore>().operator->();
   dbe_->setCurrentFolder("DT/DtCalib/NoiseDBValidation");
 
   // Get dataBase label
   labelDBRef_ = pset.getParameter<string>("labelDBRef");
   labelDB_ = pset.getParameter<string>("labelDB");
 
   diffTestName_ = "noiseDifferenceInRange";
   if (pset.exists("diffTestName"))
     diffTestName_ = pset.getParameter<string>("diffTestName");
 
   wheelTestName_ = "noiseWheelOccInRange";
   if (pset.exists("wheelTestName"))
     wheelTestName_ = pset.getParameter<string>("wheelTestName");
 
   stationTestName_ = "noiseStationOccInRange";
   if (pset.exists("stationTestName"))
     stationTestName_ = pset.getParameter<string>("stationTestName");
 
   sectorTestName_ = "noiseSectorOccInRange";
   if (pset.exists("sectorTestName"))
     sectorTestName_ = pset.getParameter<string>("sectorTestName");
 
   layerTestName_ = "noiseLayerOccInRange";
   if (pset.exists("layerTestName"))
     layerTestName_ = pset.getParameter<string>("layerTestName");
 
   outputMEsInRootFile_ = false;
   if (pset.exists("OutputFileName")) {
     outputMEsInRootFile_ = true;
     outputFileName_ = pset.getParameter<std::string>("OutputFileName");
   }
 }

DTnoiseDBValidation::~DTnoiseDBValidation ( )

override

Destructor.

Definition at line 77 of file DTnoiseDBValidation.cc.

77 {}

Member Function Documentation

void DTnoiseDBValidation::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

inlineoverride

Definition at line 42 of file DTnoiseDBValidation.h.

References bookHisto().

42 {}

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

override

Operations.

Definition at line 79 of file DTnoiseDBValidation.cc.

References funct::abs(), plotT0FromHistos::binNumber(), dqm::dqmstoreimpl::DQMStore::book1D(), DTSuperLayerId::chamberId(), dbe_, edm::EventSetup::get(), DTTopology::lastChannel(), DTLayerId::layer(), DTWireId::layerId(), hgcalDigitizer_cfi::noise, dtNoiseDBValidation_cfg::noiseRef, Scenarios_cff::scale, dqm::impl::MonitorElement::setBinLabel(), relativeConstraints::station, DTSuperLayerId::superLayer(), DTLayerId::superlayerId(), makeMuonMisalignmentScenario::wheel, and DTWireId::wire().

                                                                              {
   ESHandle<DTStatusFlag> noiseRef;
   setup.get<DTStatusFlagRcd>().get(labelDBRef_, noiseRef);
   noiseRefMap_ = &*noiseRef;
 
   ESHandle<DTStatusFlag> noiseValid;
   setup.get<DTStatusFlagRcd>().get(labelDB_, noiseValid);
   noiseMap_ = &*noiseValid;
 
   // Get the geometry
   setup.get<MuonGeometryRecord>().get(dtGeom_);
 
   LogVerbatim("NoiseDBValidation") << "[DTnoiseDBValidation] Parameters initialization";
 
   noisyCellsRef_ = 0;
   noisyCellsValid_ = 0;
 
   // Histo booking
   diffHisto_ =
       dbe_->book1D("noisyCellDiff", "percentual (wrt the previous db) total number of noisy cells", 1, 0.5, 1.5);
   diffHisto_->setBinLabel(1, "Diff");
   wheelHisto_ = dbe_->book1D("wheelOccupancy", "percentual noisy cells occupancy per wheel", 5, -2.5, 2.5);
   wheelHisto_->setBinLabel(1, "Wh-2");
   wheelHisto_->setBinLabel(2, "Wh-1");
   wheelHisto_->setBinLabel(3, "Wh0");
   wheelHisto_->setBinLabel(4, "Wh1");
   wheelHisto_->setBinLabel(5, "Wh2");
   stationHisto_ = dbe_->book1D("stationOccupancy", "percentual noisy cells occupancy per station", 4, 0.5, 4.5);
   stationHisto_->setBinLabel(1, "St1");
   stationHisto_->setBinLabel(2, "St2");
   stationHisto_->setBinLabel(3, "St3");
   stationHisto_->setBinLabel(4, "St4");
   sectorHisto_ = dbe_->book1D("sectorOccupancy", "percentual noisy cells occupancy per sector", 12, 0.5, 12.5);
   sectorHisto_->setBinLabel(1, "Sect1");
   sectorHisto_->setBinLabel(2, "Sect2");
   sectorHisto_->setBinLabel(3, "Sect3");
   sectorHisto_->setBinLabel(4, "Sect4");
   sectorHisto_->setBinLabel(5, "Sect5");
   sectorHisto_->setBinLabel(6, "Sect6");
   sectorHisto_->setBinLabel(7, "Sect7");
   sectorHisto_->setBinLabel(8, "Sect8");
   sectorHisto_->setBinLabel(9, "Sect9");
   sectorHisto_->setBinLabel(10, "Sect10");
   sectorHisto_->setBinLabel(11, "Sect11");
   sectorHisto_->setBinLabel(12, "Sect12");
   layerHisto_ = dbe_->book1D("layerOccupancy", "percentual noisy cells occupancy per layer", 3, 0.5, 3.5);
   layerHisto_->setBinLabel(1, "First 10 bins");
   layerHisto_->setBinLabel(2, "Middle bins");
   layerHisto_->setBinLabel(3, "Last 10 bins");
 
   // map initialization
   map<int, int> whMap;
   whMap.clear();
   map<int, int> stMap;
   stMap.clear();
   map<int, int> sectMap;
   sectMap.clear();
   map<int, int> layerMap;
   layerMap.clear();
 
   // Loop over reference DB entries
   for (DTStatusFlag::const_iterator noise = noiseRefMap_->begin(); noise != noiseRefMap_->end(); noise++) {
     DTWireId wireId((*noise).first.wheelId,
                     (*noise).first.stationId,
                     (*noise).first.sectorId,
                     (*noise).first.slId,
                     (*noise).first.layerId,
                     (*noise).first.cellId);
     LogVerbatim("NoiseDBValidation") << "Ref. noisy wire: " << wireId;
     ++noisyCellsRef_;
   }
 
   // Loop over validation DB entries
   for (DTStatusFlag::const_iterator noise = noiseMap_->begin(); noise != noiseMap_->end(); noise++) {
     DTWireId wireId((*noise).first.wheelId,
                     (*noise).first.stationId,
                     (*noise).first.sectorId,
                     (*noise).first.slId,
                     (*noise).first.layerId,
                     (*noise).first.cellId);
     LogVerbatim("NoiseDBValidation") << "Valid. noisy wire: " << wireId;
     ++noisyCellsValid_;
 
     whMap[(*noise).first.wheelId]++;
     stMap[(*noise).first.stationId]++;
     sectMap[(*noise).first.sectorId]++;
 
     const DTTopology &dtTopo = dtGeom_->layer(wireId.layerId())->specificTopology();
     const int lastWire = dtTopo.lastChannel();
     if ((*noise).first.cellId <= 10)
       layerMap[1]++;
     if ((*noise).first.cellId > 10 && (*noise).first.cellId < (lastWire - 10))
       layerMap[2]++;
     if ((*noise).first.cellId >= (lastWire - 10))
       layerMap[3]++;
 
     const DTChamberId chId = wireId.layerId().superlayerId().chamberId();
     if (noiseHistoMap_.find(chId) == noiseHistoMap_.end())
       bookHisto(chId);
     int binNumber = 4 * (wireId.superLayer() - 1) + wireId.layer();
     noiseHistoMap_[chId]->Fill(wireId.wire(), binNumber);
   }
 
   // histo filling
   double scale = 1 / double(noisyCellsRef_);
   diffHisto_->Fill(1, abs(noisyCellsRef_ - noisyCellsValid_) * scale);
 
   scale = 1 / double(noisyCellsValid_);
   for (map<int, int>::const_iterator wheel = whMap.begin(); wheel != whMap.end(); wheel++) {
     wheelHisto_->Fill((*wheel).first, ((*wheel).second) * scale);
   }
   for (map<int, int>::const_iterator station = stMap.begin(); station != stMap.end(); station++) {
     stationHisto_->Fill((*station).first, ((*station).second) * scale);
   }
   for (map<int, int>::const_iterator sector = sectMap.begin(); sector != sectMap.end(); sector++) {
     sectorHisto_->Fill((*sector).first, ((*sector).second) * scale);
   }
   for (map<int, int>::const_iterator layer = layerMap.begin(); layer != layerMap.end(); layer++) {
     layerHisto_->Fill((*layer).first, ((*layer).second) * scale);
   }
 }

void DTnoiseDBValidation::bookHisto ( const DTChamberId & chId )

private

Definition at line 275 of file DTnoiseDBValidation.cc.

References plotT0FromHistos::binNumber(), dqm::dqmstoreimpl::DQMStore::book2D(), dbe_, HltBtagPostValidation_cff::histoName, label, hgcalTopologyTester_cfi::layers, DTChamberId::sector(), DTChamberId::station(), DTChamber::superLayers(), and DTChamberId::wheel().

Referenced by analyze().

                                                            {
   stringstream histoName;
   histoName << "NoiseOccupancy"
             << "_W" << chId.wheel() << "_St" << chId.station() << "_Sec" << chId.sector();
 
   if (noiseHistoMap_.find(chId) == noiseHistoMap_.end()) {  // Redundant check
     // Get the chamber from the geometry
     int nWiresMax = 0;
     const DTChamber *dtchamber = dtGeom_->chamber(chId);
     const vector<const DTSuperLayer *> &superlayers = dtchamber->superLayers();
 
     // Loop over layers and find the max # of wires
     for (vector<const DTSuperLayer *>::const_iterator sl = superlayers.begin(); sl != superlayers.end();
          ++sl) {  // loop over SLs
       vector<const DTLayer *> layers = (*sl)->layers();
       for (vector<const DTLayer *>::const_iterator lay = layers.begin(); lay != layers.end();
            ++lay) {  // loop over layers
         int nWires = (*lay)->specificTopology().channels();
         if (nWires > nWiresMax)
           nWiresMax = nWires;
       }
     }
 
     noiseHistoMap_[chId] = dbe_->book2D(histoName.str(), "Noise occupancy", nWiresMax, 1, (nWiresMax + 1), 12, 1, 13);
     for (int i_sl = 1; i_sl <= 3; ++i_sl) {
       for (int i_lay = 1; i_lay <= 4; ++i_lay) {
         int binNumber = 4 * (i_sl - 1) + i_lay;
         stringstream label;
         label << "SL" << i_sl << ": L" << i_lay;
         noiseHistoMap_[chId]->setBinLabel(binNumber, label.str(), 2);
       }
     }
   }
 }

void DTnoiseDBValidation::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 269 of file DTnoiseDBValidation.cc.

References dbe_, and dqm::dqmstoreimpl::DQMStore::save().

                                  {
   // Write the histos in a ROOT file
   if (outputMEsInRootFile_)
     dbe_->save(outputFileName_);
 }

void DTnoiseDBValidation::endRun	(	edm::Run const &	run,
		edm::EventSetup const &	setup
	)

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 201 of file DTnoiseDBValidation.cc.

References QReport::getBadChannels(), and makeMuonMisalignmentScenario::wheel.

                                                                               {
   // test on difference histo
   // string testCriterionName;
   // testCriterionName =
   // parameters.getUntrackedParameter<string>("diffTestName","noiseDifferenceInRange");
   const QReport *theDiffQReport = diffHisto_->getQReport(diffTestName_);
   if (theDiffQReport) {
     vector<dqm::me_util::Channel> badChannels = theDiffQReport->getBadChannels();
     for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); channel != badChannels.end();
          channel++) {
       LogWarning("NoiseDBValidation") << " Bad partial difference of noisy channels! Contents : "
                                       << (*channel).getContents();
     }
   }
   // testCriterionName =
   // parameters.getUntrackedParameter<string>("wheelTestName","noiseWheelOccInRange");
   const QReport *theDiffQReport2 = wheelHisto_->getQReport(wheelTestName_);
   if (theDiffQReport2) {
     vector<dqm::me_util::Channel> badChannels = theDiffQReport2->getBadChannels();
     for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); channel != badChannels.end();
          channel++) {
       int wheel = (*channel).getBin() - 3;
       LogWarning("NoiseDBValidation") << " Bad percentual occupancy for wheel : " << wheel
                                       << "  Contents : " << (*channel).getContents();
     }
   }
   // testCriterionName =
   // parameters.getUntrackedParameter<string>("stationTestName","noiseStationOccInRange");
   const QReport *theDiffQReport3 = stationHisto_->getQReport(stationTestName_);
   if (theDiffQReport3) {
     vector<dqm::me_util::Channel> badChannels = theDiffQReport3->getBadChannels();
     for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); channel != badChannels.end();
          channel++) {
       LogWarning("NoiseDBValidation") << " Bad percentual occupancy for station : " << (*channel).getBin()
                                       << "  Contents : " << (*channel).getContents();
     }
   }
   // testCriterionName =
   // parameters.getUntrackedParameter<string>("sectorTestName","noiseSectorOccInRange");
   const QReport *theDiffQReport4 = sectorHisto_->getQReport(sectorTestName_);
   if (theDiffQReport4) {
     vector<dqm::me_util::Channel> badChannels = theDiffQReport4->getBadChannels();
     for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); channel != badChannels.end();
          channel++) {
       LogWarning("NoiseDBValidation") << " Bad percentual occupancy for sector : " << (*channel).getBin()
                                       << "  Contents : " << (*channel).getContents();
     }
   }
   // testCriterionName =
   // parameters.getUntrackedParameter<string>("layerTestName","noiseLayerOccInRange");
   const QReport *theDiffQReport5 = layerHisto_->getQReport(layerTestName_);
   if (theDiffQReport5) {
     vector<dqm::me_util::Channel> badChannels = theDiffQReport5->getBadChannels();
     for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); channel != badChannels.end();
          channel++) {
       if ((*channel).getBin() == 1)
         LogWarning("NoiseDBValidation") << " Bad percentual occupancy for the first 10 wires! Contents : "
                                         << (*channel).getContents();
       if ((*channel).getBin() == 2)
         LogWarning("NoiseDBValidation") << " Bad percentual occupancy for the middle wires! Contents : "
                                         << (*channel).getContents();
       if ((*channel).getBin() == 3)
         LogWarning("NoiseDBValidation") << " Bad percentual occupancy for the last 10 wires! Contents : "
                                         << (*channel).getContents();
     }
   }
 }

Member Data Documentation

DQMStore* DTnoiseDBValidation::dbe_

private

Definition at line 48 of file DTnoiseDBValidation.h.

MonitorElement* DTnoiseDBValidation::diffHisto_

private

Definition at line 68 of file DTnoiseDBValidation.h.

std::string DTnoiseDBValidation::diffTestName_

private

Definition at line 52 of file DTnoiseDBValidation.h.

edm::ESHandle<DTGeometry> DTnoiseDBValidation::dtGeom_

private

Definition at line 58 of file DTnoiseDBValidation.h.

std::string DTnoiseDBValidation::labelDB_

private

Definition at line 51 of file DTnoiseDBValidation.h.

std::string DTnoiseDBValidation::labelDBRef_

private

Definition at line 50 of file DTnoiseDBValidation.h.

MonitorElement* DTnoiseDBValidation::layerHisto_

private

Definition at line 72 of file DTnoiseDBValidation.h.

std::string DTnoiseDBValidation::layerTestName_

private

Definition at line 52 of file DTnoiseDBValidation.h.

std::map<DTChamberId, MonitorElement *> DTnoiseDBValidation::noiseHistoMap_

private

Definition at line 73 of file DTnoiseDBValidation.h.

const DTStatusFlag* DTnoiseDBValidation::noiseMap_

private

Definition at line 61 of file DTnoiseDBValidation.h.

const DTStatusFlag* DTnoiseDBValidation::noiseRefMap_

private

Definition at line 62 of file DTnoiseDBValidation.h.

int DTnoiseDBValidation::noisyCellsRef_

private

Definition at line 65 of file DTnoiseDBValidation.h.

int DTnoiseDBValidation::noisyCellsValid_

private

Definition at line 66 of file DTnoiseDBValidation.h.

std::string DTnoiseDBValidation::outputFileName_

private

Definition at line 55 of file DTnoiseDBValidation.h.

bool DTnoiseDBValidation::outputMEsInRootFile_

private

Definition at line 54 of file DTnoiseDBValidation.h.

MonitorElement* DTnoiseDBValidation::sectorHisto_

private

Definition at line 71 of file DTnoiseDBValidation.h.

std::string DTnoiseDBValidation::sectorTestName_

private

Definition at line 52 of file DTnoiseDBValidation.h.

MonitorElement* DTnoiseDBValidation::stationHisto_

private

Definition at line 70 of file DTnoiseDBValidation.h.

std::string DTnoiseDBValidation::stationTestName_

private

Definition at line 52 of file DTnoiseDBValidation.h.

MonitorElement* DTnoiseDBValidation::wheelHisto_

private

Definition at line 69 of file DTnoiseDBValidation.h.

std::string DTnoiseDBValidation::wheelTestName_

private

Definition at line 52 of file DTnoiseDBValidation.h.

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation