ESPedestalClient Class Reference

#include <ESPedestalClient.h>

Inheritance diagram for ESPedestalClient:

Public Member Functions
void	endJobAnalyze (DQMStore::IGetter &) override
	ESPedestalClient (const edm::ParameterSet &)
	~ESPedestalClient () override
Public Member Functions inherited from ESClient
virtual void	endLumiAnalyze (DQMStore::IGetter &)
	ESClient (edm::ParameterSet const &)
template<typename T >
T *	getHisto (MonitorElement *, bool=false, T *=0) const
void	setup (DQMStore::IBooker &)
virtual	~ESClient ()

Private Member Functions
void	book (DQMStore::IBooker &) override

Private Attributes
TF1 *	fg_
bool	fitPedestal_
MonitorElement *	hPed_ [2][2][40][40]
MonitorElement *	hTotN_ [2][2][40][40]
std::vector< int >	senP_
std::vector< int >	senX_
std::vector< int >	senY_
std::vector< int >	senZ_

Additional Inherited Members
Protected Attributes inherited from ESClient
bool	cloneME_
bool	debug_
bool	initialized_
std::string	prefixME_
bool	verbose_

Detailed Description

Definition at line 14 of file ESPedestalClient.h.

Constructor & Destructor Documentation

ESPedestalClient::ESPedestalClient

(

const edm::ParameterSet &

ps

)

Definition at line 13 of file ESPedestalClient.cc.

References gather_cfg::cout, FrontierConditions_GlobalTag_cff::file, edm::FileInPath::fullPath(), edm::ParameterSet::getUntrackedParameter(), hPed_, hTotN_, mps_fire::i, gen::k, funct::m, senP_, senX_, senY_, senZ_, and AlCaHLTBitMon_QueryRunRegistry::string.

    : ESClient(ps),
      fitPedestal_(ps.getUntrackedParameter<bool>("fitPedestal")),
      fg_(new TF1("fg", "gaus")),
      senZ_(),
      senP_(),
      senX_(),
      senY_() {
  for (int i = 0; i < 2; i++)
    for (int j = 0; j < 2; j++)
      for (int k = 0; k < 40; k++)
        for (int m = 0; m < 40; m++) {
          hPed_[i][j][k][m] = nullptr;
          hTotN_[i][j][k][m] = nullptr;
        }

  std::string lutPath(ps.getUntrackedParameter<edm::FileInPath>("LookupTable").fullPath());

  // read in look-up table
  int nLines, iz, ip, ix, iy, fed, kchip, pace, bundle, fiber, optorx;
  ifstream file(lutPath);

  if (file.is_open()) {
    file >> nLines;

    for (int i = 0; i < nLines; ++i) {
      file >> iz >> ip >> ix >> iy >> fed >> kchip >> pace >> bundle >> fiber >> optorx;

      senZ_.push_back(iz);
      senP_.push_back(ip);
      senX_.push_back(ix);
      senY_.push_back(iy);
    }

    file.close();
  } else {
    cout << "ESPedestalClient : Look up table file can not be found in " << lutPath << endl;
  }
}

ESPedestalClient::~ESPedestalClient ( )

override

Definition at line 53 of file ESPedestalClient.cc.

References fg_.

{ delete fg_; }

Member Function Documentation

void ESPedestalClient::book ( DQMStore::IBooker &  _ibooker )

overrideprivatevirtual

Reimplemented from ESClient.

Definition at line 112 of file ESPedestalClient.cc.

References DQMStore::IBooker::book1D(), hPed_, hTotN_, mps_fire::i, ESClient::prefixME_, senP_, senX_, senY_, senZ_, and DQMStore::IBooker::setCurrentFolder().

                                                     {
  // define histograms
  _ibooker.setCurrentFolder(prefixME_ + "/ESPedestalClient");

  char hname[300];
  for (unsigned i = 0; i < senZ_.size(); ++i) {
    int iz = (senZ_[i] == 1) ? 0 : 1;

    sprintf(hname, "Ped Z %d P %d X %d Y %d", senZ_[i], senP_[i], senX_[i], senY_[i]);
    hPed_[iz][senP_[i] - 1][senX_[i] - 1][senY_[i] - 1] = _ibooker.book1D(hname, hname, 32, 0, 32);

    sprintf(hname, "Total Noise Z %d P %d X %d Y %d", senZ_[i], senP_[i], senX_[i], senY_[i]);
    hTotN_[iz][senP_[i] - 1][senX_[i] - 1][senY_[i] - 1] = _ibooker.book1D(hname, hname, 32, 0, 32);
  }
}

void ESPedestalClient::endJobAnalyze ( DQMStore::IGetter &  _igetter )

overridevirtual

Reimplemented from ESClient.

Definition at line 55 of file ESPedestalClient.cc.

References gather_cfg::cout, ESClient::debug_, compare_using_db::dirname, fg_, fitPedestal_, DQMStore::IGetter::get(), MonitorElement::getTH1F(), hPed_, hTotN_, mps_fire::i, ESClient::prefixME_, senP_, senX_, senY_, senZ_, MonitorElement::setBinContent(), and ESClient::verbose_.

                                                              {
  if (debug_)
    cout << "ESPedestalClient: endJob" << endl;

  // Preform pedestal fit
  char hname[300];
  int iz = 0;
  if (fitPedestal_) {
    if (verbose_)
      cout << "ESPedestalClient: Fit Pedestal" << endl;

    for (unsigned i = 0; i < senZ_.size(); ++i) {
      iz = (senZ_[i] == 1) ? 0 : 1;

      for (int is = 0; is < 32; ++is) {
        string dirname = prefixME_ + "/ESPedestalTask/";
        sprintf(hname, "ADC Z %d P %d X %d Y %d Str %d", senZ_[i], senP_[i], senX_[i], senY_[i], is + 1);
        MonitorElement *meFit = _igetter.get(dirname + hname);

        if (meFit == nullptr)
          continue;
        TH1F *rootHisto = meFit->getTH1F();

        rootHisto->Fit(fg_, "Q", "", 500, 1800);
        rootHisto->Fit(fg_,
                       "RQ",
                       "",
                       fg_->GetParameter(1) - 2. * fg_->GetParameter(2),
                       fg_->GetParameter(1) + 2. * fg_->GetParameter(2));
        hPed_[iz][senP_[i] - 1][senX_[i] - 1][senY_[i] - 1]->setBinContent(is + 1, (int)(fg_->GetParameter(1) + 0.5));
        hTotN_[iz][senP_[i] - 1][senX_[i] - 1][senY_[i] - 1]->setBinContent(is + 1, fg_->GetParameter(2));
      }
    }

  } else {
    if (verbose_)
      cout << "ESPedestalClient: Use Histogram Mean" << endl;

    for (unsigned i = 0; i < senZ_.size(); ++i) {
      iz = (senZ_[i] == 1) ? 0 : 1;

      for (int is = 0; is < 32; ++is) {
        string dirname = prefixME_ + "/ESPedestalTask/";
        sprintf(hname, "ADC Z %d P %d X %d Y %d Str %d", senZ_[i], senP_[i], senX_[i], senY_[i], is + 1);
        MonitorElement *meMean = _igetter.get(dirname + hname);

        if (meMean == nullptr)
          continue;
        TH1F *rootHisto = meMean->getTH1F();

        hPed_[iz][senP_[i] - 1][senX_[i] - 1][senY_[i] - 1]->setBinContent(is + 1, (int)(rootHisto->GetMean() + 0.5));
        hTotN_[iz][senP_[i] - 1][senX_[i] - 1][senY_[i] - 1]->setBinContent(is + 1, rootHisto->GetRMS());
      }
    }
  }
}

Member Data Documentation

TF1* ESPedestalClient::fg_

private

Definition at line 28 of file ESPedestalClient.h.

Referenced by endJobAnalyze(), and ~ESPedestalClient().

bool ESPedestalClient::fitPedestal_

private

Definition at line 23 of file ESPedestalClient.h.

Referenced by endJobAnalyze().

MonitorElement* ESPedestalClient::hPed_[2][2][40][40]

private

Definition at line 25 of file ESPedestalClient.h.

Referenced by book(), endJobAnalyze(), and ESPedestalClient().

MonitorElement* ESPedestalClient::hTotN_[2][2][40][40]

private

Definition at line 26 of file ESPedestalClient.h.

Referenced by book(), endJobAnalyze(), and ESPedestalClient().

std::vector<int> ESPedestalClient::senP_

private

Definition at line 31 of file ESPedestalClient.h.

Referenced by book(), endJobAnalyze(), and ESPedestalClient().

std::vector<int> ESPedestalClient::senX_

private

Definition at line 32 of file ESPedestalClient.h.

Referenced by book(), endJobAnalyze(), and ESPedestalClient().

std::vector<int> ESPedestalClient::senY_

private

Definition at line 33 of file ESPedestalClient.h.

Referenced by book(), endJobAnalyze(), and ESPedestalClient().

std::vector<int> ESPedestalClient::senZ_

private

Definition at line 30 of file ESPedestalClient.h.

Referenced by book(), endJobAnalyze(), and ESPedestalClient().