#include <DQM/SiStripMonitorSummary/plugins/SiStripCorrelateBadStripAndNoise.cc>

Inheritance diagram for SiStripCorrelateBadStripAndNoise:

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

	~SiStripCorrelateBadStripAndNoise () 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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)

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

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)

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

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

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

void	iterateOnDets (const TrackerTopology *tTopo)

Private Attributes
unsigned long long	cacheID_noise

unsigned long long	cacheID_quality

TFile *	file

SiStripDetInfoFileReader *	fr

edm::ESHandle< SiStripNoises >	noiseHandle_

edm::ESHandle< SiStripQuality >	qualityHandle_

TrackerMap *	tkmap

std::vector< TH2F * >	vTH2

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	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 ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, 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

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

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

Definition at line 49 of file SiStripCorrelateBadStripAndNoise.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 8 of file SiStripCorrelateBadStripAndNoise.cc.

References file, fr, Utilities::operator, and tkmap.

                                                                                                 :
   cacheID_quality(0xFFFFFFFF),cacheID_noise(0xFFFFFFFF)
 {
   //now do what ever initialization is needed
   if(!edm::Service<SiStripDetInfoFileReader>().isAvailable()){
     edm::LogError("TkLayerMap") << 
       "\n------------------------------------------"
       "\nUnAvailable Service SiStripDetInfoFileReader: please insert in the configuration file an instance like"
       "\n\tprocess.SiStripDetInfoFileReader = cms.Service(\"SiStripDetInfoFileReader\")"
       "\n------------------------------------------";
   }
  
   fr=edm::Service<SiStripDetInfoFileReader>().operator->();
   file = new TFile("correlTest.root","RECREATE");
   tkmap = new TrackerMap();
 }

SiStripCorrelateBadStripAndNoise::~SiStripCorrelateBadStripAndNoise ( )

override

Definition at line 26 of file SiStripCorrelateBadStripAndNoise.cc.

27 {}

Member Function Documentation

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

inlineoverrideprivate

Definition at line 57 of file SiStripCorrelateBadStripAndNoise.h.

References DoAnalysis(), endJob(), getHisto(), getHistos(), and plotFactory::histos.

57 {};

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

overrideprivate

Definition at line 32 of file SiStripCorrelateBadStripAndNoise.cc.

References cacheID_noise, cacheID_quality, DoAnalysis(), edm::EventSetup::get(), getNoiseCache(), getQualityCache(), noiseHandle_, and qualityHandle_.

                                                                                     {
 
   if(getNoiseCache(es)==cacheID_noise && getQualityCache(es)==cacheID_quality)
     return;
   cacheID_noise=getNoiseCache(es);
   cacheID_quality=getQualityCache(es);
   
   edm::LogInfo("") << "[SiStripCorrelateBadStripAndNoise::beginRun] cacheID_quality " << cacheID_quality << " cacheID_noise " << cacheID_noise << std::endl; 
   
   es.get<SiStripQualityRcd>().get(qualityHandle_);
   es.get<SiStripNoisesRcd>().get(noiseHandle_);
 
   DoAnalysis(es);
 }

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

private

Definition at line 105 of file SiStripCorrelateBadStripAndNoise.cc.

References HcalObjRepresent::Fill(), TrackerMap::fillc(), getHistos(), getMeanNoise(), SiStripNoises::getRange(), plotFactory::histos, mps_fire::i, noiseHandle_, and tkmap.

Referenced by getQualityCache(), and iterateOnBadStrips().

                                                                                                                                                               {
 
   std::vector<TH2F *>histos;
 
   SiStripNoises::Range noiseRange = noiseHandle_->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);
   }
 }

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

private

Definition at line 48 of file SiStripCorrelateBadStripAndNoise.cc.

References edm::EventSetup::get(), iterateOnDets(), and edm::ESHandle< T >::product().

Referenced by analyze(), and beginRun().

                                                                    {
 
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   //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(tTopo);
 
 }

void SiStripCorrelateBadStripAndNoise::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 186 of file SiStripCorrelateBadStripAndNoise.cc.

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

Referenced by analyze().

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

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

private

Definition at line 171 of file SiStripCorrelateBadStripAndNoise.cc.

References dataset::name, and vTH2.

Referenced by analyze(), and 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];
 }

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

private

Definition at line 142 of file SiStripCorrelateBadStripAndNoise.cc.

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

Referenced by analyze(), and 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));
   
 }

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

private

Definition at line 132 of file SiStripCorrelateBadStripAndNoise.cc.

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

Referenced by correlateWithNoise(), and getQualityCache().

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

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

inlineprivate

Definition at line 64 of file SiStripCorrelateBadStripAndNoise.h.

References edm::EventSetup::get().

Referenced by beginRun().

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

SiStripNoisesRcd

Definition: SiStripCondDataRecords.h:34

edm::EventSetup::get

T get() const

Definition: EventSetup.h:68

unsigned long long SiStripCorrelateBadStripAndNoise::getQualityCache ( const edm::EventSetup & eSetup )

inlineprivate

Definition at line 65 of file SiStripCorrelateBadStripAndNoise.h.

References correlateWithNoise(), plotBeamSpotDB::first, edm::EventSetup::get(), getMeanNoise(), iterateOnBadStrips(), and iterateOnDets().

Referenced by beginRun().

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

SiStripQualityRcd

Definition: SiStripDependentRecords.h:40

edm::EventSetup::get

T get() const

Definition: EventSetup.h:68

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

private

Definition at line 82 of file SiStripCorrelateBadStripAndNoise.cc.

References correlateWithNoise(), SiStripBadStrip::decode(), SiStripBadStrip::data::firstStrip, fr, SiStripDetInfoFileReader::getNumberOfApvsAndStripLength(), qualityHandle_, and SiStripBadStrip::data::range.

Referenced by getQualityCache(), and iterateOnDets().

                                                                                                                                       {
 
   float percentage=0;
   for(int it=0;it<sqrange.second-sqrange.first;it++){
     unsigned int firstStrip=qualityHandle_->decode( *(sqrange.first+it) ).firstStrip;
     unsigned int range=qualityHandle_->decode( *(sqrange.first+it) ).range;
     
     correlateWithNoise(detid,tTopo,firstStrip,range);
 
     edm::LogInfo("range")<< range;
     percentage+=range;
   }
   if(percentage!=0)
     percentage/=128.*fr->getNumberOfApvsAndStripLength(detid).first;
   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;
 }

void SiStripCorrelateBadStripAndNoise::iterateOnDets ( const TrackerTopology * tTopo )

private

Definition at line 67 of file SiStripCorrelateBadStripAndNoise.cc.

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

Referenced by DoAnalysis(), and getQualityCache().

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