#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 9 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 23 of file SiStripLorentzAngleDQM.cc.

23 {}

Member Function Documentation

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

inlineoverridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 21 of file SiStripLorentzAngleDQM.h.

References fillMEsForLayer(), and fillSummaryMEs().

21 {};

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

overridevirtual

Implements SiStripBaseCondObjDQM.

Definition at line 124 of file SiStripLorentzAngleDQM.cc.

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 20 of file SiStripLorentzAngleDQM.h.

20 {};

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

overridevirtual

Reimplemented from SiStripBaseCondObjDQM.

Definition at line 40 of file SiStripLorentzAngleDQM.cc.

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 27 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 27 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 31 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 37 of file SiStripLorentzAngleDQM.h.

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

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation