#include <ESIntegrityTask.h>

Inheritance diagram for ESIntegrityTask:

Public Member Functions
	ESIntegrityTask (const edm::ParameterSet &ps)
virtual	~ESIntegrityTask ()
Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
	Analyze.
void	beginJob (void)
	BeginJob.
void	beginLuminosityBlock (const edm::LuminosityBlock &lumi, const edm::EventSetup &c)
	Begin Lumi.
void	beginRun (const edm::Run &r, const edm::EventSetup &c)
	BeginRun.
void	calculateDIFraction (void)
	Calculate Data Integrity Fraction.
void	cleanup (void)
	Cleanup.
void	endJob (void)
	EndJob.
void	endLuminosityBlock (const edm::LuminosityBlock &lumi, const edm::EventSetup &c)
	End Lumi.
void	endRun (const edm::Run &r, const edm::EventSetup &c)
	EndRun.
void	reset (void)
	Reset.
void	setup (void)
	Setup.
Private Attributes
int	dac_
edm::InputTag	dccCollections_
bool	doLumiAnalysis_
DQMStore *	dqmStore_
int	eCount_
bool	enableCleanup_
int	fed_ [2][2][40][40]
int	fiber_ [2][2][40][40]
int	firstDAC_
int	gain_
int	ievt_
bool	init_
int	isPed_
int	kchip_ [2][2][40][4]
edm::InputTag	kchipCollections_
edm::FileInPath	lookup_
MonitorElement *	meDCCCRCErr_
MonitorElement *	meDCCErr_
MonitorElement *	meDIErrors_ [2][2]
MonitorElement *	meDIErrorsLS_ [2][2]
MonitorElement *	meDIFraction_
MonitorElement *	meEVDR_
MonitorElement *	meFED_
MonitorElement *	meFiberBadStatus_
MonitorElement *	meFiberErrCode_
MonitorElement *	meFiberOff_
MonitorElement *	meGain_
MonitorElement *	meKBC_
MonitorElement *	meKEC_
MonitorElement *	meKF1_
MonitorElement *	meKF2_
MonitorElement *	meOptoBC_
MonitorElement *	meOptoRX_
bool	mergeRuns_
MonitorElement *	meSLinkCRCErr_
int	nDAC_
int	precision_
std::string	prefixME_
int	runNum_
int	runtype_
int	seqtype_
int	vDAC_ [5]

Detailed Description

Definition at line 11 of file ESIntegrityTask.h.

Constructor & Destructor Documentation

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

Definition at line 27 of file ESIntegrityTask.cc.

References gather_cfg::cout, dbtoweb::file, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, j, gen::k, m, cmsCodeRules::cppFunctionSkipper::operator, and z.

                                                       {

   init_ = false;

   dqmStore_ = Service<DQMStore>().operator->();

   prefixME_      = ps.getUntrackedParameter<string>("prefixME", "");
   enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
   mergeRuns_     = ps.getUntrackedParameter<bool>("mergeRuns", false);
   lookup_        = ps.getUntrackedParameter<FileInPath>("LookupTable");

   dccCollections_   = ps.getParameter<InputTag>("ESDCCCollections");
   kchipCollections_ = ps.getParameter<InputTag>("ESKChipCollections");

   doLumiAnalysis_   = ps.getParameter<bool>("DoLumiAnalysis");

   // 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;
         }

   int nLines_, z, iz, ip, ix, iy, fed, kchip, pace, bundle, fiber, optorx;
   ifstream file;
   file.open(lookup_.fullPath().c_str());
   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;
     }
   } 
   else {
     cout<<"ESIntegrityTask : Look up table file can not be found in "<<lookup_.fullPath().c_str()<<endl;
   }

}

ESIntegrityTask::~ESIntegrityTask ( ) [virtual]

Definition at line 75 of file ESIntegrityTask.cc.

{
}

Member Function Documentation

void ESIntegrityTask::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)		`[protected, virtual]`

Analyze.

Implements edm::EDAnalyzer.

Definition at line 267 of file ESIntegrityTask.cc.

                                                                {

   if ( ! init_ ) this->setup();

   ievt_++;

   Handle<ESRawDataCollection> dccs;
   Handle<ESLocalRawDataCollection> kchips;

   // Fill # of events
   meDCCErr_->Fill(575, 2, 1);
   meDCCCRCErr_->Fill(575, 2, 1);
   meOptoRX_->Fill(575, 2, 1);
   meOptoBC_->Fill(575, 2, 1);
   meFiberBadStatus_->Fill(575, 36, 1);
   meFiberOff_->Fill(575, 36, 1);
   meEVDR_->Fill(575, 36, 1);

   // # of DI errors
   Double_t nDIErr[56][36];
   for (int i=0; i<56; ++i) 
     for (int j=0; j<36; ++j)
       nDIErr[i][j] = 0;

   // DCC 
   vector<int> fiberStatus;
   if ( e.getByLabel(dccCollections_, dccs) ) {
     
     for (ESRawDataCollection::const_iterator dccItr = dccs->begin(); dccItr != dccs->end(); ++dccItr) {
       ESDCCHeaderBlock dcc = (*dccItr);
       
       meFED_->Fill(dcc.fedId());
       
       meDCCErr_->Fill(dcc.fedId(), dcc.getDCCErrors());
       
       // SLink CRC error
       if (dcc.getDCCErrors() == 101) {
         meSLinkCRCErr_->Fill(dcc.fedId());
         for (int j=0; j<36; ++j) nDIErr[dcc.fedId()-520][j]++;
       }
       
       if (dcc.getOptoRX0() == 129) {
         meOptoRX_->Fill(dcc.fedId(), 0);
         if (((dcc.getOptoBC0()-15) & 0x0fff) != dcc.getBX()) meOptoBC_->Fill(dcc.fedId(), 0);
       }
       if (dcc.getOptoRX1() == 129) {
         meOptoRX_->Fill(dcc.fedId(), 1);
         if (((dcc.getOptoBC1()-15) & 0x0fff) != dcc.getBX()) meOptoBC_->Fill(dcc.fedId(), 1);
       }
       if (dcc.getOptoRX2() == 129) {
         meOptoRX_->Fill(dcc.fedId(), 2);
         if (((dcc.getOptoBC2()-15) & 0x0fff) != dcc.getBX()) meOptoBC_->Fill(dcc.fedId(), 2);
       }
       
       if (dcc.getOptoRX0() == 128) {
         meDCCCRCErr_->Fill(dcc.fedId(), 0);
         for (int j=0; j<12; ++j) nDIErr[dcc.fedId()-520][j]++;
       }
       if (dcc.getOptoRX1() == 128) {
         meDCCCRCErr_->Fill(dcc.fedId(), 1);
         for (int j=12; j<24; ++j) nDIErr[dcc.fedId()-520][j]++;
       }
       if (dcc.getOptoRX2() == 128) {
         meDCCCRCErr_->Fill(dcc.fedId(), 2);
         for (int j=24; j<36; ++j) nDIErr[dcc.fedId()-520][j]++;
       }
       
       fiberStatus = dcc.getFEChannelStatus();
       
       for (unsigned int i=0; i<fiberStatus.size(); ++i) {
         if (fiberStatus[i]==4 || fiberStatus[i]==8 || fiberStatus[i]==10 || fiberStatus[i]==11 || fiberStatus[i]==12) {
           meFiberBadStatus_->Fill(dcc.fedId(), i+1, 1);
           meFiberErrCode_->Fill(fiberStatus[i]);
           nDIErr[dcc.fedId()-520][i]++;
         }
         if (fiberStatus[i]==7)
           meFiberOff_->Fill(dcc.fedId(), i+1, 1);
         if (fiberStatus[i]==6) {
           meFiberErrCode_->Fill(fiberStatus[i]);
           meEVDR_->Fill(dcc.fedId(), i+1, 1);
         }
       }
       
       runtype_   = dcc.getRunType();
       seqtype_   = dcc.getSeqType();
       dac_       = dcc.getDAC();
       gain_      = dcc.getGain();
       precision_ = dcc.getPrecision();
       
       meGain_->Fill(gain_);
     }
   } else {
     LogWarning("ESIntegrityTask") << dccCollections_  << " not available";
   }

   // KCHIP's
   if ( e.getByLabel(kchipCollections_, kchips) ) {
     
     for (ESLocalRawDataCollection::const_iterator kItr = kchips->begin(); kItr != kchips->end(); ++kItr) {
       
       ESKCHIPBlock kchip = (*kItr);
       
       meKF1_->Fill(kchip.id(), kchip.getFlag1());
       meKF2_->Fill(kchip.id(), kchip.getFlag2());
       if (kchip.getBC() != kchip.getOptoBC()) meKBC_->Fill(kchip.id());
       if (kchip.getEC() != kchip.getOptoEC()) meKEC_->Fill(kchip.id());
     }
   } else {
     LogWarning("ESIntegrityTask") << kchipCollections_  << " not available";
   }

   // Fill # of DI errors
   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 (nDIErr[fed_[iz][ip][ix][iy]-520][fiber_[iz][ip][ix][iy]] > 0)
             meDIErrors_[iz][ip]->Fill(ix+1, iy+1, 1);

         }

}

void ESIntegrityTask::beginJob ( void ) [protected, virtual]

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 78 of file ESIntegrityTask.cc.

                                   {

   ievt_ = 0;

   if ( dqmStore_ ) {
      dqmStore_->setCurrentFolder(prefixME_ + "/ESIntegrityTask");
      dqmStore_->rmdir(prefixME_ + "/ESIntegrityTask");
   }

}

void ESIntegrityTask::beginLuminosityBlock	(	const edm::LuminosityBlock &	lumi,
		const edm::EventSetup &	c
	)		`[protected, virtual]`

Begin Lumi.

Reimplemented from edm::EDAnalyzer.

Definition at line 109 of file ESIntegrityTask.cc.

References i, and j.

                                                                                                  {

  LogInfo("ESIntegrityTask") << "analyzed " << ievt_ << " events";
  // In case of Lumi based analysis SoftReset the Integrity histogram
  if (doLumiAnalysis_) {
    for (int i=0; i<2; ++i) {
      for (int j=0; j<2; ++j) {
        if (meDIErrors_[i][j]) {
          dqmStore_->softReset(meDIErrors_[i][j]);
        }
      }
    }
    ievt_ = 0;
  }
}

void ESIntegrityTask::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)		`[protected, virtual]`

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 89 of file ESIntegrityTask.cc.

References reset().

                                                                {

  if ( ! init_ ) this->setup();
  if ( ! mergeRuns_ ) this->reset();

}

void ESIntegrityTask::calculateDIFraction ( void ) [protected]

Calculate Data Integrity Fraction.

Definition at line 395 of file ESIntegrityTask.cc.

References i, j, ExpressReco_HICollisions_FallBack::x, and ExpressReco_HICollisions_FallBack::y.

                                             {

  float nValidChannels=0; 
  float nGlobalErrors=0;

  for (int i=0; i<2; ++i) {
    for (int j=0; j<2; ++j) {
      float nValidChannelsES=0; 
      float nGlobalErrorsES = 0;
      float reportSummaryES = -1;
      if (!meDIErrors_[i][j]) continue;
      for (int x=0; x<40; ++x) {
        for (int y=0; y<40; ++y) {
          float val = meDIErrors_[i][j]->getBinContent(x+1, y+1);
          if (fed_[i][j][x][y] == -1) continue;
          if (ievt_ != 0) nGlobalErrors += val/ievt_;
          nValidChannels++;
          if (ievt_ != 0) nGlobalErrorsES += val/ievt_;
          nValidChannelsES++;
        }
      }
      if (nValidChannelsES != 0) reportSummaryES = 1 - nGlobalErrorsES/nValidChannelsES;
      meDIFraction_->setBinContent(i+1, j+1, reportSummaryES);
    }
  }
  float reportSummary = -1.0;
  if (nValidChannels != 0) reportSummary = 1.0 - nGlobalErrors/nValidChannels;
  meDIFraction_->setBinContent(3,3, reportSummary);
}

void ESIntegrityTask::cleanup ( void ) [protected]

Cleanup.

Definition at line 247 of file ESIntegrityTask.cc.

                                 {

   if ( ! init_ ) return;

   if ( dqmStore_ ) {

   }

   init_ = false;

}

void ESIntegrityTask::endJob ( void ) [protected, virtual]

EndJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 259 of file ESIntegrityTask.cc.

References edm::cleanup().

                                {

   LogInfo("ESIntegrityTask") << "analyzed " << ievt_ << " events";

   if ( enableCleanup_ ) this->cleanup();

}

void ESIntegrityTask::endLuminosityBlock	(	const edm::LuminosityBlock &	lumi,
		const edm::EventSetup &	c
	)		`[protected, virtual]`

End Lumi.

Reimplemented from edm::EDAnalyzer.

Definition at line 125 of file ESIntegrityTask.cc.

                                                                                               {
  if (doLumiAnalysis_) calculateDIFraction();
}

void ESIntegrityTask::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)		`[protected, virtual]`

EndRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 96 of file ESIntegrityTask.cc.

References i, and j.

                                                              {
  // In case of Lumi based analysis Disable SoftReset from Integrity histogram to get full statistics 
  if (doLumiAnalysis_) {
    for (int i=0; i<2; ++i) {
      for (int j=0; j<2; ++j) {
        if (meDIErrors_[i][j]) {
          dqmStore_->disableSoftReset(meDIErrors_[i][j]);
        }
      }
    }
  }
}

void ESIntegrityTask::reset ( void ) [protected]

Reset.

Definition at line 129 of file ESIntegrityTask.cc.

                                {
   /*
      if ( meFED_ ) meFED_->Reset();
      if ( meGain_ ) meGain_->Reset();
      if ( meDCCErr_ ) meDCCErr_->Reset();
      if ( meOptoRX_ ) meOptoRX_->Reset();
      if ( meOptoBC_ ) meOptoBC_->Reset();
      if ( meFiberStatus_ ) meFiberStatus_->Reset();
      if ( meKF1_ ) meKF1_->Reset();
      if ( meKF2_ ) meKF2_->Reset();
      if ( meKBC_ ) meKBC_->Reset();
      if ( meKEC_ ) meKEC_->Reset();
    */
}

void ESIntegrityTask::setup ( void ) [protected]

Setup.

Definition at line 144 of file ESIntegrityTask.cc.

References trackerHits::histo, i, and j.

                               {

   init_ = true;

   char histo[200];

   if (dqmStore_) {
      dqmStore_->setCurrentFolder(prefixME_ + "/ESIntegrityTask");

      sprintf(histo, "ES FEDs used for data taking");
      meFED_ = dqmStore_->book1D(histo, histo, 56, 519.5, 575.5); 
      meFED_->setAxisTitle("ES FED", 1);
      meFED_->setAxisTitle("Num of Events", 2);

      sprintf(histo, "ES Gain used for data taking");
      meGain_ = dqmStore_->book1D(histo, histo, 2, -0.5, 1.5); 
      meGain_->setAxisTitle("ES Gain", 1);
      meGain_->setAxisTitle("Num of Events", 2);

      sprintf(histo, "ES DCC Error codes");
      meDCCErr_ = dqmStore_->book2D(histo, histo, 56, 519.5, 575.5, 6, -0.5, 5.5); 
      meDCCErr_->setAxisTitle("ES FED", 1);
      meDCCErr_->setAxisTitle("ES DCC Error code", 2);

      sprintf(histo, "ES SLink CRC Errors");
      meSLinkCRCErr_ = dqmStore_->book1D(histo, histo, 56, 519.5, 575.5); 
      meSLinkCRCErr_->setAxisTitle("ES FED", 1);
      meSLinkCRCErr_->setAxisTitle("Num of Events", 2);

      sprintf(histo, "ES DCC CRC Errors");
      meDCCCRCErr_ = dqmStore_->book2D(histo, histo, 56, 519.5, 575.5, 3, -0.5, 2.5); 
      meDCCCRCErr_->setAxisTitle("ES FED", 1);
      meDCCCRCErr_->setAxisTitle("ES OptoRX", 2);

      sprintf(histo, "ES OptoRX used for data taking");
      meOptoRX_ = dqmStore_->book2D(histo, histo, 56, 519.5, 575.5, 3, -0.5, 2.5); 
      meOptoRX_->setAxisTitle("ES FED", 1);
      meOptoRX_->setAxisTitle("ES OptoRX", 2);

      sprintf(histo, "ES OptoRX BC mismatch");
      meOptoBC_ = dqmStore_->book2D(histo, histo, 56, 519.5, 575.5, 3, -0.5, 2.5); 
      meOptoBC_->setAxisTitle("ES FED", 1);
      meOptoBC_->setAxisTitle("ES OptoRX", 2);

      sprintf(histo, "ES Fiber Bad Status");
      meFiberBadStatus_ = dqmStore_->book2D(histo, histo, 56, 519.5, 575.5, 36, 0.5, 36.5);
      meFiberBadStatus_->setAxisTitle("ES FED", 1);
      meFiberBadStatus_->setAxisTitle("Fiber Number", 2);

      sprintf(histo, "ES Fiber Error Code");
      meFiberErrCode_ = dqmStore_->book1D(histo, histo, 17, -0.5, 16.5);
      meFiberErrCode_->setAxisTitle("Fiber Error Code", 1);

      sprintf(histo, "ES Fiber Off");

      meFiberOff_ = dqmStore_->book2D(histo, histo, 56, 519.5, 575.5, 36, 0.5, 36.5);
      meFiberOff_->setAxisTitle("ES FED", 1);
      meFiberOff_->setAxisTitle("Fiber Number", 2);

      sprintf(histo, "ES Warning Event Dropped");
      meEVDR_ = dqmStore_->book2D(histo, histo, 56, 519.5, 575.5, 36, 0.5, 36.5);
      meEVDR_->setAxisTitle("ES FED", 1);
      meEVDR_->setAxisTitle("Fiber Number", 2);

      sprintf(histo, "ES KChip Flag 1 Error codes");
      meKF1_ = dqmStore_->book2D(histo, histo, 1550, -0.5, 1549.5, 16, -0.5, 15.5);
      meKF1_->setAxisTitle("ES KChip", 1);
      meKF1_->setAxisTitle("ES KChip F1 Error Code ", 2);

      sprintf(histo, "ES KChip Flag 2 Error codes");
      meKF2_ = dqmStore_->book2D(histo, histo, 1550, -0.5, 1549.5, 16, -0.5, 15.5);
      meKF2_->setAxisTitle("ES KChip", 1);
      meKF2_->setAxisTitle("ES KChip F1 Error Code ", 2);

      sprintf(histo, "ES KChip BC mismatch with OptoRX");
      meKBC_ = dqmStore_->book1D(histo, histo, 1550, -0.5, 1549.5);
      meKBC_->setAxisTitle("ES KChip", 1);
      meKBC_->setAxisTitle("Num of BC mismatch", 2);

      sprintf(histo, "ES KChip EC mismatch with OptoRX");
      meKEC_ = dqmStore_->book1D(histo, histo, 1550, -0.5, 1549.5);
      meKEC_->setAxisTitle("ES KChip", 1);
      meKEC_->setAxisTitle("Num of EC mismatch", 2);
      
      for (int i=0; i<2; ++i) 
        for (int j=0; j<2; ++j) {
          int iz = (i==0)? 1:-1;
          sprintf(histo, "ES Integrity Errors Z %d P %d", iz, j+1);
          meDIErrors_[i][j] = dqmStore_->book2D(histo, histo, 40, 0.5, 40.5, 40, 0.5, 40.5);
          meDIErrors_[i][j]->setAxisTitle("Si X", 1);
          meDIErrors_[i][j]->setAxisTitle("Si Y", 2);
        }

      if (doLumiAnalysis_) {
        sprintf(histo, "ES Good Channel Fraction");
        meDIFraction_ =  dqmStore_->book2D(histo, histo, 3, 1.0, 3.0, 3, 1.0, 3.0);
        meDIFraction_->setLumiFlag();
      }
   }
   

}