#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_, const 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 (const ModMEs &selModME_, uint32_t selDetId_, const TrackerTopology *tTopo)=0

virtual void	fillMEsForLayer (uint32_t selDetId_, const TrackerTopology *tTopo)=0

virtual void	fillModMEs (const std::vector< uint32_t > &selectedDetIds, const edm::EventSetup &es)

virtual void	fillSummaryMEs (const std::vector< uint32_t > &selectedDetIds, const edm::EventSetup &es)

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

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

void	getModMEs (ModMEs &CondObj_ME, const uint32_t &detId_, const TrackerTopology *tTopo)

std::vector< uint32_t >	GetSameLayerDetId (const std::vector< uint32_t > &activeDetIds, uint32_t selDetId, const TrackerTopology *tTopo)

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

void	getSummaryMEs (ModMEs &CondObj_ME, const uint32_t &detId_, const TrackerTopology *tTopo)

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_, const TrackerTopology *tTopo)

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

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

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

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

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

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_, contentValuesFiles::fullPath, edm::ParameterSet::getParameter(), GrandSummary_On_, HistoMaps_On_, hPSet_, maxValue, minValue, Mod_On_, reader, OfflineSiStripConfigDb_cff::SiStripDetInfoFileReader, AlCaHLTBitMon_QueryRunRegistry::string, 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(), AlCaHLTBitMon_QueryRunRegistry::string, 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, eSetup_); }
   if(SummaryOnLayerLevel_On_ || SummaryOnStringLevel_On_ ){ fillSummaryMEs(activeDetIds, eSetup_); }
 
   if(fPSet_.getParameter<bool>("TkMap_On") || hPSet_.getParameter<bool>("TkMap_On")) {
     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 119 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_, eSetup_); 
   
 }

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

Definition at line 136 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, eSetup_); 
   
 }

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

Definition at line 86 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_,
		const TrackerTopology *	tTopo
	)

private

Definition at line 503 of file SiStripBaseCondObjDQM.cc.

References DQMStore::book1D(), SiStripHistoId::createHistoId(), SiStripBaseCondObjDQM::ModMEs::CumulDistr, dqmStore_, folder_organizer, edm::ParameterSet::getParameter(), hidmanager, hPSet_, MonitorElement::setAxisTitle(), SiStripFolderOrganizer::setDetectorFolder(), AlCaHLTBitMon_QueryRunRegistry::string, 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_,tTopo); 
       
   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_,
		const TrackerTopology *	tTopo
	)

private

Definition at line 454 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(), AlCaHLTBitMon_QueryRunRegistry::string, 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_,tTopo); 
       
   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_,
		const TrackerTopology *	tTopo
	)

private

Definition at line 771 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(), AlCaHLTBitMon_QueryRunRegistry::string, 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_, tTopo).second;
     
   folder_organizer.setLayerFolder(detId_,tTopo,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 and BP Histos are plotted for each string:
   if((CondObj_name_ == "lorentzangle" || CondObj_name_ == "bpcorrection") && SummaryOnStringLevel_On_) { 
     hSummaryOfCumul_name = hidmanager.createHistoLayer(hSummaryOfCumul_description, "layer" , getStringNameAndId(detId_,tTopo).first, "") ;
   }
   else {  
     hSummaryOfCumul_name = hidmanager.createHistoLayer(hSummaryOfCumul_description, "layer" , getLayerNameAndId(detId_,tTopo).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_,
		const TrackerTopology *	tTopo
	)

private

Definition at line 837 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(), AlCaHLTBitMon_QueryRunRegistry::string, 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_,tTopo);
 
   hSummary_NchX           = sameLayerDetIds_.size(); 
   hSummary_LowX           = 0.5;
   hSummary_HighX          = sameLayerDetIds_.size()+0.5;
   
   uint32_t layer_=0;
       
   layer_ = folder_organizer.GetSubDetAndLayer(detId_, tTopo).second;
       
   folder_organizer.setLayerFolder(detId_,tTopo,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_,tTopo).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, 
                            0., 
                            0.);
   //                           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_,
		const TrackerTopology *	tTopo
	)

private

Definition at line 544 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_, reader, MonitorElement::setAxisRange(), MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), SiStripFolderOrganizer::setLayerFolder(), AlCaHLTBitMon_QueryRunRegistry::string, SiStripBaseCondObjDQM::ModMEs::SummaryOfProfileDistr, SummaryOnLayerLevel_On_, SummaryOnStringLevel_On_, DQMStore::tag(), TrackerTopology::tibIsExternalString(), TrackerTopology::tibIsInternalString(), TrackerTopology::tibIsZMinusSide(), TrackerTopology::tibIsZPlusSide(), TrackerTopology::tibLayer(), TrackerTopology::tibModule(), TrackerTopology::tibString(), TrackerTopology::tobIsZMinusSide(), TrackerTopology::tobIsZPlusSide(), TrackerTopology::tobLayer(), TrackerTopology::tobModule(), and TrackerTopology::tobRod().

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_,tTopo).second;}
   else if(CondObj_name_ == "bpcorrection" && SummaryOnStringLevel_On_) { layerId_= getStringNameAndId(detId_,tTopo).second;}
   else                                                            { layerId_= getLayerNameAndId(detId_,tTopo).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" || CondObj_name_ == "bpcorrection" ) && 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_,tTopo->tibLayer(detId_),0,0,tTopo->tibString(detId_));  
     }
     else if(subDetId_==4){  // TID
       substructure_.getTIDDetectors(activeDetIds, sameLayerDetIds_,0,0,0,0);
     }
     else if(subDetId_==5){  // TOB
       substructure_.getTOBDetectors(activeDetIds, sameLayerDetIds_, tTopo->tobLayer(detId_),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" || CondObj_name_ == "bpcorrection") && 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(tTopo->tibIsInternalString(detId_)){
         substructure_.getTIBDetectors(activeDetIds, sameLayerDetIds_, tTopo->tibLayer(detId_),0,1,tTopo->tibString(detId_)); }
       else if(tTopo->tibIsExternalString(detId_)){
     substructure_.getTIBDetectors(activeDetIds, sameLayerDetIds_, tTopo->tibLayer(detId_),0,2,tTopo->tibString(detId_)); }
     }
     else if(subDetId_==4){  // TID
       substructure_.getTIDDetectors(activeDetIds, sameLayerDetIds_, 0,0,0,0);
     }
     else if(subDetId_==5){  // TOB
       substructure_.getTOBDetectors(activeDetIds, sameLayerDetIds_, tTopo->tobLayer(detId_),0,tTopo->tobRod(detId_));
     }
     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_, tTopo).second;
       
   folder_organizer.setLayerFolder(detId_,tTopo,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" || CondObj_name_ == "bpcorrection") && SummaryOnStringLevel_On_) { 
     hSummaryOfProfile_name = hidmanager.createHistoLayer(hSummaryOfProfile_description, "layer" , getStringNameAndId(detId_,tTopo).first,"") ;
   }
   else {
     hSummaryOfProfile_name = hidmanager.createHistoLayer(hSummaryOfProfile_description, "layer" , getLayerNameAndId(detId_,tTopo).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, 
                                     0., 
                                     0.);
   //                                    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" || CondObj_name_ == "bpcorrection"){
 
     // 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(tTopo->tibIsZPlusSide(sameLayerDetIds_[i])){
            sprintf(sameLayerDetIds_Name,"%i",tTopo->tibModule(sameLayerDetIds_[i]));}
          else if(tTopo->tibIsZMinusSide(sameLayerDetIds_[i])){
            sprintf(sameLayerDetIds_Name,"%i",-tTopo->tibModule(sameLayerDetIds_[i]));}
          CondObj_ME.SummaryOfProfileDistr->setBinLabel(iBin, sameLayerDetIds_Name);
        }
        else if(subdetectorId_==5){//re-abelling for TOB
          if(tTopo->tobIsZPlusSide(sameLayerDetIds_[i]))      { sprintf(sameLayerDetIds_Name,"%i",tTopo->tobModule(sameLayerDetIds_[i]));}
          else if(tTopo->tobIsZMinusSide(sameLayerDetIds_[i])) { sprintf(sameLayerDetIds_Name,"%i",-tTopo->tobModule(sameLayerDetIds_[i]));}
          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 1187 of file SiStripBaseCondObjDQM.cc.

References tkMap.

Referenced by SiStripBaseCondObjDQM().

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

void SiStripBaseCondObjDQM::end ( void )

virtual

Definition at line 1244 of file SiStripBaseCondObjDQM.cc.

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

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

virtual void SiStripBaseCondObjDQM::fillMEsForDet	(	const ModMEs &	selModME_,
		uint32_t	selDetId_,
		const TrackerTopology *	tTopo
	)

protectedpure virtual

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

Referenced by fillModMEs().

virtual void SiStripBaseCondObjDQM::fillMEsForLayer	(	uint32_t	selDetId_,
		const TrackerTopology *	tTopo
	)

protectedpure virtual

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

Referenced by fillSummaryMEs().

void SiStripBaseCondObjDQM::fillModMEs	(	const std::vector< uint32_t > &	selectedDetIds,
		const edm::EventSetup &	es
	)

protectedvirtual

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

Definition at line 1249 of file SiStripBaseCondObjDQM.cc.

References fillMEsForDet(), edm::EventSetup::get(), and edm::ESHandle< class >::product().

Referenced by analysis(), and analysisOnDemand().

                                                                                                          {
 
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<IdealGeometryRecord>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   ModMEs CondObj_ME;
  
   for(std::vector<uint32_t>::const_iterator detIter_=selectedDetIds.begin();
       detIter_!=selectedDetIds.end();++detIter_){
     fillMEsForDet(CondObj_ME,*detIter_,tTopo);
   }
 }

void SiStripBaseCondObjDQM::fillSummaryMEs	(	const std::vector< uint32_t > &	selectedDetIds,
		const edm::EventSetup &	es
	)

protectedvirtual

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

Definition at line 1265 of file SiStripBaseCondObjDQM.cc.

References alignmentValidation::c1, CondObj_fillId_, fillMEsForLayer(), fPSet_, edm::EventSetup::get(), edm::ParameterSet::getParameter(), MonitorElement::getTH1(), MonitorElement::getTProfile(), hPSet_, getDQMSummary::iter, mergeVDriftHistosByStation::name, edm::ESHandle< class >::product(), AlCaHLTBitMon_QueryRunRegistry::string, SiStripBaseCondObjDQM::ModMEs::SummaryDistr, SummaryMEsMap_, SiStripBaseCondObjDQM::ModMEs::SummaryOfCumulDistr, and SiStripBaseCondObjDQM::ModMEs::SummaryOfProfileDistr.

Referenced by analysis(), and analysisOnDemand().

                                                                                                              {
 
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<IdealGeometryRecord>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
   
   for(std::vector<uint32_t>::const_iterator detIter_ = selectedDetIds.begin();
       detIter_!= selectedDetIds.end();detIter_++){
     fillMEsForLayer(/*SummaryMEsMap_,*/ *detIter_,tTopo);    
   }
 
   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 1192 of file SiStripBaseCondObjDQM.cc.

References TrackerMap::fill(), and tkMap.

Referenced by SiStripApvGainsDQM::fillMEsForDet(), SiStripNoisesDQM::fillMEsForLayer(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripPedestalsDQM::fillMEsForLayer(), SiStripThresholdDQM::fillMEsForLayer(), SiStripBackPlaneCorrectionDQM::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, SiStripApvGainsDQM, SiStripNoisesDQM, SiStripBackPlaneCorrectionDQM, SiStripLorentzAngleDQM, SiStripPedestalsDQM, and SiStripThresholdDQM.

Referenced by analysis(), and analysisOnDemand().

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

Definition at line 151 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, SiStripBackPlaneCorrectionDQM, 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, SiStripBackPlaneCorrectionDQM, SiStripLorentzAngleDQM, SiStripPedestalsDQM, SiStripThresholdDQM, and SiStripNoisesDQM.

Referenced by analysis(), and analysisOnDemand().

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

protected

Definition at line 936 of file SiStripBaseCondObjDQM.cc.

References i, TrackerTopology::tecSide(), TrackerTopology::tecWheel(), TrackerTopology::tibLayer(), TrackerTopology::tidSide(), TrackerTopology::tidWheel(), and TrackerTopology::tobLayer().

Referenced by bookSummaryCumulMEs(), bookSummaryMEs(), bookSummaryProfileMEs(), SiStripNoisesDQM::fillMEsForLayer(), SiStripPedestalsDQM::fillMEsForLayer(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripBackPlaneCorrectionDQM::fillMEsForLayer(), SiStripThresholdDQM::fillMEsForLayer(), SiStripApvGainsDQM::fillMEsForLayer(), SiStripQualityDQM::fillMEsForLayer(), SiStripBackPlaneCorrectionDQM::fillSummaryMEs(), 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( tTopo->tibLayer( detId_ ) ==i ){ 
     layerName << "TIB__layer__" << i;
     layerId_ = 300 + i;
       }
     }
     
   }
 
   else if( subdetectorId_ == 4 ){ //TIDD
     
     if( tTopo->tidSide( detId_ ) == 1 ) { // TIDD side 1
       
       for( unsigned int i = 1; i < 4; i++ ){
     if(tTopo->tidWheel(detId_)==i){ 
       layerName << "TID__side__1__wheel__" << i;
       layerId_ = 410 + i;
     }
       }
       
     }
 
     else if( tTopo->tidSide( detId_ ) == 2 ) { // TIDD side 2
       
       for( unsigned int i = 1; i < 4; i++ ) {
     if(tTopo->tidWheel(detId_)==i){ 
       layerName << "TID__side__2__wheel__" << i;
       layerId_ = 420 + i;
     }
       }
 
     }
 
   }
 
 
   else if( subdetectorId_ == 5 ){ // TOB
     
     for( unsigned int i = 1; i < 7; i++ ) {
       if( tTopo->tobLayer( detId_ ) == i ) { 
     layerName << "TOB__layer__" << i;
     layerId_ = 500 + i;
       }
     }
     
   }
 
   else if( subdetectorId_ == 6 ){ // TEC
     
     if( tTopo->tecSide( detId_ ) == 1) { // TEC side 1
       
       for( unsigned int i = 1; i < 10; i++ ) {
     if( tTopo->tecWheel( detId_ ) == i ) { 
       layerName << "TEC__side__1__wheel__" << i;
       layerId_ = 610 + i;
     }
       }
 
     }
 
     else if( tTopo->tecSide( detId_ ) == 2 ) { // TEC side 2
       
       for( unsigned int i = 1; i < 10; i++ ) {
     if( tTopo->tecWheel( detId_ ) == 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_,
		const TrackerTopology *	tTopo
	)

protected

Definition at line 342 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(), SiStripPedestalsDQM::fillMEsForDet(), SiStripApvGainsDQM::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_,tTopo);
   }  
   
   // --> 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_,tTopo);
   
  
   ModMEsMap_.insert( std::make_pair(detId_,CondObj_ME) );
   
 }

std::vector< uint32_t > SiStripBaseCondObjDQM::GetSameLayerDetId	(	const std::vector< uint32_t > &	activeDetIds,
		uint32_t	selDetId,
		const TrackerTopology *	tTopo
	)

protected

Definition at line 1159 of file SiStripBaseCondObjDQM.cc.

References SiStripSubStructure::getTECDetectors(), SiStripSubStructure::getTIBDetectors(), SiStripSubStructure::getTIDDetectors(), SiStripSubStructure::getTOBDetectors(), TrackerTopology::tecSide(), TrackerTopology::tecWheel(), TrackerTopology::tibLayer(), TrackerTopology::tidSide(), TrackerTopology::tidWheel(), and TrackerTopology::tobLayer().

Referenced by bookSummaryMEs(), SiStripNoisesDQM::fillMEsForLayer(), SiStripThresholdDQM::fillMEsForLayer(), SiStripPedestalsDQM::fillMEsForLayer(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripBackPlaneCorrectionDQM::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, tTopo->tibLayer(selDetId),0,0,0);  
   }
   else if(subselDetId_==4){  // TID
     substructure_.getTIDDetectors(activeDetIds, sameLayerDetIds, tTopo->tidSide(selDetId),tTopo->tidWheel(selDetId),0,0);
   }
   else if(subselDetId_==5){  // TOB
     substructure_.getTOBDetectors(activeDetIds, sameLayerDetIds, tTopo->tobLayer(selDetId),0,0);
   }
   else if(subselDetId_==6){  // TEC
     substructure_.getTECDetectors(activeDetIds, sameLayerDetIds, tTopo->tecSide(selDetId),tTopo->tecWheel(selDetId),0,0,0,0);
   }
 
   return sameLayerDetIds;
   
 }

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

protected

Definition at line 1025 of file SiStripBaseCondObjDQM.cc.

References i, TrackerTopology::tibIsExternalString(), TrackerTopology::tibIsInternalString(), TrackerTopology::tibLayer(), TrackerTopology::tibString(), TrackerTopology::tobLayer(), and TrackerTopology::tobRod().

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

                                                                                                                          {
 
   int subdetectorId_ = ((detId_>>25)&0x7);
   int layerStringId_=0;
   
   std::stringstream layerStringName;
   
   if( subdetectorId_==3 ){ //TIB
     if(tTopo->tibLayer(detId_)==1 && tTopo->tibIsInternalString(detId_)){ //1st layer int
       for( unsigned int i=1; i < 27 ;i++){
     if(tTopo->tibString(detId_)==i){  
       layerStringName << "TIB_L1_Int_Str_" << i;
       layerStringId_ = 30110+i; 
     }
       }      
     }
     else  if(tTopo->tibLayer(detId_)==1 && tTopo->tibIsExternalString(detId_)){ //1st layer ext
       for( unsigned int i=1; i < 31 ;i++){
     if(tTopo->tibString(detId_)==i){  
       layerStringName << "TIB_L1_Ext_Str_" << i;
       layerStringId_ = 301200+i; 
     }
       }      
     }
     else if(tTopo->tibLayer(detId_)==2 && tTopo->tibIsInternalString(detId_)){ //2nd layer int
       for( unsigned int i=1; i < 35 ;i++){
     if(tTopo->tibString(detId_)==i){  
       layerStringName << "TIB_L2_Int_Str_" << i;
       layerStringId_ = 302100+i; 
     }
       }      
     }
     else if(tTopo->tibLayer(detId_)==2 && tTopo->tibIsExternalString(detId_)){ //2nd layer ext
       for( unsigned int i=1; i < 39 ;i++){
     if(tTopo->tibString(detId_)==i){  
       layerStringName << "TIB_L2_Ext_Str_" << i;
       layerStringId_ = 302200+i; 
     }
       }      
     }
     else if(tTopo->tibLayer(detId_)==3 && tTopo->tibIsInternalString(detId_)){ //3rd layer int
       for( unsigned int i=1; i < 45 ;i++){
     if(tTopo->tibString(detId_)==i){  
       layerStringName << "TIB_L3_Int_Str_" << i;
       layerStringId_ = 303100+i; 
     }
       }      
     }
     else if(tTopo->tibLayer(detId_)==3 && tTopo->tibIsExternalString(detId_)){ //3rd layer ext
       for( unsigned int i=1; i < 47 ;i++){
     if(tTopo->tibString(detId_)==i){  
       layerStringName << "TIB_L3_Ext_Str_" << i;
       layerStringId_ = 303200+i; 
     }
       }      
     }
     else if(tTopo->tibLayer(detId_)==4 && tTopo->tibIsInternalString(detId_)){ //4th layer int
       for( unsigned int i=1; i < 53 ;i++){
     if(tTopo->tibString(detId_)==i){  
       layerStringName << "TIB_L4_Int_Str_" << i;
       layerStringId_ = 304100+i; 
     }
       }      
     }
     else if(tTopo->tibLayer(detId_)==4 && tTopo->tibIsExternalString(detId_)){ //4th layer ext
       for( unsigned int i=1; i < 57 ;i++){
     if(tTopo->tibString(detId_)==i){  
       layerStringName << "TIB_L4_Ext_Str_" << i;
       layerStringId_ = 304200+i; 
     }
       }      
     }
   } //TIB
 
 
   else if( subdetectorId_==5 ){ // TOB
     if(tTopo->tobLayer(detId_)==1){ //1st layer
       for( unsigned int i=1; i < 43 ;i++){
     if(tTopo->tobRod(detId_)==i){  
       layerStringName << "TOB_L1_Rod_" << i;
       layerStringId_ = 50100+i;
     }
       }      
     }
     else if(tTopo->tobLayer(detId_)==2){ //2nd layer
       for( unsigned int i=1; i < 49 ;i++){
     if(tTopo->tobRod(detId_)==i){  
       layerStringName << "TOB_L2_Rod_" << i;
       layerStringId_ = 50200+i; 
     }
       }      
     }
     else if(tTopo->tobLayer(detId_)==3){ //3rd layer
       for( unsigned int i=1; i < 55 ;i++){
     if(tTopo->tobRod(detId_)==i){  
       layerStringName << "TOB_L3_Rod_" << i;
       layerStringId_ = 50300+i; 
     }
       }      
     }
     else if(tTopo->tobLayer(detId_)==4){ //4th layer
       for( unsigned int i=1; i < 61 ;i++){
     if(tTopo->tobRod(detId_)==i){  
       layerStringName << "TOB_L4_Rod_" << i;
       layerStringId_ = 50400+i; 
     }
       }      
     }
     else if(tTopo->tobLayer(detId_)==5){ //5th layer
       for( unsigned int i=1; i < 67 ;i++){
     if(tTopo->tobRod(detId_)==i){  
       layerStringName << "TOB_L5_Rod_" << i;
       layerStringId_ = 50500+i; 
     }
       }      
     }
     else if(tTopo->tobLayer(detId_)==6){ //6st layer
       for( unsigned int i=1; i < 75 ;i++){
     if(tTopo->tobRod(detId_)==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_,
		const TrackerTopology *	tTopo
	)

protected

Definition at line 383 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(), SiStripBackPlaneCorrectionDQM::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_,tTopo).second);
   }
   else if(CondObj_name_ == "bpcorrection" && SummaryOnStringLevel_On_ ){
     SummaryMEsMap_iter = SummaryMEsMap_.find(getStringNameAndId(detId_,tTopo).second);
   }
   else {
     SummaryMEsMap_iter = SummaryMEsMap_.find(getLayerNameAndId(detId_,tTopo).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_ == "bpcorrection"  || 
       CondObj_name_ == "lorentzangle") ) {
     if(hPSet_.getParameter<bool>("FillSummaryProfileAtLayerLevel")) 
       if (!CondObj_ME.SummaryOfProfileDistr) { bookSummaryProfileMEs(CondObj_ME,detId_,tTopo); }
   }    
     
   // --> currently only genuine cumul LA
   if(   (CondObj_fillId_ =="ProfileAndCumul" || CondObj_fillId_ =="onlyCumul" ) &&
     (
      CondObj_name_ == "lorentzangle" ||  
      CondObj_name_ == "bpcorrection" ||  
      CondObj_name_ == "noise")  ) {
     if(hPSet_.getParameter<bool>("FillCumulativeSummaryAtLayerLevel"))
       if (!CondObj_ME.SummaryOfCumulDistr) { bookSummaryCumulMEs(CondObj_ME,detId_,tTopo); } 
   } 
                           
   // --> 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_,tTopo); } 
     
   } 
                           
   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_,tTopo).second,CondObj_ME) );
   }
   else if(CondObj_name_ == "bpcorrection" && 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_,tTopo).second,CondObj_ME) );
   }
   else {
     SummaryMEsMap_.insert( std::make_pair(getLayerNameAndId(detId_,tTopo).second,CondObj_ME) );
   }
 
 }

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

private

Definition at line 1197 of file SiStripBaseCondObjDQM.cc.

References python.tagInventory::entries, edm::ParameterSet::getParameter(), hPSet_, i, j, bookConverter::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("SiStripBaseCondObjDQM") << "Resetting TkMap maxValue from " << maxValue << " to " << imax;
       maxValue=imax;
     }
     //reset minValue;
     if(minValue>imin){
         edm::LogInfo("SiStripBaseCondObjDQM") << "Resetting TkMap minValue from " << minValue << " to " << imin;
       minValue=imin;
     }
   }
   
   tkMap->save(false, minValue, maxValue, TkMapname.c_str(),4500,2400);
   tkMap->setPalette(1); tkMap->showPalette(true);
   
 }

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

Definition at line 168 of file SiStripBaseCondObjDQM.cc.

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

Referenced by analysis().

                                                                       {
   
   edm::LogInfo("SiStripBaseCondObjDQM") << "[SiStripBaseCondObjDQM::selectModules] input detIds_: " << detIds_.size() << std::endl;
   
   if( fPSet_.getParameter<bool>("restrictModules")){
 
     std::vector<DetIdSelector> included_subdetsels;
     std::vector<std::string> included_subdets = fPSet_.getParameter<std::vector<std::string> >("ModulesToBeIncluded_DetIdSelector");
     for(std::vector<std::string>::const_iterator wsdps = included_subdets.begin();wsdps!=included_subdets.end();++wsdps) {
       included_subdetsels.push_back(DetIdSelector(*wsdps));
     }
     
     std::vector<uint32_t> modulesToBeIncluded;
     for(std::vector<uint32_t>::const_iterator detid=detIds_.begin();detid!=detIds_.end();++detid) {
       for(std::vector<DetIdSelector>::const_iterator detidsel=included_subdetsels.begin();detidsel!=included_subdetsels.end();++detidsel) {
     if(detidsel->isSelected(*detid)) {
       modulesToBeIncluded.push_back(*detid);
       break;
       //      std::cout << "detId: " << *detid << " is selected" << std::endl;
     }
       }
     }
     
     // -----
     // *** exclude modules ***
     std::vector<DetIdSelector> excluded_subdetsels;
     std::vector<std::string> excluded_subdets = fPSet_.getParameter<std::vector<std::string> >("ModulesToBeExcluded_DetIdSelector");
     for(std::vector<std::string>::const_iterator wsdps = excluded_subdets.begin();wsdps!=excluded_subdets.end();++wsdps) {
       excluded_subdetsels.push_back(DetIdSelector(*wsdps));
     }
     
     std::vector<uint32_t> modulesToBeExcluded;
     for(std::vector<uint32_t>::const_iterator detid=detIds_.begin();detid!=detIds_.end();++detid) {
       for(std::vector<DetIdSelector>::const_iterator detidsel=excluded_subdetsels.begin();detidsel!=excluded_subdetsels.end();++detidsel) {
     if(detidsel->isSelected(*detid)) {
       modulesToBeExcluded.push_back(*detid);
       break;
     }
       }
     }
     
     ModulesToBeExcluded_     = fPSet_.getParameter< std::vector<unsigned int> >("ModulesToBeExcluded");
     ModulesToBeIncluded_     = fPSet_.getParameter< std::vector<unsigned int> >("ModulesToBeIncluded");
     SubDetectorsToBeExcluded_= fPSet_.getParameter< std::vector<std::string> >("SubDetectorsToBeExcluded");  
     
     // vectors to be sorted otherwise the intersection is non computed properly
     
     std::sort(ModulesToBeExcluded_.begin(),ModulesToBeExcluded_.end());
     std::sort(ModulesToBeIncluded_.begin(),ModulesToBeIncluded_.end());
     
     if (modulesToBeExcluded.size()==0  && modulesToBeIncluded.size()==0 &&
     ModulesToBeExcluded_.size()==0 && ModulesToBeIncluded_.size()==0 )
       edm::LogWarning("SiStripBaseCondObjDQM") 
     << "[SiStripBaseCondObjDQM::selectModules] PLEASE CHECK : no modules to be exclude/included in your cfg"
     << std::endl; 
     
     
     modulesToBeIncluded.insert(modulesToBeIncluded.end(),ModulesToBeIncluded_.begin(),ModulesToBeIncluded_.end());
     edm::LogInfo("SiStripBaseCondObjDQM") << "[SiStripBaseCondObjDQM::selectModules] modulesToBeIncluded: " << modulesToBeIncluded.size() << std::endl;
     modulesToBeExcluded.insert(modulesToBeExcluded.end(),ModulesToBeExcluded_.begin(),ModulesToBeExcluded_.end());
     edm::LogInfo("SiStripBaseCondObjDQM") << "[SiStripBaseCondObjDQM::selectModules] modulesToBeExcluded: " << modulesToBeExcluded.size() << std::endl;
     
     // apply modules selection
     if( modulesToBeIncluded.size()>0 ){
       std::vector<uint32_t> tmp;
       // The intersection of two sets is formed only by the elements that are present in both sets
       set_intersection( detIds_.begin(), detIds_.end(),  
             modulesToBeIncluded.begin(), modulesToBeIncluded.end(),
             inserter(tmp,tmp.begin()));
       swap(detIds_,tmp);
     }
     
     std::sort(detIds_.begin(),detIds_.end());
     if(modulesToBeExcluded.size()>0) {
       for( std::vector<uint32_t>::const_iterator mod = modulesToBeExcluded.begin(); 
        mod != modulesToBeExcluded.end(); mod++){
     
     std::vector<uint32_t>::iterator detid=std::lower_bound(detIds_.begin(),detIds_.end(),*mod);
     if (detid!=detIds_.end())
       detIds_.erase(detid);
     detid--;
       }
     }
     
     // -----
     // *** 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_
     }
     
   }
   edm::LogInfo("SiStripBaseCondObjDQM") << "[SiStripBaseCondObjDQM::selectModules] output detIds_: " << detIds_.size() << std::endl;
   
   // -----
   // *** 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