#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 ()
	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
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_, FrontierConditions_GlobalTag_cff::file, edm::FileInPath::fullPath(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, gen::k, kchip_, funct::m, meFED_, meKCHIP_, slinkCRCStatus_, AlCaHLTBitMon_QueryRunRegistry::string, syncStatus_, and z.

                                                               :
   ESClient(ps),
   hFED_(0),
   hFiberOff_(0),
   hFiberBadStatus_(0),
   hKF1_(0),
   hKF2_(0),
   hKBC_(0),
   hKEC_(0), 
   hL1ADiff_(0),
   hBXDiff_(0),
   hOrbitNumberDiff_(0),
   hSLinkCRCErr_(0)
 {
   // 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] = 0;
       meKCHIP_[i][j] = 0;
     }
 
   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 ( )

Destructor.

Definition at line 76 of file ESIntegrityClient.cc.

76 {

77 }

Member Function Documentation

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

overrideprivatevirtual

Reimplemented from ESClient.

Definition at line 79 of file ESIntegrityClient.cc.

References DQMStore::IBooker::book2D(), trackerHits::histo, mps_fire::i, meFED_, meKCHIP_, ESClient::prefixME_, MonitorElement::setAxisTitle(), and DQMStore::IBooker::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 100 of file ESIntegrityClient.cc.

References ESClient::cloneME_, fed_, fedStatus_, fiber_, fiberStatus_, DQMStore::IGetter::get(), ESClient::getHisto(), hBXDiff_, hFED_, hFiberBadStatus_, hFiberOff_, hKBC_, hKEC_, hKF1_, hKF2_, hL1ADiff_, hOrbitNumberDiff_, hSLinkCRCErr_, mps_fire::i, createfilelist::int, kchip_, meFED_, meKCHIP_, ESClient::prefixME_, MonitorElement::setBinContent(), slinkCRCStatus_, and syncStatus_.

                                                                 {
 
   double nDI_FedErr[56];
   for (int i=0; i<56; ++i) nDI_FedErr[i] = 0;
 
   MonitorElement* me = 0;
   
   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); 
     }
   
 }