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_, const std::vector< uint32_t > &detIdsOnDemand)
void	analysisOnDemand (const edm::EventSetup &eSetup_, std::string requestedSubDetector, uint32_t requestedSide, uint32_t requestedLayer)
void	analysisOnDemand (const edm::EventSetup &eSetup_, uint32_t detIdOnDemand)
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
Public Types inherited from SiStripBaseCondObjDQM
typedef dqm::legacy::DQMStore	DQMStore
typedef dqm::legacy::MonitorElement	MonitorElement
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 9 of file SiStripLorentzAngleDQM.h.

Constructor & Destructor Documentation

SiStripLorentzAngleDQM()

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

Definition at line 9 of file SiStripLorentzAngleDQM.cc.

    : 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.);
  }
}

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

~SiStripLorentzAngleDQM()

SiStripLorentzAngleDQM::~SiStripLorentzAngleDQM ( )

override

Definition at line 24 of file SiStripLorentzAngleDQM.cc.

{}

Member Function Documentation

fillMEsForDet()

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

inlineoverridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 21 of file SiStripLorentzAngleDQM.h.

{};

fillMEsForLayer()

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

overridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 125 of file SiStripLorentzAngleDQM.cc.

                                                        {
  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
}

References SiStripBaseCondObjDQM::activeDetIds, gather_cfg::cout, SiStripHistoId::createHistoLayer(), dqm::impl::MonitorElement::Fill(), SiStripBaseCondObjDQM::fillTkMap(), dqmdumpme::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 fillSummaryMEs().

fillModMEs()

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

inlineoverridevirtual

Reimplemented from SiStripBaseCondObjDQM.

Definition at line 20 of file SiStripLorentzAngleDQM.h.

{};

fillSummaryMEs()

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

overridevirtual

Reimplemented from SiStripBaseCondObjDQM.

Definition at line 41 of file SiStripLorentzAngleDQM.cc.

                                                                                                              {
  // 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());
      }
    }
  }
}

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

getActiveDetIds()

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

overridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 28 of file SiStripLorentzAngleDQM.cc.

                                                                        {
  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);
  }
}

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

getCache()

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

inlineoverridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 27 of file SiStripLorentzAngleDQM.h.

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

References edm::EventSetup::get().

getConditionObject()

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

inlineoverridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 31 of file SiStripLorentzAngleDQM.h.

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

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

Referenced by getActiveDetIds().

Member Data Documentation

lorentzangleHandle_

edm::ESHandle<SiStripLorentzAngle> SiStripLorentzAngleDQM::lorentzangleHandle_

private

Definition at line 37 of file SiStripLorentzAngleDQM.h.

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