#include <HcalDigiClient.h>

Inheritance diagram for HcalDigiClient:

Public Member Functions
void	analyze (void)

void	beginJob (void)

void	beginRun (void)

void	calculateProblems (void)

void	cleanup (void)

void	endJob (void)

void	endRun (void)

bool	hasErrors_Temp (void)

bool	hasOther_Temp (void)

bool	hasWarnings_Temp (void)

	HcalDigiClient ()
	Constructors. More...

	HcalDigiClient (std::string myname)

	HcalDigiClient (std::string myname, const edm::ParameterSet &ps)

void	setup (void)

bool	test_enabled (void)

void	updateChannelStatus (std::map< HcalDetId, unsigned int > &myqual)

	~HcalDigiClient ()
	Destructor. More...

Public Member Functions inherited from HcalBaseDQClient
	HcalBaseDQClient ()

	HcalBaseDQClient (std::string s, const edm::ParameterSet &ps)

virtual void	htmlOutput (std::string htmlDir)

std::string	name ()

void	setEventSetup (const edm::EventSetup &es)

virtual void	setStatusMap (std::map< HcalDetId, unsigned int > &map)

virtual bool	validHtmlOutput ()

virtual	~HcalBaseDQClient (void)

Private Attributes
MonitorElement *	HFTiming_averageTime

int	nevts_

Additional Inherited Members
Public Attributes inherited from HcalBaseDQClient
int	badChannelStatusMask_

std::map< HcalDetId, unsigned int >	badstatusmap

const edm::EventSetup *	c

bool	cloneME_

int	debug_

DQMStore *	dqmStore_

bool	enableCleanup_

bool	enoughevents_

double	minerrorrate_

int	minevents_

std::string	name_

bool	Online_

std::string	prefixME_

MonitorElement *	ProblemCells

EtaPhiHists *	ProblemCellsByDepth

std::vector< std::string >	problemnames_

std::string	subdir_

bool	testenabled_

bool	validHtmlOutput_

Detailed Description

Definition at line 8 of file HcalDigiClient.h.

Constructor & Destructor Documentation

HcalDigiClient::HcalDigiClient ( )

inline

Constructors.

Definition at line 13 of file HcalDigiClient.h.

References HcalBaseDQClient::name_.

13 {name_="";};

HcalBaseDQClient::name_

std::string name_

Definition: HcalBaseDQClient.h:56

HcalDigiClient::HcalDigiClient ( std::string myname )

Definition at line 20 of file HcalDigiClient.cc.

References HcalBaseDQClient::name_.

 {
   name_=myname;
 }

HcalDigiClient::HcalDigiClient	(	std::string	myname,
		const edm::ParameterSet &	ps
	)

Definition at line 25 of file HcalDigiClient.cc.

 {
   name_=myname;
   enableCleanup_         = ps.getUntrackedParameter<bool>("enableCleanup",false);
   debug_                 = ps.getUntrackedParameter<int>("debug",0);
   prefixME_              = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
   if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
     prefixME_.append("/");
   subdir_                = ps.getUntrackedParameter<std::string>("DigiFolder","DigiMonitor_Hcal/"); // DigiMonitor_Hcal  
   if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
     subdir_.append("/");
   subdir_=prefixME_+subdir_;
 
   validHtmlOutput_       = ps.getUntrackedParameter<bool>("Digi_validHtmlOutput",true);
   cloneME_ = ps.getUntrackedParameter<bool>("cloneME", true);
   badChannelStatusMask_   = ps.getUntrackedParameter<int>("Digi_BadChannelStatusMask",
                                                           ps.getUntrackedParameter<int>("BadChannelStatusMask",
                                             (1<<HcalChannelStatus::HcalCellDead))); // identify channel status values to mask
   
   minerrorrate_ = ps.getUntrackedParameter<double>("Digi_minerrorrate",
                            ps.getUntrackedParameter<double>("minerrorrate",0.05));
   minevents_    = ps.getUntrackedParameter<int>("Digi_minevents",
                         ps.getUntrackedParameter<int>("minevents",1));
   Online_                = ps.getUntrackedParameter<bool>("online",false);
 
   ProblemCellsByDepth=0;
   ProblemCells=0;
 
   HFTiming_averageTime=0;
 }

HcalDigiClient::~HcalDigiClient ( )

Destructor.

Definition at line 307 of file HcalDigiClient.cc.

308 {}

Member Function Documentation

void HcalDigiClient::analyze ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 56 of file HcalDigiClient.cc.

