#include <SiStripCorrelateNoise.h>

Inheritance diagram for SiStripCorrelateNoise:

Classes
struct	Data

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

	~SiStripCorrelateNoise () override

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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

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

unsigned long long	getGainCache (const edm::EventSetup &eSetup)

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)

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

unsigned long long	getNoiseCache (const edm::EventSetup &eSetup)

Private Attributes
unsigned long long	cacheID_gain

unsigned long long	cacheID_noise

bool	equalGain

TFile *	file

SiStripDetInfoFileReader *	fr

edm::ESHandle< SiStripApvGain >	gainHandle_

SiStripApvGain *	newGain

edm::ESHandle< SiStripNoises >	noiseHandle_

SiStripApvGain *	oldGain

SiStripNoises *	refNoise

uint32_t	theRun

TrackerMap *	tkmap

std::vector< TH1F * >	vTH1

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

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

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

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

DQM/SiStripMonitorSummary/plugins/SiStripCorrelateNoise.cc

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

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

Definition at line 49 of file SiStripCorrelateNoise.h.

Constructor & Destructor Documentation

◆ SiStripCorrelateNoise()

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

explicit

Definition at line 9 of file SiStripCorrelateNoise.cc.

     : refNoise(nullptr), oldGain(nullptr), newGain(nullptr), cacheID_noise(0xFFFFFFFF), cacheID_gain(0xFFFFFFFF) {
   // now do what ever initialization is needed
   file = new TFile("correlTest.root", "RECREATE");
  
   file->cd();
   tkmap = new TrackerMap();
 }

References file, and tkmap.

◆ ~SiStripCorrelateNoise()

SiStripCorrelateNoise::~SiStripCorrelateNoise ( )

override

Definition at line 18 of file SiStripCorrelateNoise.cc.

18 {}

Member Function Documentation

◆ analyze()

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

inlineoverrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 56 of file SiStripCorrelateNoise.h.

56 {};

◆ beginRun()

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

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 22 of file SiStripCorrelateNoise.cc.

                                                                                {
   if (getNoiseCache(es) == cacheID_noise)
     return;
  
   edm::LogInfo("") << "[SiStripCorrelateNoise::beginRun]  cacheID_noise " << cacheID_noise << std::endl;
  
   es.get<SiStripNoisesRcd>().get(noiseHandle_);
   SiStripNoises *aNoise = new SiStripNoises(*noiseHandle_.product());
  
   // 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 (cacheID_noise != 0xFFFFFFFF) {
     char dir[128];
     theRun = run.run();
     sprintf(dir, "Run_%d", theRun);
     file->cd("");
     file->mkdir(dir);
     file->cd(dir);
     DoAnalysis(es, *noiseHandle_.product(), *refNoise);
     DoPlots();
   }
  
   cacheID_noise = getNoiseCache(es);
   if (refNoise != nullptr)
     delete refNoise;
   refNoise = aNoise;
 }

References cacheID_noise, checkGainCache(), DeadROC_duringRun::dir, DoAnalysis(), DoPlots(), file, edm::EventSetup::get(), get, getNoiseCache(), noiseHandle_, edm::ESHandle< T >::product(), refNoise, writedatasetfile::run, and theRun.

◆ checkGainCache()

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

private

Definition at line 53 of file SiStripCorrelateNoise.cc.

                                                                   {
   equalGain = true;
   if (getGainCache(es) != cacheID_gain) {
     es.get<SiStripApvGainRcd>().get(gainHandle_);
     if (oldGain != nullptr)
       delete oldGain;
  
     oldGain = newGain;
     newGain = new SiStripApvGain(*gainHandle_.product());
  
     if (cacheID_gain != 0xFFFFFFFF)
       equalGain = false;
     cacheID_gain = getGainCache(es);
     edm::LogInfo("") << "[SiStripCorrelateNoise::checkGainCache]  cacheID_gain " << cacheID_gain << std::endl;
   }
 }

References cacheID_gain, equalGain, gainHandle_, edm::EventSetup::get(), get, getGainCache(), newGain, oldGain, and edm::ESHandle< T >::product().

Referenced by beginRun().

◆ DoAnalysis()

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

private

Definition at line 100 of file SiStripCorrelateNoise.cc.

                                                                 {
   SiStripNoises Noise = _Noise;
   typedef std::vector<SiStripNoises::ratioData> collection;
   collection divNoise = Noise / refNoise;
  
   edm::LogInfo("") << "[Doanalysis]";
  
   // Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology *const tTopo = tTopoHandle.product();
  
   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);
     }
   }
 }

References universalConfigTemplate::collection, equalGain, TrackerMap::fill(), HcalObjRepresent::Fill(), edm::EventSetup::get(), get, getGainRatio(), getHistos(), combine::histos, mps_fire::i, DigiDM_cff::Noise, edm::ESHandle< T >::product(), refNoise, digitizers_cfi::strip, tkmap, and relativeConstraints::value.

Referenced by beginRun().

◆ DoPlots()

void SiStripCorrelateNoise::DoPlots ( )

private

Definition at line 70 of file SiStripCorrelateNoise.cc.

                                     {
   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();
 }

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

Referenced by beginRun().

◆ endJob()

void SiStripCorrelateNoise::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 217 of file SiStripCorrelateNoise.cc.

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

References file.

◆ getGainCache()

unsigned long long SiStripCorrelateNoise::getGainCache ( const edm::EventSetup & eSetup )

inlineprivate

Definition at line 67 of file SiStripCorrelateNoise.h.

                                                                {
     return eSetup.get<SiStripApvGainRcd>().cacheIdentifier();
   }

References edm::EventSetup::get().

Referenced by checkGainCache().

◆ getGainRatio()

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

private

Definition at line 146 of file SiStripCorrelateNoise.cc.

                                                                                     {
   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);
 }

References SiStripApvGain::getApvGain(), SiStripApvGain::getRange(), newGain, and oldGain.

Referenced by DoAnalysis().

◆ getHisto()

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

private

Definition at line 194 of file SiStripCorrelateNoise.cc.

                                                                     {
   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];
 }

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

Referenced by getHistos().

◆ getHistos()

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

private

Definition at line 166 of file SiStripCorrelateNoise.cc.

                                                                  {
   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));
 }

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

Referenced by DoAnalysis().

◆ getMeanNoise()

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

private

Definition at line 156 of file SiStripCorrelateNoise.cc.

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

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

◆ getNoiseCache()

unsigned long long SiStripCorrelateNoise::getNoiseCache ( const edm::EventSetup & eSetup )

inlineprivate

Definition at line 64 of file SiStripCorrelateNoise.h.

                                                                 {
     return eSetup.get<SiStripNoisesRcd>().cacheIdentifier();
   }

References edm::EventSetup::get().

Referenced by beginRun().

Member Data Documentation

◆ cacheID_gain

unsigned long long SiStripCorrelateNoise::cacheID_gain

private

Definition at line 98 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache().

◆ cacheID_noise

unsigned long long SiStripCorrelateNoise::cacheID_noise

private

Definition at line 97 of file SiStripCorrelateNoise.h.

Referenced by beginRun().

◆ equalGain

bool SiStripCorrelateNoise::equalGain

private

Definition at line 90 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache(), and DoAnalysis().

◆ file

TFile* SiStripCorrelateNoise::file

private

Definition at line 92 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(), SiStripCorrelateNoise(), tfile.TFileService::stop(), core.VertexHistograms.VertexHistograms::write(), core.TreeAnalyzerNumpy.TreeAnalyzerNumpy::write(), and ztee.GZipLog::write_block().

◆ fr

SiStripDetInfoFileReader* SiStripCorrelateNoise::fr

private

Definition at line 82 of file SiStripCorrelateNoise.h.

◆ gainHandle_

edm::ESHandle<SiStripApvGain> SiStripCorrelateNoise::gainHandle_

private

Definition at line 83 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache().

◆ newGain

SiStripApvGain * SiStripCorrelateNoise::newGain

private

Definition at line 89 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache(), and getGainRatio().

◆ noiseHandle_

edm::ESHandle<SiStripNoises> SiStripCorrelateNoise::noiseHandle_

private

Definition at line 84 of file SiStripCorrelateNoise.h.

Referenced by beginRun(), and getMeanNoise().

◆ oldGain

SiStripApvGain* SiStripCorrelateNoise::oldGain

private

Definition at line 89 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache(), and getGainRatio().

◆ refNoise

SiStripNoises* SiStripCorrelateNoise::refNoise

private

Definition at line 87 of file SiStripCorrelateNoise.h.

Referenced by beginRun(), and DoAnalysis().

◆ theRun

uint32_t SiStripCorrelateNoise::theRun

private

Definition at line 86 of file SiStripCorrelateNoise.h.

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

◆ tkmap

TrackerMap* SiStripCorrelateNoise::tkmap

private

Definition at line 95 of file SiStripCorrelateNoise.h.

Referenced by DoAnalysis(), DoPlots(), and SiStripCorrelateNoise().

◆ vTH1

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

private

Definition at line 93 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()

◆ getGainCache()

◆ getGainRatio()

◆ getHisto()

◆ getHistos()

◆ getMeanNoise()

◆ getNoiseCache()

Member Data Documentation

◆ cacheID_gain

◆ cacheID_noise

◆ equalGain

◆ file

◆ fr

◆ gainHandle_

◆ newGain

◆ noiseHandle_

◆ oldGain

◆ refNoise

◆ theRun

◆ tkmap

◆ vTH1