#include <SiStripCorrelateNoise.h>

Inheritance diagram for SiStripCorrelateNoise:

Classes
struct	Data

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

	~SiStripCorrelateNoise () 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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

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

bool	wantsStreamRuns () const noexcept

	~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

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)

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

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

void	checkGainCache (const edm::EventSetup &es)

void	DoAnalysis (const edm::EventSetup &, const SiStripNoises &, SiStripNoises &)

void	DoPlots ()

void	endJob () override

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

float	getGainRatio (const uint32_t &detid, const uint16_t &apv)

TH1F *	getHisto (const long unsigned int &index)

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

Private Attributes
bool	equalGain

TFile *	file

edm::ESGetToken< SiStripApvGain, SiStripApvGainRcd >	gainToken_

edm::ESWatcher< SiStripApvGainRcd >	gainWatcher_

std::unique_ptr< SiStripApvGain >	newGain

edm::ESGetToken< SiStripNoises, SiStripNoisesRcd >	noiseToken_

edm::ESWatcher< SiStripNoisesRcd >	noiseWatcher_

std::unique_ptr< SiStripApvGain >	oldGain

std::unique_ptr< SiStripNoises >	refNoise

uint32_t	theRun

TrackerMap *	tkmap

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

std::vector< TH1F * >	vTH1

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 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/SiStripCorrelateNoise.cc

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

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

Definition at line 46 of file SiStripCorrelateNoise.h.

Constructor & Destructor Documentation

◆ SiStripCorrelateNoise()

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

explicit

Definition at line 8 of file SiStripCorrelateNoise.cc.

     : noiseToken_{esConsumes<edm::Transition::BeginRun>()},
       gainToken_{esConsumes<edm::Transition::BeginRun>()},
       tTopoToken_{esConsumes<edm::Transition::BeginRun>()},
       refNoise(nullptr),
       oldGain(nullptr),
       newGain(nullptr) {
   // now do what ever initialization is needed
   file = new TFile("correlTest.root", "RECREATE");
 
   file->cd();
   tkmap = new TrackerMap();
 }

◆ ~SiStripCorrelateNoise()

SiStripCorrelateNoise::~SiStripCorrelateNoise ( )

overridedefault

Member Function Documentation

◆ analyze()

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

inlineoverrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 53 of file SiStripCorrelateNoise.h.

53 {};

◆ beginRun()

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

overrideprivate

Definition at line 22 of file SiStripCorrelateNoise.cc.

References edm::ESWatcher< T >::check(), checkGainCache(), DeadROC_duringRun::dir, DoAnalysis(), DoPlots(), file, edm::EventSetup::getData(), eostools::move(), noiseToken_, noiseWatcher_, refNoise, writedatasetfile::run, and theRun.

                                                                                {
   if (noiseWatcher_.check(es)) {
     auto aNoise = std::make_unique<SiStripNoises>(es.getData(noiseToken_));
     // Check if gain is the same from one noise iov to the other, otherwise cache
     // the new gain (and the old one) to rescale
     checkGainCache(es);
     if (refNoise) {
       char dir[128];
       theRun = run.run();
       sprintf(dir, "Run_%d", theRun);
       file->cd("");
       file->mkdir(dir);
       file->cd(dir);
       DoAnalysis(es, *aNoise, *refNoise);
       DoPlots();
     }
     refNoise = std::move(aNoise);
   }
 }

◆ checkGainCache()

void SiStripCorrelateNoise::checkGainCache ( const edm::EventSetup & es )

private

Definition at line 42 of file SiStripCorrelateNoise.cc.

References edm::ESWatcher< T >::check(), equalGain, gainToken_, gainWatcher_, edm::EventSetup::getData(), eostools::move(), newGain, and oldGain.

Referenced by beginRun().

                                                                   {
   equalGain = true;
   if (gainWatcher_.check(es)) {
     if (oldGain) {
       equalGain = false;
     }
     oldGain = std::move(newGain);
     newGain = std::make_unique<SiStripApvGain>(es.getData(gainToken_));
   }
 }

◆ DoAnalysis()

void SiStripCorrelateNoise::DoAnalysis	(	const edm::EventSetup &	es,
		const SiStripNoises &	_Noise,
		SiStripNoises &	refNoise
	)

private

Definition at line 83 of file SiStripCorrelateNoise.cc.

