#include <CastorChannelQualityMonitor.h>

Inheritance diagram for CastorChannelQualityMonitor:

Public Member Functions
	CastorChannelQualityMonitor ()

void	processEvent (const CastorRecHitCollection &castorHits)

void	reset ()

void	setup (const edm::ParameterSet &ps, DQMStore *dbe)

	~CastorChannelQualityMonitor ()

Public Member Functions inherited from CastorBaseMonitor
	CastorBaseMonitor ()

virtual void	clearME ()

virtual void	done ()

bool	getDiagnostics () const

int	getVerbosity () const

void	setDiagnostics (bool myval)

void	setVerbosity (int verb)

bool	vetoCell (HcalCastorDetId id)

virtual	~CastorBaseMonitor ()

Public Attributes
struct {
std::vector< float > energy

int module

int sector

}	ChannelStatus

Private Types
typedef std::map< uint32_t, float, std::less< uint32_t > >	MapType

Private Attributes
int	aboveNoisyThreshold [14][16]

int	aboveThr [14][16]

double	averageEnergy

bool	averageEnergyMethod_

int	belowDThreshold [14][16]

double	castorEnergy

int	castorModule

int	castorSector

int	counter1

int	counter2

double	dThreshold_

double	energyArray [14][16]

double	fraction

TH2F *	h_reportSummaryMap

std::string	histo

int	ievt_

bool	iRecHit

double	nThreshold_

int	numOK

bool	offline_

MonitorElement *	overallStatus

MonitorElement *	reportSummary

MonitorElement *	reportSummaryMap

int	status

double	wcounter1

double	wcounter2

Additional Inherited Members
Protected Attributes inherited from CastorBaseMonitor
std::string	baseFolder_

edm::CPUTimer	cpu_timer

int	fVerbosity

DQMStore *	m_dbe

bool	makeDiagnostics

std::string	rootFolder_

bool	showTiming

Detailed Description

Definition at line 28 of file CastorChannelQualityMonitor.h.

Member Typedef Documentation

typedef std::map<uint32_t,float,std::less<uint32_t> > CastorChannelQualityMonitor::MapType

private

Definition at line 30 of file CastorChannelQualityMonitor.h.

Constructor & Destructor Documentation

CastorChannelQualityMonitor::CastorChannelQualityMonitor ( )

Definition at line 17 of file CastorChannelQualityMonitor.cc.

References counter1, counter2, and ievt_.

                                                         {
  ievt_=0;
  counter1=0;
  counter2=0;
 }

CastorChannelQualityMonitor::~CastorChannelQualityMonitor ( )

Definition at line 27 of file CastorChannelQualityMonitor.cc.

27 {

28 }

Member Function Documentation

void CastorChannelQualityMonitor::processEvent ( const CastorRecHitCollection & castorHits )

Definition at line 103 of file CastorChannelQualityMonitor.cc.

References aboveNoisyThreshold, aboveThr, averageEnergy, averageEnergyMethod_, edm::SortedCollection< T, SORT >::begin(), belowDThreshold, castorEnergy, castorModule, castorSector, counter1, counter2, gather_cfg::cout, dThreshold_, edm::SortedCollection< T, SORT >::end(), energyArray, MonitorElement::Fill(), edm::SortedCollection< T, SORT >::find(), fraction, CastorBaseMonitor::fVerbosity, MonitorElement::getTH2F(), ievt_, iRecHit, CastorBaseMonitor::m_dbe, module, HcalCastorDetId::module(), nThreshold_, numOK, overallStatus, reportSummary, reportSummaryMap, sector, HcalCastorDetId::sector(), edm::SortedCollection< T, SORT >::size(), status, wcounter1, and wcounter2.

Referenced by CastorMonitorModule::analyze().

                                                                                       {
 
   if(fVerbosity>0) std::cout << "CastorChannelQualityMonitor::processEvent !!!" << std::endl;
 
   if(!m_dbe) { 
     if(fVerbosity>0) std::cout <<"CastorChannelQualityMonitor::processEvent => DQMStore is not instantiated !!!"<<std::endl;  
     return; 
   }
 
   if(fVerbosity>0){
      std::cout << "CastorChannelQualityMonitor: Noisy Threshold is set to: "<< nThreshold_ << std::endl;
      std::cout << "CastorChannelQualityMonitor: Dead Threshold is set to: " << dThreshold_ << std::endl;
     }
 
   castorModule = -1;  castorSector = -1; castorEnergy = -1.; 
   
   if (castorHits.size()>0) iRecHit=true;
   else iRecHit=false;
  
  
  for(CastorRecHitCollection::const_iterator recHit = castorHits.begin(); recHit != castorHits.end(); ++recHit){
 
     HcalCastorDetId CastorID = HcalCastorDetId(recHit->id());
     if(fVerbosity>0) std::cout << "Castor ID = " << CastorID << std::endl;
     CastorRecHitCollection::const_iterator rh = castorHits.find(CastorID);
     castorModule = CastorID.module();
     castorSector = CastorID.sector();
     castorEnergy = rh->energy();
     //if(ievt_ % 1000 == 0) std::cout << "==> module=" << module << " sector=" << sector << " energy= "<< energy << std::endl;
     iRecHit=true; 
    
     if(castorEnergy > nThreshold_) ++aboveNoisyThreshold[castorModule-1][castorSector-1];
     if(castorEnergy < dThreshold_) ++belowDThreshold[castorModule-1][castorSector-1];
     if(castorEnergy>1) { //if(castorEnergy>0) {
       ++aboveThr[castorModule-1][castorSector-1]; 
       energyArray[castorModule-1][castorSector-1]=energyArray[castorModule-1][castorSector-1]+castorEnergy;
     }
 
    }
 
 
 
 
   ievt_++;
 
 
 
   // if( (ievt_ == 25 || ievt_ % 500 == 0) && iRecHit ) {
   // no particular event selection done 
    if( iRecHit ) {
    
    status = -99;  numOK = 0; 
 
     // if(offline_) reportSummaryMap->Reset();
 
 
  for (int sector=0; sector<16; sector++){
   for (int module=0; module<14; module++){
   
 
    if(averageEnergyMethod_){
 
      if(aboveThr[module][sector] >0)
      averageEnergy= energyArray[module][sector]/double(aboveThr[module][sector]); // calculate the average energy in each channel
      else averageEnergy=0;
  
      if( averageEnergy >  nThreshold_ )  status= 0;   
      if( averageEnergy < dThreshold_  ) { status= -1;  
           if(fVerbosity>0) std::cout << "!!! dChannels ===> module="<< module+1 << " sector="<< sector+1 << std::endl;
       }
      if( averageEnergy <  nThreshold_ && averageEnergy > dThreshold_ )  status= 1; 
 
      if(fVerbosity>0)
       std::cout << "===> module="<< module+1 << " sector="<< sector+1 <<" *** average Energy=" 
        << averageEnergy << " => energy=" << energyArray[module][sector] << " events="
            << ievt_ << " aboveThr="<< double(aboveThr[module][sector]) <<std::endl;
     }
 
 
    else{
       //-- look at the arrays
       counter1= aboveNoisyThreshold[module][sector]; //counter1 defines how many times the energy was above a noisythreshold
       counter2= belowDThreshold[module][sector];  //counter2 defines how many times the energy was below a dthreshold
       wcounter1= double(counter1)/double(ievt_);
       wcounter2= double(counter2)/double(ievt_);
       if( wcounter1 > 0.85 )  status= 0; 
       if( wcounter2 > 0.85 ) {status= -1;  
          if(fVerbosity>0) std::cout << "!!! dChannels ===> module="<< module+1 << " sector="<< sector+1 << std::endl; 
     }
       if( wcounter1 < 0.85 && wcounter2 <  0.85 ) status= 1; 
 
       if(fVerbosity>0)
         std::cout << "===> module="<< module+1 << " sector="<< sector+1 <<" *** counter1=" 
         << counter1 << " => counter2=" << counter2 << " events="<< ievt_ 
         << " wcounter1=" << wcounter1  << " wcounter2=" <<  wcounter2
         << " *** ==> STATUS=" << status <<std::endl;      
     }
 
    // reportSummaryMap->Fill(module,sector,status);
    reportSummaryMap->getTH2F()->SetBinContent(module+1,sector+1,double(status));
    if (status == 1) numOK++;
 
       }
     }
   fraction=double(numOK)/224;
   overallStatus->Fill(fraction); reportSummary->Fill(fraction); 
   }
   
   
    return;
 } 

void CastorChannelQualityMonitor::reset ( void )

Definition at line 31 of file CastorChannelQualityMonitor.cc.

Referenced by CastorMonitorModule::reset().

31 {

32 }

void CastorChannelQualityMonitor::setup	(	const edm::ParameterSet &	ps,
		DQMStore *	dbe
	)

virtual

Reimplemented from CastorBaseMonitor.

Definition at line 39 of file CastorChannelQualityMonitor.cc.

References aboveNoisyThreshold, aboveThr, averageEnergy, averageEnergyMethod_, belowDThreshold, DQMStore::book2D(), DQMStore::bookFloat(), cuy::col, counter1, counter2, gather_cfg::cout, dThreshold_, energyArray, MonitorElement::Fill(), fraction, CastorBaseMonitor::fVerbosity, MonitorElement::getTH2F(), edm::ParameterSet::getUntrackedParameter(), h_reportSummaryMap, ievt_, iRecHit, CastorBaseMonitor::m_dbe, nThreshold_, NULL, numOK, offline_, overallStatus, reportSummary, reportSummaryMap, CastorBaseMonitor::rootFolder_, DQMStore::setCurrentFolder(), CastorBaseMonitor::setup(), status, wcounter1, and wcounter2.

Referenced by CastorMonitorModule::CastorMonitorModule().

                                                                                {
 
   CastorBaseMonitor::setup(ps,dbe);
    averageEnergyMethod_ = ps.getUntrackedParameter<bool>("averageEnergyMethod", true); 
    offline_             = ps.getUntrackedParameter<bool>("OfflineMode", false); 
    nThreshold_          = ps.getUntrackedParameter<double>("nThreshold", 0);
    dThreshold_          = ps.getUntrackedParameter<double>("dThreshold", 0);
 
   if ( m_dbe !=NULL ) {    
     m_dbe->setCurrentFolder(rootFolder_+"CastorChannelQuality");
     reportSummary    = m_dbe->bookFloat("RecHit Energy based reportSummary");
     reportSummaryMap = m_dbe->book2D("RecHitEnergyBasedSummaryMap","RecHitEnergyBasedSummaryMap",14,0.0,14.0,16,0.0,16.0);
     if(offline_){
       h_reportSummaryMap =reportSummaryMap->getTH2F();
       h_reportSummaryMap->SetOption("textcolz");
       h_reportSummaryMap->GetXaxis()->SetTitle("module");
       h_reportSummaryMap->GetYaxis()->SetTitle("sector");
     }
     
     overallStatus = m_dbe->bookFloat("RecHit Energy based fraction of good channels");
    } 
 
   else{
   if(fVerbosity>0) std::cout << "CastorChannelQualityMonitor::setup - NO DQMStore service" << std::endl; 
  }
 
   // baseFolder_ = rootFolder_+"CastorChannelQualityMonitor";
 
   if(fVerbosity>0) std::cout << "CastorChannelQualityMonitor::setup (start)" << std::endl;
 
   for (int row=0; row<14; row++) {
     for (int col=0; col<16; col++){
         aboveNoisyThreshold[row][col] = 0;
         belowDThreshold[row][col]     = 0;
         energyArray[row][col]         = 0;
         aboveThr[row][col]         = 0;
     }
   }
 
   ievt_=0; counter1=0; counter2=0;  wcounter1=0.; wcounter2=0.; 
   fraction=0.; 
   overallStatus->Fill(fraction); reportSummary->Fill(fraction); 
   status = -99; numOK = 0; 
   iRecHit=false;
   averageEnergy=0.;
 
   if(fVerbosity>0) std::cout << "CastorChannelQualityMonitor::setup (end)" << std::endl;
 
   return;
 }