#include <DTnoiseDBValidation.h>

Inheritance diagram for DTnoiseDBValidation:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup)
void	beginRun (const edm::Run &run, const edm::EventSetup &setup)
	Operations.
	DTnoiseDBValidation (const edm::ParameterSet &pset)
	Constructor.
void	endJob ()
void	endRun (edm::Run const &, edm::EventSetup const &)
virtual	~DTnoiseDBValidation ()
	Destructor.
Private Member Functions
void	bookHisto (const DTChamberId &)
Private Attributes
DQMStore *	dbe_
MonitorElement *	diffHisto_
std::string	diffTestName_
edm::ESHandle< DTGeometry >	dtGeom_
std::string	labelDB_
std::string	labelDBRef_
MonitorElement *	layerHisto_
std::string	layerTestName_
std::map< DTChamberId, MonitorElement * >	noiseHistoMap_
const DTStatusFlag *	noiseMap_
const DTStatusFlag *	noiseRefMap_
int	noisyCellsRef_
int	noisyCellsValid_
std::string	outputFileName_
bool	outputMEsInRootFile_
MonitorElement *	sectorHisto_
std::string	sectorTestName_
MonitorElement *	stationHisto_
std::string	stationTestName_
MonitorElement *	wheelHisto_
std::string	wheelTestName_

Detailed Description

Plot the noise from the DB comparaison

Date:: 2011/08/02 14:23:50

Revision:: 1.4

Author:: G. Mila - INFN Torino

Definition at line 31 of file DTnoiseDBValidation.h.

Constructor & Destructor Documentation

DTnoiseDBValidation::DTnoiseDBValidation ( const edm::ParameterSet & pset )

Constructor.

Definition at line 42 of file DTnoiseDBValidation.cc.

References dbe_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), cppFunctionSkipper::operator, and DQMStore::setCurrentFolder().

                                                                 {

  LogVerbatim("NoiseDBValidation") << "[DTnoiseDBValidation] Constructor called!";

  // Get the DQM needed services
  dbe_ = edm::Service<DQMStore>().operator->();
  dbe_->setCurrentFolder("DT/DtCalib/NoiseDBValidation");

  // Get dataBase label
  labelDBRef_ = pset.getParameter<string>("labelDBRef");
  labelDB_ = pset.getParameter<string>("labelDB");

  diffTestName_ = "noiseDifferenceInRange";
  if( pset.exists("diffTestName") ) diffTestName_ = pset.getParameter<string>("diffTestName");
 
  wheelTestName_ = "noiseWheelOccInRange";
  if( pset.exists("wheelTestName") ) wheelTestName_ = pset.getParameter<string>("wheelTestName");

  stationTestName_ = "noiseStationOccInRange";
  if( pset.exists("stationTestName") ) stationTestName_ = pset.getParameter<string>("stationTestName");

  sectorTestName_ = "noiseSectorOccInRange";
  if( pset.exists("sectorTestName") ) sectorTestName_ = pset.getParameter<string>("sectorTestName");

  layerTestName_ = "noiseLayerOccInRange";
  if( pset.exists("layerTestName") ) layerTestName_ = pset.getParameter<string>("layerTestName");
  
  outputMEsInRootFile_ = false;
  if( pset.exists("OutputFileName") ){
     outputMEsInRootFile_ = true;
     outputFileName_ = pset.getParameter<std::string>("OutputFileName");
  }
}

DTnoiseDBValidation::~DTnoiseDBValidation ( ) [virtual]

Destructor.

Definition at line 77 of file DTnoiseDBValidation.cc.

{}

Member Function Documentation

void DTnoiseDBValidation::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)		`[inline, virtual]`

Implements edm::EDAnalyzer.

Definition at line 43 of file DTnoiseDBValidation.h.

{}

void DTnoiseDBValidation::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)		`[virtual]`

Operations.

Reimplemented from edm::EDAnalyzer.

Definition at line 79 of file DTnoiseDBValidation.cc.

References abs, plotTTrigFromHistos::binNumber(), DQMStore::book1D(), DTSuperLayerId::chamberId(), dbe_, edm::EventSetup::get(), DTTopology::lastChannel(), DTLayerId::layer(), DTWireId::layerId(), pileupReCalc_HLTpaths::scale, MonitorElement::setBinLabel(), relativeConstraints::station, DTSuperLayerId::superLayer(), DTLayerId::superlayerId(), and DTWireId::wire().

                                                                             {
  ESHandle<DTStatusFlag> noiseRef;
  setup.get<DTStatusFlagRcd>().get(labelDBRef_, noiseRef);
  noiseRefMap_ = &*noiseRef;
 
  ESHandle<DTStatusFlag> noiseValid;
  setup.get<DTStatusFlagRcd>().get(labelDB_, noiseValid);
  noiseMap_ = &*noiseValid;

  // Get the geometry
  setup.get<MuonGeometryRecord>().get(dtGeom_);

  LogVerbatim("NoiseDBValidation")<<"[DTnoiseDBValidation] Parameters initialization";
 
  noisyCellsRef_ = 0;
  noisyCellsValid_ = 0;

  // Histo booking
  diffHisto_ = dbe_->book1D("noisyCellDiff", "percentual (wrt the previous db) total number of noisy cells",1, 0.5, 1.5);
  diffHisto_->setBinLabel(1,"Diff");
  wheelHisto_ = dbe_->book1D("wheelOccupancy", "percentual noisy cells occupancy per wheel",5, -2.5, 2.5);
  wheelHisto_->setBinLabel(1,"Wh-2");
  wheelHisto_->setBinLabel(2,"Wh-1");
  wheelHisto_->setBinLabel(3,"Wh0");
  wheelHisto_->setBinLabel(4,"Wh1");
  wheelHisto_->setBinLabel(5,"Wh2");
  stationHisto_ = dbe_->book1D("stationOccupancy", "percentual noisy cells occupancy per station",4, 0.5, 4.5);
  stationHisto_->setBinLabel(1,"St1");
  stationHisto_->setBinLabel(2,"St2");
  stationHisto_->setBinLabel(3,"St3");
  stationHisto_->setBinLabel(4,"St4");
  sectorHisto_ = dbe_->book1D("sectorOccupancy", "percentual noisy cells occupancy per sector",12, 0.5, 12.5);
  sectorHisto_->setBinLabel(1,"Sect1");
  sectorHisto_->setBinLabel(2,"Sect2");
  sectorHisto_->setBinLabel(3,"Sect3");
  sectorHisto_->setBinLabel(4,"Sect4");
  sectorHisto_->setBinLabel(5,"Sect5");
  sectorHisto_->setBinLabel(6,"Sect6");
  sectorHisto_->setBinLabel(7,"Sect7");
  sectorHisto_->setBinLabel(8,"Sect8");
  sectorHisto_->setBinLabel(9,"Sect9");
  sectorHisto_->setBinLabel(10,"Sect10");
  sectorHisto_->setBinLabel(11,"Sect11");
  sectorHisto_->setBinLabel(12,"Sect12");
  layerHisto_ = dbe_->book1D("layerOccupancy", "percentual noisy cells occupancy per layer",3, 0.5, 3.5);
  layerHisto_->setBinLabel(1,"First 10 bins");
  layerHisto_->setBinLabel(2,"Middle bins");
  layerHisto_->setBinLabel(3,"Last 10 bins");

  // map initialization
  map<int, int> whMap;
  whMap.clear();
  map<int, int> stMap;
  stMap.clear();
  map<int, int> sectMap;
  sectMap.clear();
  map<int, int> layerMap;
  layerMap.clear();

  // Loop over reference DB entries
  for(DTStatusFlag::const_iterator noise = noiseRefMap_->begin();
      noise != noiseRefMap_->end(); noise++) {
    DTWireId wireId((*noise).first.wheelId,
                    (*noise).first.stationId,
                    (*noise).first.sectorId,
                    (*noise).first.slId,
                    (*noise).first.layerId,
                    (*noise).first.cellId);
    LogVerbatim("NoiseDBValidation") << "Ref. noisy wire: " << wireId;
    ++noisyCellsRef_;
  }

  // Loop over validation DB entries
  for(DTStatusFlag::const_iterator noise = noiseMap_->begin();
      noise != noiseMap_->end(); noise++) {
    DTWireId wireId((*noise).first.wheelId,
                    (*noise).first.stationId,
                    (*noise).first.sectorId,
                    (*noise).first.slId,
                    (*noise).first.layerId,
                    (*noise).first.cellId);
    LogVerbatim("NoiseDBValidation") << "Valid. noisy wire: " <<  wireId;
    ++noisyCellsValid_;

    whMap[(*noise).first.wheelId]++;
    stMap[(*noise).first.stationId]++;
    sectMap[(*noise).first.sectorId]++;
    
    const DTTopology& dtTopo = dtGeom_->layer(wireId.layerId())->specificTopology();
    const int lastWire = dtTopo.lastChannel();
    if((*noise).first.cellId<=10)
      layerMap[1]++;
    if((*noise).first.cellId>10 && (*noise).first.cellId<(lastWire-10))
      layerMap[2]++;
    if((*noise).first.cellId>=(lastWire-10))
      layerMap[3]++;

    const DTChamberId chId = wireId.layerId().superlayerId().chamberId();
    if( noiseHistoMap_.find(chId) == noiseHistoMap_.end() ) bookHisto(chId);
    int binNumber = 4*(wireId.superLayer() - 1) + wireId.layer(); 
    noiseHistoMap_[chId]->Fill(wireId.wire(),binNumber);
  }

  //histo filling
  double scale = 1/double(noisyCellsRef_);
  diffHisto_->Fill(1,abs(noisyCellsRef_ - noisyCellsValid_)*scale);

  scale = 1/double(noisyCellsValid_);
  for(map<int, int >::const_iterator wheel = whMap.begin();
        wheel != whMap.end();
        wheel++) {
    wheelHisto_->Fill((*wheel).first, ((*wheel).second)*scale);
  }
  for(map<int, int >::const_iterator station = stMap.begin();
        station != stMap.end();
        station++) {
    stationHisto_->Fill((*station).first, ((*station).second)*scale);
  }
  for(map<int, int >::const_iterator sector = sectMap.begin();
        sector != sectMap.end();
        sector++) {
    sectorHisto_->Fill((*sector).first, ((*sector).second)*scale);
  }
  for(map<int, int >::const_iterator layer = layerMap.begin();
        layer != layerMap.end();
        layer++) {
    layerHisto_->Fill((*layer).first, ((*layer).second)*scale);
  }

}

void DTnoiseDBValidation::bookHisto ( const DTChamberId & chId ) [private]

Definition at line 273 of file DTnoiseDBValidation.cc.

References plotTTrigFromHistos::binNumber(), DQMStore::book2D(), dbe_, label, DTChamberId::sector(), DTChamberId::station(), DTChamber::superLayers(), and DTChamberId::wheel().

                                                           {
  stringstream histoName;
  histoName << "NoiseOccupancy" 
            << "_W" << chId.wheel() 
            <<"_St" << chId.station() 
            << "_Sec" << chId.sector();

  if( noiseHistoMap_.find(chId) == noiseHistoMap_.end() ){ // Redundant check
     // Get the chamber from the geometry
     int nWiresMax = 0;
     const DTChamber* dtchamber = dtGeom_->chamber(chId);
     const vector<const DTSuperLayer*> superlayers = dtchamber->superLayers();

     // Loop over layers and find the max # of wires
     for(vector<const DTSuperLayer*>::const_iterator sl = superlayers.begin();
                                                     sl != superlayers.end(); ++sl) { // loop over SLs
        vector<const DTLayer*> layers = (*sl)->layers();
        for(vector<const DTLayer*>::const_iterator lay = layers.begin();
                                                   lay != layers.end(); ++lay) { // loop over layers
           int nWires = (*lay)->specificTopology().channels();
           if(nWires > nWiresMax) nWiresMax = nWires;
        }
     }

     noiseHistoMap_[chId] = dbe_->book2D(histoName.str(),"Noise occupancy",nWiresMax,1,(nWiresMax + 1),12,1,13);
     for(int i_sl = 1; i_sl <= 3; ++i_sl) {
        for(int i_lay = 1; i_lay <= 4; ++i_lay) {
           int binNumber = 4*(i_sl - 1) + i_lay;
           stringstream label;
           label << "SL" << i_sl << ": L" << i_lay;
           noiseHistoMap_[chId]->setBinLabel(binNumber,label.str(),2);
        }
     }
  }
  
}

void DTnoiseDBValidation::endJob ( void ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 268 of file DTnoiseDBValidation.cc.

References dbe_, and DQMStore::save().

                                 {
  // Write the histos in a ROOT file
  if(outputMEsInRootFile_) dbe_->save(outputFileName_);
}

void DTnoiseDBValidation::endRun	(	edm::Run const &	run,
		edm::EventSetup const &	setup
	)		`[virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 210 of file DTnoiseDBValidation.cc.

References QReport::getBadChannels().

                                                                              {

  // test on difference histo
  //string testCriterionName;
  //testCriterionName = parameters.getUntrackedParameter<string>("diffTestName","noiseDifferenceInRange");
  const QReport * theDiffQReport = diffHisto_->getQReport(diffTestName_);
  if(theDiffQReport) {
      vector<dqm::me_util::Channel> badChannels = theDiffQReport->getBadChannels();
      for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); 
           channel != badChannels.end(); channel++) {
        LogWarning("NoiseDBValidation") << " Bad partial difference of noisy channels! Contents : " << (*channel).getContents();
      }
  }
  //testCriterionName = parameters.getUntrackedParameter<string>("wheelTestName","noiseWheelOccInRange");
  const QReport * theDiffQReport2 = wheelHisto_->getQReport(wheelTestName_);
  if(theDiffQReport2) {
      vector<dqm::me_util::Channel> badChannels = theDiffQReport2->getBadChannels();
      for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); 
           channel != badChannels.end(); channel++) {
        int wheel = (*channel).getBin()-3;
        LogWarning("NoiseDBValidation") << " Bad percentual occupancy for wheel : " << wheel << "  Contents : " << (*channel).getContents();
      }
  }
  //testCriterionName = parameters.getUntrackedParameter<string>("stationTestName","noiseStationOccInRange");
  const QReport * theDiffQReport3 = stationHisto_->getQReport(stationTestName_);
  if(theDiffQReport3) {
      vector<dqm::me_util::Channel> badChannels = theDiffQReport3->getBadChannels();
      for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); 
           channel != badChannels.end(); channel++) {
        LogWarning("NoiseDBValidation") << " Bad percentual occupancy for station : " << (*channel).getBin() << "  Contents : " << (*channel).getContents();
      }
  }
  //testCriterionName = parameters.getUntrackedParameter<string>("sectorTestName","noiseSectorOccInRange");
  const QReport * theDiffQReport4 = sectorHisto_->getQReport(sectorTestName_);
  if(theDiffQReport4) {
      vector<dqm::me_util::Channel> badChannels = theDiffQReport4->getBadChannels();
      for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); 
           channel != badChannels.end(); channel++) {
        LogWarning("NoiseDBValidation") << " Bad percentual occupancy for sector : " << (*channel).getBin() << "  Contents : " << (*channel).getContents();
      }
  }
  //testCriterionName = parameters.getUntrackedParameter<string>("layerTestName","noiseLayerOccInRange");
  const QReport * theDiffQReport5 = layerHisto_->getQReport(layerTestName_);
  if(theDiffQReport5) {
      vector<dqm::me_util::Channel> badChannels = theDiffQReport5->getBadChannels();
      for (vector<dqm::me_util::Channel>::iterator channel = badChannels.begin(); 
           channel != badChannels.end(); channel++) {
        if((*channel).getBin()==1)
          LogWarning("NoiseDBValidation") << " Bad percentual occupancy for the first 10 wires! Contents : " << (*channel).getContents();
        if((*channel).getBin()==2)
          LogWarning("NoiseDBValidation") << " Bad percentual occupancy for the middle wires! Contents : " << (*channel).getContents();
        if((*channel).getBin()==3)
          LogWarning("NoiseDBValidation") << " Bad percentual occupancy for the last 10 wires! Contents : "<<(*channel).getContents();
      }
  }

}