#include <DTNoiseTest.h>

Inheritance diagram for DTNoiseTest:

Public Member Functions
	DTNoiseTest (const edm::ParameterSet &ps)
	Constructor. More...

virtual	~DTNoiseTest ()
	Destructor. More...

Public Member Functions inherited from DQMEDHarvester
virtual void	analyze (edm::Event const &, edm::EventSetup const &) final

virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

virtual void	beginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDHarvester (void)

virtual void	endJob () final

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

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

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

Public Member Functions inherited from edm::one::EDAnalyzerBase
	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
void	bookHistos (DQMStore::IBooker &, const DTChamberId &ch, std::string folder, std::string histoTag)
	book the new ME More...

void	bookHistos (DQMStore::IBooker &, const DTLayerId &ch, int nWire, std::string folder, std::string histoTag)

void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override

void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)
	DQM Client Diagnostic. More...

std::string	getMEName (const DTChamberId &ch)
	Get the ME name. More...

std::string	getMEName (const DTLayerId &ly)

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

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

Private Attributes
bool	bookingdone

bool	debug

std::map< std::string, std::map< uint32_t, MonitorElement * > >	histos

edm::ESHandle< DTGeometry >	muonGeom

unsigned int	nLumiSegs

edm::ParameterSet	parameters

int	prescaleFactor

int	run

std::vector< DTWireId >	theNoisyChannels

edm::ESHandle< DTTtrig >	tTrigMap

int	updates

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

DQM Test Client

A. Gresele - INFN Trento G. Mila - INFN Torino M. Zanetti - CERN PH

threadsafe version (//-) oct/nov 2014 - WATWanAbdullah -ncpp-um-my

Definition at line 55 of file DTNoiseTest.h.

Constructor & Destructor Documentation

DTNoiseTest::DTNoiseTest ( const edm::ParameterSet & ps )

Constructor.

Definition at line 39 of file DTNoiseTest.cc.

References Parameters::parameters, and dtDQMClient_cfg::prescaleFactor.

                                                  {
 
   edm::LogVerbatim ("noise") <<"[DTNoiseTest]: Constructor";  
 
   parameters = ps;
   
   prescaleFactor = parameters.getUntrackedParameter<int>("diagnosticPrescale", 1);
 
   bookingdone = 0;
 
 }

DTNoiseTest::~DTNoiseTest ( )

virtual

Destructor.

Definition at line 53 of file DTNoiseTest.cc.

                          {
 
   edm::LogVerbatim ("noise") <<"DTNoiseTest: analyzed " << updates << " events";
 
 }

Member Function Documentation

void DTNoiseTest::bookHistos	(	DQMStore::IBooker &	,
		const DTChamberId &	ch,
		std::string	folder,
		std::string	histoTag
	)

protected

book the new ME

void DTNoiseTest::bookHistos	(	DQMStore::IBooker &	,
		const DTLayerId &	ch,
		int	nWire,
		std::string	folder,
		std::string	histoTag
	)

protected

void DTNoiseTest::dqmEndJob	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

overrideprotectedvirtual

Implements DQMEDHarvester.

Definition at line 221 of file DTNoiseTest.cc.

                                                                                 {
 
   edm::LogVerbatim ("noise") <<"[DTNoiseTest] endjob called!";
 }

void DTNoiseTest::dqmEndLuminosityBlock	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	context
	)

protectedvirtual

DQM Client Diagnostic.

Reimplemented from DQMEDHarvester.

Definition at line 59 of file DTNoiseTest.cc.

References PDRates::average, begin, bookHistos(), DTTimeUnits::counts, end, spr::find(), edm::EventSetup::get(), DQMStore::IGetter::get(), QReport::getBadChannels(), MonitorElement::getTH2F(), mergeVDriftHistosByStation::histos, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), Mean, cmsHarvester::nevents, Parameters::parameters, dtDQMClient_cfg::prescaleFactor, edm::second(), DQMStore::IBooker::setCurrentFolder(), crabStatusFromReport::statusMap, AlCaHLTBitMon_QueryRunRegistry::string, DTSuperLayerId::superLayer(), and w.

                                                                                                 {
 
   if (!bookingdone) {
     ibooker.setCurrentFolder("DT/Tests/Noise");
   updates = 0;
   // Get the geometry
   context.get<MuonGeometryRecord>().get(muonGeom);
   }
   bookingdone = 1; 
 
   edm::LogVerbatim ("noise") <<"[DTNoiseTest]: End of LS transition, performing the DQM client operation";
 
   // counts number of lumiSegs 
   nLumiSegs = lumiSeg.id().luminosityBlock();
 
   // prescale factor
   if ( nLumiSegs%prescaleFactor != 0 ) return;
 
   edm::LogVerbatim ("noise") <<"[DTNoiseTest]: "<<nLumiSegs<<" updates";
 
   ESHandle<DTStatusFlag> statusMap;
   context.get<DTStatusFlagRcd>().get(statusMap);
   
   context.get<DTTtrigRcd>().get(tTrigMap);
   float tTrig, tTrigRMS, kFactor;
 
   string histoTag;
   // loop over chambers
   vector<const DTChamber*>::const_iterator ch_it = muonGeom->chambers().begin();
   vector<const DTChamber*>::const_iterator ch_end = muonGeom->chambers().end();
 
   for (; ch_it != ch_end; ++ch_it) {
     DTChamberId ch = (*ch_it)->id();
     vector<const DTSuperLayer*>::const_iterator sl_it = (*ch_it)->superLayers().begin(); 
     vector<const DTSuperLayer*>::const_iterator sl_end = (*ch_it)->superLayers().end();
         
     MonitorElement * noiseME = igetter.get(getMEName(ch));
     if (noiseME) {
       TH2F * noiseHisto = noiseME->getTH2F();
 
       // WARNING uncorrect normalization!! TO BE PROVIDED CENTRALLY
       double nevents = (int) noiseHisto->GetEntries();  
     
       double normalization =0;
 
       float average=0;
       float nOfChannels=0;
       float noiseStatistics=0;
       int newNoiseChannels=0;
 
       for(; sl_it != sl_end; ++sl_it) {
     const DTSuperLayerId & slID = (*sl_it)->id();
         
         // ttrig and rms are counts
     tTrigMap->get(slID, tTrig, tTrigRMS, kFactor, DTTimeUnits::counts);
     if (tTrig==0) tTrig=1;
     const double ns_s = 1e9*(32/25);
     normalization = ns_s/float(tTrig*nevents);
         
     noiseHisto->Scale(normalization);
         
     // loop over layers
         
     for (int binY=(slID.superLayer()-1)*4+1 ; binY <= (slID.superLayer()-1)*4+4; binY++) {
           
       int Y = binY - 4*(slID.superLayer()-1);
           
       // the layer
           
       const DTLayerId theLayer(slID,Y);
          
       // loop over channels 
       for (int binX=1; binX <= noiseHisto->GetNbinsX(); binX++) {
         
         if (noiseHisto->GetBinContent(binX,binY) > parameters.getUntrackedParameter<int>("HzThreshold", 300))
           theNoisyChannels.push_back(DTWireId(theLayer, binX));
           
         // get rid of the dead channels
         else {
           average += noiseHisto->GetBinContent(binX,binY); 
           nOfChannels++; 
         }
       }
     }
         
     if (nOfChannels) noiseStatistics = average/nOfChannels;
     histoTag = "NoiseAverage";
 
     if (histos[histoTag].find((*ch_it)->id().rawId()) == histos[histoTag].end()) bookHistos(ibooker,(*ch_it)->id(),string("NoiseAverage"), histoTag );
     histos[histoTag].find((*ch_it)->id().rawId())->second->setBinContent(slID.superLayer(),noiseStatistics); 
 
     for ( vector<DTWireId>::const_iterator nb_it = theNoisyChannels.begin();
           nb_it != theNoisyChannels.end(); ++nb_it) {
           
       bool isNoisy = false;
       bool isFEMasked = false;
       bool isTDCMasked = false;
       bool isTrigMask = false;
       bool isDead = false;
       bool isNohv = false;
       statusMap->cellStatus((*nb_it), isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
                   
       if (!isNoisy) newNoiseChannels++;
     }
     theNoisyChannels.clear();
     histoTag = "NewNoisyChannels";
 
     if (histos[histoTag].find((*ch_it)->id().rawId()) == histos[histoTag].end()) 
                bookHistos(ibooker,(*ch_it)->id(),string("NewNoisyChannels"), histoTag );
     histos[histoTag].find((*ch_it)->id().rawId())->second->setBinContent(slID.superLayer(), newNoiseChannels);   
       }
     }
     //To compute the Noise Mean test
     vector<const DTSuperLayer*>::const_iterator sl2_it = (*ch_it)->superLayers().begin(); 
     vector<const DTSuperLayer*>::const_iterator sl2_end = (*ch_it)->superLayers().end();
     for(; sl2_it != sl2_end; ++sl2_it) {
       vector<const DTLayer*>::const_iterator l_it = (*sl2_it)->layers().begin(); 
       vector<const DTLayer*>::const_iterator l_end = (*sl2_it)->layers().end();
       for(; l_it != l_end; ++l_it) {
     
     DTLayerId lID = (*l_it)->id();
     MonitorElement * noisePerEventME = igetter.get(getMEName(lID));
 
     if (noisePerEventME) {
       TH2F * noiseHistoPerEvent = noisePerEventME->getTH2F();
       int nWires = muonGeom->layer(lID)->specificTopology().channels();
       double MeanNumerator=0, MeanDenominator=0;
       histoTag = "MeanDigiPerEvent";
       for (int w=1; w<=nWires; w++){
         for(int numDigi=1; numDigi<=10; numDigi++){
           MeanNumerator+=(noiseHistoPerEvent->GetBinContent(w,numDigi)*(numDigi-1));
           MeanDenominator+=noiseHistoPerEvent->GetBinContent(w,numDigi);
         }
         double Mean=MeanNumerator/MeanDenominator;
         if (histos[histoTag].find((*l_it)->id().rawId()) == histos[histoTag].end()) bookHistos(ibooker,(*l_it)->id(),nWires, string("MeanDigiPerEvent"), histoTag );
         histos[histoTag].find((*l_it)->id().rawId())->second->setBinContent(w, Mean);   
       } 
     }
       }
     }
   }
   
   // Noise Mean test 
   histoTag = "MeanDigiPerEvent";
   string MeanCriterionName = parameters.getUntrackedParameter<string>("meanTestName","NoiseMeanInRange");
   for(map<uint32_t, MonitorElement*>::const_iterator hMean = histos[histoTag].begin();
       hMean != histos[histoTag].end();
       hMean++) {
     const QReport * theMeanQReport = (*hMean).second->getQReport(MeanCriterionName);
     if(theMeanQReport) {
       vector<dqm::me_util::Channel> badChannels = theMeanQReport->getBadChannels();
       for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); 
        channel != badChannels.end(); channel++) {
     LogVerbatim ("tTrigCalibration")<<"LayerId : "<<(*hMean).first<<" Bad mean channels: "
                                         <<(*channel).getBin()<<"  Contents : "<<(*channel).getContents();
       }
     }
   }
   
 }

string DTNoiseTest::getMEName ( const DTChamberId & ch )

protected

Get the ME name.

Definition at line 227 of file DTNoiseTest.cc.

References Parameters::parameters, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, and DTChamberId::wheel().

                                                     {
   
   stringstream wheel; wheel << ch.wheel();  
   stringstream station; station << ch.station();    
   stringstream sector; sector << ch.sector();   
   
   string folderRoot = parameters.getUntrackedParameter<string>("folderRoot", "Collector/FU0/");
   string folderTag = parameters.getUntrackedParameter<string>("folderTag", "Occupancies");
   string folderName = 
     folderRoot + "DT/DTDigiTask/Wheel" +  wheel.str() +
     "/Station" + station.str() +
     "/Sector" + sector.str() + "/" + folderTag + "/";
 
   string histoTag = parameters.getUntrackedParameter<string>("histoTag", "OccupancyNoise_perCh");
   string histoname = folderName + histoTag  
     + "_W" + wheel.str() 
     + "_St" + station.str() 
     + "_Sec" + sector.str(); 
     
     
   return histoname;
   
 }

string DTNoiseTest::getMEName ( const DTLayerId & ly )

protected

Definition at line 251 of file DTNoiseTest.cc.

References DTLayerId::layer(), Parameters::parameters, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superlayer(), and DTChamberId::wheel().

                                                   {
   
   stringstream wheel; wheel << ly.wheel();  
   stringstream station; station << ly.station();    
   stringstream sector; sector << ly.sector();
   stringstream superLayer; superLayer << ly.superlayer();
   stringstream layer; layer << ly.layer();
   
   string folderRoot = parameters.getUntrackedParameter<string>("folderRoot", "Collector/FU0/");
   string folderTag = parameters.getUntrackedParameter<string>("folderTagForDigiPerEventTest", "DigiPerEvent");
   string folderName = 
     folderRoot + "DT/DTDigiForNoiseTask/Wheel" +  wheel.str() +
     "/Station" + station.str() +
     "/Sector" + sector.str() + "/" + folderTag + "/";
   
   string histoTag = parameters.getUntrackedParameter<string>("histoTagForDigiPerEventTest", "DigiPerEvent");
   string histoname = folderName + histoTag  
     + "_W" + wheel.str() 
     + "_St" + station.str() 
     + "_Sec" + sector.str()
     + "_SL" + superLayer.str()
     + "_L" + layer.str();
     
     
   return histoname;
 
 }

Member Data Documentation

bool DTNoiseTest::bookingdone

private

Definition at line 90 of file DTNoiseTest.h.

bool DTNoiseTest::debug

private

Definition at line 84 of file DTNoiseTest.h.

Referenced by rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), and pkg.AbstractPkg::write().

std::map<std::string, std::map<uint32_t, MonitorElement*> > DTNoiseTest::histos

private

Definition at line 100 of file DTNoiseTest.h.

Referenced by cuy.ValElement::__init__(), cuy.additionElement::__init__(), cuy.superimposeElement::__init__(), and cuy.graphElement::__init__().

edm::ESHandle<DTGeometry> DTNoiseTest::muonGeom

private

Definition at line 93 of file DTNoiseTest.h.

unsigned int DTNoiseTest::nLumiSegs

private

Definition at line 86 of file DTNoiseTest.h.

edm::ParameterSet DTNoiseTest::parameters

private

Definition at line 92 of file DTNoiseTest.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

int DTNoiseTest::prescaleFactor

private

Definition at line 87 of file DTNoiseTest.h.

int DTNoiseTest::run

private

Definition at line 88 of file DTNoiseTest.h.

Referenced by generateEDF.LumiInfo::__str__(), Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), generateEDF.LumiInfo::fixXingInfo(), and dqm_interfaces.DirWalkerFile::walk().

std::vector<DTWireId> DTNoiseTest::theNoisyChannels

private

Definition at line 97 of file DTNoiseTest.h.

edm::ESHandle<DTTtrig> DTNoiseTest::tTrigMap

private

Definition at line 94 of file DTNoiseTest.h.

int DTNoiseTest::updates

private

Definition at line 85 of file DTNoiseTest.h.

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation