SiStripPopConHistoryDQM.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"

#include "DQMOffline/CalibTracker/plugins/SiStripPopConHistoryDQMBase.h"

class SiStripPopConHistoryDQM : public SiStripPopConHistoryDQMBase
{
public:
  explicit SiStripPopConHistoryDQM(const edm::ParameterSet& pset)
    : SiStripPopConHistoryDQMBase(pset)
  {}

  void initES(const edm::EventSetup&) override;

  ~SiStripPopConHistoryDQM() override;
private:
  uint32_t returnDetComponent(const MonitorElement* ME) const override;
  bool setDBLabelsForUser(const std::string& keyName, std::vector<std::string>& userDBContent, const std::string& quantity) const override;
  bool setDBValuesForUser(const MonitorElement* me, HDQMSummary::InputVector& values, const std::string& quantity) const override;
private:
  const TrackerTopology* trackerTopo_;
};

SiStripPopConHistoryDQM::~SiStripPopConHistoryDQM() {}

void SiStripPopConHistoryDQM::initES(const edm::EventSetup& setup)
{
  edm::ESHandle<TrackerTopology> tTopo;
  setup.get<TrackerTopologyRcd>().get(tTopo);
  trackerTopo_ = tTopo.product();
}

uint32_t SiStripPopConHistoryDQM::returnDetComponent(const MonitorElement* ME) const
{
  LogTrace("SiStripHistoryDQMService") <<  "[SiStripHistoryDQMService::returnDetComponent]";
  const std::string& str{ME->getName()};
  const size_t __key_length__=7;
  const size_t __detid_length__=9;

  uint32_t layer=0,side=0;

  if ( str.find("__det__") != std::string::npos ) {
    return atoi(str.substr(str.find("__det__")+__key_length__,__detid_length__).c_str());
  }
  //TIB
  else if ( str.find("TIB") != std::string::npos ) {
    if ( str.find("layer")!= std::string::npos )
      layer = atoi(str.substr(str.find("layer__")+__key_length__,1).c_str());
    return trackerTopo_->tibDetId(layer,0,0,0,0,0).rawId();
  }
  //TOB
  else if ( str.find("TOB") != std::string::npos ) {
    if ( str.find("layer") != std::string::npos )
      layer = atoi(str.substr(str.find("layer__")+__key_length__,1).c_str());
    return trackerTopo_->tobDetId(layer,0,0,0,0).rawId();
  }
  //TID
  else if ( str.find("TID") != std::string::npos ) {
    if ( str.find("side") != std::string::npos ) {
      side = atoi(str.substr(str.find("_side__")+__key_length__,1).c_str());
      if ( str.find("wheel") != std::string::npos ) {
    layer = atoi(str.substr(str.find("wheel__")+__key_length__,1).c_str());
      }
    }
    return trackerTopo_->tidDetId(side,layer,0,0,0,0).rawId();
  }
  //TEC
  else if ( str.find("TEC") != std::string::npos ) {
    if ( str.find("side") != std::string::npos ) {
      side = atoi(str.substr(str.find("_side__")+__key_length__,1).c_str());
      if ( str.find("wheel") != std::string::npos ) {
    layer = atoi(str.substr(str.find("wheel__")+__key_length__,1).c_str());
      }
    }
    return trackerTopo_->tecDetId(side,layer,0,0,0,0,0).rawId();
  }
  else
    return DetId(DetId::Tracker,0).rawId(); //Full Tracker
}

//Example on how to define an user function for the statistic extraction
bool SiStripPopConHistoryDQM::setDBLabelsForUser(const std::string& keyName, std::vector<std::string>& userDBContent, const std::string& quantity) const
{
  if (quantity == "user_2DYmean") {
    userDBContent.push_back(keyName+sep()+std::string("yMean"));
    userDBContent.push_back(keyName+sep()+std::string("yError"));
  } else {
    edm::LogError("SiStripHistoryDQMService") << "ERROR: quantity does not exist in SiStripHistoryDQMService::setDBValuesForUser(): " << quantity;
    return false;
  }
  return true;
}

bool SiStripPopConHistoryDQM::setDBValuesForUser(const MonitorElement* me, HDQMSummary::InputVector& values, const std::string& quantity) const
{
  if (quantity == "user_2DYmean") {
    TH2F* Hist = (TH2F*) me->getTH2F();
    values.push_back( Hist->GetMean(2) );
    values.push_back( Hist->GetRMS(2) );
  } else {
    edm::LogError("SiStripHistoryDQMService") << "ERROR: quantity does not exist in SiStripHistoryDQMService::setDBValuesForUser(): " << quantity;
    return false;
  }
  return true;
}

#include "DQMOffline/CalibTracker/plugins/SiStripPopConDQMEDHarvester.h"
#include "FWCore/Framework/interface/MakerMacros.h"
using SiStripDQMHistoryPopCon = SiStripPopConDQMEDHarvester<SiStripPopConHistoryDQM>;
DEFINE_FWK_MODULE(SiStripDQMHistoryPopCon);