CastorDigiMonitor Class Reference

#include <CastorDigiMonitor.h>

Public Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
	CastorDigiMonitor (const edm::ParameterSet &ps)
int	ModSecToIndex (int module, int sector)
void	processEvent (const CastorDigiCollection &cast, const CastorDbService &cond)
	~CastorDigiMonitor ()

Private Attributes
int	fVerbosity
MonitorElement *	h2digierr
MonitorElement *	h2QmeanMap
MonitorElement *	h2QmeantsvsCh
MonitorElement *	h2QrmsTSvsCh
MonitorElement *	h2QtsvsCh
MonitorElement *	h2qualityMap
MonitorElement *	h2repsum
MonitorElement *	h2status
MonitorElement *	h2TSratio
MonitorElement *	hdigisize
MonitorElement *	hModule
MonitorElement *	hQIErms [10]
MonitorElement *	hReport
MonitorElement *	hSector
MonitorElement *	hTSratio
int	ievt_
float	QIEerrThreshold = 0.0001
double	QmeanTS [224][10]
float	Qrms_DEAD
double	QrmsTS [224][10]
float	RatioThresh1 = 0.
std::string	subsystemname_
int	TS_MAX = 10
const int	TSped = 0

Detailed Description

Definition at line 16 of file CastorDigiMonitor.h.

Constructor & Destructor Documentation

CastorDigiMonitor::CastorDigiMonitor

(

const edm::ParameterSet &

ps

)

Definition at line 17 of file CastorDigiMonitor.cc.

References gather_cfg::cout, fVerbosity, edm::ParameterSet::getUntrackedParameter(), Qrms_DEAD, RatioThresh1, AlCaHLTBitMon_QueryRunRegistry::string, subsystemname_, and TS_MAX.

{
 fVerbosity = ps.getUntrackedParameter<int>("debug",0);
 if(fVerbosity) std::cout<<"CastorDigiMonitor Constructor: "<<this<<std::endl;
subsystemname_=ps.getUntrackedParameter<std::string>("subSystemFolder","Castor");
 RatioThresh1 = ps.getUntrackedParameter<double>("ratioThreshold",0.9);
 Qrms_DEAD = ps.getUntrackedParameter<double>("QrmsDead",0.01); //fC
 Qrms_DEAD = Qrms_DEAD*Qrms_DEAD;
 TS_MAX = ps.getUntrackedParameter<double>("qieTSmax",6);
}

CastorDigiMonitor::~CastorDigiMonitor

(

)

Definition at line 28 of file CastorDigiMonitor.cc.

{ }

Member Function Documentation

void CastorDigiMonitor::bookHistograms	(	DQMStore::IBooker &	,
		edm::Run const &	,
		edm::EventSetup const &
	)

Definition at line 30 of file CastorDigiMonitor.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookFloat(), DQMStore::IBooker::bookProfile2D(), gather_cfg::cout, fVerbosity, MonitorElement::getTH1F(), MonitorElement::getTH2F(), MonitorElement::getTProfile2D(), h2digierr, h2QmeanMap, h2QmeantsvsCh, h2QrmsTSvsCh, h2QtsvsCh, h2qualityMap, h2repsum, h2status, h2TSratio, hdigisize, hModule, hQIErms, hReport, hSector, hTSratio, ievt_, QIEerrThreshold, QmeanTS, QrmsTS, alignCSCRings::s, indexGen::s2, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and subsystemname_.

Referenced by CastorMonitorModule::bookHistograms().

{
  char s[60];
  if(fVerbosity>0) std::cout << "CastorDigiMonitor::beginRun (start)" << std::endl;
  char sTileIndex[50];
  sprintf(sTileIndex,"Tile(=moduleZ*16+sector#phi)");

  ievt_=0;

  ibooker.setCurrentFolder(subsystemname_ + "/CastorDigiMonitor");

  std::string s2 = "CASTOR QIE_capID+er+dv";
  h2digierr=ibooker.bookProfile2D(s2,s2,14,0.,14., 16,0.,16.,100,0,1.e10,"");
  h2digierr->getTProfile2D()->GetXaxis()->SetTitle("Module Z");
  h2digierr->getTProfile2D()->GetYaxis()->SetTitle("Sector #phi");
  h2digierr->getTProfile2D()->SetMaximum(1.);
  h2digierr->getTProfile2D()->SetMinimum(QIEerrThreshold);
  h2digierr->getTProfile2D()->SetOption("colz");

  sprintf(s,"CASTORreportSummaryMap");
   h2repsum=ibooker.bookProfile2D(s,s,14,0.,14., 16,0.,16.,100,0,1.e10,"");
   h2repsum->getTProfile2D()->GetXaxis()->SetTitle("Module Z");
   h2repsum->getTProfile2D()->GetYaxis()->SetTitle("Sector #phi");
   h2repsum->getTProfile2D()->SetMaximum(1.);
   h2repsum->getTProfile2D()->SetMinimum(QIEerrThreshold);
   h2repsum->getTProfile2D()->SetOption("colz");

  sprintf(s,"CASTOR DeadChannelsMap");
    h2status = ibooker.book2D(s,s,14,0.,14., 16,0.,16.);   
    h2status->getTH2F()->GetXaxis()->SetTitle("Module Z");
    h2status->getTH2F()->GetYaxis()->SetTitle("Sector #phi");
    h2status->getTH2F()->SetOption("colz");

  sprintf(s,"CASTOR TSmax Significance Map");
    h2TSratio = ibooker.book2D(s,s,14,0.,14., 16,0.,16.);   
    h2TSratio->getTH2F()->GetXaxis()->SetTitle("Module Z");
    h2TSratio->getTH2F()->GetYaxis()->SetTitle("Sector #phi");
    h2TSratio->getTH2F()->SetOption("colz");

  sprintf(s,"CASTOR TSmax Significance All chan");
    hTSratio = ibooker.book1D(s,s,105,0.,1.05);   

    sprintf(s,"DigiSize");
        hdigisize = ibooker.book1D(s,s,20,0.,20.);
    sprintf(s,"ModuleZ(fC)_allTS");
        hModule = ibooker.book1D(s,s,14,0.,14.);
    hModule->getTH1F()->GetXaxis()->SetTitle("ModuleZ");
    hModule->getTH1F()->GetYaxis()->SetTitle("QIE(fC)");
    sprintf(s,"Sector #phi(fC)_allTS");
        hSector = ibooker.book1D(s,s,16,0.,16.);
    hSector->getTH1F()->GetXaxis()->SetTitle("Sector #phi");
    hSector->getTH1F()->GetYaxis()->SetTitle("QIE(fC)");

    sprintf(s,"QfC=f(x=Tile y=TS) (cumulative)");
      h2QtsvsCh = ibooker.book2D(s,s,224,0.,224., 10,0.,10.);
      h2QtsvsCh->getTH2F()->GetXaxis()->SetTitle(sTileIndex);
      h2QtsvsCh->getTH2F()->GetYaxis()->SetTitle("TS");
      h2QtsvsCh->getTH2F()->SetOption("colz");

    sprintf(s,"QmeanfC=f(Tile TS)");
      h2QmeantsvsCh = ibooker.book2D(s,s,224,0.,224., 10,0.,10.);
      h2QmeantsvsCh->getTH2F()->GetXaxis()->SetTitle(sTileIndex);
      h2QmeantsvsCh->getTH2F()->GetYaxis()->SetTitle("Time Slice");
      h2QmeantsvsCh->getTH2F()->SetOption("colz");

    sprintf(s,"QrmsfC=f(Tile TS)");
      h2QrmsTSvsCh = ibooker.book2D(s,s,224,0.,224., 10,0.,10.);
      h2QrmsTSvsCh->getTH2F()->GetXaxis()->SetTitle(sTileIndex);
      h2QrmsTSvsCh->getTH2F()->GetYaxis()->SetTitle("TS");
      h2QrmsTSvsCh->getTH2F()->SetOption("colz");

 sprintf(s,"CASTOR data quality");
    h2qualityMap = ibooker.book2D(s,s,14, 0,14, 16, 0,16);
    h2qualityMap->getTH2F()->GetXaxis()->SetTitle("module Z");
    h2qualityMap->getTH2F()->GetYaxis()->SetTitle("Sector #phi");
    h2qualityMap->getTH2F()->SetOption("colz");
    
   hReport = ibooker.bookFloat("CASTOR reportSummary");

 sprintf(s,"QmeanfC_map(allTS)");
      h2QmeanMap = ibooker.book2D(s,s,14,0.,14., 16,0.,16.);
      h2QmeanMap->getTH2F()->GetXaxis()->SetTitle("Module Z");
      h2QmeanMap->getTH2F()->GetYaxis()->SetTitle("Sector #phi");
      h2QmeanMap->getTH2F()->SetOption("textcolz");

 for(int ts=0; ts<=1; ts++) {
   sprintf(s,"QIErms_TS=%d",ts);
   hQIErms[ts] = ibooker.book1D(s,s,1000,0.,100.);  
   hQIErms[ts]->getTH1F()->GetXaxis()->SetTitle("QIErms(fC)");
 }

 for(int ind=0; ind<224; ind++) for(int ts=0; ts<10; ts++) 
   QrmsTS[ind][ts] = QmeanTS[ind][ts]= 0.;

 if(fVerbosity>0) 
std::cout<<"CastorDigiMonitor::bookingHist(end)"<<std::endl;
 return;
}

int CastorDigiMonitor::ModSecToIndex	(	int	module,
		int	sector
	)

Definition at line 255 of file CastorDigiMonitor.cc.

Referenced by processEvent().

                                                           {
  int ind = sector + module*16;
  if(ind>223) ind=223;
  return(ind);
}

void CastorDigiMonitor::processEvent	(	const CastorDigiCollection &	cast,
		const CastorDbService &	cond
	)

Definition at line 131 of file CastorDigiMonitor.cc.

References a, HcalQIESample::adc(), edm::SortedCollection< T, SORT >::begin(), HcalQIESample::capid(), gather_cfg::cout, HcalQIESample::dv(), edm::SortedCollection< T, SORT >::empty(), edm::SortedCollection< T, SORT >::end(), HcalQIESample::er(), MonitorElement::Fill(), objects.autophobj::float, fVerbosity, MonitorElement::getTH1F(), MonitorElement::getTH2F(), h2digierr, h2QmeanMap, h2QmeantsvsCh, h2QrmsTSvsCh, h2QtsvsCh, h2qualityMap, h2repsum, h2status, h2TSratio, hdigisize, hModule, hQIErms, hReport, hSector, hTSratio, mps_fire::i, CastorDataFrame::id(), ievt_, isum, LedMonAdc2fc, funct::m, mod(), ModSecToIndex(), HcalCastorDetId::module(), lumiQueryAPI::q, QmeanTS, Qrms_DEAD, QrmsTS, alignCSCRings::r, RatioThresh1, alignCSCRings::s, CastorDataFrame::sample(), HcalCastorDetId::sector(), CastorDataFrame::size(), mathSSE::sqrt(), TS_MAX, and TSped.

Referenced by CastorMonitorModule::analyze().

                                 {
  if(fVerbosity>0) std::cout << "CastorDigiMonitor::processEvent (begin)"<< std::endl;

 if(castorDigis.empty()) {
  for(int mod=0; mod<14; mod++) for(int sec=0; sec<16; sec++) 
        h2repsum->Fill(mod,sec,0.);
  if(fVerbosity>0) std::cout<<"CastorPSMonitor::processEvent NO Castor Digis"<<std::endl;
  return;
 }

 for(CastorDigiCollection::const_iterator j=castorDigis.begin();
     j!=castorDigis.end(); j++)
 {
   const CastorDataFrame digi = (const CastorDataFrame)(*j);    
 
   int capid1 = digi.sample(0).capid();
   hdigisize->Fill(digi.size());
   double sum = 0.;
   int err=0, err2=0;
   int module = digi.id().module()-1;
   int sector = digi.id().sector()-1;
   for (int i=0; i<digi.size(); i++) {
     int capid = digi.sample(i).capid();
     int dv = digi.sample(i).dv();
     int er = digi.sample(i).er();
     int rawd = digi.sample(i).adc();
     rawd = rawd&0x7F;
     err |= (capid != capid1) | er<<1 | (!dv)<<2; // =0
     err2 += (capid != capid1) | er | (!dv); // =0
//     if(err !=0) continue;
     int ind = ModSecToIndex(module,sector);
     h2QtsvsCh->Fill(ind,i,LedMonAdc2fc[rawd]);
     float q = LedMonAdc2fc[rawd];
     sum += q;  //     sum += LedMonAdc2fc[rawd];
        QrmsTS[ind][i] += (q*q);
        QmeanTS[ind][i] += q;
     if(err != 0 && fVerbosity>0)
   std::cout<<"event/idigi=" <<ievt_<<"/"<<i<<" cap=cap1_dv_er_err: "<<
    capid <<"="<< capid1 <<" "<< dv <<" "<< er<<" "<< err << std::endl;
     if(capid1 < 3) capid1 = capid+1;
     else capid1 = 0;
   }
   h2digierr->Fill(module,sector,err);
   h2repsum->Fill(module,sector,1.-err2/digi.size());
//   hBunchOcc->Fill(iBunch,sum);
 } //end for(CastorDigiCollection::const_iterator ...

 ievt_++;

 const float repChanBAD = 0.9;
 const float repChanWarning = 0.95;
 if(ievt_ %100 != 0) return;
   float ModuleSum[14], SectorSum[16];
   for(int m=0; m<14; m++) ModuleSum[m]=0.;
   for(int s=0; s<16; s++) SectorSum[s]=0.;
   for(int mod=0; mod<14; mod++) for(int sec=0; sec<16; sec++) {
     for(int ts=0; ts<=1; ts++) {
     int ind = ModSecToIndex(mod,sec);
       double Qmean = QmeanTS[ind][ts]/ievt_;
       double Qrms = sqrt(QrmsTS[ind][ts]/ievt_ - Qmean*Qmean);
       hQIErms[ts]->Fill(Qrms);
     }

     double sum=0.;
     for(int ts=1; ts<=TS_MAX; ts++) {
     int ind = ModSecToIndex(mod,sec) + 1;
       double a=h2QtsvsCh->getTH2F()->GetBinContent(ind,ts);
  h2QmeantsvsCh->getTH2F()->SetBinContent(ind,ts,a/double(ievt_));
       sum += a;
       double Qmean = QmeanTS[ind-1][ts-1]/ievt_;
       double Qrms = QrmsTS[ind-1][ts-1]/ievt_ - Qmean*Qmean;
  h2QrmsTSvsCh->getTH2F()->SetBinContent(ind,ts,sqrt(Qrms));
     }
     sum /= double(ievt_);
     ModuleSum[mod] += sum;
     SectorSum[sec] += sum;
     float isum = float(int(sum*10.+0.5))/10.;
     h2QmeanMap->getTH2F()->SetBinContent(mod+1,sec+1,isum);
   } // end for(int mod=0; mod<14; mod++) for(int sec=0; sec<16; sec++) 

   for(int mod=0; mod<14; mod++) 
    hModule->getTH1F()->SetBinContent(mod+1,ModuleSum[mod]);
   for(int sec=0; sec<16; sec++) 
    hSector->getTH1F()->SetBinContent(sec+1,SectorSum[sec]);

  int nGoodCh = 0;
  for(int mod=0; mod<14; mod++) for(int sec=0; sec<16;sec++) {
    int ind = ModSecToIndex(mod,sec);
    double Qmean = QmeanTS[ind][TSped]/ievt_;
    double Qrms = QrmsTS[ind][TSped]/ievt_ - Qmean*Qmean;
    float ChanStatus = 0.;
    if(Qrms < Qrms_DEAD) ChanStatus = 1.;
    h2status->getTH2F()->SetBinContent(mod+1,sec+1,ChanStatus);

    int tsm = 0;
    float am=0., amin = 10000.;
    for(int ts=0; ts<TS_MAX; ts++) {
      float a = h2QmeantsvsCh->getTH2F()->GetBinContent(ind+1,ts+1);
      if(am < a) {am = a; tsm = ts;}
      if(a < amin) amin = a;
    }

    double sum = 0.;
    for(int ts=0; ts<TS_MAX; ts++) if(ts != tsm) 
      sum += h2QmeantsvsCh->getTH2F()->GetBinContent(ind+1,ts+1);
    float r = 1.; // worth case - no peak
//    if(am > 0.) r = sum/(TS_MAX-1)/am;
    if(am > amin) r = (sum/(TS_MAX-1)-amin)/(am - amin);
    h2TSratio->getTH2F()->SetBinContent(mod+1,sec+1,r);
    hTSratio->Fill(r);
    float statusTS = 1.0;
    if(r > RatioThresh1) statusTS = repChanWarning;
    else if(r > 0.99) statusTS = repChanBAD;
    float gChanStatus = statusTS;
    if(ChanStatus > 0.) gChanStatus = repChanBAD; // RMS
// if(h2digierr->getTProfile2D()->GetBinContent(mod+1,sec+1)>QIEerrThreshold) gChanStatus = repChanBAD;
    h2qualityMap->getTH2F()->SetBinContent(mod+1,sec+1,gChanStatus);
    if(gChanStatus > repChanBAD) ++nGoodCh;
  }
  hReport->Fill(float(nGoodCh)/224.);
  return;
 }

Member Data Documentation

int CastorDigiMonitor::fVerbosity

private

Definition at line 27 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), CastorDigiMonitor(), and processEvent().

MonitorElement* CastorDigiMonitor::h2digierr

private

Definition at line 36 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::h2QmeanMap

private

Definition at line 42 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::h2QmeantsvsCh

private

Definition at line 41 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::h2QrmsTSvsCh

private

Definition at line 31 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::h2QtsvsCh

private

Definition at line 40 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::h2qualityMap

private

Definition at line 38 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::h2repsum

private

Definition at line 37 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::h2status

private

Definition at line 35 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::h2TSratio

private

Definition at line 34 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::hdigisize

private

Definition at line 45 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::hModule

private

Definition at line 43 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::hQIErms[10]

private

Definition at line 32 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::hReport

private

Definition at line 39 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::hSector

private

Definition at line 44 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorDigiMonitor::hTSratio

private

Definition at line 33 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

int CastorDigiMonitor::ievt_

private

Definition at line 28 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

float CastorDigiMonitor::QIEerrThreshold = 0.0001

private

Definition at line 50 of file CastorDigiMonitor.h.

Referenced by bookHistograms().

double CastorDigiMonitor::QmeanTS[224][10]

private

Definition at line 51 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

float CastorDigiMonitor::Qrms_DEAD

private

Definition at line 29 of file CastorDigiMonitor.h.

Referenced by CastorDigiMonitor(), and processEvent().

double CastorDigiMonitor::QrmsTS[224][10]

private

Definition at line 51 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and processEvent().

float CastorDigiMonitor::RatioThresh1 = 0.

private

Definition at line 49 of file CastorDigiMonitor.h.

Referenced by CastorDigiMonitor(), and processEvent().

std::string CastorDigiMonitor::subsystemname_

private

Definition at line 26 of file CastorDigiMonitor.h.

Referenced by bookHistograms(), and CastorDigiMonitor().

int CastorDigiMonitor::TS_MAX = 10

private

Definition at line 48 of file CastorDigiMonitor.h.

Referenced by CastorDigiMonitor(), and processEvent().

const int CastorDigiMonitor::TSped = 0

private

Definition at line 52 of file CastorDigiMonitor.h.

Referenced by processEvent().