#include <ESPedestalClient.h>

Inheritance diagram for ESPedestalClient:

Public Member Functions
void	endJobAnalyze (DQMStore::IGetter &) override

	ESPedestalClient (const edm::ParameterSet &)

	~ESPedestalClient ()

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, mergeVDriftHistosByStation::file, edm::FileInPath::fullPath(), edm::ParameterSet::getUntrackedParameter(), hPed_, hTotN_, i, j, roll_playback::k, m, submitDQMOfflineCAF::nLines, 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] = 0;
       hTotN_[i][j][k][m] = 0;
     }
 
    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 ( )

Definition at line 55 of file ESPedestalClient.cc.

References fg_.

                                     {
   delete fg_;
 }

Member Function Documentation

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

overrideprivatevirtual

Reimplemented from ESClient.

Definition at line 115 of file ESPedestalClient.cc.

References DQMStore::IBooker::book1D(), hPed_, hTotN_, 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 59 of file ESPedestalClient.cc.

References gather_cfg::cout, ESClient::debug_, fg_, fitPedestal_, DQMStore::IGetter::get(), MonitorElement::getTH1F(), hPed_, hTotN_, 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==0) 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==0) 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());
 
      } 
       }
    }
 }