#include <SiPixelInformationExtractor.h>

Public Member Functions
void	bookNoisyPixels (DQMStore *bei, float noiseRate, bool Tier0Flag)

void	computeStatus (MonitorElement *mE, double &colorValue, std::pair< double, double > &norm)
	(Documentation under construction). More...

void	createImages (DQMStore *bei)

void	findNoisyPixels (DQMStore *bei, bool init, float noiseRate, int noiseRateDenominator, edm::EventSetup const &eSetup)

int	getDetId (MonitorElement *mE)
	(Documentation under construction). More...

void	getMEList (DQMStore *bei, std::map< std::string, int > &mEHash)
	(Documentation under construction). More...

std::string	getMEType (MonitorElement *mE)

void	getNormalization (MonitorElement *mE, std::pair< double, double > &norm, std::string theMEType)
	(Documentation under construction). More...

void	getNormalization2D (MonitorElement *mE, std::pair< double, double > &normX, std::pair< double, double > &normY, std::string theMEType)
	(Documentation under construction). More...

void	readConfiguration ()
	Read Configuration File. More...

bool	readConfiguration (std::map< std::string, std::vector< std::string > > &layoutMap, std::map< std::string, std::map< std::string, std::string > > &qtestsMap, std::map< std::string, std::vector< std::string > > &meQTestsMap)

void	selectMEList (DQMStore bei, std::string &name, std::vector< MonitorElement > &mes)
	(Documentation under construction). More...

	SiPixelInformationExtractor (bool offlineXMLfile)
	Constructor of the SiPixelInformationExtractor class. More...

	~SiPixelInformationExtractor ()
	Destructor of the SiPixelInformationExtractor class. More...

Private Member Functions
void	createDummiesFromLayout ()

void	fillModuleAndHistoList (DQMStore *bei, std::vector< std::string > &modules, std::map< std::string, std::string > &histos)
	(Documentation under construction). More...

void	getItemList (const std::multimap< std::string, std::string > &req_map, std::string item_name, std::vector< std::string > &items)
	(Documentation under construction). More...

std::string	getItemValue (const std::multimap< std::string, std::string > &req_map, std::string item_name)
	(Documentation under construction). More...

std::string	getItemValue (std::multimap< std::string, std::string > &req_map, std::string item_name)

bool	goToDir (DQMStore *bei, std::string &sname)
	(Documentation under construction). More...

bool	hasItem (std::multimap< std::string, std::string > &req_map, std::string item_name)
	(Documentation under construction). More...

void	printAlarmList (DQMStore *bei, std::ostringstream &str_val)
	(Documentation under construction). More...

void	printModuleHistoList (DQMStore *bei, std::ostringstream &str_val)
	(Documentation under construction). More...

void	printSummaryHistoList (DQMStore *bei, std::ostringstream &str_val)
	(Documentation under construction). More...

void	selectColor (std::string &col, int status)

void	selectColor (std::string &col, std::vector< QReport * > &reports)

void	selectImage (std::string &name, int status)

void	selectImage (std::string &name, std::vector< QReport * > &reports)

Private Attributes
SiPixelActionExecutor *	actionExecutor_

int	alarmCounter_

MonitorElement *	BarrelNdigisFREQProjection

SiPixelConfigParser *	configParser_

SiPixelConfigWriter *	configWriter_

int	count

MonitorElement *	EndcapNdigisFREQProjection

bool	endOfModules_

int	errcount

MonitorElement *	EventRateBarrelPixels

MonitorElement *	EventRateEndcapPixels

bool	flagHotModule_

bool	gotDigis

SiPixelHistoPlotter *	histoPlotter_

std::map< std::string, std::vector< std::string > >	layoutMap

SiPixelLayoutParser *	layoutParser_

std::map< std::string, std::vector< std::string > >	meQTestsMap

std::ofstream	myfile_

int	nevents_

std::map< uint32_t, std::vector< std::pair < std::pair< int, int >, float > > >	noisyDetIds_

bool	offlineXMLfile_

std::map< std::string, std::map< std::string, std::string > >	qtestsMap

bool	readMeMap_

bool	readQTestMap_

bool	readReference_

edm::ESHandle < SiPixelFedCablingMap >	theCablingMap

Detailed Description

Definition at line 41 of file SiPixelInformationExtractor.h.

Constructor & Destructor Documentation

SiPixelInformationExtractor::SiPixelInformationExtractor ( bool offlineXMLfile )

Constructor of the SiPixelInformationExtractor class.

Definition at line 67 of file SiPixelInformationExtractor.cc.

References histoPlotter_, and readReference_.

                                                                             : offlineXMLfile_(offlineXMLfile) {
   edm::LogInfo("SiPixelInformationExtractor") << 
     " Creating SiPixelInformationExtractor " << "\n" ;
   
   readReference_ = false;
   histoPlotter_=0;
   histoPlotter_ = new SiPixelHistoPlotter();
 }

SiPixelInformationExtractor::~SiPixelInformationExtractor ( )

Destructor of the SiPixelInformationExtractor class.

Definition at line 80 of file SiPixelInformationExtractor.cc.

References histoPlotter_.

                                                           {
   edm::LogInfo("SiPixelInformationExtractor") << 
     " Deleting SiPixelInformationExtractor " << "\n" ;
   
   if (histoPlotter_) delete histoPlotter_;
 }

Member Function Documentation

void SiPixelInformationExtractor::bookNoisyPixels	(	DQMStore *	bei,
		float	noiseRate,
		bool	Tier0Flag
	)

Definition at line 1346 of file SiPixelInformationExtractor.cc.

References DQMStore::book1D(), DQMStore::cd(), EventRateBarrelPixels, EventRateEndcapPixels, MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().

                                                                                                  {
 //std::cout<<"BOOK NOISY PIXEL MEs!"<<std::endl;
   bei->cd();
   if(noiseRate_>=0.){
     bei->setCurrentFolder("Pixel/Barrel");
     EventRateBarrelPixels = bei->book1D("barrelEventRate","Digi event rate for all Barrel pixels",1000,0.,0.01);
     EventRateBarrelPixels->setAxisTitle("Event Rate",1);
     EventRateBarrelPixels->setAxisTitle("Number of Pixels",2);
     bei->cd();  
     bei->setCurrentFolder("Pixel/Endcap");
     EventRateEndcapPixels = bei->book1D("endcapEventRate","Digi event rate for all Endcap pixels",1000,0.,0.01);
     EventRateEndcapPixels->setAxisTitle("Event Rate",1);
     EventRateEndcapPixels->setAxisTitle("Number of Pixels",2);
   }
 }

void SiPixelInformationExtractor::computeStatus	(	MonitorElement *	theME,
		double &	colorValue,
		std::pair< double, double > &	norm
	)

(Documentation under construction).

Definition at line 912 of file SiPixelInformationExtractor.cc.

References MonitorElement::getMean(), getMEType(), getNormalization(), getNormalization2D(), reco::return(), and mathSSE::sqrt().

 {
   double normalizationX = 1 ;
   double normalizationY = 1 ;
   double meanX          = 0 ;
   double meanY          = 0 ;
   
   colorValue = 0 ;
 
   pair<double,double> normX ;
   pair<double,double> normY ;
 
   string theMEType = getMEType(theME) ;
 
 //   cout << ACRed << ACReverse
 //        << "[SiPixelInformationExtractor::computeStatus()]"
 //        << ACPlain
 //        << " Computing average for "
 //        << theME->getName()
 //        << endl ;
 
   if( theMEType.find("TH1") != string::npos)
   {
    meanX = (double)theME->getMean();
    getNormalization(theME, normX, "TH1") ;
    normalizationX = fabs( normX.second - normX.first) ;
    if( normalizationX == 0 ) {normalizationX=1.E-20;}
    colorValue  = meanX / normalizationX ;
    norm.first  = normX.first ;
    norm.second = normX.second ;
   }
   
   if( theMEType.find("TH2") != string::npos)
   {
    meanX = (double)theME->getMean(1);
    meanY = (double)theME->getMean(2);
    getNormalization2D(theME, normX, normY, "TH2") ;
    normalizationX = fabs( normX.second - normX.first) ;
    normalizationY = fabs( normY.second - normY.first) ;
    if( normalizationX == 0 ) {normalizationX=1.E-20;}
    if( normalizationY == 0 ) {normalizationY=1.E-20;}
    double cVX = meanX / normalizationX ;
    double cVY = meanY / normalizationY ;
    colorValue = sqrt(cVX*cVX + cVY*cVY) ;
    if( normalizationX >= normalizationY )
    { 
     norm.first  = normX.first;
     norm.second = normX.second ;
    } else { 
     norm.first  = normY.first;
     norm.second = normY.second ;
    }
 //   cout << ACBlue << ACBold << ACReverse
 //        << "[SiPixelInformationExtractor::computeStatus()]"
 //  << ACPlain << "    "
 //  << theME->getName()
 //  << " meanX:Y "
 //  << meanX << ":" << meanY
 //  << " normX:Y " 
 //  << norm.first << ":" << norm.second
 //  << endl ;
   } 
  
   return ;
 }

void SiPixelInformationExtractor::createDummiesFromLayout ( )

private

void SiPixelInformationExtractor::createImages ( DQMStore * bei )

Definition at line 1604 of file SiPixelInformationExtractor.cc.

References SiPixelHistoPlotter::createPlots(), and histoPlotter_.

                                                            {
   histoPlotter_->createPlots(bei);
 }

void SiPixelInformationExtractor::fillModuleAndHistoList	(	DQMStore *	bei,
		std::vector< std::string > &	modules,
		std::map< std::string, std::string > &	histos
	)

private

(Documentation under construction).

This method

(Documentation under construction).

This method

Definition at line 318 of file SiPixelInformationExtractor.cc.

References DQMStore::cd(), relmon_rootfiles_spy::contents, spr::find(), DQMStore::get(), DQMStore::getMEs(), MonitorElement::getRootObject(), DQMStore::getSubdirs(), DQMStore::goUp(), and DQMStore::pwd().

                                                              {
 //cout<<"entering SiPixelInformationExtractor::fillModuleAndHistoList"<<endl;
   string currDir = bei->pwd();
   //cout<<"currDir= "<<currDir<<endl;
   if(currDir.find("Module_") != string::npos){
     if(histos.size() == 0){
       vector<string> contents = bei->getMEs();
       for (vector<string>::const_iterator it = contents.begin(); it != contents.end(); it++) {
     string hname          = (*it).substr(0, (*it).find("_siPixel"));
     if(hname==" ") hname = (*it).substr(0, (*it).find("_generalTracks"));
         string fullpathname   = bei->pwd() + "/" + (*it); 
        // cout<<"fullpathname="<<fullpathname<<endl;
         MonitorElement * me   = bei->get(fullpathname);
         string htype          = "undefined" ;
         if(me) htype = me->getRootObject()->IsA()->GetName() ;
     //cout<<"hname="<<hname<<endl;
     //if(htype=="TH1F" || htype=="TH1D"){
         histos[hname] = htype ;
         string mId=" ";
     if(hname.find("ndigis")            !=string::npos) mId = (*it).substr((*it).find("ndigis_siPixelDigis_")+20, 9);
     if(mId==" " && hname.find("nclusters") !=string::npos) mId = (*it).substr((*it).find("nclusters_siPixelClusters_")+26, 9);
         if(mId==" " && hname.find("residualX") !=string::npos) mId = (*it).substr((*it).find("residualX_ctfWithMaterialTracks_")+32, 9);
         if(mId==" " && hname.find("NErrors") !=string::npos) mId = (*it).substr((*it).find("NErrors_siPixelDigis_")+21, 9);
         if(mId==" " && hname.find("ClustX") !=string::npos) mId = (*it).substr((*it).find("ClustX_siPixelRecHit_")+21, 9);
         if(mId==" " && hname.find("pixelAlive") !=string::npos) mId = (*it).substr((*it).find("pixelAlive_siPixelCalibDigis_")+29, 9);
         if(mId==" " && hname.find("Gain1d") !=string::npos) mId = (*it).substr((*it).find("Gain1d_siPixelCalibDigis_")+25, 9);
         if(mId!=" ") modules.push_back(mId);
         //cout<<"mId="<<mId<<endl;
     //}
       }    
     }
   } else {  
     vector<string> subdirs = bei->getSubdirs();
     for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
       if((bei->pwd()).find("Barrel")==string::npos && (bei->pwd()).find("Endcap")==string::npos) bei->goUp();
       bei->cd(*it);
       fillModuleAndHistoList(bei, modules, histos);
       bei->goUp();
     }
   }
 //  fillBarrelList(bei, modules, histos);
 //cout<<"leaving SiPixelInformationExtractor::fillModuleAndHistoList"<<endl;
 }

void SiPixelInformationExtractor::findNoisyPixels	(	DQMStore *	bei,
		bool	init,
		float	noiseRate,
		int	noiseRateDenominator,
		edm::EventSetup const &	eSetup
	)

Definition at line 1365 of file SiPixelInformationExtractor.cc.

References DQMStore::cd(), counter, sipixelobjects::ElectronicIndex::dcol, sipixelobjects::LocalPixel::DcolPxid::dcol, cond::rpcobgas::detid, RefreshWebPage::dname, endOfModules_, EventRateBarrelPixels, EventRateEndcapPixels, NewTree::fid, MonitorElement::Fill(), edm::EventSetup::get(), DQMStore::get(), getDetId(), DQMStore::getMEs(), MonitorElement::getPathname(), DQMStore::getSubdirs(), MonitorElement::getTH2F(), DQMStore::goUp(), i, sipixelobjects::PixelROC::idInDetUnit(), PixelBarrelName::isHalfModule(), j, sipixelobjects::ElectronicIndex::link, python.multivaluedict::map(), myfile_, PixelBarrelName::name(), PixelEndcapName::name(), nevents_, noisyDetIds_, getHLTPrescaleColumns::path, edm::ESHandle< class >::product(), DQMStore::pwd(), sipixelobjects::ElectronicIndex::pxid, sipixelobjects::LocalPixel::DcolPxid::pxid, sipixelobjects::ElectronicIndex::roc, sipixelobjects::LocalPixel::rocCol(), sipixelobjects::LocalPixel::rocRow(), DQMStore::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), theCablingMap, and relativeConstraints::value.

 {
 //cout<<"Entering SiPixelInformationExtractor::findNoisyPixels with noiseRate set to "<<noiseRate_<<endl;
 
   
   if(init){
     endOfModules_=false;
     nevents_=noiseRateDenominator_;
     if(nevents_ == -1){
       bei->cd();
       bei->setCurrentFolder("Pixel/EventInfo");
       nevents_ = (bei->get("Pixel/EventInfo/processedEvents"))->getIntValue();
     }
     bei->cd();  
     myfile_.open ("NoisyPixelList.txt", ios::app);
     myfile_ << "Noise summary, ran over " << nevents_ << " events, threshold was set to " << noiseRate_ <<  std::endl;
   }
   string currDir = bei->pwd();
   string dname = currDir.substr(currDir.find_last_of("/")+1);
 
 
   if(dname.find("Module_")!=string::npos){
     vector<string> meVec = bei->getMEs();
     for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
       string full_path = currDir + "/" + (*it);
       if(full_path.find("hitmap_siPixelDigis")!=string::npos){
         //broken HV bond:
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_01/Panel_2/Module_2")!=string::npos) continue;
         //?noisy?
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_12/Panel_1/Module_4")!=string::npos) continue;
         //ROG with HV problem (short?):
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_10/Panel_1/Module_3")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_10/Panel_1/Module_4")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_10/Panel_2/Module_2")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_10/Panel_2/Module_3")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_11/Panel_1/Module_3")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_11/Panel_1/Module_4")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_11/Panel_2/Module_2")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_11/Panel_2/Module_3")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_12/Panel_1/Module_3")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_12/Panel_1/Module_4")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_12/Panel_2/Module_2")!=string::npos) continue;
     //if(currDir.find("HalfCylinder_mI/Disk_1/Blade_12/Panel_2/Module_3")!=string::npos) continue;
         MonitorElement * me = bei->get(full_path);
         if (!me) continue;
     int detid=getDetId(me); int pixcol=-1; int pixrow=-1; 
 
     //cout<<"detid= "<<detid<<endl;
     std::vector<std::pair<std::pair<int, int>, float> > noisyPixelsInModule;
     TH2F * hothisto = me->getTH2F();
     if(hothisto){
       for(int i=1; i!=hothisto->GetNbinsX()+1; i++){
         for(int j=1; j!=hothisto->GetNbinsY()+1; j++){
           float value = (hothisto->GetBinContent(i,j))/float(nevents_);
           if(me->getPathname().find("Barrel")!=string::npos){
             EventRateBarrelPixels = bei->get("Pixel/Barrel/barrelEventRate");
             if(EventRateBarrelPixels) EventRateBarrelPixels->Fill(value);
           }else if(me->getPathname().find("Endcap")!=string::npos){
             EventRateEndcapPixels = bei->get("Pixel/Endcap/endcapEventRate");
             if(EventRateEndcapPixels) EventRateEndcapPixels->Fill(value);
           }
           if(value > noiseRate_){
             pixcol = i-1;
             pixrow = j-1;
         //cout<<"pixcol= "<<pixcol<<" , pixrow= "<<pixrow<<" , value= "<<value<<endl;
  
             std::pair<int, int> address(pixcol, pixrow);
             std::pair<std::pair<int, int>, float>  PixelStats(address, value);
             noisyPixelsInModule.push_back(PixelStats);
           }
             }
       }
     }
     noisyDetIds_[detid] = noisyPixelsInModule;
     //if(noisyPixelsInModule.size()>=20) cout<<"This module has 20 or more hot pixels: "<<detid<<","<<bei->pwd()<<","<<noisyPixelsInModule.size()<<endl;
       }
     }
   }
   vector<string> subDirVec = bei->getSubdirs();  
   for (vector<string>::const_iterator ic = subDirVec.begin();
        ic != subDirVec.end(); ic++) {
     if((*ic).find("AdditionalPixelErrors")!=string::npos) continue;
     bei->cd(*ic);
     init=false;
     findNoisyPixels(bei,init,noiseRate_,noiseRateDenominator_,eSetup);
     bei->goUp();
   }
 
   if(bei->pwd().find("EventInfo")!=string::npos) endOfModules_ = true;
   
   if(!endOfModules_) return;
   // myfile_ <<"am in "<<bei->pwd()<<" now!"<<endl;
   if(currDir == "Pixel/EventInfo/reportSummaryContents"){
     eSetup.get<SiPixelFedCablingMapRcd>().get(theCablingMap);
     std::vector<std::pair<sipixelobjects::DetectorIndex,double> > pixelvec;
     std::map<uint32_t,int> myfedmap;
     std::map<uint32_t,std::string> mynamemap;
     int realfedID = -1;
     //int Nnoisies = noisyDetIds_.size();
     //cout<<"Number of noisy modules: "<<Nnoisies<<endl;
     int counter = 0;
     int n_noisyrocs_all = 0;
     int n_noisyrocs_barrel = 0;
     int n_noisyrocs_endcap = 0;
     int n_verynoisyrocs_all = 0;
     int n_verynoisyrocs_barrel = 0;
     int n_verynoisyrocs_endcap = 0;
 
     for(int fid = 0; fid < 40; fid++){
     for(std::map<uint32_t, std::vector< std::pair<std::pair<int, int>, float> > >::const_iterator it = noisyDetIds_.begin(); 
         it != noisyDetIds_.end(); it++){
       uint32_t detid = (*it).first;
       std::vector< std::pair<std::pair<int, int>, float> > noisyPixels = (*it).second;
       //cout<<noisyPixels.size()<<" noisy pixels in a module: "<<detid<<endl;
       // now convert into online conventions:
       for(int fedid=0; fedid<=40; ++fedid){
     SiPixelFrameConverter converter(theCablingMap.product(),fedid);
     uint32_t newDetId = detid;
     if(converter.hasDetUnit(newDetId)){
       realfedID=fedid;
       break;   
     }
       }
       if(fid == realfedID){
       //cout<<"FED ID is = "<<realfedID<<endl;
       if(realfedID==-1) continue; 
       DetId detId(detid);
       uint32_t detSubId = detId.subdetId();
       std::string outputname;
       bool HalfModule = false;
       if (detSubId == 2){   //FPIX
     PixelEndcapName nameworker(detid);
     outputname = nameworker.name();
       } else if(detSubId == 1){   //BPIX
     PixelBarrelName nameworker(detid);
     outputname = nameworker.name();
     HalfModule = nameworker.isHalfModule();
 
       } else{
     continue;
       } 
       std::map<int,int> myrocmap;
       myfedmap[detid]=realfedID;
       mynamemap[detid]=outputname;
       
       for(std::vector< std::pair< std::pair<int,int>, float> >::const_iterator pxl = noisyPixels.begin(); 
           pxl != noisyPixels.end(); pxl++){
         std::pair<int,int> offlineaddress = (*pxl).first;
     float Noise_frac = (*pxl).second;
     int offlineColumn = offlineaddress.first;
         int offlineRow = offlineaddress.second;
         counter++;
         //cout<<"noisy pixel counter: "<<counter<<endl;
 
         sipixelobjects::ElectronicIndex cabling; 
         SiPixelFrameConverter formatter(theCablingMap.product(),realfedID);
         sipixelobjects::DetectorIndex detector = {detid, offlineRow, offlineColumn};      
     formatter.toCabling(cabling,detector);
         // cabling should now contain cabling.roc and cabling.dcol  and cabling.pxid
         // however, the coordinates now need to be converted from dcl,pxid to the row,col coordinates used in the calibration info 
         sipixelobjects::LocalPixel::DcolPxid loc;
         loc.dcol = cabling.dcol;
         loc.pxid = cabling.pxid;
     
     
     // OLD version, not 31X compatible:
 //        const sipixelobjects::PixelFEDCabling *theFed= theCablingMap.product()->fed(realfedID);
 //  const sipixelobjects::PixelFEDLink * link = theFed->link(cabling.link);
 //  const sipixelobjects::PixelROC *theRoc = link->roc(cabling.roc);
 //  sipixelobjects::LocalPixel locpixel(loc);
     
     
     // FIX to adhere to new cabling map. To be replaced with CalibTracker/SiPixelTools detid - > hardware id classes ASAP.
     //        const sipixelobjects::PixelFEDCabling *theFed= theCablingMap.product()->fed(realfedID);
     //        const sipixelobjects::PixelFEDLink * link = theFed->link(cabling.link);
     //        const sipixelobjects::PixelROC *theRoc = link->roc(cabling.roc);
         sipixelobjects::LocalPixel locpixel(loc);
         assert(realfedID >= 0);
         assert(cabling.link >= 0);
         assert(cabling.roc >= 0);
     sipixelobjects::CablingPathToDetUnit path = {static_cast<unsigned int>(realfedID), 
                                                      static_cast<unsigned int>(cabling.link),
                                                      static_cast<unsigned int>(cabling.roc)};  
     const sipixelobjects::PixelROC *theRoc = theCablingMap->findItem(path);
     // END of FIX
     
         int onlineColumn = locpixel.rocCol();
         int onlineRow= locpixel.rocRow();
     myrocmap[(theRoc->idInDetUnit())]++;
 
     // ROC numbers in the barrel go from 8 to 15 instead of 0 to 7 in half modules.  This is a 
     // fix to get the roc number, and add 8 to it if:
     // it's a Barrel module AND on the minus side AND a Half module
 
     int rocnumber = -1;
 
     if((detSubId == 1) && (outputname.find("mO")!=string::npos || outputname.find("mI")!=string::npos) && (HalfModule)){
       rocnumber = theRoc->idInDetUnit() + 8;
     }
     else{
       rocnumber = theRoc->idInDetUnit();
     }
 
         //cout<<counter<<" : \t detid= "<<detid<<" , OFF col,row= "<<offlineColumn<<","<<offlineRow<<" , ON roc,col,row= "<<theRoc->idInDetUnit()<<","<<onlineColumn<<","<<onlineRow<<endl;
         myfile_ <<"NAME: "<<outputname<<" , DETID: "<<detid<<" , OFFLINE: col,row: "<<offlineColumn<<","<<offlineRow<<"  \t , ONLINE: roc,col,row: "<<rocnumber<<","<<onlineColumn<<","<<onlineRow<< "  \t , fed,dcol,pixid,link: "<<realfedID<<","<<loc.dcol<<","<<loc.pxid<<","<<cabling.link << ", Noise fraction: " << Noise_frac << std::endl;
       }
       for(std::map<int, int>::const_iterator nrc = myrocmap.begin(); nrc != myrocmap.end(); nrc++){
     if((*nrc).second > 0){
       n_noisyrocs_all++;
       if(detSubId == 2){
         n_noisyrocs_endcap++;
       } else if(detSubId == 1){
         n_noisyrocs_barrel++;}
     }
     if((*nrc).second > 40){
       n_verynoisyrocs_all++;
       if(detSubId == 2){
         n_verynoisyrocs_endcap++;
       } else if(detSubId == 1){
         n_verynoisyrocs_barrel++;}
     }
       }
       }
 
     }
     }
     myfile_ << "There are " << n_noisyrocs_all << " noisy ROCs (ROCs with at least 1 noisy pixel) in the entire detector. " << n_noisyrocs_endcap << " are in the FPIX and " << n_noisyrocs_barrel << " are in the BPIX. " << endl;
     myfile_ << "There are " << n_verynoisyrocs_all << " highly noisy ROCs (ROCs with at least 10% of all pixels passing the noise threshold) in the entire detector. " << n_verynoisyrocs_endcap << " are in the FPIX and " << n_verynoisyrocs_barrel << " are in the BPIX. " << endl;
 
   }
   myfile_.close();
 //cout<<"...leaving SiPixelInformationExtractor::findNoisyPixels!"<<endl;
   return;
 }

int SiPixelInformationExtractor::getDetId ( MonitorElement * mE )

(Documentation under construction).

Given a pointer to ME returns the associated detId

Definition at line 1224 of file SiPixelInformationExtractor.cc.

References MonitorElement::getName().

Referenced by findNoisyPixels().

 {
 //cout<<"In SiPixelInformationExtractor::getDetId: for ME= "<<mE->getName()<<endl;
  string mEName = mE->getName();
 
  int detId = 0;
  
  if( mEName.find("_3") != string::npos )
  {
   string detIdString = mEName.substr((mEName.find_last_of("_"))+1,9);
   //cout<<"string: "<<detIdString<<endl;
   std::istringstream isst;
   isst.str(detIdString);
   isst>>detId;
 // } else {
 //  cout << ACYellow << ACBold
 //       << "[SiPixelInformationExtractor::getDetId()] "
 //       << ACPlain
 //       << "Could not extract detId from "
 //       << mEName
 //       << endl ;
  }
   //cout<<"returning with: "<<detId<<endl;
   return detId ;
   
 }

void SiPixelInformationExtractor::getItemList	(	const std::multimap< std::string, std::string > &	req_map,
		std::string	item_name,
		std::vector< std::string > &	items
	)

private

(Documentation under construction).

This method

Definition at line 685 of file SiPixelInformationExtractor.cc.

                                                  {
 //cout<<"entering SiPixelInformationExtractor::getItemList"<<endl;
   items.clear();
   for (multimap<string, string>::const_iterator it = req_map.begin();
        it != req_map.end(); it++) {
     if (it->first == item_name) {
       items.push_back(it->second);
     }
   }
 //cout<<"leaving SiPixelInformationExtractor::getItemList"<<endl;
 }

std::string SiPixelInformationExtractor::getItemValue	(	const std::multimap< std::string, std::string > &	req_map,
		std::string	item_name
	)

private

(Documentation under construction).

This method

Definition at line 718 of file SiPixelInformationExtractor.cc.

References AlCaHLTBitMon_QueryRunRegistry::string, and relativeConstraints::value.

                                              {
 //cout<<"entering SiPixelInformationExtractor::getItemValue for item: "<<item_name<<endl;
   std::multimap<std::string,std::string>::const_iterator pos = req_map.find(item_name);
   std::string value = " ";
   if (pos != req_map.end()) {
     value = pos->second;
   }
 //  cout<<"value = "<<value<<endl;
   return value;
 //cout<<"leaving SiPixelInformationExtractor::getItemValue"<<endl;
 }

std::string SiPixelInformationExtractor::getItemValue	(	std::multimap< std::string, std::string > &	req_map,
		std::string	item_name
	)

private

Definition at line 730 of file SiPixelInformationExtractor.cc.

References AlCaHLTBitMon_QueryRunRegistry::string, and relativeConstraints::value.

                                              {
 //cout<<"entering SiPixelInformationExtractor::getItemValue for item: "<<item_name<<endl;
   std::multimap<std::string,std::string>::iterator pos = req_map.find(item_name);
   std::string value = " ";
   if (pos != req_map.end()) {
 //  cout<<"item found!"<<endl;
     value = pos->second;
   }
 //  cout<<"value = "<<value<<endl;
   return value;
 //cout<<"leaving SiPixelInformationExtractor::getItemValue"<<endl;
 }

void SiPixelInformationExtractor::getMEList	(	DQMStore *	bei,
		std::map< std::string, int > &	mEHash
	)

(Documentation under construction).

Definition at line 1255 of file SiPixelInformationExtractor.cc.

References ACBold, ACPlain, ACRed, DQMStore::cd(), relmon_rootfiles_spy::contents, gather_cfg::cout, DQMStore::getMEs(), DQMStore::getSubdirs(), DQMStore::goUp(), and DQMStore::pwd().

 {
   string currDir = bei->pwd();
    
 //   cout << ACRed << ACBold
 //        << "[SiPixelInformationExtractor::getMEList()]"
 //        << ACPlain
 //        << " Requesting ME list in " 
 //        << currDir
 //        << endl ;
        
   string theME ;
    
   // Get ME from Collector/FU0/Tracker/PixelEndcap/HalfCylinder_pX/Disk_X/Blade_XX/Panel_XX/Module_XX
   if (currDir.find("Module_") != string::npos ||
       currDir.find("FED_") != string::npos)  
   {
     vector<string> contents = bei->getMEs(); 
        
     for (vector<string>::const_iterator it = contents.begin(); it != contents.end(); it++) 
     {
       theME = (*it) ;
 //       cout << ACRed << ACReverse
 //            << "[SiPixelInformationExtractor::getMEList()]"
 //            << ACPlain
 //            << " ME: " 
 //            << (*it)
 //            << endl ;
       if(theME.find("siPixel")==string::npos && theME.find("ctfWithMaterialTracks")==string::npos) 
       {
        cout << ACRed << ACBold
             << "[SiPixelInformationExtractor::getMEList()]"
         << ACPlain
         << " ----> Skipping " 
         << (*it)
         << endl ;
        continue ;
       } // If the ME is not a Pixel one, skip
       string full_path = currDir + "/" + (*it);
       string mEName = theME.substr(0,theME.find_first_of("_"));
       mEHash[mEName]++ ;
     }
     
     return;
   } else {  // If not yet reached the desired level in the directory tree, recursively go down one level more
     vector<string> subdirs = bei->getSubdirs();
     for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) 
     {
       bei->cd(*it);
       getMEList(bei, mEHash);
       bei->goUp();
     }
   }
 }

std::string SiPixelInformationExtractor::getMEType ( MonitorElement * mE )

Definition at line 260 of file SiPixelInformationExtractor.cc.

References MonitorElement::getRootObject().

Referenced by computeStatus().

 {
   string qtype = theMe->getRootObject()->IsA()->GetName() ;
   if(         qtype.find("TH1") != string::npos )
   {
     return "TH1" ;
   } else if ( qtype.find("TH2") != string::npos  ) {
     return "TH2" ;
   } else if ( qtype.find("TH3") != string::npos ) {
     return "TH3" ;
   }
   return "TH1" ;
 }

void SiPixelInformationExtractor::getNormalization	(	MonitorElement *	theME,
		std::pair< double, double > &	norm,
		std::string	theMEType
	)

(Documentation under construction).

Definition at line 984 of file SiPixelInformationExtractor.cc.

References MonitorElement::getNbinsX().

Referenced by computeStatus().

 {
   double normLow  = 0 ;
   double normHigh = 0 ;
 
   if( theMEType.find("TH1") != string::npos)
   {
    normHigh    = (double)theME->getNbinsX() ;
    norm.first  = normLow  ;
    norm.second = normHigh ;
   }
 }

void SiPixelInformationExtractor::getNormalization2D	(	MonitorElement *	theME,
		std::pair< double, double > &	normX,
		std::pair< double, double > &	normY,
		std::string	theMEType
	)

(Documentation under construction).

Definition at line 1003 of file SiPixelInformationExtractor.cc.

References MonitorElement::getNbinsX(), and MonitorElement::getNbinsY().

Referenced by computeStatus().

 {
   double normLow  = 0 ;
   double normHigh = 0 ;
 
   if( theMEType.find("TH2") != string::npos )
   {
    normHigh    = (double)theME->getNbinsX() ;
    normX.first  = normLow  ;
    normX.second = normHigh ;
    normHigh    = (double)theME->getNbinsY() ;
    normY.first  = normLow  ;
    normY.second = normHigh ;
 //   cout << ACCyan << ACBold << ACReverse
 //        << "[SiPixelInformationExtractor::getNormalization2D()]"
 //  << ACPlain << " "
 //  << theME->getName()
 //  << " normX: " 
 //  << normX.first << ":" << normX.second
 //  << " normY: " 
 //  << normY.first << ":" << normY.second
 //  << endl ;
   }
 }

bool SiPixelInformationExtractor::goToDir	(	DQMStore *	bei,
		std::string &	sname
	)

private

(Documentation under construction).

This method

Definition at line 830 of file SiPixelInformationExtractor.cc.

References DQMStore::cd(), TrackerOfflineValidation_Dqm_cff::dirName, and DQMStore::pwd().

                                                         { 
 //cout<<"entering SiPixelInformationExtractor::goToDir"<<endl;
   bei->cd();
   //if(flg) bei->cd("Collector/Collated");
   bei->cd(sname);
   string dirName = bei->pwd();
   if (dirName.find(sname) != string::npos) return true;
   else return false;  
 //cout<<"leaving SiPixelInformationExtractor::goToDir"<<endl;
 }

bool SiPixelInformationExtractor::hasItem	(	std::multimap< std::string, std::string > &	req_map,
		std::string	item_name
	)

private

(Documentation under construction).

This method

Definition at line 704 of file SiPixelInformationExtractor.cc.

                                        {
 //cout<<"entering SiPixelInformationExtractor::hasItem"<<endl;
   multimap<string,string>::iterator pos = req_map.find(item_name);
   if (pos != req_map.end()) return true;
   return false;  
 //cout<<"leaving SiPixelInformationExtractor::hasItem"<<endl;
 }

void SiPixelInformationExtractor::printAlarmList	(	DQMStore *	bei,
		std::ostringstream &	str_val
	)

private

(Documentation under construction).

This method

(Documentation under construction).

Returns a stringstream containing an HTML-formatted list of alarms for the current directory. This is a recursive method.

Definition at line 601 of file SiPixelInformationExtractor.cc.

References alarmCounter_, DQMStore::cd(), RefreshWebPage::dname, DQMStore::get(), DQMStore::getMEs(), MonitorElement::getQReports(), DQMStore::getStatus(), DQMStore::getSubdirs(), DQMStore::goUp(), DQMStore::pwd(), and selectImage().

                                                                         {
 //cout<<"entering SiPixelInformationExtractor::printAlarmList"<<endl;
 //   cout << ACRed << ACBold
 //        << "[SiPixelInformationExtractor::printAlarmList()]"
 //        << ACPlain
 //        << " Enter" 
 //        << endl ;
   string currDir = bei->pwd();
   string dname = currDir.substr(currDir.find_last_of("/")+1);
   string image_name;
   selectImage(image_name,bei->getStatus(currDir));
   if(image_name!="images/LI_green.gif")
     str_val << " <li>\n"
             << "  <a href=\"#\" id=\"" << currDir << "\">\n   " 
         <<     dname 
         << "  </a>\n"
         << "  <img src=\"" 
             <<     image_name 
         << "\">" << endl;
   vector<string> subDirVec = bei->getSubdirs();
 
   vector<string> meVec = bei->getMEs();
    
   if (subDirVec.size() == 0 && meVec.size() == 0) {
     str_val <<  "</li> "<< endl;    
     return;
   }
   str_val << "<ul>" << endl;
   for (vector<string>::const_iterator it = meVec.begin();
        it != meVec.end(); it++) {
     string full_path = currDir + "/" + (*it);
 
     MonitorElement * me = bei->get(full_path);
     
     if (!me) continue;
     std::vector<QReport *> my_map = me->getQReports();
     if (my_map.size() > 0) {
       string image_name1;
       selectImage(image_name1,my_map);
       if(image_name1!="images/LI_green.gif") {
         alarmCounter_++;
         str_val << "    <li class=\"dhtmlgoodies_sheet.gif\">\n"
             << "     <input id  = \"selectedME\""
             << "        folder  = \"" << currDir << "\""
             << "        type    = \"checkbox\""
             << "        name    = \"selected\""
             << "        class   = \"smallCheckBox\""
             << "        value   = \"" << (*it) << "\""
             << "        onclick = \"javascript:IMGC.selectedIMGCItems()\" />\n"
 //          << "     <a href=\"javascript:IMGC.updateIMGC('" << currDir << "')\">\n       " 
             << "     <a href=\"javascript:IMGC.plotFromPath('" << currDir << "')\">\n       " 
             <<    (*it) << "\n"
             << "     </a>\n"
         << "     <img src=\""
         <<        image_name1 
         << "\">"
             << "    </li>" 
             << endl;
     }   
     }
   }
   for (vector<string>::const_iterator ic = subDirVec.begin();
        ic != subDirVec.end(); ic++) {
     bei->cd(*ic);
     printAlarmList(bei, str_val);
     bei->goUp();
   }
   str_val << "</ul> "<< endl;  
   str_val << "</li> "<< endl;  
 //   cout << ACGreen << ACBold
 //        << "[SiPixelInformationExtractor::printAlarmList()]"
 //        << ACPlain
 //        << " Done" 
 //        << endl ;
 //cout<<"leaving SiPixelInformationExtractor::printAlarmList"<<endl;
 }

void SiPixelInformationExtractor::printModuleHistoList	(	DQMStore *	bei,
		std::ostringstream &	str_val
	)

private

(Documentation under construction).

This method

(Documentation under construction).

This method

Definition at line 407 of file SiPixelInformationExtractor.cc.

References DQMStore::cd(), RefreshWebPage::dname, DQMStore::getMEs(), DQMStore::getSubdirs(), DQMStore::goUp(), and DQMStore::pwd().

                                                                               {
 //cout<<"entering SiPixelInformationExtractor::printModuleHistoList"<<endl;
   string currDir = bei->pwd();
   string dname = currDir.substr(currDir.find_last_of("/")+1);
   str_val << " <li>\n"
       << "  <a href=\"#\" id=\"" << currDir << "\">\n   " 
       <<     dname << "\n"
       << "  </a>\n"
       << endl << endl;
 
   vector<string> meVec     = bei->getMEs(); 
   
   vector<string> subDirVec = bei->getSubdirs();
   if ( meVec.size()== 0  && subDirVec.size() == 0 ) {
     str_val << " </li>" << endl;    
     return;
   }
   str_val << "\n   <ul>" << endl; 
   for (vector<string>::const_iterator it  = meVec.begin();
                                       it != meVec.end(); it++) {
     if ((*it).find("_siPixel")!=string::npos || 
         (*it).find("_ctfWithMaterialTracks")!=string::npos) {
       string qit = (*it) ;
       string temp_s;
       if(qit.find("siPixel")!=string::npos || qit.find("ctfWithMaterialTracks")!=string::npos) { temp_s = qit.substr(0,qit.find_first_of("_")); }
       str_val << "    <li class=\"dhtmlgoodies_sheet.gif\">\n"
           << "     <input id      = \"selectedME\""
           << "            folder  = \"" << currDir << "\""
           << "            type    = \"checkbox\""
           << "            name    = \"selected\""
           << "            class   = \"smallCheckBox\""
           << "            value   = \"" << (*it) << "\""
           << "            onclick = \"javascript:IMGC.selectedIMGCItems()\" />\n"
 //        << "     <a href=\"javascript:IMGC.updateIMGC('" << currDir << "')\">\n       " 
           << "     <a href=\"javascript:IMGC.plotFromPath('" << currDir << "')\">\n       " 
 //        <<        temp_s << "\n"
           <<        (*it) << "\n"
           << "     </a>\n"
           << "    </li>" 
           << endl;
     }
   }
   for (vector<string>::const_iterator ic  = subDirVec.begin();
                                       ic != subDirVec.end(); ic++) {
     bei->cd(*ic);
     printModuleHistoList(bei, str_val);
     bei->goUp();
   }
   str_val << "   </ul>" << endl;  
   str_val << "  </li>"  << endl;  
 //cout<<"leaving SiPixelInformationExtractor::printModuleHistoList"<<endl;
 }

void SiPixelInformationExtractor::printSummaryHistoList	(	DQMStore *	bei,
		std::ostringstream &	str_val
	)

private

(Documentation under construction).

This method

(Documentation under construction).

Returns a stringstream containing an HTML-formatted list of ME in the current directory. This is a recursive method.

Definition at line 502 of file SiPixelInformationExtractor.cc.

References DQMStore::cd(), RefreshWebPage::dname, DQMStore::getMEs(), DQMStore::getSubdirs(), DQMStore::goUp(), and DQMStore::pwd().

                                                                                {
 //cout<<"entering SiPixelInformationExtractor::printSummaryHistoList"<<endl;
   string currDir = bei->pwd();
   string dname = currDir.substr(currDir.find_last_of("/")+1);
   if (dname.find("Module_") ==0 || dname.find("FED_")==0) return;
   str_val << " <li>\n"
           << "  <a href=\"#\" id=\"" << currDir << "\">\n   " 
       <<     dname 
       << "  </a>" 
       << endl;
 
   vector<string> meVec     = bei->getMEs(); 
   
   vector<string> subDirVec = bei->getSubdirs();
   if ( meVec.size()== 0  && subDirVec.size() == 0 ) {
     str_val << " </li> "<< endl;    
     return;
   }
   str_val << "\n   <ul>" << endl;      
   for (vector<string>::const_iterator it = meVec.begin();
        it != meVec.end(); it++) {
     if ((*it).find("SUM") == 0) {
       str_val << "    <li class=\"dhtmlgoodies_sheet.gif\">\n"
           << "     <input id      = \"selectedME\""
           << "            folder  = \"" << currDir << "\""
           << "            type    = \"checkbox\""
           << "            name    = \"selected\""
           << "            class   = \"smallCheckBox\""
           << "            value   = \"" << (*it) << "\""
           << "            onclick = \"javascript:IMGC.selectedIMGCItems()\" />\n"
 //              << "     <a href=\"javascript:IMGC.updateIMGC('" << currDir << "')\">\n       " 
               << "     <a href=\"javascript:IMGC.plotFromPath('" << currDir << "')\">\n       " 
           <<       (*it) << "\n"
           << "     </a>\n"
           << "    </li>" 
           << endl;
     }
   }
 
   for (vector<string>::const_iterator ic = subDirVec.begin();
        ic != subDirVec.end(); ic++) {
     bei->cd(*ic);
     printSummaryHistoList(bei, str_val);
     bei->goUp();
   }
   str_val << "   </ul> "<< endl;  
   str_val << "  </li> "<< endl;  
 //cout<<"leaving SiPixelInformationExtractor::printSummaryHistoList"<<endl;
 }

void SiPixelInformationExtractor::readConfiguration ( )

Read Configuration File.

Definition at line 91 of file SiPixelInformationExtractor.cc.

91 { }

bool SiPixelInformationExtractor::readConfiguration	(	std::map< std::string, std::vector< std::string > > &	layoutMap,
		std::map< std::string, std::map< std::string, std::string > > &	qtestsMap,
		std::map< std::string, std::vector< std::string > > &	meQTestsMap
	)

void SiPixelInformationExtractor::selectColor	(	std::string &	col,
		int	status
	)

private

void SiPixelInformationExtractor::selectColor	(	std::string &	col,
		std::vector< QReport * > &	reports
	)

private

void SiPixelInformationExtractor::selectImage	(	std::string &	name,
		int	status
	)

private

Referenced by printAlarmList().

void SiPixelInformationExtractor::selectImage	(	std::string &	name,
		std::vector< QReport * > &	reports
	)

private

void SiPixelInformationExtractor::selectMEList	(	DQMStore *	bei,
		std::string &	name,
		std::vector< MonitorElement * > &	mes
	)

(Documentation under construction).

Definition at line 1036 of file SiPixelInformationExtractor.cc.

References DQMStore::cd(), relmon_rootfiles_spy::contents, DQMStore::get(), DQMStore::getMEs(), DQMStore::getSubdirs(), DQMStore::goUp(), and DQMStore::pwd().

 {  
 //  cout<<"In SiPixelInformationExtractor::selectMEList: "<<endl;
   string currDir = bei->pwd();
    
   string theME ;
    
   // Get ME from Collector/FU0/Tracker/PixelEndcap/HalfCylinder_pX/Disk_X/Blade_XX/Panel_XX/Module_XX
   if (currDir.find("Module_") != string::npos ||
       currDir.find("FED_") != string::npos)  
   {
     vector<string> contents = bei->getMEs(); 
        
     for (vector<string>::const_iterator it = contents.begin(); it != contents.end(); it++) 
     {
       theME = (*it) ;
       if(theME.find("siPixel")==string::npos && theME.find("ctfWithMaterialTracks")==string::npos) {continue ;} // If the ME is not a siPixel or ctfWithMaterialTrack one, skip
       string temp_s = theME.substr(0,theME.find_first_of("_"));
       //cout<<"should be the variable name: temp_s= "<<temp_s<<endl;
       if (temp_s == theMEName)  
       {
         string full_path = currDir + "/" + (*it);
 
         MonitorElement * me = bei->get(full_path.c_str());
     
         if (me) {mes.push_back(me);}
       }
     }
     return;
   } else {  // If not yet reached the desired level in the directory tree, recursively go down one level more
     vector<string> subdirs = bei->getSubdirs();
     for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) 
     {
       bei->cd(*it);
       selectMEList(bei, theMEName, mes);
       bei->goUp();
     }
   }
 }