#include <SiStripCorrelateBadStripAndNoise.h>

Inheritance diagram for SiStripCorrelateBadStripAndNoise:

Public Member Functions
	SiStripCorrelateBadStripAndNoise (const edm::ParameterSet &)

	~SiStripCorrelateBadStripAndNoise () override=default

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

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

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

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

void	correlateWithNoise (const uint32_t &detid, const TrackerTopology &tTopo, const uint32_t &firstStrip, const uint32_t &range)

void	DoAnalysis (const edm::EventSetup &)

void	endJob () override

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

TH2F *	getHisto (const long unsigned int &index)

void	getHistos (const uint32_t &detid, const TrackerTopology &tTopo, std::vector< TH2F *> &histos)

float	getMeanNoise (const SiStripNoises::Range &noiseRange, const uint32_t &first, const uint32_t &range)

void	iterateOnBadStrips (const uint32_t &detid, const TrackerTopology &tTopo, const TrackerGeometry &tGeom, SiStripQuality::Range &sqrange)

void	iterateOnDets (const TrackerTopology &tTopo, const TrackerGeometry &tGeom)

Private Attributes
TFile *	file

const SiStripNoises *	noises_

edm::ESGetToken< SiStripNoises, SiStripNoisesRcd >	noiseToken_

edm::ESWatcher< SiStripNoisesRcd >	noiseWatcher_

const SiStripQuality *	quality_

edm::ESGetToken< SiStripQuality, SiStripQualityRcd >	qualityToken_

edm::ESWatcher< SiStripQualityRcd >	qualityWatcher_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

TrackerMap *	tkmap

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

std::vector< TH2F * >	vTH2

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)

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 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<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)

Detailed Description

DQM/SiStripMonitorSummary/plugins/SiStripCorrelateBadStripAndNoise.cc

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 46 of file SiStripCorrelateBadStripAndNoise.h.

Constructor & Destructor Documentation

◆ SiStripCorrelateBadStripAndNoise()

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

explicit

Definition at line 8 of file SiStripCorrelateBadStripAndNoise.cc.

     : qualityToken_{esConsumes<edm::Transition::BeginRun>()},
       noiseToken_{esConsumes<edm::Transition::BeginRun>()},
       tTopoToken_{esConsumes<edm::Transition::BeginRun>()},
       tkGeomToken_{esConsumes<edm::Transition::BeginRun>()} {
   // now do what ever initialization is needed
   file = new TFile("correlTest.root", "RECREATE");
   tkmap = new TrackerMap();
 }

◆ ~SiStripCorrelateBadStripAndNoise()

SiStripCorrelateBadStripAndNoise::~SiStripCorrelateBadStripAndNoise ( )

overridedefault

Member Function Documentation

◆ analyze()

void SiStripCorrelateBadStripAndNoise::analyze	(	const edm::Event &	,
		const edm::EventSetup &
	)

inlineoverrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 53 of file SiStripCorrelateBadStripAndNoise.h.

53 {};

◆ beginRun()

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

overrideprivate

Definition at line 18 of file SiStripCorrelateBadStripAndNoise.cc.

References edm::ESWatcher< T >::check(), DoAnalysis(), edm::EventSetup::getData(), HLT_2022v15_cff::newQuality, noises_, noiseToken_, noiseWatcher_, quality_, qualityToken_, and qualityWatcher_.

                                                                                           {
   auto newNoise = noiseWatcher_.check(es);
   auto newQuality = qualityWatcher_.check(es);
   if (newNoise || newQuality) {
     edm::LogInfo("") << "[SiStripCorrelateBadStripAndNoise::beginRun]" << std::endl;
 
     quality_ = &es.getData(qualityToken_);
     noises_ = &es.getData(noiseToken_);
 
     DoAnalysis(es);
   }
 }

◆ correlateWithNoise()

void SiStripCorrelateBadStripAndNoise::correlateWithNoise	(	const uint32_t &	detid,
		const TrackerTopology &	tTopo,
		const uint32_t &	firstStrip,
		const uint32_t &	range
	)

private

Definition at line 80 of file SiStripCorrelateBadStripAndNoise.cc.