References universalConfigTemplate::collection, equalGain, TrackerMap::fill(), HcalObjRepresent::Fill(), edm::EventSetup::getData(), getGainRatio(), getHistos(), combine::histos, mps_fire::i, SiStripCondTypes::Noise, refNoise, nano_mu_digi_cff::strip, tkmap, tTopoToken_, and relativeConstraints::value.

Referenced by beginRun().

                                                                 {
   SiStripNoises Noise = _Noise;
   typedef std::vector<SiStripNoises::ratioData> collection;
   collection divNoise = Noise / refNoise;
 
   edm::LogInfo("") << "[Doanalysis]";
 
   const auto tTopo = &es.getData(tTopoToken_);
 
   std::vector<TH1F *> histos;
 
   collection::const_iterator iter = divNoise.begin();
   collection::const_iterator iterE = divNoise.end();
 
   float value;
   float gainRatio = 1.;
   // Divide result by d
   for (; iter != iterE; ++iter) {
     getHistos(iter->detid, tTopo, histos);
 
     size_t strip = 0, stripE = iter->values.size();
     size_t apvNb = 7;
 
     for (; strip < stripE; ++strip) {
       if (!equalGain && strip / 128 != apvNb) {
         apvNb = strip / 128;
         if (apvNb < 6)
           gainRatio = getGainRatio(iter->detid, apvNb);
         else
           edm::LogInfo("") << "[Doanalysis] detid " << iter->detid << " strip " << strip << " apvNb " << apvNb;
       }
       // edm::LogInfo("") << "[Doanalysis] detid " << iter->detid << " strip "
       // << strip << " value " << iter->values[strip];
       value = iter->values[strip] * gainRatio;
       tkmap->fill(iter->detid, value);
       for (size_t i = 0; i < histos.size(); ++i)
         histos[i]->Fill(value);
     }
   }
 }

◆ DoPlots()

void SiStripCorrelateNoise::DoPlots ( )

private

Definition at line 53 of file SiStripCorrelateNoise.cc.

References correctionTermsCaloMet_cff::C, DeadROC_duringRun::dir, file, mps_fire::i, TrackerMap::save(), theRun, tkmap, and vTH1.

Referenced by beginRun().

                                     {
   TCanvas *C = new TCanvas();
   C->Divide(2, 2);
 
   char outName[128];
   sprintf(outName, "Run_%d.png", theRun);
   for (size_t i = 0; i < vTH1.size(); i++)
     if (vTH1[i] != nullptr) {
       if (i % 100 == 0) {
         C->cd(i / 100);
         vTH1[i]->SetLineColor(i / 100);
         vTH1[i]->Draw();
         C->cd(i / 100)->SetLogy();
       }
       vTH1[i]->Write();
     }
 
   C->Print(outName);
   delete C;
 
   vTH1.clear();
   file->cd("");
 
   char dir[128];
   sprintf(dir, "Run_%d_TkMap.png", theRun);
   tkmap->save(false, 0, 5, dir);
   delete tkmap;
   tkmap = new TrackerMap();
 }

◆ endJob()

void SiStripCorrelateNoise::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 187 of file SiStripCorrelateNoise.cc.

References file.

                                    {
   file->Write();
   file->Close();
 }

◆ endRun()

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

inlineoverrideprivate

Definition at line 54 of file SiStripCorrelateNoise.h.

54 {};

◆ getGainRatio()

float SiStripCorrelateNoise::getGainRatio	(	const uint32_t &	detid,
		const uint16_t &	apv
	)

private

Definition at line 126 of file SiStripCorrelateNoise.cc.

References newGain, and oldGain.

Referenced by DoAnalysis().

                                                                                     {
   SiStripApvGain::Range oldRange = oldGain->getRange(detid);
   SiStripApvGain::Range newRange = newGain->getRange(detid);
 
   if (oldRange.first == oldRange.second || newRange.first == newRange.second)
     return 1.;
 
   return oldGain->getApvGain(apv, oldRange) / newGain->getApvGain(apv, newRange);
 }

◆ getHisto()

TH1F * SiStripCorrelateNoise::getHisto ( const long unsigned int & index )

private

Definition at line 164 of file SiStripCorrelateNoise.cc.

References Skims_PA_cff::name, AlCaHLTBitMon_QueryRunRegistry::string, theRun, and vTH1.

Referenced by getHistos().

                                                                     {
   if (vTH1.size() < index + 1)
     vTH1.resize(index + 1, nullptr);
 
   if (vTH1[index] == nullptr) {
     char name[128];
     std::string SubD;
     if (index < 200)
       SubD = "TIB";
     else if (index < 300)
       SubD = "TID";
     else if (index < 400)
       SubD = "TOB";
     else
       SubD = "TEC";
     sprintf(name, "%d_%lu__%s", theRun, index, SubD.c_str());
     edm::LogInfo("") << "[getHisto] creating index " << index << std::endl;
     vTH1[index] = new TH1F(name, name, 200, -0.5, 10.5);
   }
 
   return vTH1[index];
 }

◆ getHistos()

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

private

Definition at line 136 of file SiStripCorrelateNoise.cc.

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

Referenced by DoAnalysis().

                                                                  {
   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(100 + 100 * subdet));
   histos.push_back(getHisto(index));
 }

Member Data Documentation

◆ equalGain

bool SiStripCorrelateNoise::equalGain

private

Definition at line 82 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache(), and DoAnalysis().

◆ file

TFile* SiStripCorrelateNoise::file

private

Definition at line 84 of file SiStripCorrelateNoise.h.

Referenced by beginRun(), python.rootplot.utilities.RootFile::cd(), DoPlots(), endJob(), ztee.GZipLog::finish(), ztee.GZipLog::flush_block(), python.rootplot.utilities.RootFile::get(), tfile.TFileService::stop(), core.VertexHistograms.VertexHistograms::write(), core.TreeAnalyzerNumpy.TreeAnalyzerNumpy::write(), and ztee.GZipLog::write_block().

◆ gainToken_

edm::ESGetToken<SiStripApvGain, SiStripApvGainRcd> SiStripCorrelateNoise::gainToken_

private

Definition at line 76 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache().

◆ gainWatcher_

edm::ESWatcher<SiStripApvGainRcd> SiStripCorrelateNoise::gainWatcher_

private

Definition at line 74 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache().

◆ newGain

std::unique_ptr<SiStripApvGain> SiStripCorrelateNoise::newGain

private

Definition at line 81 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache(), and getGainRatio().

◆ noiseToken_

edm::ESGetToken<SiStripNoises, SiStripNoisesRcd> SiStripCorrelateNoise::noiseToken_

private

Definition at line 75 of file SiStripCorrelateNoise.h.

Referenced by beginRun().

◆ noiseWatcher_

edm::ESWatcher<SiStripNoisesRcd> SiStripCorrelateNoise::noiseWatcher_

private

Definition at line 73 of file SiStripCorrelateNoise.h.

Referenced by beginRun().

◆ oldGain

std::unique_ptr<SiStripApvGain> SiStripCorrelateNoise::oldGain

private

Definition at line 81 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache(), and getGainRatio().

◆ refNoise

std::unique_ptr<SiStripNoises> SiStripCorrelateNoise::refNoise

private

Definition at line 80 of file SiStripCorrelateNoise.h.

Referenced by beginRun(), and DoAnalysis().

◆ theRun

uint32_t SiStripCorrelateNoise::theRun

private

Definition at line 79 of file SiStripCorrelateNoise.h.

Referenced by beginRun(), DoPlots(), and getHisto().

◆ tkmap

TrackerMap* SiStripCorrelateNoise::tkmap

private

Definition at line 87 of file SiStripCorrelateNoise.h.

Referenced by DoAnalysis(), and DoPlots().

◆ tTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiStripCorrelateNoise::tTopoToken_

private

Definition at line 77 of file SiStripCorrelateNoise.h.

Referenced by DoAnalysis().

◆ vTH1

std::vector<TH1F *> SiStripCorrelateNoise::vTH1

private

Definition at line 85 of file SiStripCorrelateNoise.h.

Referenced by DoPlots(), and getHisto().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SiStripCorrelateNoise()

◆ ~SiStripCorrelateNoise()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ checkGainCache()

◆ DoAnalysis()

◆ DoPlots()

◆ endJob()

◆ endRun()

◆ getGainRatio()

◆ getHisto()

◆ getHistos()

Member Data Documentation

◆ equalGain

◆ file

◆ gainToken_

◆ gainWatcher_

◆ newGain

◆ noiseToken_

◆ noiseWatcher_

◆ oldGain

◆ refNoise

◆ theRun

◆ tkmap

◆ tTopoToken_

◆ vTH1