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 (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, 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 ( ) [inline, virtual]

Definition at line 53 of file SiStripBaseCondObjDQM.h.

{};

Member Function Documentation

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

Definition at line 51 of file SiStripBaseCondObjDQM.cc.

References activeDetIds, cacheID_current, cacheID_memory, lut2db_cfg::filename, fillModMEs(), fillSummaryMEs(), fPSet_, getActiveDetIds(), SiStripDetInfoFileReader::getAllDetIds(), getCache(), getConditionObject(), edm::ParameterSet::getParameter(), hPSet_, edm::EventSetup::iovSyncValue(), maxValue, minValue, Mod_On_, reader, 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_,
		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::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::bookCumulMEs	(	SiStripBaseCondObjDQM::ModMEs &	CondObj_ME,
		const uint32_t &	detId_
	)		`[private]`

Definition at line 466 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 417 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 733 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 799 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, 
                                                   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_
	)		`[private]`

Definition at line 507 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(), python::rootplot::argparse::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, 
                                                            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"){

    // 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 1149 of file SiStripBaseCondObjDQM.cc.

References tkMap.

Referenced by SiStripBaseCondObjDQM().

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

void SiStripBaseCondObjDQM::end ( void ) [virtual]

Definition at line 1205 of file SiStripBaseCondObjDQM.cc.

Referenced by SiStripMonitorCondData::endRun().

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

virtual void SiStripBaseCondObjDQM::fillMEsForDet	(	ModMEs	selModME_,
		uint32_t	selDetId_
	)		`[protected, pure virtual]`

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

Referenced by fillModMEs().

virtual void SiStripBaseCondObjDQM::fillMEsForLayer ( uint32_t selDetId_ ) [protected, pure virtual]

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

Referenced by fillSummaryMEs().

void SiStripBaseCondObjDQM::fillModMEs ( const std::vector< uint32_t > & selectedDetIds ) [protected, virtual]

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

Definition at line 1212 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 ) [protected, virtual]

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

Definition at line 1222 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 1154 of file SiStripBaseCondObjDQM.cc.

References TrackerMap::fill(), and tkMap.

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

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

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

Implemented in SiStripApvGainsDQM, SiStripCablingDQM, SiStripLorentzAngleDQM, SiStripNoisesDQM, SiStripPedestalsDQM, SiStripQualityDQM, 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 SiStripApvGainsDQM, SiStripCablingDQM, SiStripLorentzAngleDQM, SiStripNoisesDQM, SiStripPedestalsDQM, SiStripQualityDQM, and SiStripThresholdDQM.

Referenced by analysis(), and analysisOnDemand().

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

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

Referenced by analysis(), and analysisOnDemand().

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

Definition at line 898 of file SiStripBaseCondObjDQM.cc.

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

Referenced by bookSummaryCumulMEs(), bookSummaryMEs(), bookSummaryProfileMEs(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripQualityDQM::fillMEsForLayer(), SiStripApvGainsDQM::fillMEsForLayer(), SiStripThresholdDQM::fillMEsForLayer(), SiStripPedestalsDQM::fillMEsForLayer(), SiStripNoisesDQM::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 314 of file SiStripBaseCondObjDQM.cc.

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

Referenced by SiStripThresholdDQM::fillMEsForDet(), SiStripQualityDQM::fillMEsForDet(), SiStripPedestalsDQM::fillMEsForDet(), SiStripNoisesDQM::fillMEsForDet(), and SiStripApvGainsDQM::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 1121 of file SiStripBaseCondObjDQM.cc.

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

Referenced by bookSummaryMEs(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripApvGainsDQM::fillMEsForLayer(), SiStripQualityDQM::fillMEsForLayer(), SiStripThresholdDQM::fillMEsForLayer(), SiStripNoisesDQM::fillMEsForLayer(), and SiStripPedestalsDQM::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 987 of file SiStripBaseCondObjDQM.cc.

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

Referenced by bookSummaryCumulMEs(), bookSummaryProfileMEs(), SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripNoisesDQM::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 355 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 SiStripLorentzAngleDQM::fillMEsForLayer(), SiStripApvGainsDQM::fillMEsForLayer(), SiStripQualityDQM::fillMEsForLayer(), SiStripThresholdDQM::fillMEsForLayer(), SiStripNoisesDQM::fillMEsForLayer(), and SiStripPedestalsDQM::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 1159 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(),4500,2400);
  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");  

  // vectors to be sorted otherwise the intersection is non computed properly

  std::sort(ModulesToBeExcluded_.begin(),ModulesToBeExcluded_.end());
  std::sort(ModulesToBeIncluded_.begin(),ModulesToBeIncluded_.end());

  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