#include <ESIntegrityClient.h>

Inheritance diagram for ESIntegrityClient:

Public Member Functions
void	endLumiAnalyze (DQMStore::IGetter &) override
	Analyze. More...

	ESIntegrityClient (const edm::ParameterSet &ps)
	Constructor. More...

	~ESIntegrityClient () override
	Destructor. More...

Public Member Functions inherited from ESClient
virtual void	endJobAnalyze (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
int	fed_ [2][2][40][40]

int	fedStatus_ [56]

int	fiber_ [2][2][40][40]

int	fiberStatus_ [56][36]

TH1F *	hBXDiff_

TH1F *	hFED_

TH2F *	hFiberBadStatus_

TH2F *	hFiberOff_

TH1F *	hKBC_

TH1F *	hKEC_

TH2F *	hKF1_

TH2F *	hKF2_

TH1F *	hL1ADiff_

TH1F *	hOrbitNumberDiff_

TH1F *	hSLinkCRCErr_

int	kchip_ [2][2][40][40]

MonitorElement *	meFED_ [2][2]

MonitorElement *	meKCHIP_ [2][2]

int	slinkCRCStatus_ [56]

int	syncStatus_ [56]

Additional Inherited Members
Public Types inherited from ESClient
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Protected Attributes inherited from ESClient
bool	cloneME_

bool	debug_

bool	initialized_

std::string	prefixME_

bool	verbose_

Detailed Description

Definition at line 9 of file ESIntegrityClient.h.

Constructor & Destructor Documentation

ESIntegrityClient::ESIntegrityClient ( const edm::ParameterSet & ps )

Constructor.

Definition at line 13 of file ESIntegrityClient.cc.

References gather_cfg::cout, fed_, fedStatus_, fiber_, fiberStatus_, mergeVDriftHistosByStation::file, edm::FileInPath::fullPath(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, dqmiolumiharvest::j, isotrackApplyRegressor::k, kchip_, visualization-live-secondInstance_cfg::m, meFED_, meKCHIP_, submitDQMOfflineCAF::nLines, slinkCRCStatus_, AlCaHLTBitMon_QueryRunRegistry::string, syncStatus_, and z.

     : ESClient(ps),
       hFED_(nullptr),
       hFiberOff_(nullptr),
       hFiberBadStatus_(nullptr),
       hKF1_(nullptr),
       hKF2_(nullptr),
       hKBC_(nullptr),
       hKEC_(nullptr),
       hL1ADiff_(nullptr),
       hBXDiff_(nullptr),
       hOrbitNumberDiff_(nullptr),
       hSLinkCRCErr_(nullptr) {
   // read in look-up table
   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) {
           fed_[i][j][k][m] = -1;
           kchip_[i][j][k][m] = -1;
           fiber_[i][j][k][m] = -1;
         }
 
   for (int i = 0; i < 56; ++i) {
     fedStatus_[i] = -1;
     syncStatus_[i] = -1;
     slinkCRCStatus_[i] = -1;
     for (int j = 0; j < 36; ++j)
       fiberStatus_[i][j] = -1;
   }
 
   for (int i = 0; i < 2; ++i)
     for (int j = 0; j < 2; ++j) {
       meFED_[i][j] = nullptr;
       meKCHIP_[i][j] = nullptr;
     }
 
   std::string lutPath(ps.getUntrackedParameter<edm::FileInPath>("LookupTable").fullPath());
 
   int nLines, z, 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;
 
       z = (iz == -1) ? 2 : iz;
       fed_[z - 1][ip - 1][ix - 1][iy - 1] = fed;
       kchip_[z - 1][ip - 1][ix - 1][iy - 1] = kchip;
       fiber_[z - 1][ip - 1][ix - 1][iy - 1] = (fiber - 1) + (optorx - 1) * 12;
     }
 
     file.close();
   } else {
     cout << "ESIntegrityClient : Look up table file can not be found in " << lutPath << endl;
   }
 }

ESIntegrityClient::~ESIntegrityClient ( )

override

Destructor.

Definition at line 73 of file ESIntegrityClient.cc.

73 {}

Member Function Documentation

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

overrideprivatevirtual

Reimplemented from ESClient.

Definition at line 75 of file ESIntegrityClient.cc.

References dqm::implementation::IBooker::book2D(), timingPdfMaker::histo, mps_fire::i, dqmiolumiharvest::j, meFED_, meKCHIP_, ESClient::prefixME_, dqm::impl::MonitorElement::setAxisTitle(), and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                       {
   char histo[200];
 
   _ibooker.setCurrentFolder(prefixME_ + "/ESIntegrityClient");
 
   for (int i = 0; i < 2; ++i)
     for (int j = 0; j < 2; ++j) {
       int iz = (i == 0) ? 1 : -1;
       snprintf(histo, 200, "ES Integrity Summary 1 Z %d P %d", iz, j + 1);
       meFED_[i][j] = _ibooker.book2D(histo, histo, 40, 0.5, 40.5, 40, 0.5, 40.5);
       meFED_[i][j]->setAxisTitle("Si X", 1);
       meFED_[i][j]->setAxisTitle("Si Y", 2);
 
       snprintf(histo, 200, "ES Integrity Summary 2 Z %d P %d", iz, j + 1);
       meKCHIP_[i][j] = _ibooker.book2D(histo, histo, 40, 0.5, 40.5, 40, 0.5, 40.5);
       meKCHIP_[i][j]->setAxisTitle("Si X", 1);
       meKCHIP_[i][j]->setAxisTitle("Si Y", 2);
     }
 }

void ESIntegrityClient::endLumiAnalyze ( DQMStore::IGetter & _igetter )

overridevirtual

Analyze.

Reimplemented from ESClient.

Definition at line 95 of file ESIntegrityClient.cc.

References ESClient::cloneME_, fed_, fedStatus_, fiber_, fiberStatus_, dqm::implementation::IGetter::get(), ESClient::getHisto(), hBXDiff_, hFED_, hFiberBadStatus_, hFiberOff_, hKBC_, hKEC_, hKF1_, hKF2_, hL1ADiff_, hOrbitNumberDiff_, hSLinkCRCErr_, mps_fire::i, dqmiolumiharvest::j, kchip_, hlt_dqm_clientPB-live_cfg::me, meFED_, meKCHIP_, ESClient::prefixME_, dqm::impl::MonitorElement::setBinContent(), slinkCRCStatus_, and syncStatus_.

                                                                 {
   double nDI_FedErr[56];
   for (int i = 0; i < 56; ++i)
     nDI_FedErr[i] = 0;
 
   MonitorElement *me = nullptr;
 
   me = _igetter.get(prefixME_ + "/ESIntegrityTask/ES FEDs used for data taking");
   hFED_ = getHisto(me, cloneME_, hFED_);
 
   me = _igetter.get(prefixME_ + "/ESIntegrityTask/ES Fiber Off");
   hFiberOff_ = getHisto(me, cloneME_, hFiberOff_);
 
   me = _igetter.get(prefixME_ + "/ESIntegrityTask/ES Fiber Bad Status");
   hFiberBadStatus_ = getHisto(me, cloneME_, hFiberBadStatus_);
 
   me = _igetter.get(prefixME_ + "/ESIntegrityTask/ES SLink CRC Errors");
   hSLinkCRCErr_ = getHisto(me, cloneME_, hSLinkCRCErr_);
 
   int xval = 0;
   int nevFEDs = 0;
   if (hFED_) {
     for (int i = 1; i <= 56; ++i)
       if (nevFEDs < hFED_->GetBinContent(i))
         nevFEDs = (int)hFED_->GetBinContent(i);
   }
 
   // FED integrity
   for (int i = 1; i <= 56; ++i) {
     if (hFED_) {
       if (hFED_->GetBinContent(i) > 0)
         fedStatus_[i - 1] = 1;
     }
     if (hSLinkCRCErr_) {
       if (hSLinkCRCErr_->GetBinContent(i) > 0)
         slinkCRCStatus_[i - 1] = 1;
     }
     for (int j = 0; j < 36; ++j) {
       if (hFiberBadStatus_) {
         if (hFiberBadStatus_->GetBinContent(i, j + 1) > 0)
           fiberStatus_[i - 1][j] = 2;  // bad
         else
           fiberStatus_[i - 1][j] = 1;  // good
       }
       if (hFiberOff_)
         if (hFiberOff_->GetBinContent(i, j + 1) > 0) {
           fiberStatus_[i - 1][j] = 0;  // off
         }
     }
   }
 
   // obtain MEs from ESRawDataTask
   me = _igetter.get(prefixME_ + "/ESRawDataTask/ES L1A DCC errors");
   hL1ADiff_ = getHisto(me, cloneME_, hL1ADiff_);
 
   me = _igetter.get(prefixME_ + "/ESRawDataTask/ES BX DCC errors");
   hBXDiff_ = getHisto(me, cloneME_, hBXDiff_);
 
   me = _igetter.get(prefixME_ + "/ESRawDataTask/ES Orbit Number DCC errors");
   hOrbitNumberDiff_ = getHisto(me, cloneME_, hOrbitNumberDiff_);
 
   for (int i = 1; i <= 56; ++i) {
     if (hL1ADiff_ && hBXDiff_) {
       if (hL1ADiff_->GetBinContent(i) > 0 || hBXDiff_->GetBinContent(i) > 0) {
         syncStatus_[i - 1] = 1;
         if (hL1ADiff_->GetBinContent(i) > nDI_FedErr[i - 1])
           nDI_FedErr[i - 1] = hL1ADiff_->GetBinContent(i);
         if (hBXDiff_->GetBinContent(i) > nDI_FedErr[i - 1])
           nDI_FedErr[i - 1] = hBXDiff_->GetBinContent(i);
       }
     }
     // if (hOrbitNumberDiff_->GetBinContent(i) > 0) syncStatus_[i-1] = 1;
   }
 
   // KCHIP integrity
   me = _igetter.get(prefixME_ + "/ESIntegrityTask/ES KChip Flag 1 Error codes");
   hKF1_ = getHisto(me, cloneME_, hKF1_);
 
   me = _igetter.get(prefixME_ + "/ESIntegrityTask/ES KChip Flag 2 Error codes");
   hKF2_ = getHisto(me, cloneME_, hKF2_);
 
   me = _igetter.get(prefixME_ + "/ESIntegrityTask/ES KChip BC mismatch with OptoRX");
   hKBC_ = getHisto(me, cloneME_, hKBC_);
 
   me = _igetter.get(prefixME_ + "/ESIntegrityTask/ES KChip EC mismatch with OptoRX");
   hKEC_ = getHisto(me, cloneME_, hKEC_);
 
   Int_t kchip_xval[1550];
   for (int i = 0; i < 1550; ++i) {
     xval = 3;
     Int_t kErr = 0;
     for (int j = 1; j < 16; ++j) {
       if (hKF1_) {
         if (hKF1_->GetBinContent(i, j + 1) > 0) {
           xval = 2;
           kErr++;
         }
       }
       if (hKF2_) {
         if (hKF2_->GetBinContent(i, j + 1) > 0) {
           xval = 4;
           kErr++;
         }
       }
     }
     if (hKBC_) {
       if (hKBC_->GetBinContent(i) > 0) {
         xval = 5;
         kErr++;
       }
     }
     if (hKEC_) {
       if (hKEC_->GetBinContent(i) > 0) {
         xval = 6;
         kErr++;
       }
     }
     if (kErr > 1)
       xval = 7;
     kchip_xval[i] = xval;
   }
 
   // detailed KCHIP summary (i.e. summary 2)
   for (int iz = 0; iz < 2; ++iz)
     for (int ip = 0; ip < 2; ++ip)
       for (int ix = 0; ix < 40; ++ix)
         for (int iy = 0; iy < 40; ++iy) {
           if (fed_[iz][ip][ix][iy] == -1)
             continue;
           if (fedStatus_[fed_[iz][ip][ix][iy] - 520] == -1 ||
               fiberStatus_[fed_[iz][ip][ix][iy] - 520][fiber_[iz][ip][ix][iy]] == 0)
             kchip_xval[kchip_[iz][ip][ix][iy] - 1] = 0;
           if ((kchip_[iz][ip][ix][iy] - 2) >= 0)
             meKCHIP_[iz][ip]->setBinContent(ix + 1, iy + 1, kchip_xval[kchip_[iz][ip][ix][iy] - 2]);
         }
 
   // FED, Fiber, KCHIP summary (i.e. summary 1)
   Int_t nErr = 0;
   for (int iz = 0; iz < 2; ++iz)
     for (int ip = 0; ip < 2; ++ip)
       for (int ix = 0; ix < 40; ++ix)
         for (int iy = 0; iy < 40; ++iy) {
           if (fed_[iz][ip][ix][iy] == -1)
             continue;
           nErr = 0;
 
           if (fedStatus_[fed_[iz][ip][ix][iy] - 520] == 1) {
             if (hFED_) {
               if (hFED_->GetBinContent(fed_[iz][ip][ix][iy] - 520 + 1) == nevFEDs)
                 xval = 3;
               else {
                 xval = 4;  // FED problem
                 nErr++;
               }
             }
 
             if (fiberStatus_[fed_[iz][ip][ix][iy] - 520][fiber_[iz][ip][ix][iy]] == 2) {
               xval = 2;  // fiber problem
               nErr++;
             }
             if (fiberStatus_[fed_[iz][ip][ix][iy] - 520][fiber_[iz][ip][ix][iy]] == 0) {
               xval = 0;  // fiber off
             }
             if (kchip_xval[kchip_[iz][ip][ix][iy] - 1] != 3 && kchip_xval[kchip_[iz][ip][ix][iy] - 1] != 0) {
               xval = 5;  // kchip problem
               nErr++;
             }
             if (syncStatus_[fed_[iz][ip][ix][iy] - 520] == 1) {
               xval = 6;
               nErr++;
             }
             if (slinkCRCStatus_[fed_[iz][ip][ix][iy] - 520] == 1) {
               xval = 8;
               nErr++;
             }
             if (nErr > 1)
               xval = 7;
           } else if (fedStatus_[fed_[iz][ip][ix][iy] - 520] == -1) {
             xval = 0;
           }
 
           meFED_[iz][ip]->setBinContent(ix + 1, iy + 1, xval);
         }
 }