SiStripLorentzAngleDQM Class Reference

#include <SiStripLorentzAngleDQM.h>

Inheritance diagram for SiStripLorentzAngleDQM:

Public Member Functions
void	fillMEsForDet (const ModMEs &selModME_, uint32_t selDetId_, const TrackerTopology *tTopo) override
void	fillMEsForLayer (uint32_t selDetId_, const TrackerTopology *tTopo) override
void	fillModMEs (const std::vector< uint32_t > &selectedDetIds, const edm::EventSetup &es) override
void	fillSummaryMEs (const std::vector< uint32_t > &selectedDetIds, const edm::EventSetup &es) override
void	getActiveDetIds (const edm::EventSetup &eSetup) override
unsigned long long	getCache (const edm::EventSetup &eSetup) override
void	getConditionObject (const edm::EventSetup &eSetup) override
	SiStripLorentzAngleDQM (const edm::EventSetup &eSetup, edm::RunNumber_t iRun, edm::ParameterSet const &hPSet, edm::ParameterSet const &fPSet)
	~SiStripLorentzAngleDQM () override
Public Member Functions inherited from SiStripBaseCondObjDQM
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 ()
std::vector< uint32_t >	getCabledModules ()
void	selectModules (std::vector< uint32_t > &detIds_, const TrackerTopology *tTopo)
	SiStripBaseCondObjDQM (const edm::EventSetup &eSetup, edm::RunNumber_t iRun, edm::ParameterSet const &hPSet, edm::ParameterSet const &fPSet)
virtual	~SiStripBaseCondObjDQM ()

Private Attributes
edm::ESHandle< SiStripLorentzAngle >	lorentzangleHandle_

Additional Inherited Members
Protected Member Functions inherited from SiStripBaseCondObjDQM
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 inherited from SiStripBaseCondObjDQM
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_
std::unique_ptr< TkHistoMap >	Tk_HM_
std::unique_ptr< TkHistoMap >	Tk_HM_H
std::unique_ptr< TkHistoMap >	Tk_HM_L
TrackerMap *	tkMap
std::vector< int >	tkMapScaler

Detailed Description

Definition at line 11 of file SiStripLorentzAngleDQM.h.

Constructor & Destructor Documentation

SiStripLorentzAngleDQM::SiStripLorentzAngleDQM	(	const edm::EventSetup &	eSetup,
		edm::RunNumber_t	iRun,
		edm::ParameterSet const &	hPSet,
		edm::ParameterSet const &	fPSet
	)

Definition at line 8 of file SiStripLorentzAngleDQM.cc.

References edm::EventSetup::get(), SiStripBaseCondObjDQM::HistoMaps_On_, edm::ESHandle< T >::product(), and SiStripBaseCondObjDQM::Tk_HM_.

                                                        :SiStripBaseCondObjDQM(eSetup, iRun, hPSet, fPSet){

  // Build the Histo_TkMap:
  if ( HistoMaps_On_ ) {
    edm::ESHandle<TkDetMap> tkDetMapHandle;
    eSetup.get<TrackerTopologyRcd>().get(tkDetMapHandle);
    Tk_HM_ = std::make_unique<TkHistoMap>(tkDetMapHandle.product(), "SiStrip/Histo_Map","LA_TkMap",0.);
  }
}

SiStripLorentzAngleDQM::~SiStripLorentzAngleDQM ( )

override

Definition at line 25 of file SiStripLorentzAngleDQM.cc.

{}

Member Function Documentation

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

inlineoverridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 25 of file SiStripLorentzAngleDQM.h.

References fillMEsForLayer(), and fillSummaryMEs().

{};

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

overridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 133 of file SiStripLorentzAngleDQM.cc.

References SiStripBaseCondObjDQM::activeDetIds, gather_cfg::cout, SiStripHistoId::createHistoLayer(), MonitorElement::Fill(), SiStripBaseCondObjDQM::fillTkMap(), plotBeamSpotDB::first, SiStripBaseCondObjDQM::fPSet_, SiStripBaseCondObjDQM::getLayerNameAndId(), SiStripLorentzAngle::getLorentzAngle(), edm::ParameterSet::getParameter(), SiStripBaseCondObjDQM::GetSameLayerDetId(), SiStripBaseCondObjDQM::getStringNameAndId(), SiStripBaseCondObjDQM::getSummaryMEs(), SiStripSubStructure::getTECDetectors(), SiStripSubStructure::getTIBDetectors(), SiStripSubStructure::getTIDDetectors(), SiStripSubStructure::getTOBDetectors(), SiStripBaseCondObjDQM::hidmanager, SiStripBaseCondObjDQM::HistoMaps_On_, SiStripBaseCondObjDQM::hPSet_, mps_fire::i, lorentzangleHandle_, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), SiStripBaseCondObjDQM::SummaryMEsMap_, SiStripBaseCondObjDQM::ModMEs::SummaryOfCumulDistr, SiStripBaseCondObjDQM::ModMEs::SummaryOfProfileDistr, SiStripBaseCondObjDQM::SummaryOnStringLevel_On_, StripSubdetector::TEC, StripSubdetector::TIB, TrackerTopology::tibIsExternalString(), TrackerTopology::tibIsInternalString(), TrackerTopology::tibLayer(), TrackerTopology::tibString(), StripSubdetector::TID, SiStripBaseCondObjDQM::Tk_HM_, StripSubdetector::TOB, TrackerTopology::tobLayer(), and TrackerTopology::tobRod().

Referenced by fillMEsForDet(), and fillSummaryMEs().

                                                                                               {

  SiStripHistoId hidmanager;

      
  std::string hSummaryOfProfile_description;
  hSummaryOfProfile_description  = hPSet_.getParameter<std::string>("SummaryOfProfile_description");

  std::string hSummary_name; 

  int subDetId_ = DetId(selDetId_).subdetId();
  
  if( subDetId_<3 || subDetId_>6 ){ 
    edm::LogError("SiStripLorentzAngleDQM")
      << "[SiStripLorentzAngleDQM::fillMEsForLayer] WRONG INPUT : no such subdetector type : "
      << subDetId_ << " no folder set!" 
      << std::endl;
    return;
  }

  std::vector<uint32_t> sameLayerDetIds_;
  sameLayerDetIds_.clear();

  if (SummaryOnStringLevel_On_) {  //FILLING FOR STRING LEVEL

    hSummary_name = hidmanager.createHistoLayer(hSummaryOfProfile_description, "layer", getStringNameAndId(selDetId_,tTopo).first, "") ;
    std::map<uint32_t, ModMEs>::iterator selMEsMapIter_ = SummaryMEsMap_.find(getStringNameAndId(selDetId_,tTopo).second);
    
    ModMEs selME_;
    if ( selMEsMapIter_ != SummaryMEsMap_.end())
      selME_ =selMEsMapIter_->second;

    getSummaryMEs(selME_,selDetId_,tTopo);
  
    // -----                       
    sameLayerDetIds_.clear();

    switch ( DetId(selDetId_).subdetId() ) {
      case StripSubdetector::TIB:
        if(tTopo->tibIsInternalString(selDetId_)){
          SiStripSubStructure::getTIBDetectors(activeDetIds, sameLayerDetIds_, tTopo, tTopo->tibLayer(selDetId_),0,1,tTopo->tibString(selDetId_));
        }
        if(tTopo->tibIsExternalString(selDetId_)){
          SiStripSubStructure::getTIBDetectors(activeDetIds, sameLayerDetIds_, tTopo, tTopo->tibLayer(selDetId_),0,2,tTopo->tibString(selDetId_));
        }
        break;
      case StripSubdetector::TID:
        SiStripSubStructure::getTIDDetectors(activeDetIds, sameLayerDetIds_, tTopo, 0,0,0,0);
        break;
      case StripSubdetector::TOB:
        SiStripSubStructure::getTOBDetectors(activeDetIds, sameLayerDetIds_, tTopo, tTopo->tobLayer(selDetId_),0,tTopo->tobRod(selDetId_));
        break;
      case StripSubdetector::TEC:
        SiStripSubStructure::getTECDetectors(activeDetIds, sameLayerDetIds_, tTopo, 0,0,0,0,0,0);
        break;
    }

    // -----

    for(unsigned int i=0;i< sameLayerDetIds_.size(); i++){
    selME_.SummaryOfProfileDistr->Fill(i+1,lorentzangleHandle_->getLorentzAngle(sameLayerDetIds_[i]));

    // Fill the Histo_TkMap+TkMap with the LA:
        if(HistoMaps_On_ ) Tk_HM_->fill(sameLayerDetIds_[i], lorentzangleHandle_->getLorentzAngle(sameLayerDetIds_[i]));

      std::cout<<sameLayerDetIds_[i]<<"\t"<<lorentzangleHandle_->getLorentzAngle(sameLayerDetIds_[i])<<std::endl;

    if(fPSet_.getParameter<bool>("TkMap_On") || hPSet_.getParameter<bool>("TkMap_On")){
      fillTkMap(sameLayerDetIds_[i], lorentzangleHandle_->getLorentzAngle(sameLayerDetIds_[i]));
    }
    
    } 
    
    std::string hSummaryOfCumul_description;
    hSummaryOfCumul_description  = hPSet_.getParameter<std::string>("SummaryOfCumul_description");
    
    std::string hSummaryOfCumul_name; 
    
    if( subDetId_<3 || subDetId_>6 ){ 
      edm::LogError("SiStripLorentzAngleDQM")
    << "[SiStripLorentzAngleDQM::fillMEsForLayer] WRONG INPUT : no such subdetector type : "
    << subDetId_ << " no folder set!" 
    << std::endl;
      return;
    }
    
    hSummaryOfCumul_name = hidmanager.createHistoLayer(hSummaryOfCumul_description, "layer", getStringNameAndId(selDetId_,tTopo).first, "") ;
    
    for(unsigned int i=0;i< sameLayerDetIds_.size(); i++){
    selME_.SummaryOfCumulDistr->Fill(lorentzangleHandle_->getLorentzAngle(sameLayerDetIds_[i]));
    } 
  } //FILLING FOR STRING LEVEL
  
  
  else { //FILLING FOR LAYER LEVEL

    std::map<uint32_t, ModMEs>::iterator selMEsMapIter_ = SummaryMEsMap_.find(getLayerNameAndId(selDetId_,tTopo).second);
    
    ModMEs selME_;
    if ( selMEsMapIter_ != SummaryMEsMap_.end())
      selME_ =selMEsMapIter_->second;
    
    getSummaryMEs(selME_,selDetId_,tTopo);
    
    if(hPSet_.getParameter<bool>("FillSummaryProfileAtLayerLevel")){

      hSummary_name = hidmanager.createHistoLayer(hSummaryOfProfile_description, "layer", getLayerNameAndId(selDetId_,tTopo).first, "") ;    
    
      // -----                         
      sameLayerDetIds_.clear();

      sameLayerDetIds_=GetSameLayerDetId(activeDetIds,selDetId_,tTopo);     

      for(unsigned int i=0;i< sameLayerDetIds_.size(); i++){
      selME_.SummaryOfProfileDistr->Fill(i+1,lorentzangleHandle_->getLorentzAngle(sameLayerDetIds_[i]));

      // Fill the Histo_TkMap with LA:
      if(HistoMaps_On_ ) Tk_HM_->fill(sameLayerDetIds_[i], lorentzangleHandle_->getLorentzAngle(sameLayerDetIds_[i]));

    if(fPSet_.getParameter<bool>("TkMap_On") || hPSet_.getParameter<bool>("TkMap_On")){
      fillTkMap(sameLayerDetIds_[i], lorentzangleHandle_->getLorentzAngle(sameLayerDetIds_[i]));

    }

      } 
    }//if Fill ...

    if(hPSet_.getParameter<bool>("FillCumulativeSummaryAtLayerLevel")){

      std::string hSummaryOfCumul_description;
      hSummaryOfCumul_description  = hPSet_.getParameter<std::string>("SummaryOfCumul_description");
    
      std::string hSummaryOfCumul_name; 
    
      if( subDetId_<3 || subDetId_>6 ){ 
    edm::LogError("SiStripLorentzAngleDQM")
      << "[SiStripLorentzAngleDQM::fillMEsForLayer] WRONG INPUT : no such subdetector type : "
      << subDetId_ << " no folder set!" 
      << std::endl;
    return;
      }
    
      hSummaryOfCumul_name = hidmanager.createHistoLayer(hSummaryOfCumul_description, "layer", getLayerNameAndId(selDetId_,tTopo).first, "") ;
    
      for(unsigned int i=0;i< sameLayerDetIds_.size(); i++){
      selME_.SummaryOfCumulDistr->Fill(lorentzangleHandle_->getLorentzAngle(sameLayerDetIds_[i]));
      }
    }//if Fill ... 
  } //FILLING FOR LAYER LEVEL
  
}  

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

inlineoverridevirtual

Reimplemented from SiStripBaseCondObjDQM.

Definition at line 24 of file SiStripLorentzAngleDQM.h.

{};

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

overridevirtual

Reimplemented from SiStripBaseCondObjDQM.

Definition at line 48 of file SiStripLorentzAngleDQM.cc.

References alignmentValidation::c1, fillMEsForLayer(), SiStripBaseCondObjDQM::fPSet_, edm::EventSetup::get(), SiStripBaseCondObjDQM::getLayerNameAndId(), edm::ParameterSet::getParameter(), SiStripBaseCondObjDQM::getStringNameAndId(), MonitorElement::getTH1(), MonitorElement::getTProfile(), SiStripBaseCondObjDQM::hPSet_, mps_fire::i, dataset::name, edm::ESHandle< T >::product(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), SiStripBaseCondObjDQM::SummaryMEsMap_, SiStripBaseCondObjDQM::ModMEs::SummaryOfCumulDistr, SiStripBaseCondObjDQM::ModMEs::SummaryOfProfileDistr, SiStripBaseCondObjDQM::SummaryOnLayerLevel_On_, and SiStripBaseCondObjDQM::SummaryOnStringLevel_On_.

Referenced by fillMEsForDet().

                                                                                                              {

  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle;
  es.get<TrackerTopologyRcd>().get(tTopoHandle);
  const TrackerTopology* const tTopo = tTopoHandle.product();

  // -----
  // LA on layer-level : fill at once all detIds belonging to same layer when encountering first detID in the layer 

  bool fillNext = true; 
  for(unsigned int i=0;i<selectedDetIds.size();i++){                        
    int subDetId_ = DetId(selectedDetIds[i]).subdetId();
    if( subDetId_<3 ||subDetId_>6 ){ 
      edm::LogError("SiStripLorentzAngle")
    << "[SiStripLorentzAngle::fillSummaryMEs] WRONG INPUT : no such subdetector type : "
    << subDetId_ << " and detId " << selectedDetIds[i] << " therefore no filling!" 
    << std::endl;
    }    
    else if (SummaryOnLayerLevel_On_) {    
      if( fillNext) { fillMEsForLayer(/*SummaryMEsMap_,*/ selectedDetIds[i],tTopo);} 
      if( getLayerNameAndId(selectedDetIds[i+1],tTopo)==getLayerNameAndId(selectedDetIds[i],tTopo)){ fillNext=false;}
      else { fillNext=true;}
    } 
    else if (SummaryOnStringLevel_On_) {
      if( fillNext) { fillMEsForLayer(/*SummaryMEsMap_,*/ selectedDetIds[i],tTopo);} 
      if( getStringNameAndId(selectedDetIds[i+1],tTopo)==getStringNameAndId(selectedDetIds[i],tTopo)){ fillNext=false;}
      else { fillNext=true;}
    } 
  }

  for (std::map<uint32_t, ModMEs>::iterator iter=SummaryMEsMap_.begin(); iter!=SummaryMEsMap_.end(); iter++){

    ModMEs selME;
    selME = iter->second;

    if(SummaryOnStringLevel_On_){

      if (fPSet_.getParameter<bool>("OutputSummaryProfileAtLayerLevelAsImage")){

    TCanvas c1("c1");
    selME.SummaryOfProfileDistr->getTProfile()->Draw();
    std::string name (selME.SummaryOfProfileDistr->getTProfile()->GetTitle());
    name+=".png";
    c1.Print(name.c_str());
      }

      if (fPSet_.getParameter<bool>("OutputCumulativeSummaryAtLayerLevelAsImage")){

    TCanvas c2("c2");
    selME.SummaryOfCumulDistr->getTH1()->Draw();
    std::string name2 (selME.SummaryOfCumulDistr->getTH1()->GetTitle());
    name2+=".png";
    c2.Print(name2.c_str());
      }

    }
    else{
      if(hPSet_.getParameter<bool>("FillSummaryProfileAtLayerLevel") && fPSet_.getParameter<bool>("OutputSummaryProfileAtLayerLevelAsImage")){

    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>("FillCumulativeSummaryAtLayerLevel") && fPSet_.getParameter<bool>("OutputCumulativeSummaryAtLayerLevelAsImage")){

    TCanvas c1("c1");
    selME.SummaryOfCumulDistr->getTH1()->Draw();
    std::string name (selME.SummaryOfCumulDistr->getTH1()->GetTitle());
    name+=".png";
    c1.Print(name.c_str());
      }
    }

  }

}

void SiStripLorentzAngleDQM::getActiveDetIds ( const edm::EventSetup &  eSetup )

overridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 30 of file SiStripLorentzAngleDQM.cc.

References SiStripBaseCondObjDQM::activeDetIds, getConditionObject(), SiStripLorentzAngle::getLorentzAngles(), and lorentzangleHandle_.

                                                                        {
  
  getConditionObject(eSetup);

  std::map<uint32_t,float>::const_iterator LAMapIter_;
  std::map<uint32_t,float> LAMap_ = lorentzangleHandle_->getLorentzAngles();

  for( LAMapIter_ = LAMap_.begin();
       LAMapIter_!= LAMap_.end(); LAMapIter_++){

    activeDetIds.push_back((*LAMapIter_).first);
  }

}

unsigned long long SiStripLorentzAngleDQM::getCache ( const edm::EventSetup &  eSetup )

inlineoverridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 30 of file SiStripLorentzAngleDQM.h.

References edm::EventSetup::get().

{ return eSetup.get<SiStripLorentzAngleRcd>().cacheIdentifier();}

void SiStripLorentzAngleDQM::getConditionObject ( const edm::EventSetup &  eSetup )

inlineoverridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 32 of file SiStripLorentzAngleDQM.h.

References SiStripBaseCondObjDQM::cacheID_current, SiStripBaseCondObjDQM::cacheID_memory, edm::EventSetup::get(), and lorentzangleHandle_.

Referenced by getActiveDetIds().

                                                                {
    eSetup.get<SiStripLorentzAngleRcd>().get(lorentzangleHandle_);
    cacheID_memory = cacheID_current;
  }

Member Data Documentation

edm::ESHandle<SiStripLorentzAngle> SiStripLorentzAngleDQM::lorentzangleHandle_

private

Definition at line 38 of file SiStripLorentzAngleDQM.h.

Referenced by fillMEsForLayer(), getActiveDetIds(), and getConditionObject().