References HcalObjRepresent::Fill(), TrackerMap::fillc(), getHistos(), getMeanNoise(), SiStripNoises::getRange(), combine::histos, mps_fire::i, noises_, FastTimerService_cff::range, and tkmap.

Referenced by iterateOnBadStrips().

                                                                                  {
   std::vector<TH2F *> histos;
 
   SiStripNoises::Range noiseRange = noises_->getRange(detid);
   edm::LogInfo("Domenico") << "detid " << detid << " first " << firstStrip << " range " << range;
   float meanAPVNoise = getMeanNoise(noiseRange, firstStrip / 128, 128);
 
   // float meanNoiseHotStrips=getMeanNoise(noiseRange,firstStrip,range);
   for (size_t theStrip = firstStrip; theStrip < firstStrip + range; theStrip++) {
     float meanNoiseHotStrips = getMeanNoise(noiseRange, theStrip, 1);
 
     // Get the histogram for this detid
     getHistos(detid, tTopo, histos);
     float yvalue = range < 21 ? 1. * range : 21;
 
     for (size_t i = 0; i < histos.size(); ++i)
       histos[i]->Fill(meanNoiseHotStrips / meanAPVNoise - 1., yvalue);
 
     if (meanNoiseHotStrips / meanAPVNoise - 1. < -0.3)
       tkmap->fillc(detid, 0xFF0000);
     else
       tkmap->fillc(detid, 0x0000FF);
   }
 }

◆ DoAnalysis()

void SiStripCorrelateBadStripAndNoise::DoAnalysis ( const edm::EventSetup & es )

private

Definition at line 31 of file SiStripCorrelateBadStripAndNoise.cc.

References edm::EventSetup::getData(), iterateOnDets(), tkGeomToken_, and tTopoToken_.

Referenced by beginRun().

                                                                          {
   // Loop on quality bad stirps
   // for each strip, look at the noise
   // evalaute the mean apv noise and the ratio among strip noise and
   // meanApvNoise put the value in the histo in terms of ratio Vs percentage of
   // badStrips per APV
 
   // Fill an histo per subdet and layer (and plus && minus for TEC/TID)
   edm::LogInfo("") << "[Doanalysis]";
   iterateOnDets(es.getData(tTopoToken_), es.getData(tkGeomToken_));
 }

◆ endJob()

void SiStripCorrelateBadStripAndNoise::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 160 of file SiStripCorrelateBadStripAndNoise.cc.

References file, mps_fire::i, TrackerMap::save(), tkmap, and vTH2.

                                               {
   for (size_t i = 0; i < vTH2.size(); i++)
     if (vTH2[i] != nullptr)
       vTH2[i]->Write();
 
   file->Write();
   file->Close();
 
   tkmap->save(true, 0, 0, "testTkMap.png");
 }

◆ endRun()

void SiStripCorrelateBadStripAndNoise::endRun	(	const edm::Run &	run,
		const edm::EventSetup &	es
	)

inlineoverrideprivate

Definition at line 54 of file SiStripCorrelateBadStripAndNoise.h.

54 {};

◆ getHisto()

TH2F * SiStripCorrelateBadStripAndNoise::getHisto ( const long unsigned int & index )

private

Definition at line 146 of file SiStripCorrelateBadStripAndNoise.cc.

References Skims_PA_cff::name, and vTH2.

Referenced by getHistos().

                                                                                {
   if (vTH2.size() < index + 1)
     vTH2.resize(index + 1, nullptr);
 
   if (vTH2[index] == nullptr) {
     char name[128];
     sprintf(name, "%lu", index);
     edm::LogInfo("") << "[getHisto] creating index " << index << std::endl;
     vTH2[index] = new TH2F(name, name, 50, -2., 2., 21, 0.5, 21.5);
   }
 
   return vTH2[index];
 }

◆ getHistos()

void SiStripCorrelateBadStripAndNoise::getHistos	(	const uint32_t &	detid,
		const TrackerTopology &	tTopo,
		std::vector< TH2F *> &	histos
	)

private

Definition at line 118 of file SiStripCorrelateBadStripAndNoise.cc.

References a, getHisto(), combine::histos, TrackerTopology::tecSide(), TrackerTopology::tecWheel(), TrackerTopology::tibLayer(), TrackerTopology::tidSide(), TrackerTopology::tidWheel(), and TrackerTopology::tobLayer().

Referenced by correlateWithNoise().

                                                                             {
   histos.clear();
 
   int subdet = -999;
   int component = -999;
   SiStripDetId a(detid);
   if (a.subdetId() == 3) {
     subdet = 0;
     component = tTopo.tibLayer(detid);
   } else if (a.subdetId() == 4) {
     subdet = 1;
     component = tTopo.tidSide(detid) == 2 ? tTopo.tidWheel(detid) : tTopo.tidWheel(detid) + 3;
   } else if (a.subdetId() == 5) {
     subdet = 2;
     component = tTopo.tobLayer(detid);
   } else if (a.subdetId() == 6) {
     subdet = 3;
     component = tTopo.tecSide(detid) == 2 ? tTopo.tecWheel(detid) : tTopo.tecWheel(detid) + 9;
   }
 
   int index = 100 + subdet * 100 + component;
 
   histos.push_back(getHisto(subdet));
   histos.push_back(getHisto(index));
 }

◆ getMeanNoise()

float SiStripCorrelateBadStripAndNoise::getMeanNoise	(	const SiStripNoises::Range &	noiseRange,
		const uint32_t &	first,
		const uint32_t &	range
	)

private

Definition at line 108 of file SiStripCorrelateBadStripAndNoise.cc.

References SiStripNoises::getNoise(), SiStripPI::mean, noises_, and FastTimerService_cff::range.

Referenced by correlateWithNoise().

                                                                             {
   float mean = 0;
   for (size_t istrip = firstStrip; istrip < firstStrip + range; istrip++) {
     mean += noises_->getNoise(istrip, noiseRange);
   }
   return mean / (1. * range);
 }

◆ iterateOnBadStrips()

void SiStripCorrelateBadStripAndNoise::iterateOnBadStrips	(	const uint32_t &	detid,
		const TrackerTopology &	tTopo,
		const TrackerGeometry &	tGeom,
		SiStripQuality::Range &	sqrange
	)

private

Definition at line 55 of file SiStripCorrelateBadStripAndNoise.cc.

References correlateWithNoise(), SiStripBadStrip::decode(), SiStripBadStrip::data::firstStrip, TrackerGeometry::idToDet(), quality_, SiStripBadStrip::data::range, and FastTimerService_cff::range.

Referenced by iterateOnDets().

                                                                                         {
   float percentage = 0;
   for (int it = 0; it < sqrange.second - sqrange.first; it++) {
     unsigned int firstStrip = quality_->decode(*(sqrange.first + it)).firstStrip;
     unsigned int range = quality_->decode(*(sqrange.first + it)).range;
 
     correlateWithNoise(detid, tTopo, firstStrip, range);
 
     edm::LogInfo("range") << range;
     percentage += range;
   }
   if (percentage != 0)
     percentage /= dynamic_cast<const StripGeomDetUnit *>(tGeom.idToDet(detid))->specificTopology().nstrips();
   if (percentage > 1)
     edm::LogError("SiStripQualityStatistics") << "PROBLEM detid " << detid << " value " << percentage << std::endl;
 
   //------- Global Statistics on percentage of bad components along the IOVs
   //------//
   if (percentage != 0)
     edm::LogInfo("") << "percentage " << detid << " " << percentage;
 }

◆ iterateOnDets()

void SiStripCorrelateBadStripAndNoise::iterateOnDets	(	const TrackerTopology &	tTopo,
		const TrackerGeometry &	tGeom
	)

private

Definition at line 43 of file SiStripCorrelateBadStripAndNoise.cc.

References SiStripBadStrip::getDataVectorBegin(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), iterateOnBadStrips(), and quality_.

Referenced by DoAnalysis().

                                                                                                                {
   const auto rbegin = quality_->getRegistryVectorBegin();
   const auto rend = quality_->getRegistryVectorEnd();
   for (auto rp = rbegin; rp != rend; ++rp) {
     const uint32_t detid = rp->detid;
 
     auto sqrange =
         SiStripQuality::Range(quality_->getDataVectorBegin() + rp->ibegin, quality_->getDataVectorBegin() + rp->iend);
     iterateOnBadStrips(detid, tTopo, tGeom, sqrange);
   }
 }