#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 edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
	Analyze. More...

void	beginJob ()
	BeginJob. More...

void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)

void	beginRun (const edm::Run &r, const edm::EventSetup &c)
	BeginRun. More...

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

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

void	endJob ()
	Endjob. More...

void	endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &c)
	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::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Attributes
DQMStore *	dbe

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

DQM Test Client

Date:
2010/01/05 10:15:46

Revision:
1.8

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

Definition at line 52 of file DTNoiseTest.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 38 of file DTNoiseTest.cc.

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

                                                  {
 
   edm::LogVerbatim ("noise") <<"[DTNoiseTest]: Constructor";  
 
   parameters = ps;
   
   dbe = edm::Service<DQMStore>().operator->();
   dbe->setCurrentFolder("DT/Tests/Noise");
 
   prescaleFactor = parameters.getUntrackedParameter<int>("diagnosticPrescale", 1);
 
 }

DTNoiseTest::~DTNoiseTest ( )

virtual

Destructor.

Definition at line 53 of file DTNoiseTest.cc.

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

Member Function Documentation

void DTNoiseTest::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 91 of file DTNoiseTest.cc.

                                                                         {
 
   updates++;
   edm::LogVerbatim ("noise") << "[DTNoiseTest]: "<<updates<<" events";
 
 }

void DTNoiseTest::beginJob ( void )

protectedvirtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 61 of file DTNoiseTest.cc.

                           {
 
   edm::LogVerbatim ("noise") <<"[DTNoiseTest]: BeginJob";
 
   updates = 0;
 
 }

void DTNoiseTest::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	context
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 81 of file DTNoiseTest.cc.

References edm::LuminosityBlockBase::run(), and DTTTrigCorrFirst::run.

                                                                                                 {
 
   edm::LogVerbatim ("noise") <<"[DTNoiseTest]: Begin of LS transition";
 
   // Get the run number
   run = lumiSeg.run();
 
 }

void DTNoiseTest::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 70 of file DTNoiseTest.cc.

References edm::EventSetup::get().

                                                                          {
 
   edm::LogVerbatim ("noise") <<"[DTNoiseTest]: BeginRun";
 
   // Get the geometry
   context.get<MuonGeometryRecord>().get(muonGeom);
 
 }

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

protected

book the new ME

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

protected

void DTNoiseTest::endJob ( void )

protectedvirtual

Endjob.

Reimplemented from edm::EDAnalyzer.

Definition at line 262 of file DTNoiseTest.cc.

                         {
 
   edm::LogVerbatim ("noise") <<"[DTNoiseTest] endjob called!";
   
   //if ( parameters.getUntrackedParameter<bool>("writeHisto", true) ) 
   //  dbe->save(parameters.getUntrackedParameter<string>("outputFile", "DTNoiseTest.root"));
   
   dbe->rmdir("DT/Tests/Noise");
 }

void DTNoiseTest::endLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	c
	)

protectedvirtual

DQM Client Diagnostic.

Reimplemented from edm::EDAnalyzer.

Definition at line 99 of file DTNoiseTest.cc.

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

                                                                                               {
 
   // counts number of updats (online mode) or number of events (standalone mode)
   //updates++;
   // if running in standalone perform diagnostic only after a reasonalbe amount of events
   //if ( parameters.getUntrackedParameter<bool>("runningStandalone", false) &&
   //     updates%parameters.getUntrackedParameter<int>("diagnosticPrescale", 1000) != 0 ) return;
 
   //edm::LogVerbatim ("noise") <<"[DTNoiseTest]: "<<updates<<" updates";
 
 
   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<DTChamber*>::const_iterator ch_it = muonGeom->chambers().begin();
   vector<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 = dbe->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((*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((*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 = dbe->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((*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();
     // FIXME: getMessage() sometimes returns and invalid string (null pointer inside QReport data member)
     // LogVerbatim ("tTrigCalibration") << "-------- LayerId : "<<(*hMean).first<<"  "<<theMeanQReport->getMessage()<<" ------- "<<theMeanQReport->getStatus(); 
       }
     }
   }
   
 }

string DTNoiseTest::getMEName ( const DTChamberId & ch )

protected

Get the ME name.

Definition at line 273 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 297 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

DQMStore* DTNoiseTest::dbe

private

Definition at line 102 of file DTNoiseTest.h.

bool DTNoiseTest::debug

private

Definition at line 96 of file DTNoiseTest.h.

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

private

Definition at line 115 of file DTNoiseTest.h.

edm::ESHandle<DTGeometry> DTNoiseTest::muonGeom

private

Definition at line 105 of file DTNoiseTest.h.

unsigned int DTNoiseTest::nLumiSegs

private

Definition at line 98 of file DTNoiseTest.h.

edm::ParameterSet DTNoiseTest::parameters

private

Definition at line 104 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 99 of file DTNoiseTest.h.

int DTNoiseTest::run

private

Definition at line 100 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 111 of file DTNoiseTest.h.

edm::ESHandle<DTTtrig> DTNoiseTest::tTrigMap

private

Definition at line 106 of file DTNoiseTest.h.

int DTNoiseTest::updates

private

Definition at line 97 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