HcalDigiClient Class Reference

#include <HcalDigiClient.h>

Inheritance diagram for HcalDigiClient:

Public Member Functions
void	analyze (DQMStore::IBooker &, DQMStore::IGetter &)
void	beginRun (void)
void	calculateProblems (DQMStore::IBooker &, DQMStore::IGetter &)
void	cleanup (void)
void	endJob (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
virtual void	calculateProblems (void)
void	getLogicalMap (const edm::EventSetup &es)
	HcalBaseDQClient ()
	HcalBaseDQClient (std::string s, const edm::ParameterSet &ps)
virtual void	htmlOutput (DQMStore::IBooker &, DQMStore::IGetter &, std::string htmlDir)
std::string	name ()
void	setEventSetup (const edm::EventSetup &es)
virtual void	setStatusMap (std::map< HcalDetId, unsigned int > &map)
virtual bool	validHtmlOutput (DQMStore::IBooker &, DQMStore::IGetter &)
virtual	~HcalBaseDQClient (void)

Private Member Functions
void	setupProblemCells (DQMStore::IBooker &, DQMStore::IGetter &)

Private Attributes
bool	doProblemCellSetup_
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_
bool	enableCleanup_
bool	enoughevents_
HcalLogicalMap *	logicalMap_
double	minerrorrate_
int	minevents_
std::string	name_
bool	needLogicalMap_
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_.

{name_="";};

HcalDigiClient::HcalDigiClient

(

std::string

myname

)

Definition at line 18 of file HcalDigiClient.cc.

References HcalBaseDQClient::name_.

{
  name_=myname;
}

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

Definition at line 23 of file HcalDigiClient.cc.

References HcalBaseDQClient::badChannelStatusMask_, HcalBaseDQClient::cloneME_, HcalBaseDQClient::debug_, doProblemCellSetup_, HcalBaseDQClient::enableCleanup_, edm::ParameterSet::getUntrackedParameter(), HcalChannelStatus::HcalCellDead, HFTiming_averageTime, HcalBaseDQClient::minerrorrate_, HcalBaseDQClient::minevents_, HcalBaseDQClient::name_, HcalBaseDQClient::Online_, HcalBaseDQClient::prefixME_, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQClient::subdir_, and HcalBaseDQClient::validHtmlOutput_.

{
  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;

  doProblemCellSetup_ = true;
}

HcalDigiClient::~HcalDigiClient

(

)

Destructor.

Definition at line 299 of file HcalDigiClient.cc.

References HcalBaseDQClient::ProblemCellsByDepth.

{
  if ( ProblemCellsByDepth ) delete ProblemCellsByDepth;
}

Member Function Documentation

void HcalDigiClient::analyze ( DQMStore::IBooker &  ib, DQMStore::IGetter &  ig  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 56 of file HcalDigiClient.cc.

References HcalBaseDQClient::calculateProblems(), HcalBaseDQClient::cloneME_, gather_cfg::cout, HcalBaseDQClient::debug_, doProblemCellSetup_, DQMStore::IGetter::get(), MonitorElement::getTH2F(), HFTiming_averageTime, MonitorElement::Reset(), alignCSCRings::s, MonitorElement::setBinContent(), setupProblemCells(), AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQClient::subdir_, x, and detailsBasic3DVector::y.

{
  if (debug_>2) std::cout <<"\tHcalDigiClient::analyze()"<<std::endl;
  if ( doProblemCellSetup_ ) setupProblemCells(ib,ig);
  calculateProblems(ib,ig);

  // 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=ig.get(s.c_str());
  if (me!=0)
    TimingStudyTime=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_,TimingStudyTime, debug_);

  s=subdir_+"HFTimingStudy/sumplots/HFTiming_Occupancy";
  me=ig.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::beginRun ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 245 of file HcalDigiClient.cc.

References HcalBaseDQClient::enoughevents_, and nevts_.

{
  enoughevents_=false;
  nevts_=0;
}

void HcalDigiClient::calculateProblems	(	DQMStore::IBooker &	ib,
		DQMStore::IGetter &	ig
	)

Definition at line 92 of file HcalDigiClient.cc.

References HcalBaseDQClient::badstatusmap, HcalObjRepresent::CalcIeta(), HcalBaseDQClient::cloneME_, gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, HcalBaseDQClient::enoughevents_, eta(), MonitorElement::Fill(), HcalObjRepresent::FillUnphysicalHEHFBins(), DQMStore::IGetter::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(), AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQClient::subdir_, and ecaldqm::zside().

{
 if (debug_>2) std::cout <<"\t\tHcalDigiClient::calculateProblems()"<<std::endl;
  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=ig.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=ig.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 254 of file HcalDigiClient.cc.

{}

void HcalDigiClient::endJob ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 219 of file HcalDigiClient.cc.

{}

bool HcalDigiClient::hasErrors_Temp ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 256 of file HcalDigiClient.cc.

References HcalObjRepresent::CalcIeta(), gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, HLT_25ns14e33_v1_cff::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 289 of file HcalDigiClient.cc.

{return false;}

bool HcalDigiClient::hasWarnings_Temp ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 288 of file HcalDigiClient.cc.

{return false;}

void HcalDigiClient::setup ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 253 of file HcalDigiClient.cc.

{}

void HcalDigiClient::setupProblemCells ( DQMStore::IBooker &  ib, DQMStore::IGetter &  ig  )

private

Definition at line 221 of file HcalDigiClient.cc.

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

Referenced by analyze().

{
  ib.setCurrentFolder(subdir_);
  problemnames_.clear();
  ProblemCells=ib.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;
  ib.setCurrentFolder(subdir_+"problem_digis");
  ProblemCellsByDepth = new EtaPhiHists();
  ProblemCellsByDepth->setup(ib," Problem Digi Rate");
  for (unsigned int i=0; i<ProblemCellsByDepth->depth.size();++i)
    problemnames_.push_back(ProblemCellsByDepth->depth[i]->getName());

  ib.setCurrentFolder(subdir_+"HFTimingStudy");
  HFTiming_averageTime=ib.book2D("HFTimingStudy_Average_Time","HFTimingStudy Average Time (time sample)",83,-41.5,41.5,72,0.5,72.5);

  doProblemCellSetup_ = false;

}

bool HcalDigiClient::test_enabled ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 290 of file HcalDigiClient.cc.

{return true;}

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

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 293 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

bool HcalDigiClient::doProblemCellSetup_

private

Definition at line 39 of file HcalDigiClient.h.

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

MonitorElement* HcalDigiClient::HFTiming_averageTime

private

Definition at line 36 of file HcalDigiClient.h.

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

int HcalDigiClient::nevts_

private

Definition at line 35 of file HcalDigiClient.h.

Referenced by beginRun().