References calculateProblems(), HcalBaseDQClient::cloneME_, gather_cfg::cout, HcalBaseDQClient::debug_, HcalBaseDQClient::dqmStore_, DQMStore::get(), MonitorElement::getTH2F(), HFTiming_averageTime, MonitorElement::Reset(), alignCSCRings::s, MonitorElement::setBinContent(), HcalBaseDQClient::subdir_, vdt::x, and detailsBasic3DVector::y.

Referenced by endRun().

 {
   if (debug_>2) std::cout <<"\tHcalDigiClient::analyze()"<<std::endl;
   calculateProblems();
 
   // Get Pawel's timing plots to form averages
   TH2F* TimingStudyTime=0;
   TH2F* TimingStudyOcc=0;
   std::string s=subdir_+"HFTimingStudy/sumplots/HFTiming_Total_Time";
   
   MonitorElement* me=dqmStore_->get(s.c_str());
   if (me!=0)
     TimingStudyTime=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_,TimingStudyTime, debug_);
 
   s=subdir_+"HFTimingStudy/sumplots/HFTiming_Occupancy";
   me=dqmStore_->get(s.c_str());
   if (me!=0)
     TimingStudyOcc=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_,TimingStudyOcc, debug_);
 
   if (HFTiming_averageTime!=0)
     {
       HFTiming_averageTime->Reset();
       if (TimingStudyTime!=0 && TimingStudyOcc!=0)
         {
           int etabins=(HFTiming_averageTime->getTH2F())->GetNbinsX();
           int phibins=(HFTiming_averageTime->getTH2F())->GetNbinsY();
           for (int x=1;x<=etabins;++x)
         for (int y=1;y<=phibins;++y)
           if (TimingStudyOcc->GetBinContent(x,y)!=0)
             HFTiming_averageTime->setBinContent(x,y,TimingStudyTime->GetBinContent(x,y)*1./TimingStudyOcc->GetBinContent(x,y));
         }
       HFTiming_averageTime->getTH2F()->SetMinimum(0);
     }
 }

void HcalDigiClient::beginJob ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 218 of file HcalDigiClient.cc.

References gather_cfg::cout, HcalBaseDQClient::debug_, HcalBaseDQClient::dqmStore_, cppFunctionSkipper::operator, and DQMStore::showDirStructure().

 {
   dqmStore_ = edm::Service<DQMStore>().operator->();
   if (debug_>0) 
     {
       std::cout <<"<HcalDigiClient::beginJob()>  Displaying dqmStore directory structure:"<<std::endl;
       dqmStore_->showDirStructure();
     }
 }

void HcalDigiClient::beginRun ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 230 of file HcalDigiClient.cc.

References DQMStore::book2D(), gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, HcalBaseDQClient::dqmStore_, HcalBaseDQClient::enoughevents_, MonitorElement::getName(), HFTiming_averageTime, i, nevts_, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, HcalBaseDQClient::problemnames_, DQMStore::setCurrentFolder(), EtaPhiHists::setup(), and HcalBaseDQClient::subdir_.

 {
   enoughevents_=false;
   if (!dqmStore_) 
     {
       if (debug_>0) std::cout <<"<HcalDigiClient::beginRun> dqmStore does not exist!"<<std::endl;
       return;
     }
   dqmStore_->setCurrentFolder(subdir_);
   problemnames_.clear();
   ProblemCells=dqmStore_->book2D(" ProblemDigis",
                  " Problem Digi Rate for all HCAL;ieta;iphi",
                  85,-42.5,42.5,
                  72,0.5,72.5);
   problemnames_.push_back(ProblemCells->getName());
   if (debug_>1)
     std::cout << "Tried to create ProblemCells Monitor Element in directory "<<subdir_<<"  \t  Failed?  "<<(ProblemCells==0)<<std::endl;
   dqmStore_->setCurrentFolder(subdir_+"problem_digis");
   ProblemCellsByDepth = new EtaPhiHists();
   ProblemCellsByDepth->setup(dqmStore_," Problem Digi Rate");
   for (unsigned int i=0; i<ProblemCellsByDepth->depth.size();++i)
     problemnames_.push_back(ProblemCellsByDepth->depth[i]->getName());
 
   nevts_=0;
 
   dqmStore_->setCurrentFolder(subdir_+"HFTimingStudy");
   HFTiming_averageTime=dqmStore_->book2D("HFTimingStudy_Average_Time","HFTimingStudy Average Time (time sample)",83,-41.5,41.5,72,0.5,72.5);
 }

void HcalDigiClient::calculateProblems ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 91 of file HcalDigiClient.cc.

References HcalBaseDQClient::badstatusmap, HcalObjRepresent::CalcIeta(), HcalBaseDQClient::cloneME_, gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, HcalBaseDQClient::dqmStore_, HcalBaseDQClient::enoughevents_, eta(), MonitorElement::Fill(), HcalObjRepresent::FillUnphysicalHEHFBins(), DQMStore::get(), MonitorElement::getBinContent(), MonitorElement::getTH2F(), HcalBarrel, HcalEmpty, HcalEndcap, HcalEtaPhiHistNames(), HcalForward, HcalOuter, i, isHB(), isHE(), isHF(), isHO(), HcalBaseDQClient::minerrorrate_, HcalBaseDQClient::minevents_, HcalBaseDQClient::name(), phi, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, MonitorElement::Reset(), alignCSCRings::s, MonitorElement::setBinContent(), and HcalBaseDQClient::subdir_.

Referenced by analyze().

 {
  if (debug_>2) std::cout <<"\t\tHcalDigiClient::calculateProblems()"<<std::endl;
   if(!dqmStore_) return;
   int totalevents=0;
   int etabins=0, phibins=0, zside=0;
   double problemvalue=0;
 
   // Clear away old problems
   if (ProblemCells!=0)
     {
       ProblemCells->Reset();
       (ProblemCells->getTH2F())->SetMaximum(1.05);
       (ProblemCells->getTH2F())->SetMinimum(0.);
     }
   for  (unsigned int d=0;ProblemCellsByDepth!=0 && d<ProblemCellsByDepth->depth.size();++d)
     {
       if (ProblemCellsByDepth->depth[d]!=0) 
     {
       ProblemCellsByDepth->depth[d]->Reset();
       (ProblemCellsByDepth->depth[d]->getTH2F())->SetMaximum(1.05);
       (ProblemCellsByDepth->depth[d]->getTH2F())->SetMinimum(0.);
     }
     }
 
   // Get histograms that are used in testing
   TH2F* BadDigisByDepth[4];
   TH2F* GoodDigisByDepth[4];
 
   std::vector<std::string> name = HcalEtaPhiHistNames();
 
   bool gothistos=true;
 
   MonitorElement* me;
   for (int i=0;i<4;++i)
     {
       std::string s=subdir_+"bad_digis/bad_digi_occupancy/"+name[i]+"Bad Digi Map";
       me=dqmStore_->get(s.c_str());
       if (me==0) 
     {
       gothistos=false;
       if (debug_>0) std::cout <<"<HcalDigiClient::calculateProblems> Could not get histogram with name "<<s<<std::endl;
     }
       BadDigisByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, BadDigisByDepth[i], debug_);
 
       s=subdir_+"good_digis/digi_occupancy/"+name[i]+" Digi Eta-Phi Occupancy Map";
       me=dqmStore_->get(s.c_str());
       if (me==0) 
     {
       gothistos=false;
       if (debug_>0) std::cout <<"<HcalDigiClient::calculateProblems> Could not get histogram with name "<<s<<std::endl;
     }
       GoodDigisByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, BadDigisByDepth[i], debug_);
     }
 
   if (gothistos==false)
     {
       if (debug_>0) std::cout <<"<HcalDigiClient::calculateProblems> Unable to get all necessary histograms to evaluate problem rate"<<std::endl;
       return;
     }
 
   for (unsigned int d=0;ProblemCellsByDepth!=0 && d<ProblemCellsByDepth->depth.size();++d)
     {
       if (ProblemCellsByDepth->depth[d]==0) continue;
 
       if (BadDigisByDepth[d]==0 || GoodDigisByDepth[d]==0) continue;
       totalevents=(int)GoodDigisByDepth[d]->GetBinContent(0,0);
       if (totalevents<minevents_ ) continue;
       enoughevents_=true;
       etabins=(ProblemCellsByDepth->depth[d]->getTH2F())->GetNbinsX();
       phibins=(ProblemCellsByDepth->depth[d]->getTH2F())->GetNbinsY();
       for (int eta=0;eta<etabins;++eta)
     {
       int ieta=CalcIeta(eta,d+1);
       if (ieta==-9999) continue;
       for (int phi=0;phi<phibins;++phi)
         {
           problemvalue=0; // problem fraction sums over all three tests
           if (BadDigisByDepth[d]->GetBinContent(eta+1,phi+1) > 0) // bad cells found
         problemvalue=(BadDigisByDepth[d]->GetBinContent(eta+1,phi+1)*1./(BadDigisByDepth[d]->GetBinContent(eta+1,phi+1)+GoodDigisByDepth[d]->GetBinContent(eta+1,phi+1)));
           
           zside=0;
           if (isHF(eta,d+1)) // shift ieta by 1 for HF
         ieta<0 ? zside = -1 : zside = 1;
           
           // For problem cells that exceed our allowed rate,
           // set the values to -1 if the cells are already marked in the status database
           if (problemvalue>minerrorrate_)
         {
           HcalSubdetector subdet=HcalEmpty;
           if (isHB(eta,d+1))subdet=HcalBarrel;
           else if (isHE(eta,d+1)) subdet=HcalEndcap;
           else if (isHF(eta,d+1)) subdet=HcalForward;
           else if (isHO(eta,d+1)) subdet=HcalOuter;
           HcalDetId hcalid(subdet, ieta, phi+1, (int)(d+1));
           if (badstatusmap.find(hcalid)!=badstatusmap.end())
             problemvalue=999;
         }
           ProblemCellsByDepth->depth[d]->setBinContent(eta+1,phi+1,problemvalue);
           if (ProblemCells!=0) ProblemCells->Fill(ieta+zside,phi+1,problemvalue);
         } // loop on phi
     } // loop on eta
     } // loop on depth
   
   if (ProblemCells==0)
     {
       if (debug_>0) std::cout <<"<HcalDigiClient::analyze> ProblemCells histogram does not exist!"<<std::endl;
       return;
     }
 
   // Normalization of ProblemCell plot, in the case where there are errors in multiple depths
   etabins=(ProblemCells->getTH2F())->GetNbinsX();
   phibins=(ProblemCells->getTH2F())->GetNbinsY();
   for (int eta=0;eta<etabins;++eta)
     {
       for (int phi=0;phi<phibins;++phi)
     {
       if (ProblemCells->getBinContent(eta+1,phi+1)>1. && ProblemCells->getBinContent(eta+1,phi+1)<999)
         ProblemCells->setBinContent(eta+1,phi+1,1.);
     }
     }
   FillUnphysicalHEHFBins(*ProblemCellsByDepth);
   FillUnphysicalHEHFBins(ProblemCells);
   return;
 }

