#include <DTOccupancyTest.h>

Inheritance diagram for DTOccupancyTest:

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

virtual	~DTOccupancyTest ()
	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

	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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &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	beginRun (edm::Run const &run, edm::EventSetup const &context)
	BeginRun. More...

void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &)
	Endjob. More...

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

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 Member Functions
void	bookHistos (DQMStore::IBooker &, const int wheelId, std::string folder, std::string histoTag)
	book the summary histograms More...

int	getIntegral (TH2F *histo, int, int, int, int, bool)

std::string	getMEName (std::string histoTag, const DTChamberId &chId)
	Get the ME name. More...

int	runOccupancyTest (TH2F *histo, const DTChamberId &chId, float &chamberPercentage)

std::string	topFolder () const

Private Attributes
bool	bookingdone

MonitorElement *	glbSummaryHisto

int	lsCounter

std::set< DTLayerId >	monitoredLayers

edm::ESHandle< DTGeometry >	muonGeom

std::string	nameMonitoredHisto

int	nChannelDead

int	nChannelTotal

int	nevents

int	nMinEvts

int	nMinEvtsPC

TNtuple *	ntuple

int	nZeroEvtsPC

TFile *	rootFile

bool	runOnAllHitsOccupancies

bool	runOnInTimeOccupancies

bool	runOnNoiseOccupancies

MonitorElement *	summaryHisto

bool	tpMode

std::map< int, MonitorElement * >	wheelHistos

bool	writeRootFile

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)

Detailed Description

DQM Test Client

Author
G. Cerminara - University and INFN Torino

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

Definition at line 39 of file DTOccupancyTest.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 31 of file DTOccupancyTest.cc.

References edm::ParameterSet::getUntrackedParameter(), and cmsHarvester::nevents.

                                                          {
   LogVerbatim ("DTDQM|DTMonitorClient|DTOccupancyTest") << "[DTOccupancyTest]: Constructor";
 
   // Get the DQM service
 
   lsCounter = 0;
 
   writeRootFile  = ps.getUntrackedParameter<bool>("writeRootFile", false);
   if(writeRootFile) {
     rootFile = new TFile("DTOccupancyTest.root","RECREATE");
     ntuple = new TNtuple("OccupancyNtuple", "OccupancyNtuple", "ls:wh:st:se:lay1MeanCell:lay1RMS:lay2MeanCell:lay2RMS:lay3MeanCell:lay3RMS:lay4MeanCell:lay4RMS:lay5MeanCell:lay5RMS:lay6MeanCell:lay6RMS:lay7MeanCell:lay7RMS:lay8MeanCell:lay8RMS:lay9MeanCell:lay9RMS:lay10MeanCell:lay10RMS:lay11MeanCell:lay11RMS:lay12MeanCell:lay12RMS");
   }
   
   // switch on the mode for running on test pulses (different top folder)
   tpMode = ps.getUntrackedParameter<bool>("testPulseMode", false);
   
   runOnAllHitsOccupancies =  ps.getUntrackedParameter<bool>("runOnAllHitsOccupancies", true);
   runOnNoiseOccupancies =  ps.getUntrackedParameter<bool>("runOnNoiseOccupancies", false);
   runOnInTimeOccupancies = ps.getUntrackedParameter<bool>("runOnInTimeOccupancies", false);
   nMinEvts  = ps.getUntrackedParameter<int>("nEventsCert", 5000);
   nMinEvtsPC  = ps.getUntrackedParameter<int>("nEventsMinPC", 2200);
   nZeroEvtsPC  = ps.getUntrackedParameter<int>("nEventsZeroPC", 30);
 
   bookingdone = 0;
 
   // Event counter
   nevents = 0;
 
 }

DTOccupancyTest::~DTOccupancyTest ( )

virtual

Destructor.

Definition at line 61 of file DTOccupancyTest.cc.

                                  {
   LogVerbatim ("DTDQM|DTMonitorClient|DTOccupancyTest") << " destructor called" << endl;
 }

Member Function Documentation

void DTOccupancyTest::beginRun	(	edm::Run const &	run,
		edm::EventSetup const &	context
	)

protectedvirtual

BeginRun.

Reimplemented from DQMEDHarvester.

Definition at line 66 of file DTOccupancyTest.cc.

References edm::EventSetup::get().

                                                                           {
 
   LogVerbatim ("DTDQM|DTMonitorClient|DTOccupancyTest") << "[DTOccupancyTest]: BeginRun";
 
   // Get the geometry
   context.get<MuonGeometryRecord>().get(muonGeom);
 
 }

void DTOccupancyTest::bookHistos	(	DQMStore::IBooker &	,
		const int	wheelId,
		std::string	folder,
		std::string	histoTag
	)

private

book the summary histograms

Definition at line 213 of file DTOccupancyTest.cc.

References DQMStore::IBooker::book2D(), MonitorElement::setBinLabel(), and DQMStore::IBooker::setCurrentFolder().

                                                                                       {
   // Set the current folder
   stringstream wheel; wheel << wheelId; 
 
   ibooker.setCurrentFolder(topFolder());
 
   // build the histo name
   string histoName = histoTag + "_W" + wheel.str(); 
   
   
   LogVerbatim ("DTDQM|DTMonitorClient|DTOccupancyTest") <<"[DTOccupancyTest]: booking wheel histo:"
                             << histoName 
                             << " (tag "
                             << histoTag
                             << ") in: "
                             << topFolder() + "Wheel"+ wheel.str() + "/" + folder << endl;
   
   string histoTitle = "Occupancy summary WHEEL: "+wheel.str();
   if(tpMode) {
     histoTitle = "TP Occupancy summary WHEEL: "+wheel.str();
   }
 
   wheelHistos[wheelId] = ibooker.book2D(histoName,histoTitle,12,1,13,4,1,5);
   wheelHistos[wheelId]->setBinLabel(1,"MB1",2);
   wheelHistos[wheelId]->setBinLabel(2,"MB2",2);
   wheelHistos[wheelId]->setBinLabel(3,"MB3",2);
   wheelHistos[wheelId]->setBinLabel(4,"MB4",2);
   wheelHistos[wheelId]->setAxisTitle("sector",1);
 }

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

protectedvirtual

Endjob.

Implements DQMEDHarvester.

Definition at line 199 of file DTOccupancyTest.cc.

                                                                                     {
 
   LogVerbatim ("DTDQM|DTMonitorClient|DTOccupancyTest") << "[DTOccupancyTest] endjob called!";
   if(writeRootFile) {
     rootFile->cd();
     ntuple->Write();
     rootFile->Close();
   }
 }

void DTOccupancyTest::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 75 of file DTOccupancyTest.cc.

References DQMStore::IBooker::book2D(), bookHistos(), chambers, DQMStore::IGetter::get(), MonitorElement::getFloatValue(), MonitorElement::getTH2F(), timingPdfMaker::histo, mps_fire::result, DTChamberId::sector(), MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), indexGen::title, and DTChamberId::wheel().

                                                                                                    {
   if (!bookingdone) {
 
   // Book the summary histos
   //   - one summary per wheel
   for(int wh = -2; wh <= 2; ++wh) { // loop over wheels
     bookHistos(ibooker,wh, string("Occupancies"), "OccupancySummary");
     }
 
   ibooker.setCurrentFolder(topFolder());
   string title = "Occupancy Summary";
   if(tpMode) {
     title = "Test Pulse Occupancy Summary";
     }
   //   - global summary with alarms
   summaryHisto = ibooker.book2D("OccupancySummary",title.c_str(),12,1,13,5,-2,3);
   summaryHisto->setAxisTitle("sector",1);
   summaryHisto->setAxisTitle("wheel",2);
   
   //   - global summary with percentages
   glbSummaryHisto = ibooker.book2D("OccupancyGlbSummary",title.c_str(),12,1,13,5,-2,3);
   glbSummaryHisto->setAxisTitle("sector",1);
   glbSummaryHisto->setAxisTitle("wheel",2);
 
 
   // assign the name of the input histogram
   if(runOnAllHitsOccupancies) {
     nameMonitoredHisto = "OccupancyAllHits_perCh";
     } else if(runOnNoiseOccupancies) {
     nameMonitoredHisto = "OccupancyNoise_perCh";
     } else if(runOnInTimeOccupancies) {
     nameMonitoredHisto = "OccupancyInTimeHits_perCh";
     } else { // default is AllHits histo
     nameMonitoredHisto = "OccupancyAllHits_perCh";
     }
 
   }
   bookingdone = 1; 
 
 
   LogVerbatim ("DTDQM|DTMonitorClient|DTOccupancyTest")
     <<"[DTOccupancyTest]: End of LS transition, performing the DQM client operation";
   lsCounter++;
 
 
 
   // Reset the global summary
   summaryHisto->Reset();
   glbSummaryHisto->Reset();
 
   nChannelTotal = 0;
   nChannelDead = 0;
 
   // Get all the DT chambers
   vector<const DTChamber*> chambers = muonGeom->chambers();
 
   for(vector<const DTChamber*>::const_iterator chamber = chambers.begin();
       chamber != chambers.end(); ++chamber) {  // Loop over all chambers
     DTChamberId chId = (*chamber)->id();
 
     MonitorElement * chamberOccupancyHisto = igetter.get(getMEName(nameMonitoredHisto, chId));  
 
     // Run the tests on the plot for the various granularities
     if(chamberOccupancyHisto != 0) {
       // Get the 2D histo
       TH2F* histo = chamberOccupancyHisto->getTH2F();
       float chamberPercentage = 1.;
       int result = runOccupancyTest(histo, chId, chamberPercentage);
       int sector = chId.sector();
 
       if(sector == 13) {
     sector = 4;
     float resultSect4 = wheelHistos[chId.wheel()]->getBinContent(sector, chId.station());
     if(resultSect4 > result) {
       result = (int)resultSect4;
     }
       } else if(sector == 14) {
     sector = 10;
     float resultSect10 = wheelHistos[chId.wheel()]->getBinContent(sector, chId.station());
     if(resultSect10 > result) {
       result = (int)resultSect10;
     }
       }
       
       // the 2 MB4 of Sect 4 and 10 count as half a chamber
       if((sector == 4 || sector == 10) && chId.station() == 4) 
     chamberPercentage = chamberPercentage/2.;
 
       wheelHistos[chId.wheel()]->setBinContent(sector, chId.station(),result);
       if(result > summaryHisto->getBinContent(sector, chId.wheel()+3)) {
     summaryHisto->setBinContent(sector, chId.wheel()+3, result);
       }
       glbSummaryHisto->Fill(sector, chId.wheel(), chamberPercentage*1./4.);
     } else {
       LogVerbatim ("DTDQM|DTMonitorClient|DTOccupancyTest") << "[DTOccupancyTest] ME: "
                       << getMEName(nameMonitoredHisto, chId) << " not found!" << endl;
     }
 
   }
 
   string nEvtsName = "DT/EventInfo/Counters/nProcessedEventsDigi";
 
   MonitorElement * meProcEvts = igetter.get(nEvtsName);
 
   if (meProcEvts) {
     int nProcEvts = meProcEvts->getFloatValue();
     glbSummaryHisto->setEntries(nProcEvts < nMinEvts ? 10. : nProcEvts);
     summaryHisto->setEntries(nProcEvts < nMinEvts ? 10. : nProcEvts);
   } else {
     glbSummaryHisto->setEntries(nMinEvts +1);
     summaryHisto->setEntries(nMinEvts + 1);
     LogVerbatim ("DTDQM|DTMonitorClient|DTOccupancyTest") << "[DTOccupancyTest] ME: "
                <<  nEvtsName << " not found!" << endl;
   }
 
   // Fill the global summary
   // Check for entire sectors off and report them on the global summary
   //FIXME: TODO
 
   if(writeRootFile) ntuple->AutoSave("SaveSelf");
 
 }

int DTOccupancyTest::getIntegral	(	TH2F *	histo,
		int	firstBinX,
		int	lastBinX,
		int	firstBinY,
		int	lastBinY,
		bool	doall
	)

private

Definition at line 271 of file DTOccupancyTest.cc.

References i, and j.

                                                                                                                   {
 
   int sum = 0;
     for (Int_t i = firstBinX; i < lastBinX+1; i++) {
       for (Int_t j = firstBinY; j < lastBinY+1; j++) {
         
         if (histo->GetBinContent(i,j) >0){
           if (!doall) return 1;
       sum += histo->GetBinContent(i,j);
     }
       }
     }
 
     return sum;
 }

string DTOccupancyTest::getMEName	(	std::string	histoTag,
		const DTChamberId &	chId
	)

private

Get the ME name.

Definition at line 246 of file DTOccupancyTest.cc.

References hlt_dqm_clientPB-live_cfg::folder, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, and DTChamberId::wheel().

                                                                           {
 
   stringstream wheel; wheel << chId.wheel();
   stringstream station; station << chId.station();
   stringstream sector; sector << chId.sector();
 
 
   string folderRoot = topFolder() + "Wheel" + wheel.str() +
     "/Sector" + sector.str() +
     "/Station" + station.str() + "/";
 
   string folder = "Occupancies/";
   
   // build the histo name
   string histoName = histoTag 
     + "_W" + wheel.str() 
     + "_St" + station.str() 
     + "_Sec" + sector.str();
 
   string histoname = folderRoot + histoName;
 
   return histoname;
 }

int DTOccupancyTest::runOccupancyTest	(	TH2F *	histo,
		const DTChamberId &	chId,
		float &	chamberPercentage
	)

private

Definition at line 297 of file DTOccupancyTest.cc.

References LogTrace, DTChamberId::sector(), and DTChamberId::station().

                                                 {
   int nBinsX = histo->GetNbinsX();
 
   LogTrace("DTDQM|DTMonitorClient|DTOccupancyTest") << "--- Occupancy test for chamber: " << chId << endl;
 
   int compDeadCell = 0;
   int totCell = 0;
   int totOccup = 0;
 
   for(int slay = 1; slay <= 3; ++slay) { // loop over SLs
     int binYlow = ((slay-1)*4)+1;
 
     if(chId.station() == 4 && slay == 2) continue;
     for(int lay = 1; lay <= 4; ++lay) { // loop over layers
       DTLayerId layID(chId,slay,lay);
       int firstWire = muonGeom->layer(layID)->specificTopology().firstChannel();
       int nWires = muonGeom->layer(layID)->specificTopology().channels();
       int binY = binYlow+(lay-1);      
       int totalDeadCells = 0;
       int nDeadCellsInARow = 1;
       int nDeadCellsInARowMax = 0;
       int nCellsZeroCount = 0;
       bool previousIsDead = false;
 
       int interDeadCells = 0;
 
       totCell += nWires;
 
       for(int cell = firstWire; cell != (nWires+firstWire); ++cell) { // loop over cells
     double cellOccup = histo->GetBinContent(cell,binY);
         LogTrace("DTDQM|DTMonitorClient|DTOccupancyTest") << "       cell occup: " << cellOccup;
     totOccup += cellOccup;
 
     if(cellOccup == 0) {
           nCellsZeroCount++;
           totalDeadCells++;
           if(previousIsDead) {
         nDeadCellsInARow++;
       }
       else {
         //      if(interDeadCells > 3) {
         if(nDeadCellsInARow > nDeadCellsInARowMax) nDeadCellsInARowMax = nDeadCellsInARow;
         nDeadCellsInARow = 1; 
         //      }
       }
           previousIsDead = true;
           interDeadCells = 0;
           LogTrace("DTDQM|DTMonitorClient|DTOccupancyTest") << "       below reference" << endl;
         } else {
           previousIsDead = false;
           interDeadCells++;
     }
     //   // 3 cells not dead between a group of dead cells don't break the count
     if(nDeadCellsInARow > nDeadCellsInARowMax) nDeadCellsInARowMax = nDeadCellsInARow;
       }
       compDeadCell += totalDeadCells;
       if(nDeadCellsInARowMax >= 7.) {
         histo->SetBinContent(nBinsX+1,binY,-1.);
       }
     }
   }
 
 
   nChannelTotal += totCell;
   nChannelDead += compDeadCell;
   chamberPercentage = 1.-(float(compDeadCell)/totCell);
 
   int min_occup = nZeroEvtsPC*20;
   if (chId.station() == 3) min_occup = nZeroEvtsPC*3;
   if (chId.station() == 2) min_occup = nZeroEvtsPC*8;
   if ((chId.station() == 4) && (chId.sector() == 9)) min_occup = nZeroEvtsPC*3;
   if ((chId.station() == 4) && (chId.sector() == 10)) min_occup = nZeroEvtsPC*3;
   if ((chId.station() == 4) && (chId.sector() == 11)) min_occup = nZeroEvtsPC*3;
   if ((chId.station() == 4) && (chId.sector() == 14)) min_occup = nZeroEvtsPC*3;
 
 
   if (totOccup < min_occup) return 4;
   if (totOccup < nMinEvtsPC) chamberPercentage = 1.;
 
   if (chamberPercentage < 0.2) return 4;
   if (chamberPercentage < 0.5) return 3;
   if (chamberPercentage < 0.75) return 2;
   if (chamberPercentage < 0.9) return 1;
 
   return 0;
 }

string DTOccupancyTest::topFolder ( ) const

private

Definition at line 386 of file DTOccupancyTest.cc.

References AlCaHLTBitMon_QueryRunRegistry::string.

                                         {
   if(tpMode) return string("DT/10-TestPulses/");
   return string("DT/01-Digi/");
 }

Member Data Documentation

bool DTOccupancyTest::bookingdone

private

Definition at line 102 of file DTOccupancyTest.h.

MonitorElement* DTOccupancyTest::glbSummaryHisto

private

Definition at line 83 of file DTOccupancyTest.h.

int DTOccupancyTest::lsCounter

private

Definition at line 87 of file DTOccupancyTest.h.

std::set<DTLayerId> DTOccupancyTest::monitoredLayers

private

Definition at line 85 of file DTOccupancyTest.h.

edm::ESHandle<DTGeometry> DTOccupancyTest::muonGeom

private

Definition at line 78 of file DTOccupancyTest.h.

std::string DTOccupancyTest::nameMonitoredHisto

private

Definition at line 100 of file DTOccupancyTest.h.

int DTOccupancyTest::nChannelDead

private

Definition at line 105 of file DTOccupancyTest.h.

int DTOccupancyTest::nChannelTotal

private

Definition at line 104 of file DTOccupancyTest.h.

int DTOccupancyTest::nevents

private

Definition at line 76 of file DTOccupancyTest.h.

int DTOccupancyTest::nMinEvts

private

Definition at line 88 of file DTOccupancyTest.h.

int DTOccupancyTest::nMinEvtsPC

private

Definition at line 89 of file DTOccupancyTest.h.

TNtuple* DTOccupancyTest::ntuple

private

Definition at line 94 of file DTOccupancyTest.h.

int DTOccupancyTest::nZeroEvtsPC

private

Definition at line 90 of file DTOccupancyTest.h.

TFile* DTOccupancyTest::rootFile

private

Definition at line 93 of file DTOccupancyTest.h.

bool DTOccupancyTest::runOnAllHitsOccupancies

private

Definition at line 97 of file DTOccupancyTest.h.

bool DTOccupancyTest::runOnInTimeOccupancies

private

Definition at line 99 of file DTOccupancyTest.h.

bool DTOccupancyTest::runOnNoiseOccupancies

private

Definition at line 98 of file DTOccupancyTest.h.

MonitorElement* DTOccupancyTest::summaryHisto

private

Definition at line 82 of file DTOccupancyTest.h.

bool DTOccupancyTest::tpMode

private

Definition at line 95 of file DTOccupancyTest.h.

std::map< int, MonitorElement* > DTOccupancyTest::wheelHistos

private

Definition at line 81 of file DTOccupancyTest.h.

bool DTOccupancyTest::writeRootFile

private

Definition at line 92 of file DTOccupancyTest.h.

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation