ADDITON OF TK_HISTO_MAP. More...

#include <SiStripBaseCondObjDQM.h>

Inheritance diagram for SiStripBaseCondObjDQM:

Classes
struct	ModMEs

Public Member Functions
void	analysis (const edm::EventSetup &eSetup_)

void	analysisOnDemand (const edm::EventSetup &eSetup_, uint32_t detIdOnDemand)

void	analysisOnDemand (const edm::EventSetup &eSetup_, std::vector< uint32_t > detIdsOnDemand)

void	analysisOnDemand (const edm::EventSetup &eSetup_, std::string requestedSubDetector, uint32_t requestedSide, uint32_t requestedLayer)

virtual void	end ()

virtual void	getActiveDetIds (const edm::EventSetup &eSetup)=0

std::vector< uint32_t >	getCabledModules ()

virtual unsigned long long	getCache (const edm::EventSetup &eSetup_)=0

virtual void	getConditionObject (const edm::EventSetup &eSetup_)=0

void	selectModules (std::vector< uint32_t > &detIds_)

	SiStripBaseCondObjDQM (const edm::EventSetup &eSetup, edm::ParameterSet const &hPSet, edm::ParameterSet const &fPSet)

virtual	~SiStripBaseCondObjDQM ()

Protected Member Functions
virtual void	fillMEsForDet (ModMEs selModME_, uint32_t selDetId_)=0

virtual void	fillMEsForLayer (std::map< uint32_t, ModMEs > selModMEsMap_, uint32_t selDetId_)=0

virtual void	fillModMEs (const std::vector< uint32_t > &selectedDetIds)

virtual void	fillSummaryMEs (const std::vector< uint32_t > &selectedDetIds)

void	fillTkMap (const uint32_t &detid, const float &value)

std::pair< std::string, uint32_t >	getLayerNameAndId (const uint32_t &detId_)

void	getModMEs (ModMEs &CondObj_ME, const uint32_t &detId_)

std::vector< uint32_t >	GetSameLayerDetId (std::vector< uint32_t > activeDetIds, uint32_t selDetId)

std::pair< std::string, uint32_t >	getStringNameAndId (const uint32_t &detId_)

void	getSummaryMEs (ModMEs &CondObj_ME, const uint32_t &detId_)

Protected Attributes
std::vector< uint32_t >	activeDetIds

std::vector< uint32_t >	all_DetIds

unsigned long long	cacheID_current

unsigned long long	cacheID_memory

std::string	CondObj_fillId_

std::string	CondObj_name_

const edm::EventSetup &	eSetup_

edm::ParameterSet	fPSet_

bool	GrandSummary_On_

bool	HistoMaps_On_

edm::ParameterSet	hPSet_

double	maxValue

double	minValue

bool	Mod_On_

std::map< uint32_t, ModMEs >	ModMEsMap_

SiStripDetInfoFileReader *	reader

std::map< uint32_t, ModMEs >	SummaryMEsMap_

bool	SummaryOnLayerLevel_On_

bool	SummaryOnStringLevel_On_

TkHistoMap *	Tk_HM_

TkHistoMap *	Tk_HM_H

TkHistoMap *	Tk_HM_L

TrackerMap *	tkMap

std::vector< int >	tkMapScaler

Private Member Functions
void	bookCumulMEs (SiStripBaseCondObjDQM::ModMEs &CondObj_ME, const uint32_t &detId_)

void	bookProfileMEs (SiStripBaseCondObjDQM::ModMEs &CondObj_ME, const uint32_t &detId_)

void	bookSummaryCumulMEs (SiStripBaseCondObjDQM::ModMEs &CondObj_ME, const uint32_t &detId_)

void	bookSummaryMEs (SiStripBaseCondObjDQM::ModMEs &CondObj_ME, const uint32_t &detId_)

void	bookSummaryProfileMEs (SiStripBaseCondObjDQM::ModMEs &CondObj_ME, const uint32_t &detId_)

void	bookTkMap (const std::string &TkMapname)

void	saveTkMap (const std::string &TkMapname, double minValue, double maxValue)

Private Attributes
std::string	condDataMonitoringMode_

edm::ESHandle< SiStripDetCabling >	detCablingHandle_

DQMStore *	dqmStore_

SiStripFolderOrganizer	folder_organizer

SiStripHistoId	hidmanager

std::vector< uint32_t >	ModulesToBeExcluded_

std::vector< uint32_t >	ModulesToBeIncluded_

std::vector< std::string >	SubDetectorsToBeExcluded_

Detailed Description

ADDITON OF TK_HISTO_MAP.

Definition at line 41 of file SiStripBaseCondObjDQM.h.

Constructor & Destructor Documentation

SiStripBaseCondObjDQM::SiStripBaseCondObjDQM	(	const edm::EventSetup &	eSetup,
		edm::ParameterSet const &	hPSet,
		edm::ParameterSet const &	fPSet
	)

Definition at line 9 of file SiStripBaseCondObjDQM.cc.

References bookTkMap(), CondObj_fillId_, CondObj_name_, fPSet_, edm::ParameterSet::getParameter(), GrandSummary_On_, HistoMaps_On_, hPSet_, maxValue, minValue, Mod_On_, reader, python.OfflineSiStripConfigDb_cff::SiStripDetInfoFileReader, SummaryOnLayerLevel_On_, and SummaryOnStringLevel_On_.

                                                                            :
   eSetup_(eSetup),
   hPSet_(hPSet),
   fPSet_(fPSet),
   cacheID_memory(0),
   dqmStore_(edm::Service<DQMStore>().operator->()){
 
   reader = new SiStripDetInfoFileReader(edm::FileInPath(std::string("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat") ).fullPath());
   
   Mod_On_                  = fPSet_.getParameter<bool>("Mod_On");
   HistoMaps_On_            = fPSet_.getParameter<bool>("HistoMaps_On");
   SummaryOnLayerLevel_On_  = fPSet_.getParameter<bool>("SummaryOnLayerLevel_On");
   SummaryOnStringLevel_On_ = fPSet_.getParameter<bool>("SummaryOnStringLevel_On");
 
   GrandSummary_On_         = fPSet_.getParameter<bool>("GrandSummary_On");
 
   CondObj_fillId_    = hPSet_.getParameter<std::string>("CondObj_fillId");
   CondObj_name_      = hPSet_.getParameter<std::string>("CondObj_name");
  
 
   //Warning message from wrong input:
   if(SummaryOnLayerLevel_On_ && SummaryOnStringLevel_On_){
     edm::LogWarning("SiStripBaseCondObjDQM") 
        << "[SiStripBaseCondObjDQM::SiStripBaseCondObjDQMs] PLEASE CHECK : String and layer level options can not be activated together"
        << std::endl; 
   }
 
   //The OR of the two conditions allow to switch on this feature for all the components (if the FillConditions_PSet has the TkMap_On =true) or for single MEs (if the PSet for a ME has the TkMap_On =true)
   if(fPSet_.getParameter<bool>("TkMap_On") || hPSet_.getParameter<bool>("TkMap_On")) bookTkMap(hPSet_.getParameter<std::string>("TkMapName"));
 
 
   minValue=hPSet_.getParameter<double>("minValue");
   maxValue=hPSet_.getParameter<double>("maxValue");
 
 }

virtual SiStripBaseCondObjDQM::~SiStripBaseCondObjDQM ( )

inlinevirtual

Definition at line 53 of file SiStripBaseCondObjDQM.h.

53 {};

Member Function Documentation

void SiStripBaseCondObjDQM::analysis ( const edm::EventSetup & eSetup_ )

Definition at line 51 of file SiStripBaseCondObjDQM.cc.

References activeDetIds, cacheID_current, cacheID_memory, edm::IOVSyncValue::eventID(), lut2db_cfg::filename, fillModMEs(), fillSummaryMEs(), fPSet_, getActiveDetIds(), SiStripDetInfoFileReader::getAllDetIds(), getCache(), getConditionObject(), edm::ParameterSet::getParameter(), hPSet_, edm::EventSetup::iovSyncValue(), maxValue, minValue, Mod_On_, reader, edm::EventID::run(), saveTkMap(), selectModules(), SummaryOnLayerLevel_On_, and SummaryOnStringLevel_On_.

Referenced by SiStripClassToMonitorCondData::analyseCondData(), and SiStripMonitorCondData::analyze().

                                                                  {
  
   cacheID_current=  getCache(eSetup_);
   
   if (cacheID_memory == cacheID_current) return;
   
   getConditionObject(eSetup_);
 
   //The OR of the two conditions allows to switch on this feature for all the components (if the FillConditions_PSet has the ActiveDetIds_On =true) or for single MEs (if the PSet for a ME has the ActiveDetIds_On =true)
   if(fPSet_.getParameter<bool>("ActiveDetIds_On") || hPSet_.getParameter<bool>("ActiveDetIds_On"))
     getActiveDetIds(eSetup_);
   else
     activeDetIds=reader->getAllDetIds();
 
   selectModules(activeDetIds);
 
   if(Mod_On_ )                                            { fillModMEs(activeDetIds); }
   if(SummaryOnLayerLevel_On_ || SummaryOnStringLevel_On_ ){ fillSummaryMEs(activeDetIds);}
 
   std::string filename = hPSet_.getParameter<std::string>("TkMapName");
   if (filename!=""){
     char sRun[128];
     sprintf(sRun,"_Run_%d",eSetup_.iovSyncValue().eventID().run());
     filename.insert(filename.find("."),sRun);
     
     saveTkMap(filename.c_str(), minValue, maxValue);
   }
 }

void SiStripBaseCondObjDQM::analysisOnDemand	(	const edm::EventSetup &	eSetup_,
		uint32_t	detIdOnDemand
	)

Definition at line 117 of file SiStripBaseCondObjDQM.cc.

References cacheID_current, cacheID_memory, fillModMEs(), getCache(), and getConditionObject().

Referenced by analysisOnDemand(), SiStripClassToMonitorCondData::getLayerMEsOnDemand(), and SiStripClassToMonitorCondData::getModMEsOnDemand().

                                                                                                   {
  
   unsigned long long cacheID_current=  getCache(eSetup_);
   
   if (cacheID_memory == cacheID_current) return;
   
   getConditionObject(eSetup_);
   
   std::vector<uint32_t> vdetIdsOnDemand_;
   vdetIdsOnDemand_.push_back(detIdOnDemand); // fillModMEs needs a vector 
 
   fillModMEs(vdetIdsOnDemand_); 
   
 }

void SiStripBaseCondObjDQM::analysisOnDemand	(	const edm::EventSetup &	eSetup_,
		std::vector< uint32_t >	detIdsOnDemand
	)

Definition at line 134 of file SiStripBaseCondObjDQM.cc.

References cacheID_current, cacheID_memory, fillSummaryMEs(), getCache(), and getConditionObject().

                                                                                                               {
  
   unsigned long long cacheID_current=  getCache(eSetup_);
   
   if (cacheID_memory == cacheID_current) return;
   
   getConditionObject(eSetup_);
   
   fillSummaryMEs(detIdsOnDemand); 
   
 }

void SiStripBaseCondObjDQM::analysisOnDemand	(	const edm::EventSetup &	eSetup_,
		std::string	requestedSubDetector,
		uint32_t	requestedSide,
		uint32_t	requestedLayer
	)

Definition at line 84 of file SiStripBaseCondObjDQM.cc.

References activeDetIds, analysisOnDemand(), getActiveDetIds(), getConditionObject(), SiStripSubStructure::getTECDetectors(), SiStripSubStructure::getTIBDetectors(), SiStripSubStructure::getTIDDetectors(), and SiStripSubStructure::getTOBDetectors().

                                                 {
   
   getConditionObject(eSetup_);
   getActiveDetIds(eSetup_);
    
   std::vector<uint32_t> requestedDetIds_;
   requestedDetIds_.clear();
   
   SiStripSubStructure substructure_;
   
   if(requestedSubDetector=="TIB"){ 
       substructure_.getTIBDetectors( activeDetIds, requestedDetIds_, requestedLayer,0,0,0);
   }
   else if(requestedSubDetector=="TID"){ 
       substructure_.getTIDDetectors( activeDetIds, requestedDetIds_, requestedSide,requestedLayer,0,0);
   }
   else if(requestedSubDetector=="TOB"){  
       substructure_.getTOBDetectors( activeDetIds, requestedDetIds_, requestedLayer,0,0);
   }
   else if(requestedSubDetector=="TEC"){  
       substructure_.getTECDetectors( activeDetIds, requestedDetIds_, requestedSide,requestedLayer,0,0,0,0);
   } 
 
   analysisOnDemand(eSetup_,requestedDetIds_);
  
 }

void SiStripBaseCondObjDQM::bookCumulMEs	(	SiStripBaseCondObjDQM::ModMEs &	CondObj_ME,
		const uint32_t &	detId_
	)

private

Definition at line 461 of file SiStripBaseCondObjDQM.cc.

References DQMStore::book1D(), SiStripHistoId::createHistoId(), SiStripBaseCondObjDQM::ModMEs::CumulDistr, dqmStore_, folder_organizer, edm::ParameterSet::getParameter(), hidmanager, hPSet_, MonitorElement::setAxisTitle(), SiStripFolderOrganizer::setDetectorFolder(), and DQMStore::tag().

Referenced by getModMEs().

                                                                                                        {
 
   int    hCumul_NchX    = 0;
   double hCumul_LowX    = 0;
   double hCumul_HighX   = 0;
      
   std::string hCumul_description;
   hCumul_description   = hPSet_.getParameter<std::string>("Cumul_description");
       
   std::string hCumul_xTitle, hCumul_yTitle;
   hCumul_xTitle        = hPSet_.getParameter<std::string>("Cumul_xTitle");
   hCumul_yTitle        = hPSet_.getParameter<std::string>("Cumul_yTitle");
       
   hCumul_NchX          = hPSet_.getParameter<int>("Cumul_NchX");
   hCumul_LowX          = hPSet_.getParameter<double>("Cumul_LowX");
   hCumul_HighX         = hPSet_.getParameter<double>("Cumul_HighX");
       
   folder_organizer.setDetectorFolder(detId_); 
       
   std::string hCumul_name; 
   hCumul_name   = hidmanager.createHistoId(hCumul_description  , "det", detId_); ;
       
   std::string hCumul_title;
   hCumul_title   = hCumul_name ;
       
   CondObj_ME.CumulDistr = dqmStore_->book1D(hCumul_name, 
                                             hCumul_title, 
                         hCumul_NchX, 
                         hCumul_LowX, 
                         hCumul_HighX);
   CondObj_ME.CumulDistr->setAxisTitle(hCumul_xTitle,1);
   CondObj_ME.CumulDistr->setAxisTitle(hCumul_yTitle,2);
   dqmStore_->tag(CondObj_ME.CumulDistr, detId_);
       
 } 

void SiStripBaseCondObjDQM::bookProfileMEs	(	SiStripBaseCondObjDQM::ModMEs &	CondObj_ME,
		const uint32_t &	detId_
	)

private

Definition at line 412 of file SiStripBaseCondObjDQM.cc.

References DQMStore::book1D(), CondObj_name_, SiStripHistoId::createHistoId(), dqmStore_, folder_organizer, SiStripDetInfoFileReader::getNumberOfApvsAndStripLength(), edm::ParameterSet::getParameter(), hidmanager, hPSet_, SiStripBaseCondObjDQM::ModMEs::ProfileDistr, reader, MonitorElement::setAxisTitle(), SiStripFolderOrganizer::setDetectorFolder(), and DQMStore::tag().

Referenced by getModMEs().

                                                                                                          {
      
   int   hProfile_NchX  = 0;
   double hProfile_LowX  = 0;
   double hProfile_HighX = 0;
   
   std::string hProfile_description;
   hProfile_description   = hPSet_.getParameter<std::string>("Profile_description");
       
   std::string hProfile_xTitle, hProfile_yTitle;
   hProfile_xTitle          = hPSet_.getParameter<std::string>("Profile_xTitle");
   hProfile_yTitle          = hPSet_.getParameter<std::string>("Profile_yTitle");
             
   if( CondObj_name_!= "apvgain" ){
     
     int nStrip      = reader->getNumberOfApvsAndStripLength(detId_).first*128;
     
     hProfile_NchX           = nStrip;
     hProfile_LowX           = 0.5;
     hProfile_HighX          = nStrip+0.5;
   }
   else {
     
     int nApv      = reader->getNumberOfApvsAndStripLength(detId_).first;
     
     hProfile_NchX           = nApv;
     hProfile_LowX           = 0.5;
     hProfile_HighX          = nApv+0.5;
   }
       
   folder_organizer.setDetectorFolder(detId_); 
       
   std::string hProfile_Name; 
   hProfile_Name = hidmanager.createHistoId(hProfile_description, "det", detId_);
       
   std::string hProfile;
   hProfile = hProfile_Name ;
     
   CondObj_ME.ProfileDistr = dqmStore_->book1D(hProfile_Name, hProfile, hProfile_NchX, hProfile_LowX, hProfile_HighX);
   CondObj_ME.ProfileDistr->setAxisTitle(hProfile_xTitle,1);
   CondObj_ME.ProfileDistr->setAxisTitle(hProfile_yTitle,2);
   dqmStore_->tag(CondObj_ME.ProfileDistr, detId_);
   
 }

void SiStripBaseCondObjDQM::bookSummaryCumulMEs	(	SiStripBaseCondObjDQM::ModMEs &	CondObj_ME,
		const uint32_t &	detId_
	)

private

Definition at line 727 of file SiStripBaseCondObjDQM.cc.

References DQMStore::book1D(), CondObj_name_, SiStripHistoId::createHistoLayer(), dqmStore_, first, folder_organizer, getLayerNameAndId(), edm::ParameterSet::getParameter(), getStringNameAndId(), SiStripFolderOrganizer::GetSubDetAndLayer(), hidmanager, hPSet_, sistrip::layer_, MonitorElement::setAxisTitle(), SiStripFolderOrganizer::setLayerFolder(), SiStripBaseCondObjDQM::ModMEs::SummaryOfCumulDistr, SummaryOnStringLevel_On_, and DQMStore::tag().

Referenced by getSummaryMEs().

                                                                                                               {
     
   int   hSummaryOfCumul_NchX    = 0;
   double hSummaryOfCumul_LowX    = 0;
   double hSummaryOfCumul_HighX   = 0;
     
   std::string hSummaryOfCumul_description;
   hSummaryOfCumul_description  = hPSet_.getParameter<std::string>("SummaryOfCumul_description");
     
   std::string hSummaryOfCumul_xTitle, hSummaryOfCumul_yTitle;
   hSummaryOfCumul_xTitle        = hPSet_.getParameter<std::string>("SummaryOfCumul_xTitle");
   hSummaryOfCumul_yTitle        = hPSet_.getParameter<std::string>("SummaryOfCumul_yTitle");
     
   hSummaryOfCumul_NchX          = hPSet_.getParameter<int>("SummaryOfCumul_NchX");
   hSummaryOfCumul_LowX          = hPSet_.getParameter<double>("SummaryOfCumul_LowX");
   hSummaryOfCumul_HighX         = hPSet_.getParameter<double>("SummaryOfCumul_HighX");
     
   uint32_t layer_=0;
     
   layer_ = folder_organizer.GetSubDetAndLayer(detId_).second;
     
   folder_organizer.setLayerFolder(detId_,layer_); 
     
   std::string hSummaryOfCumul_name; 
   
   // ---
   int subdetectorId_ = ((detId_>>25)&0x7);
   
   if( subdetectorId_<3 || subdetectorId_>6 ){ 
     edm::LogError("SiStripBaseCondObjDQM")
        << "[SiStripBaseCondObjDQM::bookSummaryCumulMEs] WRONG INPUT : no such subdetector type : "
        << subdetectorId_ << " no folder set!" 
        << std::endl;
     return;
   }
   // ---
   
   // LA Histos are plotted for each string:
   if(CondObj_name_ == "lorentzangle" && SummaryOnStringLevel_On_) { 
     hSummaryOfCumul_name = hidmanager.createHistoLayer(hSummaryOfCumul_description, "layer" , getStringNameAndId(detId_).first, "") ;
   }
   else {  
     hSummaryOfCumul_name = hidmanager.createHistoLayer(hSummaryOfCumul_description, "layer" , getLayerNameAndId(detId_).first, "") ;
   }
 
     
   std::string hSummaryOfCumul_title;
   hSummaryOfCumul_title   = hSummaryOfCumul_name ;
     
   CondObj_ME.SummaryOfCumulDistr = dqmStore_->book1D(hSummaryOfCumul_name, 
                                                      hSummaryOfCumul_title, 
                              hSummaryOfCumul_NchX, 
                              hSummaryOfCumul_LowX, 
                              hSummaryOfCumul_HighX);
 
   CondObj_ME.SummaryOfCumulDistr->setAxisTitle(hSummaryOfCumul_xTitle,1);
   CondObj_ME.SummaryOfCumulDistr->setAxisTitle(hSummaryOfCumul_yTitle,2);
     
   dqmStore_->tag(CondObj_ME.SummaryOfCumulDistr, layer_);
     
 }

void SiStripBaseCondObjDQM::bookSummaryMEs	(	SiStripBaseCondObjDQM::ModMEs &	CondObj_ME,
		const uint32_t &	detId_
	)

private

Definition at line 793 of file SiStripBaseCondObjDQM.cc.

References activeDetIds, DQMStore::bookProfile(), SiStripHistoId::createHistoLayer(), dqmStore_, first, folder_organizer, getLayerNameAndId(), edm::ParameterSet::getParameter(), GetSameLayerDetId(), SiStripFolderOrganizer::GetSubDetAndLayer(), hidmanager, hPSet_, i, sistrip::layer_, MonitorElement::setAxisRange(), MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), SiStripFolderOrganizer::setLayerFolder(), SiStripBaseCondObjDQM::ModMEs::SummaryDistr, and DQMStore::tag().

Referenced by getSummaryMEs().

                                                                                                          {
   
   std::vector<uint32_t> sameLayerDetIds_;
 
   int   hSummary_NchX    = 0;
   double hSummary_LowX    = 0;
   double hSummary_HighX   = 0;
      
   std::string hSummary_description;
   hSummary_description  = hPSet_.getParameter<std::string>("Summary_description");
       
   std::string hSummary_xTitle, hSummary_yTitle;
   hSummary_xTitle        = hPSet_.getParameter<std::string>("Summary_xTitle");
   hSummary_yTitle        = hPSet_.getParameter<std::string>("Summary_yTitle");
   
   int hSummary_NchY;
   double hSummary_LowY, hSummary_HighY;
   hSummary_NchY          = hPSet_.getParameter<int>("Summary_NchY");
   hSummary_LowY          = hPSet_.getParameter<double>("Summary_LowY");
   hSummary_HighY         = hPSet_.getParameter<double>("Summary_HighY");
        
 
   // -----
   // get detIds belonging to same layer to fill X-axis with detId-number
                        
   
   sameLayerDetIds_.clear();
    
   sameLayerDetIds_=GetSameLayerDetId(activeDetIds,detId_);
 
   hSummary_NchX           = sameLayerDetIds_.size(); 
   hSummary_LowX           = 0.5;
   hSummary_HighX          = sameLayerDetIds_.size()+0.5;
   
   uint32_t layer_=0;
       
   layer_ = folder_organizer.GetSubDetAndLayer(detId_).second;
       
   folder_organizer.setLayerFolder(detId_,layer_); 
       
   std::string hSummary_name; 
   
   // ---
   int subdetectorId_ = ((detId_>>25)&0x7);
   
  
   if( subdetectorId_<3 ||subdetectorId_>6 ){ 
     edm::LogError("SiStripBaseCondObjDQM")
        << "[SiStripBaseCondObjDQM::bookSummaryMEs] WRONG INPUT : no such subdetector type : "
        << subdetectorId_ << " no folder set!" 
        << std::endl;
     return;
   }
   // ---
   
   hSummary_name = hidmanager.createHistoLayer(hSummary_description, 
                                                        "layer" , 
                                 getLayerNameAndId(detId_).first,
                             "") ;
       
   std::string hSummary_title;
   hSummary_title   = hSummary_name ;
            
   CondObj_ME.SummaryDistr = dqmStore_->bookProfile(hSummary_name, 
                                                    hSummary_title, 
                            hSummary_NchX, 
                            hSummary_LowX, 
                            hSummary_HighX, 
                            hSummary_NchY, 
                            hSummary_LowY, 
                            hSummary_HighY);
   CondObj_ME.SummaryDistr->setAxisTitle(hSummary_xTitle,1);
   CondObj_ME.SummaryDistr->setAxisTitle(hSummary_yTitle,2);
   CondObj_ME.SummaryDistr->setAxisRange(hSummary_LowY, hSummary_HighY,2);
  
   // -----
   // in order to get the right detId-number labelled in right bin of x-axis
   unsigned int iBin=0;
     
   for(unsigned int i=0;i< sameLayerDetIds_.size(); i++){
     
     iBin++;
     char sameLayerDetIds_Name[1024];
     sprintf(sameLayerDetIds_Name,"%u",sameLayerDetIds_[i]);
     if(iBin%100==0)
       CondObj_ME.SummaryDistr->setBinLabel(iBin, sameLayerDetIds_Name);
     
   }
   // -----
       
   dqmStore_->tag(CondObj_ME.SummaryDistr, layer_);
       
 } 

void SiStripBaseCondObjDQM::bookSummaryProfileMEs	(	SiStripBaseCondObjDQM::ModMEs &	CondObj_ME,
		const uint32_t &	detId_
	)

private

Definition at line 502 of file SiStripBaseCondObjDQM.cc.

References activeDetIds, DQMStore::bookProfile(), CondObj_name_, SiStripHistoId::createHistoLayer(), dqmStore_, first, folder_organizer, getLayerNameAndId(), SiStripDetInfoFileReader::getNumberOfApvsAndStripLength(), edm::ParameterSet::getParameter(), getStringNameAndId(), SiStripFolderOrganizer::GetSubDetAndLayer(), SiStripSubStructure::getTECDetectors(), SiStripSubStructure::getTIBDetectors(), SiStripSubStructure::getTIDDetectors(), SiStripSubStructure::getTOBDetectors(), hidmanager, hPSet_, i, sistrip::layer_, align::tib::layerNumber(), module(), reader, align::tob::rodNumber(), MonitorElement::setAxisRange(), MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), SiStripFolderOrganizer::setLayerFolder(), align::tib::stringNumber(), SiStripBaseCondObjDQM::ModMEs::SummaryOfProfileDistr, SummaryOnLayerLevel_On_, SummaryOnStringLevel_On_, and DQMStore::tag().

Referenced by getSummaryMEs().

                                                                                                                 {
   
   std::vector<uint32_t> sameLayerDetIds_;
 
   int   hSummaryOfProfile_NchX    = 0;
   double hSummaryOfProfile_LowX    = 0;
   double hSummaryOfProfile_HighX   = 0;
      
   std::string hSummaryOfProfile_description;
   hSummaryOfProfile_description  = hPSet_.getParameter<std::string>("SummaryOfProfile_description");
       
   std::string hSummaryOfProfile_xTitle, hSummaryOfProfile_yTitle;
   hSummaryOfProfile_xTitle        = hPSet_.getParameter<std::string>("SummaryOfProfile_xTitle");
   hSummaryOfProfile_yTitle        = hPSet_.getParameter<std::string>("SummaryOfProfile_yTitle");
   
   int hSummaryOfProfile_NchY;
   double hSummaryOfProfile_LowY, hSummaryOfProfile_HighY;
   hSummaryOfProfile_NchY          = hPSet_.getParameter<int>("SummaryOfProfile_NchY");
   hSummaryOfProfile_LowY          = hPSet_.getParameter<double>("SummaryOfProfile_LowY");
   hSummaryOfProfile_HighY         = hPSet_.getParameter<double>("SummaryOfProfile_HighY");
   
   int nStrip, nApv, layerId_;    
   
   if(CondObj_name_ == "lorentzangle" && SummaryOnStringLevel_On_) { layerId_= getStringNameAndId(detId_).second;}
   else                                                          { layerId_= getLayerNameAndId(detId_).second;}
 
 
      
   if( CondObj_name_ == "pedestal" || CondObj_name_ == "noise"|| CondObj_name_ == "lowthreshold" || CondObj_name_ == "highthreshold" ){ // plot in strip number
     
     if( (layerId_ > 610 && layerId_ < 620) || // TID & TEC have 768 strips at maximum
         (layerId_ > 620 && layerId_ < 630) ||
         (layerId_ > 410 && layerId_ < 414) ||
         (layerId_ > 420 && layerId_ < 424) ){ nStrip =768;} 
     else { nStrip      = reader->getNumberOfApvsAndStripLength(detId_).first*128;}
     
     hSummaryOfProfile_NchX           = nStrip;
     hSummaryOfProfile_LowX           = 0.5;
     hSummaryOfProfile_HighX          = nStrip+0.5;
   
   }  
   else if( (CondObj_name_ == "lorentzangle" && SummaryOnLayerLevel_On_) || CondObj_name_ == "quality"){ // plot in detId-number
 
     // -----
     // get detIds belonging to same layer to fill X-axis with detId-number
                        
     uint32_t subDetId_ =  ((detId_>>25)&0x7);
     SiStripSubStructure substructure_;
   
     sameLayerDetIds_.clear();
   
     if(subDetId_==3){  //  TIB
       substructure_.getTIBDetectors(activeDetIds, sameLayerDetIds_,TIBDetId(detId_).layerNumber(),0,0,TIBDetId(detId_).stringNumber());  
     }
     else if(subDetId_==4){  // TID
       substructure_.getTIDDetectors(activeDetIds, sameLayerDetIds_,0,0,0,0);
     }
     else if(subDetId_==5){  // TOB
       substructure_.getTOBDetectors(activeDetIds, sameLayerDetIds_, TOBDetId(detId_).layerNumber(),0,0);
     }
     else if(subDetId_==6){  // TEC
       substructure_.getTECDetectors(activeDetIds, sameLayerDetIds_, 0,0,0,0,0,0);
     }
 
     hSummaryOfProfile_NchX           = sameLayerDetIds_.size(); 
     hSummaryOfProfile_LowX           = 0.5;
     hSummaryOfProfile_HighX          = sameLayerDetIds_.size()+0.5;
  
   } 
   else if( CondObj_name_ == "lorentzangle" && SummaryOnStringLevel_On_){ // plot in detId-number
 
     // -----
     // get detIds belonging to same string to fill X-axis with detId-number
                        
     uint32_t subDetId_ =  ((detId_>>25)&0x7);
     SiStripSubStructure substructure_;
     
     sameLayerDetIds_.clear(); 
     
     if(subDetId_==3){  //  TIB    
       if(TIBDetId(detId_).isInternalString()){
         substructure_.getTIBDetectors(activeDetIds, sameLayerDetIds_, TIBDetId(detId_).layerNumber(),0,1,TIBDetId(detId_).stringNumber()); }
       else if(TIBDetId(detId_).isExternalString()){
     substructure_.getTIBDetectors(activeDetIds, sameLayerDetIds_, TIBDetId(detId_).layerNumber(),0,2,TIBDetId(detId_).stringNumber()); }
     }
     else if(subDetId_==4){  // TID
       substructure_.getTIDDetectors(activeDetIds, sameLayerDetIds_, 0,0,0,0);
     }
     else if(subDetId_==5){  // TOB
       substructure_.getTOBDetectors(activeDetIds, sameLayerDetIds_, TOBDetId(detId_).layerNumber(),0,TOBDetId(detId_).rodNumber());
     }
     else if(subDetId_==6){  // TEC
       substructure_.getTECDetectors(activeDetIds, sameLayerDetIds_, 0,0,0,0,0,0);
     }
 
 
     hSummaryOfProfile_NchX           = sameLayerDetIds_.size(); 
     hSummaryOfProfile_LowX           = 0.5;
     hSummaryOfProfile_HighX          = sameLayerDetIds_.size()+0.5;
  
   } 
   else if( CondObj_name_ == "apvgain"){
  
     if( (layerId_ > 610 && layerId_ < 620) || // TID & TEC have 6 apvs at maximum
         (layerId_ > 620 && layerId_ < 630) ||
         (layerId_ > 410 && layerId_ < 414) ||
         (layerId_ > 420 && layerId_ < 424) ){ nApv =6;} 
     else { nApv     = reader->getNumberOfApvsAndStripLength(detId_).first;}
     
     hSummaryOfProfile_NchX           = nApv;
     hSummaryOfProfile_LowX           = 0.5;
     hSummaryOfProfile_HighX          = nApv+0.5;
  
   }
   else {
     edm::LogWarning("SiStripBaseCondObjDQM") 
        << "[SiStripBaseCondObjDQM::bookSummaryProfileMEs] PLEASE CHECK : x-axis label in your cfg"
        << std::endl; 
   }
   
   uint32_t layer_=0;
       
   layer_ = folder_organizer.GetSubDetAndLayer(detId_).second;
       
   folder_organizer.setLayerFolder(detId_,layer_); 
       
   std::string hSummaryOfProfile_name; 
   
   // ---
   int subdetectorId_ = ((detId_>>25)&0x7);
   
  
   if( subdetectorId_<3 ||subdetectorId_>6 ){ 
     edm::LogError("SiStripBaseCondObjDQM")
        << "[SiStripBaseCondObjDQM::bookSummaryProfileMEs] WRONG INPUT : no such subdetector type : "
        << subdetectorId_ << " no folder set!" 
        << std::endl;
     return;
   }
   // ---
   
   if(CondObj_name_ == "lorentzangle" && SummaryOnStringLevel_On_) { 
     hSummaryOfProfile_name = hidmanager.createHistoLayer(hSummaryOfProfile_description, "layer" , getStringNameAndId(detId_).first,"") ;
   }
   else {
     hSummaryOfProfile_name = hidmanager.createHistoLayer(hSummaryOfProfile_description, "layer" , getLayerNameAndId(detId_).first,"") ;
   }
       
   std::string hSummaryOfProfile_title;
   hSummaryOfProfile_title   = hSummaryOfProfile_name ;
            
   CondObj_ME.SummaryOfProfileDistr = dqmStore_->bookProfile(hSummaryOfProfile_name, 
                                                             hSummaryOfProfile_title, 
                                     hSummaryOfProfile_NchX, 
                                     hSummaryOfProfile_LowX, 
                                     hSummaryOfProfile_HighX, 
                                     hSummaryOfProfile_NchY, 
                                     hSummaryOfProfile_LowY, 
                                     hSummaryOfProfile_HighY);
   CondObj_ME.SummaryOfProfileDistr->setAxisTitle(hSummaryOfProfile_xTitle,1);
   CondObj_ME.SummaryOfProfileDistr->setAxisTitle(hSummaryOfProfile_yTitle,2);
   CondObj_ME.SummaryOfProfileDistr->setAxisRange(hSummaryOfProfile_LowY, hSummaryOfProfile_HighY,2);
  
   // -----
   // in order to get the right detId-number labelled in right bin of x-axis
   
   if( CondObj_name_ == "quality" ){
     
     unsigned int iBin=0;
     
     for(unsigned int i=0;i< sameLayerDetIds_.size(); i++){
     
       iBin++;
       char sameLayerDetIds_Name[1024];
       sprintf(sameLayerDetIds_Name,"%u",sameLayerDetIds_[i]);
       CondObj_ME.SummaryOfProfileDistr->setBinLabel(iBin, sameLayerDetIds_Name);
     
     }
   } 
   if( CondObj_name_ == "lorentzangle"){
 
     // Put the detIds for the -z side as following the geometrical order:
       reverse(sameLayerDetIds_.begin(), sameLayerDetIds_.begin()+sameLayerDetIds_.size()/2);
 
       unsigned int iBin=0;
        for(unsigned int i=0;i< sameLayerDetIds_.size(); i++){ 
      iBin++;
      if (!SummaryOnStringLevel_On_){
        // remove the label for detIds:
          CondObj_ME.SummaryOfProfileDistr->setBinLabel(iBin, "");
      }  
      
      if (SummaryOnStringLevel_On_){
       // Label with module position instead of detIds:
        char sameLayerDetIds_Name[1024];
        if(subdetectorId_==3){//re-abelling for TIB
          if(TIBDetId(sameLayerDetIds_[i]).isZPlusSide()){
            sprintf(sameLayerDetIds_Name,"%i",TIBDetId(sameLayerDetIds_[i]).module());}
          else if(TIBDetId(sameLayerDetIds_[i]).isZMinusSide()){
            sprintf(sameLayerDetIds_Name,"%i",-TIBDetId(sameLayerDetIds_[i]).module());}
          CondObj_ME.SummaryOfProfileDistr->setBinLabel(iBin, sameLayerDetIds_Name);
        }
        else if(subdetectorId_==5){//re-abelling for TOB
          if(TOBDetId(sameLayerDetIds_[i]).isZPlusSide())      { sprintf(sameLayerDetIds_Name,"%i",TOBDetId(sameLayerDetIds_[i]).module());}
          else if(TOBDetId(sameLayerDetIds_[i]).isZMinusSide()) { sprintf(sameLayerDetIds_Name,"%i",-TOBDetId(sameLayerDetIds_[i]).module());}
          CondObj_ME.SummaryOfProfileDistr->setBinLabel(iBin, sameLayerDetIds_Name);
        }
      }
        } 
  
        
 
 
   // -----
       
   dqmStore_->tag(CondObj_ME.SummaryOfProfileDistr, layer_);
       
   } // if "lorentzangle"
 
 }

void SiStripBaseCondObjDQM::bookTkMap ( const std::string & TkMapname )

private

Definition at line 1141 of file SiStripBaseCondObjDQM.cc.

References tkMap.

Referenced by SiStripBaseCondObjDQM().

                                                                {
   tkMap= new TrackerMap(TkMapname.c_str());
 }

void SiStripBaseCondObjDQM::end ( void )

virtual

Definition at line 1197 of file SiStripBaseCondObjDQM.cc.

Referenced by Types.LuminosityBlockRange::cppID(), Types.EventRange::cppID(), and SiStripMonitorCondData::endRun().

                                {
   edm::LogInfo("SiStripBaseCondObjDQM") 
     << "SiStripBase::end"
     << std::endl; 
 }

virtual void SiStripBaseCondObjDQM::fillMEsForDet	(	ModMEs	selModME_,
		uint32_t	selDetId_
	)

protectedpure virtual

Implemented in SiStripQualityDQM, SiStripApvGainsDQM, SiStripCablingDQM, SiStripPedestalsDQM, SiStripThresholdDQM, SiStripLorentzAngleDQM, and SiStripNoisesDQM.

Referenced by fillModMEs().

virtual void SiStripBaseCondObjDQM::fillMEsForLayer	(	std::map< uint32_t, ModMEs >	selModMEsMap_,
		uint32_t	selDetId_
	)

protectedpure virtual

Implemented in SiStripQualityDQM, SiStripApvGainsDQM, SiStripCablingDQM, SiStripLorentzAngleDQM, SiStripPedestalsDQM, SiStripThresholdDQM, and SiStripNoisesDQM.

Referenced by fillSummaryMEs().

void SiStripBaseCondObjDQM::fillModMEs ( const std::vector< uint32_t > & selectedDetIds )

protectedvirtual

Reimplemented in SiStripQualityDQM, SiStripApvGainsDQM, SiStripCablingDQM, SiStripLorentzAngleDQM, SiStripPedestalsDQM, and SiStripThresholdDQM.

Definition at line 1204 of file SiStripBaseCondObjDQM.cc.

References fillMEsForDet().

Referenced by analysis(), and analysisOnDemand().

                                                                                 {
   ModMEs CondObj_ME;
  
   for(std::vector<uint32_t>::const_iterator detIter_=selectedDetIds.begin();
                                            detIter_!=selectedDetIds.end();++detIter_){
     fillMEsForDet(CondObj_ME,*detIter_);
   }
 }

void SiStripBaseCondObjDQM::fillSummaryMEs ( const std::vector< uint32_t > & selectedDetIds )

protectedvirtual

Reimplemented in SiStripQualityDQM, SiStripLorentzAngleDQM, SiStripApvGainsDQM, SiStripCablingDQM, SiStripPedestalsDQM, and SiStripThresholdDQM.

Definition at line 1214 of file SiStripBaseCondObjDQM.cc.

References alignmentValidation::c1, CondObj_fillId_, fillMEsForLayer(), fPSet_, edm::ParameterSet::getParameter(), MonitorElement::getTH1(), MonitorElement::getTProfile(), hPSet_, mergeVDriftHistosByStation::name, SiStripBaseCondObjDQM::ModMEs::SummaryDistr, SummaryMEsMap_, SiStripBaseCondObjDQM::ModMEs::SummaryOfCumulDistr, and SiStripBaseCondObjDQM::ModMEs::SummaryOfProfileDistr.

Referenced by analysis(), and analysisOnDemand().

                                                                                     {
   
   for(std::vector<uint32_t>::const_iterator detIter_ = selectedDetIds.begin();
       detIter_!= selectedDetIds.end();detIter_++){
     fillMEsForLayer(SummaryMEsMap_, *detIter_);    
   }
 
   for (std::map<uint32_t, ModMEs>::iterator iter=SummaryMEsMap_.begin(); iter!=SummaryMEsMap_.end(); iter++){
 
     ModMEs selME;
     selME = iter->second;
 
     if(hPSet_.getParameter<bool>("FillSummaryProfileAtLayerLevel") && fPSet_.getParameter<bool>("OutputSummaryProfileAtLayerLevelAsImage")){
 
       if( CondObj_fillId_ =="onlyProfile" || CondObj_fillId_ =="ProfileAndCumul"){
 
     TCanvas c1("c1");
     selME.SummaryOfProfileDistr->getTProfile()->Draw();
     std::string name (selME.SummaryOfProfileDistr->getTProfile()->GetTitle());
     name+=".png";
     c1.Print(name.c_str());
       }
     }
     if(hPSet_.getParameter<bool>("FillSummaryAtLayerLevel") && fPSet_.getParameter<bool>("OutputSummaryAtLayerLevelAsImage")){
 
       TCanvas c1("c1");
       selME.SummaryDistr->getTH1()->Draw();
       std::string name (selME.SummaryDistr->getTH1()->GetTitle());
       name+=".png";
       c1.Print(name.c_str());
     }
     if(hPSet_.getParameter<bool>("FillCumulativeSummaryAtLayerLevel") && fPSet_.getParameter<bool>("OutputCumulativeSummaryAtLayerLevelAsImage")){
 
       if( CondObj_fillId_ =="onlyCumul" || CondObj_fillId_ =="ProfileAndCumul"){
 
     TCanvas c1("c1");
     selME.SummaryOfCumulDistr->getTH1()->Draw();
     std::string name (selME.SummaryOfCumulDistr->getTH1()->GetTitle());
     name+=".png";
     c1.Print(name.c_str());
       }
     }
 
   }
  
 }

void SiStripBaseCondObjDQM::fillTkMap	(	const uint32_t &	detid,
		const float &	value
	)

protected

Definition at line 1146 of file SiStripBaseCondObjDQM.cc.

References TrackerMap::fill(), and tkMap.

Referenced by SiStripApvGainsDQM::fillMEsForDet(), SiStripNoisesDQM::fillMEsForLayer(), SiStripThresholdDQM::fillMEsForLayer(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripPedestalsDQM::fillMEsForLayer(), SiStripApvGainsDQM::fillMEsForLayer(), SiStripQualityDQM::fillMEsForLayer(), and SiStripCablingDQM::getActiveDetIds().

                                                                               {
   tkMap->fill(detid,value);
 }

virtual void SiStripBaseCondObjDQM::getActiveDetIds ( const edm::EventSetup & eSetup )

pure virtual

Implemented in SiStripQualityDQM, SiStripCablingDQM, SiStripNoisesDQM, SiStripApvGainsDQM, SiStripLorentzAngleDQM, SiStripPedestalsDQM, and SiStripThresholdDQM.

Referenced by analysis(), and analysisOnDemand().

std::vector< uint32_t > SiStripBaseCondObjDQM::getCabledModules ( )

Definition at line 149 of file SiStripBaseCondObjDQM.cc.

References detCablingHandle_, eSetup_, and edm::EventSetup::get().

                                                             {     
  
   std::vector<uint32_t> cabledDetIds_;  
   eSetup_.get<SiStripDetCablingRcd>().get(detCablingHandle_);
   detCablingHandle_->addActiveDetectorsRawIds(cabledDetIds_);
 
   return cabledDetIds_;  
 
 }

virtual unsigned long long SiStripBaseCondObjDQM::getCache ( const edm::EventSetup & eSetup_ )

pure virtual

Implemented in SiStripQualityDQM, SiStripApvGainsDQM, SiStripCablingDQM, SiStripLorentzAngleDQM, SiStripPedestalsDQM, SiStripThresholdDQM, and SiStripNoisesDQM.

Referenced by analysis(), and analysisOnDemand().

virtual void SiStripBaseCondObjDQM::getConditionObject ( const edm::EventSetup & eSetup_ )

pure virtual

Implemented in SiStripQualityDQM, SiStripApvGainsDQM, SiStripCablingDQM, SiStripLorentzAngleDQM, SiStripPedestalsDQM, SiStripThresholdDQM, and SiStripNoisesDQM.

Referenced by analysis(), and analysisOnDemand().

std::pair< std::string, uint32_t > SiStripBaseCondObjDQM::getLayerNameAndId ( const uint32_t & detId_ )

protected

Definition at line 890 of file SiStripBaseCondObjDQM.cc.

References i, TOBDetId::layer(), TIBDetId::layer(), TIDDetId::side(), TECDetId::side(), TIDDetId::wheel(), and TECDetId::wheel().

Referenced by bookSummaryCumulMEs(), bookSummaryMEs(), bookSummaryProfileMEs(), SiStripNoisesDQM::fillMEsForLayer(), SiStripPedestalsDQM::fillMEsForLayer(), SiStripThresholdDQM::fillMEsForLayer(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripApvGainsDQM::fillMEsForLayer(), SiStripQualityDQM::fillMEsForLayer(), SiStripLorentzAngleDQM::fillSummaryMEs(), and getSummaryMEs().

                                                                                           {
 
   int subdetectorId_ = ((detId_>>25)&0x7);
   int layerId_=0;
 
   std::stringstream layerName;
   
   if( subdetectorId_ == 3 ){ //TIB
 
     for( unsigned int i = 1; i < 5; i++ ){
       if( TIBDetId( detId_ ).layer() ==i ){ 
     layerName << "TIB__layer__" << i;
     layerId_ = 300 + i;
       }
     }
     
   }
 
   else if( subdetectorId_ == 4 ){ //TIDD
     
     if( TIDDetId( detId_ ).side() == 1 ) { // TIDD side 1
       
       for( unsigned int i = 1; i < 4; i++ ){
     if(TIDDetId(detId_).wheel()==i){ 
       layerName << "TID__side__1__wheel__" << i;
       layerId_ = 410 + i;
     }
       }
       
     }
 
     else if( TIDDetId( detId_ ).side() == 2 ) { // TIDD side 2
       
       for( unsigned int i = 1; i < 4; i++ ) {
     if(TIDDetId(detId_).wheel()==i){ 
       layerName << "TID__side__2__wheel__" << i;
       layerId_ = 420 + i;
     }
       }
 
     }
 
   }
 
 
   else if( subdetectorId_ == 5 ){ // TOB
     
     for( unsigned int i = 1; i < 7; i++ ) {
       if( TOBDetId( detId_ ).layer() == i ) { 
     layerName << "TOB__layer__" << i;
     layerId_ = 500 + i;
       }
     }
     
   }
 
   else if( subdetectorId_ == 6 ){ // TEC
     
     if( TECDetId( detId_ ).side() == 1) { // TEC side 1
       
       for( unsigned int i = 1; i < 10; i++ ) {
     if( TECDetId( detId_ ).wheel() == i ) { 
       layerName << "TEC__side__1__wheel__" << i;
       layerId_ = 610 + i;
     }
       }
 
     }
 
     else if( TECDetId( detId_ ).side() == 2 ) { // TEC side 2
       
       for( unsigned int i = 1; i < 10; i++ ) {
     if( TECDetId( detId_ ).wheel() == i ) { 
       layerName << "TEC__side__2__wheel__" << i;
       layerId_ = 620 + i;
     }
       }
 
      }
   }
   
   return std::make_pair( layerName.str(), layerId_ );
 
 }

void SiStripBaseCondObjDQM::getModMEs	(	ModMEs &	CondObj_ME,
		const uint32_t &	detId_
	)

protected

Definition at line 309 of file SiStripBaseCondObjDQM.cc.

References bookCumulMEs(), bookProfileMEs(), CondObj_fillId_, CondObj_name_, SiStripBaseCondObjDQM::ModMEs::CumulDistr, ModMEsMap_, SiStripBaseCondObjDQM::ModMEs::ProfileDistr, and MonitorElement::Reset().

Referenced by SiStripNoisesDQM::fillMEsForDet(), SiStripThresholdDQM::fillMEsForDet(), SiStripApvGainsDQM::fillMEsForDet(), SiStripPedestalsDQM::fillMEsForDet(), and SiStripQualityDQM::fillMEsForDet().

                                                                                {
   
   std::map< uint32_t, ModMEs >::const_iterator ModMEsMap_iter = ModMEsMap_.find(detId_);
 
   if (ModMEsMap_iter != ModMEsMap_.end()){ 
   
     CondObj_ME=ModMEsMap_iter->second;
   
     if( ( CondObj_fillId_ =="ProfileAndCumul" || CondObj_fillId_ =="onlyProfile") && CondObj_ME.ProfileDistr ) { 
       CondObj_ME.ProfileDistr ->Reset();    
     }
        
     if( (CondObj_fillId_ =="ProfileAndCumul" || CondObj_fillId_ =="onlyCumul" ) &&  CondObj_ME.CumulDistr ){
       CondObj_ME.CumulDistr ->Reset();
     }
     else {
       edm::LogWarning("SiStripBaseCondObjDQM") 
     << "[SiStripBaseCondObjDQM::getModMEs] PLEASE CHECK : CondObj_fillId option mispelled";
     }
     return; 
 
   }
   
   // --> profile defined for all CondData
   if ( (CondObj_fillId_ =="ProfileAndCumul" || CondObj_fillId_ =="onlyProfile")) {
     bookProfileMEs(CondObj_ME,detId_);
   }  
   
   // --> cumul currently only defined for noise and apvgain
   if( (CondObj_fillId_ =="ProfileAndCumul" || CondObj_fillId_ =="onlyCumul" )
       &&(CondObj_name_ == "noise" || CondObj_name_ == "apvgain")          ) bookCumulMEs(CondObj_ME,detId_);
   
  
   ModMEsMap_.insert( std::make_pair(detId_,CondObj_ME) );
   
 }

std::vector< uint32_t > SiStripBaseCondObjDQM::GetSameLayerDetId	(	std::vector< uint32_t >	activeDetIds,
		uint32_t	selDetId
	)

protected

Definition at line 1113 of file SiStripBaseCondObjDQM.cc.

References SiStripSubStructure::getTECDetectors(), SiStripSubStructure::getTIBDetectors(), SiStripSubStructure::getTIDDetectors(), and SiStripSubStructure::getTOBDetectors().

Referenced by bookSummaryMEs(), SiStripNoisesDQM::fillMEsForLayer(), SiStripThresholdDQM::fillMEsForLayer(), SiStripPedestalsDQM::fillMEsForLayer(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripApvGainsDQM::fillMEsForLayer(), and SiStripQualityDQM::fillMEsForLayer().

                                                                                                                {
  
   std::vector<uint32_t> sameLayerDetIds;
   sameLayerDetIds.clear();
 
   SiStripSubStructure substructure_;
   
   uint32_t subselDetId_ =  ((selDetId>>25)&0x7);
 
   if(subselDetId_==3){  //  TIB
     substructure_.getTIBDetectors(activeDetIds, sameLayerDetIds, TIBDetId(selDetId).layer(),0,0,0);  
   }
   else if(subselDetId_==4){  // TID
     substructure_.getTIDDetectors(activeDetIds, sameLayerDetIds, TIDDetId(selDetId).side(),TIDDetId(selDetId).wheel(),0,0);
   }
   else if(subselDetId_==5){  // TOB
     substructure_.getTOBDetectors(activeDetIds, sameLayerDetIds, TOBDetId(selDetId).layer(),0,0);
   }
   else if(subselDetId_==6){  // TEC
     substructure_.getTECDetectors(activeDetIds, sameLayerDetIds, TECDetId(selDetId).side(),TECDetId(selDetId).wheel(),0,0,0,0);
   }
 
   return sameLayerDetIds;
   
 }

std::pair< std::string, uint32_t > SiStripBaseCondObjDQM::getStringNameAndId ( const uint32_t & detId_ )

protected

Definition at line 979 of file SiStripBaseCondObjDQM.cc.

References i, TIBDetId::isExternalString(), TIBDetId::isInternalString(), TOBDetId::layer(), TIBDetId::layer(), TOBDetId::rodNumber(), and TIBDetId::stringNumber().

Referenced by bookSummaryCumulMEs(), bookSummaryProfileMEs(), SiStripNoisesDQM::fillMEsForLayer(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripLorentzAngleDQM::fillSummaryMEs(), and getSummaryMEs().

                                                                                            {
 
   int subdetectorId_ = ((detId_>>25)&0x7);
   int layerStringId_=0;
   
   std::stringstream layerStringName;
   
   if( subdetectorId_==3 ){ //TIB
     if(TIBDetId(detId_).layer()==1 && TIBDetId(detId_).isInternalString()){ //1st layer int
       for( unsigned int i=1; i < 27 ;i++){
     if(TIBDetId(detId_).stringNumber()==i){  
       layerStringName << "TIB_L1_Int_Str_" << i;
       layerStringId_ = 30110+i; 
     }
       }      
     }
     else  if(TIBDetId(detId_).layer()==1 && TIBDetId(detId_).isExternalString()){ //1st layer ext
       for( unsigned int i=1; i < 31 ;i++){
     if(TIBDetId(detId_).stringNumber()==i){  
       layerStringName << "TIB_L1_Ext_Str_" << i;
       layerStringId_ = 301200+i; 
     }
       }      
     }
     else if(TIBDetId(detId_).layer()==2 && TIBDetId(detId_).isInternalString()){ //2nd layer int
       for( unsigned int i=1; i < 35 ;i++){
     if(TIBDetId(detId_).stringNumber()==i){  
       layerStringName << "TIB_L2_Int_Str_" << i;
       layerStringId_ = 302100+i; 
     }
       }      
     }
     else if(TIBDetId(detId_).layer()==2 && TIBDetId(detId_).isExternalString()){ //2nd layer ext
       for( unsigned int i=1; i < 39 ;i++){
     if(TIBDetId(detId_).stringNumber()==i){  
       layerStringName << "TIB_L2_Ext_Str_" << i;
       layerStringId_ = 302200+i; 
     }
       }      
     }
     else if(TIBDetId(detId_).layer()==3 && TIBDetId(detId_).isInternalString()){ //3rd layer int
       for( unsigned int i=1; i < 45 ;i++){
     if(TIBDetId(detId_).stringNumber()==i){  
       layerStringName << "TIB_L3_Int_Str_" << i;
       layerStringId_ = 303100+i; 
     }
       }      
     }
     else if(TIBDetId(detId_).layer()==3 && TIBDetId(detId_).isExternalString()){ //3rd layer ext
       for( unsigned int i=1; i < 47 ;i++){
     if(TIBDetId(detId_).stringNumber()==i){  
       layerStringName << "TIB_L3_Ext_Str_" << i;
       layerStringId_ = 303200+i; 
     }
       }      
     }
     else if(TIBDetId(detId_).layer()==4 && TIBDetId(detId_).isInternalString()){ //4th layer int
       for( unsigned int i=1; i < 53 ;i++){
     if(TIBDetId(detId_).stringNumber()==i){  
       layerStringName << "TIB_L4_Int_Str_" << i;
       layerStringId_ = 304100+i; 
     }
       }      
     }
     else if(TIBDetId(detId_).layer()==4 && TIBDetId(detId_).isExternalString()){ //4th layer ext
       for( unsigned int i=1; i < 57 ;i++){
     if(TIBDetId(detId_).stringNumber()==i){  
       layerStringName << "TIB_L4_Ext_Str_" << i;
       layerStringId_ = 304200+i; 
     }
       }      
     }
   } //TIB
 
 
   else if( subdetectorId_==5 ){ // TOB
     if(TOBDetId(detId_).layer()==1){ //1st layer
       for( unsigned int i=1; i < 43 ;i++){
     if(TOBDetId(detId_).rodNumber()==i){  
       layerStringName << "TOB_L1_Rod_" << i;
       layerStringId_ = 50100+i;
     }
       }      
     }
     else if(TOBDetId(detId_).layer()==2){ //2nd layer
       for( unsigned int i=1; i < 49 ;i++){
     if(TOBDetId(detId_).rodNumber()==i){  
       layerStringName << "TOB_L2_Rod_" << i;
       layerStringId_ = 50200+i; 
     }
       }      
     }
     else if(TOBDetId(detId_).layer()==3){ //3rd layer
       for( unsigned int i=1; i < 55 ;i++){
     if(TOBDetId(detId_).rodNumber()==i){  
       layerStringName << "TOB_L3_Rod_" << i;
       layerStringId_ = 50300+i; 
     }
       }      
     }
     else if(TOBDetId(detId_).layer()==4){ //4th layer
       for( unsigned int i=1; i < 61 ;i++){
     if(TOBDetId(detId_).rodNumber()==i){  
       layerStringName << "TOB_L4_Rod_" << i;
       layerStringId_ = 50400+i; 
     }
       }      
     }
     else if(TOBDetId(detId_).layer()==5){ //5th layer
       for( unsigned int i=1; i < 67 ;i++){
     if(TOBDetId(detId_).rodNumber()==i){  
       layerStringName << "TOB_L5_Rod_" << i;
       layerStringId_ = 50500+i; 
     }
       }      
     }
     else if(TOBDetId(detId_).layer()==6){ //6st layer
       for( unsigned int i=1; i < 75 ;i++){
     if(TOBDetId(detId_).rodNumber()==i){  
       layerStringName << "TOB_L6_Rod_" << i;
       layerStringId_ = 50600+i; 
     }
       }      
     }
   }//TOB
 
   return std::make_pair( layerStringName.str(), layerStringId_ );
 
 }

void SiStripBaseCondObjDQM::getSummaryMEs	(	ModMEs &	CondObj_ME,
		const uint32_t &	detId_
	)

protected

Definition at line 350 of file SiStripBaseCondObjDQM.cc.

References bookSummaryCumulMEs(), bookSummaryMEs(), bookSummaryProfileMEs(), CondObj_fillId_, CondObj_name_, getLayerNameAndId(), edm::ParameterSet::getParameter(), getStringNameAndId(), hPSet_, edm::second(), SiStripBaseCondObjDQM::ModMEs::SummaryDistr, SummaryMEsMap_, SiStripBaseCondObjDQM::ModMEs::SummaryOfCumulDistr, SiStripBaseCondObjDQM::ModMEs::SummaryOfProfileDistr, and SummaryOnStringLevel_On_.

Referenced by SiStripNoisesDQM::fillMEsForLayer(), SiStripThresholdDQM::fillMEsForLayer(), SiStripPedestalsDQM::fillMEsForLayer(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripApvGainsDQM::fillMEsForLayer(), and SiStripQualityDQM::fillMEsForLayer().

                                                                                    {
 
   std::map<uint32_t, ModMEs>::const_iterator SummaryMEsMap_iter;
 
   if(CondObj_name_ == "lorentzangle" && SummaryOnStringLevel_On_ ){
     SummaryMEsMap_iter = SummaryMEsMap_.find(getStringNameAndId(detId_).second);
   }
   else {
     SummaryMEsMap_iter = SummaryMEsMap_.find(getLayerNameAndId(detId_).second);
   }
    
   if (SummaryMEsMap_iter != SummaryMEsMap_.end()){ return;}
 
   //FIXME t's not good that the base class has to know about which derived class shoudl exist.
   // please modify this part. implement virtual functions, esplicited in the derived classes
   // --> currently only profile summary defined for all condition objects except quality
   if(  (CondObj_fillId_ =="ProfileAndCumul" || CondObj_fillId_ =="onlyProfile" ) &&
      (
       CondObj_name_ == "pedestal"     || 
       CondObj_name_ == "noise"         || 
       CondObj_name_ == "lowthreshold"  || 
       CondObj_name_ == "highthreshold" || 
       CondObj_name_ == "apvgain"       || 
       CondObj_name_ == "lorentzangle") ) {
     if(hPSet_.getParameter<bool>("FillSummaryProfileAtLayerLevel")) 
            if (CondObj_ME.SummaryOfProfileDistr) { bookSummaryProfileMEs(CondObj_ME,detId_);
       }  
     
 }
     
   // --> currently only genuine cumul LA
   if(   (CondObj_fillId_ =="ProfileAndCumul" || CondObj_fillId_ =="onlyCumul" ) &&
     (CondObj_name_ == "lorentzangle" ||  CondObj_name_ == "noise")  ) {
     if(hPSet_.getParameter<bool>("FillCumulativeSummaryAtLayerLevel"))
       if (CondObj_ME.SummaryOfCumulDistr) { bookSummaryCumulMEs(CondObj_ME,detId_); } 
   } 
                           
   // --> currently only summary as a function of detId for noise, pedestal and apvgain 
   if(      CondObj_name_ == "noise"         ||
        CondObj_name_ == "lowthreshold"  || 
        CondObj_name_ == "highthreshold" || 
        CondObj_name_ == "apvgain"       || 
        CondObj_name_ == "pedestal"      || 
        CondObj_name_ == "quality"           ) {
     if(hPSet_.getParameter<bool>("FillSummaryAtLayerLevel"))          
       if (CondObj_ME.SummaryDistr) { bookSummaryMEs(CondObj_ME,detId_); } 
     
   } 
                           
   if(CondObj_name_ == "lorentzangle" && SummaryOnStringLevel_On_) {
     //FIXME getStringNameandId takes time. not need to call it every timne. put the call at the beginning of the method and caache the string 
     SummaryMEsMap_.insert( std::make_pair(getStringNameAndId(detId_).second,CondObj_ME) );
   }
   else {
     SummaryMEsMap_.insert( std::make_pair(getLayerNameAndId(detId_).second,CondObj_ME) );
   }
 
 }

void SiStripBaseCondObjDQM::saveTkMap	(	const std::string &	TkMapname,
		double	minValue,
		double	maxValue
	)

private

Definition at line 1151 of file SiStripBaseCondObjDQM.cc.

References python.tagInventory::entries, edm::ParameterSet::getParameter(), hPSet_, i, j, max(), min, TrackerMap::save(), TrackerMap::setPalette(), TrackerMap::showPalette(), tkMap, and tkMapScaler.

Referenced by analysis().

                                                                                                  {
   if(tkMapScaler.size()!=0){
     //check that saturation is below x%  below minValue and above minValue, and in case re-arrange.
     float th=hPSet_.getParameter<double>("saturatedFraction");
 
     size_t imin=0,imax=0;
     float entries=0 ;
     for(size_t i=0;i<tkMapScaler.size();++i)
       entries+=tkMapScaler[i];
 
     float min=0 ;
     for(size_t i=0;(i<tkMapScaler.size()) && (min<th);++i){
       min+=tkMapScaler[i]/entries;
       imin=i;
     }
 
     float max=0;
     // for(size_t i=tkMapScaler.size()-1;(i>=0) && (max<th);--i){ // Wrong
     // Since i is unsigned, i >= 0 is always true,
     // and the loop termination condition is never reached.
     // We offset the loop index by one to fix this.
     for(size_t j=tkMapScaler.size();(j>0) && (max<th);--j){ 
       size_t i = j - 1;
       max+=tkMapScaler[i]/entries;
       imax=i;
     }
     
     //reset maxValue;
     if(maxValue<imax){
       edm::LogInfo("")<< "Resetting TkMap maxValue from " << maxValue << " to " << imax;
       maxValue=imax;
     }
     //reset minValue;
     if(minValue>imin){
       edm::LogInfo("")<< "Resetting TkMap minValue from " << minValue << " to " << imin;
       minValue=imin;
     }
   }
 
   tkMap->save(false, minValue, maxValue, TkMapname.c_str());
   tkMap->setPalette(1); tkMap->showPalette(true);
 
 }

void SiStripBaseCondObjDQM::selectModules ( std::vector< uint32_t > & detIds_ )

Definition at line 166 of file SiStripBaseCondObjDQM.cc.

References fPSet_, edm::ParameterSet::getParameter(), SiStripSubStructure::getTECDetectors(), SiStripSubStructure::getTIBDetectors(), SiStripSubStructure::getTIDDetectors(), SiStripSubStructure::getTOBDetectors(), i, ModulesToBeExcluded_, ModulesToBeIncluded_, python.multivaluedict::sort(), SubDetectorsToBeExcluded_, swap(), and tmp.

Referenced by analysis().

                                                                       {
    
   ModulesToBeExcluded_     = fPSet_.getParameter< std::vector<unsigned int> >("ModulesToBeExcluded");
   ModulesToBeIncluded_     = fPSet_.getParameter< std::vector<unsigned int> >("ModulesToBeIncluded");
   SubDetectorsToBeExcluded_= fPSet_.getParameter< std::vector<std::string> >("SubDetectorsToBeExcluded");  
  
   if(fPSet_.getParameter<bool>("restrictModules") 
      && ModulesToBeExcluded_.size()==0 
      && ModulesToBeIncluded_.size()==0 ){
     edm::LogWarning("SiStripBaseCondObjDQM") 
        << "[SiStripBaseCondObjDQM::selectModules] PLEASE CHECK : no modules to be exclude/included in your cfg"
        << std::endl; 
   }
 
  
  
   
   // --> detIds to start with
 
   if( fPSet_.getParameter<bool>("restrictModules")){
       
     if( ModulesToBeIncluded_.size()>0 ){
      std::vector<uint32_t> tmp;
      tmp.clear();
      set_intersection( detIds_.begin(), detIds_.end(),  
                        ModulesToBeIncluded_.begin(), ModulesToBeIncluded_.end(),
                inserter(tmp,tmp.begin()));
      swap(detIds_,tmp);
    }
     
   }
 
   
 
   // -----
   // *** exclude modules ***
   
   if( fPSet_.getParameter<bool>("restrictModules") ){
     
     std::sort(detIds_.begin(),detIds_.end());
 
     for( std::vector<uint32_t>::const_iterator modIter_  = ModulesToBeExcluded_.begin(); 
                                                modIter_ != ModulesToBeExcluded_.end(); modIter_++){
       
       std::vector<uint32_t>::iterator detIter_=std::lower_bound(detIds_.begin(),detIds_.end(),*modIter_);
       detIds_.erase(detIter_);
       detIter_--;
      
     }
   
   }
   // *** exclude modules ***
   // -----
    
 
   // -----
   // *** restrict to a particular subdetector ***
    
   if( *(SubDetectorsToBeExcluded_.begin()) !="none" ){
     
     std::vector<uint32_t> tmp;
 
     SiStripSubStructure substructure_;
     
     for( std::vector<std::string>::const_iterator modIter_  = SubDetectorsToBeExcluded_.begin(); 
                                                  modIter_ != SubDetectorsToBeExcluded_.end(); modIter_++){
       tmp.clear();
 
       if (*modIter_=="TIB")     { substructure_.getTIBDetectors(detIds_, tmp, 0,0,0,0);}
       else if (*modIter_=="TOB") { substructure_.getTOBDetectors(detIds_, tmp, 0,0,0);}
       else if (*modIter_=="TID") { substructure_.getTIDDetectors(detIds_, tmp, 0,0,0,0);}
       else if (*modIter_=="TEC") { substructure_.getTECDetectors(detIds_, tmp, 0,0,0,0,0,0);}
       else {
         edm::LogWarning("SiStripBaseCondObjDQM") 
        << "[SiStripBaseCondObjDQM::selectModules] PLEASE CHECK : no correct (name) subdetector to be excluded in your cfg"
        << std::endl; 
       }
 
       std::vector<uint32_t>::iterator iterBegin_=std::lower_bound(detIds_.begin(),
                                                               detIds_.end(),
                                   *min_element(tmp.begin(), tmp.end()));
                                   
       std::vector<uint32_t>::iterator iterEnd_=std::lower_bound(detIds_.begin(),
                                                             detIds_.end(),
                                 *max_element(tmp.begin(), tmp.end()));
                                 
       for(std::vector<uint32_t>::iterator detIter_ = iterEnd_;
                                          detIter_!= iterBegin_-1;detIter_--){
       detIds_.erase(detIter_);
       } 
 
     } // loop SubDetectorsToBeExcluded_
   }
 
   
   // -----
   // *** fill only one Module per layer ***
 
   if(fPSet_.getParameter<std::string>("ModulesToBeFilled") == "onlyOneModulePerLayer"){
    
     std::vector<uint32_t> tmp;
     std::vector<uint32_t> layerDetIds;
 
     SiStripSubStructure substructure_;
     
     for(unsigned int i=1; i<5 ; i++){
       tmp.clear();
       substructure_.getTIBDetectors(detIds_, tmp, i,0,0,0);
       if(tmp.size() !=0) { layerDetIds.push_back(*(tmp.begin()));}
     }
     for(unsigned int i=1; i<7 ; i++){
       tmp.clear();
       substructure_.getTOBDetectors(detIds_, tmp, i,0,0);
       if(tmp.size() !=0) { layerDetIds.push_back(*(tmp.begin()));}
     }
     for(unsigned int i=1; i<4 ; i++){
       tmp.clear();
       substructure_.getTIDDetectors(detIds_, tmp, 1,i,0,0);
       if(tmp.size() !=0) { layerDetIds.push_back(*(tmp.begin()));}
       substructure_.getTIDDetectors(detIds_, tmp, 2,i,0,0);
       if(tmp.size() !=0) { layerDetIds.push_back(*(tmp.begin()));}
     }  
     for(unsigned int i=1; i<10 ; i++){
       tmp.clear();
       substructure_.getTECDetectors(detIds_, tmp, 1,i,0,0,0,0);
       if(tmp.size() !=0) { layerDetIds.push_back(*(tmp.begin()));}
       substructure_.getTECDetectors(detIds_, tmp, 2,i,0,0,0,0);
       if(tmp.size() !=0) { layerDetIds.push_back(*(tmp.begin()));}
     }
     
     detIds_.clear();
     detIds_=layerDetIds;
 
   }
   // -----
 
    
 } //selectModules