void HcalDigiClient::cleanup ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 262 of file HcalDigiClient.cc.

262 {}

void HcalDigiClient::endJob ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 228 of file HcalDigiClient.cc.

228 {}

void HcalDigiClient::endRun ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 259 of file HcalDigiClient.cc.

References analyze().

259 {analyze();}

HcalDigiClient::analyze

void analyze(void)

Definition: HcalDigiClient.cc:56

bool HcalDigiClient::hasErrors_Temp ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 264 of file HcalDigiClient.cc.

References HcalObjRepresent::CalcIeta(), gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, MonitorElement::getTH2F(), HcalBaseDQClient::minerrorrate_, HcalBaseDQClient::ProblemCells, and HcalBaseDQClient::ProblemCellsByDepth.

 {
   if (!ProblemCells)
     {
       if (debug_>1) std::cout <<"<HcalDigiClient::hasErrors_Temp>  ProblemCells histogram does not exist!"<<std::endl;
       return false;
     }
   int problemcount=0;
   int ieta=-9999;
 
   for (int depth=0;depth<4; ++depth)
     {
       int etabins  = (ProblemCells->getTH2F())->GetNbinsX();
       int phibins  = (ProblemCells->getTH2F())->GetNbinsY();
       for (int hist_eta=0;hist_eta<etabins;++hist_eta)
         {
           for (int hist_phi=0; hist_phi<phibins;++hist_phi)
             {
               ieta=CalcIeta(hist_eta,depth+1);
           if (ieta==-9999) continue;
           if (ProblemCellsByDepth->depth[depth]==0)
         continue;
           if (ProblemCellsByDepth->depth[depth]->getBinContent(hist_eta,hist_phi)>minerrorrate_)
         ++problemcount;
         } // for (int hist_phi=1;...)
     } // for (int hist_eta=1;...)
     } // for (int depth=0;...)
   
   if (problemcount>0) return true;
   return false;
 }

bool HcalDigiClient::hasOther_Temp ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 297 of file HcalDigiClient.cc.

297 {return false;}

bool HcalDigiClient::hasWarnings_Temp ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 296 of file HcalDigiClient.cc.

296 {return false;}

void HcalDigiClient::setup ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 261 of file HcalDigiClient.cc.

261 {}

bool HcalDigiClient::test_enabled ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 298 of file HcalDigiClient.cc.

298 {return true;}

void HcalDigiClient::updateChannelStatus ( std::map< HcalDetId, unsigned int > & myqual )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 301 of file HcalDigiClient.cc.

 {
   // digi client does not alter channel status yet;
   // look at dead cell or hot cell clients for example code
 } //void HcalDigiClient::updateChannelStatus

Member Data Documentation

MonitorElement* HcalDigiClient::HFTiming_averageTime

private

Definition at line 37 of file HcalDigiClient.h.

Referenced by analyze(), beginRun(), and HcalDigiClient().

int HcalDigiClient::nevts_

private

Definition at line 36 of file HcalDigiClient.h.

Referenced by beginRun().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation