DTOccupancyTest Class Reference

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