HcalRecHitClient Class Reference

#include <HcalRecHitClient.h>

Inheritance diagram for HcalRecHitClient:

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)
	HcalRecHitClient ()
	Constructors. More...
	HcalRecHitClient (std::string myname)
	HcalRecHitClient (std::string myname, const edm::ParameterSet &ps)
void	setup (void)
bool	test_enabled (void)
void	updateChannelStatus (std::map< HcalDetId, unsigned int > &myqual)
	~HcalRecHitClient ()
	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_
EtaPhiHists *	meEnergyByDepth
EtaPhiHists *	meEnergyThreshByDepth
MonitorElement *	meHBEnergy_1D
MonitorElement *	meHBEnergyRMS_1D
MonitorElement *	meHBEnergyRMSThresh_1D
MonitorElement *	meHBEnergyThresh_1D
MonitorElement *	meHEEnergy_1D
MonitorElement *	meHEEnergyRMS_1D
MonitorElement *	meHEEnergyRMSThresh_1D
MonitorElement *	meHEEnergyThresh_1D
MonitorElement *	meHFEnergy_1D
MonitorElement *	meHFEnergyRMS_1D
MonitorElement *	meHFEnergyRMSThresh_1D
MonitorElement *	meHFEnergyThresh_1D
MonitorElement *	meHOEnergy_1D
MonitorElement *	meHOEnergyRMS_1D
MonitorElement *	meHOEnergyRMSThresh_1D
MonitorElement *	meHOEnergyThresh_1D
EtaPhiHists *	meSqrtSumEnergy2ByDepth
EtaPhiHists *	meSqrtSumEnergy2ThreshByDepth
EtaPhiHists *	meTimeByDepth
EtaPhiHists *	meTimeThreshByDepth
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 10 of file HcalRecHitClient.h.

Constructor & Destructor Documentation

HcalRecHitClient::HcalRecHitClient ( )

inline

Constructors.

Definition at line 15 of file HcalRecHitClient.h.

References HcalBaseDQClient::name_.

{name_="";};

HcalRecHitClient::HcalRecHitClient

(

std::string

myname

)

Definition at line 18 of file HcalRecHitClient.cc.

References HcalBaseDQClient::name_.

{
  name_=myname;
}

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

Definition at line 23 of file HcalRecHitClient.cc.

References HcalBaseDQClient::badChannelStatusMask_, HcalBaseDQClient::cloneME_, HcalBaseDQClient::debug_, doProblemCellSetup_, HcalBaseDQClient::enableCleanup_, HcalBaseDQClient::enoughevents_, edm::ParameterSet::getUntrackedParameter(), meEnergyByDepth, meEnergyThreshByDepth, meSqrtSumEnergy2ByDepth, meSqrtSumEnergy2ThreshByDepth, meTimeByDepth, meTimeThreshByDepth, 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>("RecHitFolder","RecHitMonitor_Hcal/"); // RecHitMonitor_Hcal  
  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
    subdir_.append("/");
  subdir_=prefixME_+subdir_;

  validHtmlOutput_       = ps.getUntrackedParameter<bool>("RecHit_validHtmlOutput",true);
  cloneME_ = ps.getUntrackedParameter<bool>("cloneME", true);
  badChannelStatusMask_   = ps.getUntrackedParameter<int>("RecHit_BadChannelStatusMask",
                                                          ps.getUntrackedParameter<int>("BadChannelStatusMask",
                                            0)); // identify channel status values to mask

  minerrorrate_ = ps.getUntrackedParameter<double>("RecHit_minerrorrate",
                           ps.getUntrackedParameter<double>("minerrorrate",0.01));
  minevents_    = ps.getUntrackedParameter<int>("RecHit_minevents",
                        ps.getUntrackedParameter<int>("minevents",1));
  enoughevents_=false;
  Online_                = ps.getUntrackedParameter<bool>("online",false);

  ProblemCells=0;
  ProblemCellsByDepth=0;

  meEnergyByDepth=0;
  meEnergyThreshByDepth=0;
  meTimeByDepth=0;
  meTimeThreshByDepth=0;
  meSqrtSumEnergy2ByDepth=0;
  meSqrtSumEnergy2ThreshByDepth=0;

  doProblemCellSetup_ = true;
}

HcalRecHitClient::~HcalRecHitClient

(

)

Destructor.

Definition at line 505 of file HcalRecHitClient.cc.

References meEnergyByDepth, meEnergyThreshByDepth, meSqrtSumEnergy2ByDepth, meSqrtSumEnergy2ThreshByDepth, meTimeByDepth, meTimeThreshByDepth, and HcalBaseDQClient::ProblemCellsByDepth.

{
  if ( ProblemCellsByDepth ) delete ProblemCellsByDepth;
  if ( meEnergyByDepth ) delete meEnergyByDepth;
  if ( meEnergyThreshByDepth ) delete meEnergyThreshByDepth;
  if ( meTimeByDepth ) delete meTimeByDepth;
  if ( meTimeThreshByDepth ) delete meTimeThreshByDepth;
  if ( meSqrtSumEnergy2ByDepth ) delete meSqrtSumEnergy2ByDepth;
  if ( meSqrtSumEnergy2ThreshByDepth ) delete meSqrtSumEnergy2ThreshByDepth;

}

Member Function Documentation

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

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 62 of file HcalRecHitClient.cc.

References HcalObjRepresent::CalcIeta(), HcalBaseDQClient::calculateProblems(), HcalBaseDQClient::cloneME_, gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, doProblemCellSetup_, eta(), MonitorElement::Fill(), HcalObjRepresent::FillUnphysicalHEHFBins(), DQMStore::IGetter::get(), HcalBarrel, HcalEndcap, HcalEtaPhiHistNames(), HcalForward, HcalOuter, i, isHB(), isHE(), isHF(), isHO(), meEnergyByDepth, meEnergyThreshByDepth, meHBEnergy_1D, meHBEnergyRMS_1D, meHBEnergyRMSThresh_1D, meHBEnergyThresh_1D, meHEEnergy_1D, meHEEnergyRMS_1D, meHEEnergyRMSThresh_1D, meHEEnergyThresh_1D, meHFEnergy_1D, meHFEnergyRMS_1D, meHFEnergyRMSThresh_1D, meHFEnergyThresh_1D, meHOEnergy_1D, meHOEnergyRMS_1D, meHOEnergyRMSThresh_1D, meHOEnergyThresh_1D, meTimeByDepth, meTimeThreshByDepth, HcalBaseDQClient::name(), phi, funct::pow(), MonitorElement::Reset(), RMS, alignCSCRings::s, setupProblemCells(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQClient::subdir_, and validDetId().

{
  if (debug_>2) std::cout <<"\tHcalRecHitClient::analyze()"<<std::endl;

  if ( doProblemCellSetup_ ) setupProblemCells(ib,ig);
 
  TH2F* OccupancyByDepth[4];
  TH2F* SumEnergyByDepth[4];
  TH2F* SumTimeByDepth[4];
  TH2F* SqrtSumEnergy2ByDepth[4];

  TH2F* OccupancyThreshByDepth[4];
  TH2F* SumEnergyThreshByDepth[4];
  TH2F* SumTimeThreshByDepth[4];
  TH2F* SqrtSumEnergy2ThreshByDepth[4];

  std::vector<std::string> name = HcalEtaPhiHistNames();

  MonitorElement* me;
  bool gotHistos=true;

  for (int i=0;i<4;++i)
    {
      std::string s=subdir_+"Distributions_AllRecHits/"+name[i]+"RecHit Occupancy";
      me=ig.get(s.c_str());
      if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
      OccupancyByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, OccupancyByDepth[i], debug_);
      s=subdir_+"Distributions_AllRecHits/sumplots/"+name[i]+"RecHit Summed Energy GeV";
      me=ig.get(s.c_str());
      if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
      SumEnergyByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SumEnergyByDepth[i], debug_);
      s=subdir_+"Distributions_AllRecHits/sumplots/"+name[i]+"RecHit Summed Time nS";
      me=ig.get(s.c_str());
      if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
      SumTimeByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SumTimeByDepth[i], debug_);
      s=subdir_+"Distributions_AllRecHits/sumplots/"+name[i]+"RecHit Sqrt Summed Energy2 GeV";
      me=ig.get(s.c_str());
      if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
      SqrtSumEnergy2ByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SqrtSumEnergy2ByDepth[i], debug_);

      // Threshold histograms
      s=subdir_+"Distributions_PassedMinBias/"+name[i]+"Above Threshold RecHit Occupancy";
      me=ig.get(s.c_str());
      if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
      OccupancyThreshByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, OccupancyThreshByDepth[i], debug_);
      s=subdir_+"Distributions_PassedMinBias/sumplots/"+name[i]+"Above Threshold RecHit Summed Energy GeV";
      me=ig.get(s.c_str());
      if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
      SumEnergyThreshByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SumEnergyThreshByDepth[i], debug_);
      s=subdir_+"Distributions_PassedMinBias/sumplots/"+name[i]+"Above Threshold RecHit Summed Time nS";
      me=ig.get(s.c_str());
      if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
      SumTimeThreshByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SumTimeThreshByDepth[i], debug_);
      s=subdir_+"Distributions_PassedMinBias/sumplots/"+name[i]+"Above Threshold RecHit Sqrt Summed Energy2 GeV";
      me=ig.get(s.c_str());
      if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
      SqrtSumEnergy2ThreshByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SqrtSumEnergy2ThreshByDepth[i], debug_);
    }
  if (gotHistos==false)
    {
      if (debug_>0) std::cout <<"<HcalRecHitClient::calculateProblems()> Not all histograms could be found; skipping normalization"<<std::endl;
      return;
    }

  // Clear histograms before re-filling
  meHBEnergy_1D->Reset();
  meHBEnergyRMS_1D->Reset();
  meHEEnergy_1D->Reset();
  meHEEnergyRMS_1D->Reset();
  meHOEnergy_1D->Reset();
  meHOEnergyRMS_1D->Reset();
  meHFEnergy_1D->Reset();
  meHFEnergyRMS_1D->Reset();
  meHBEnergyThresh_1D->Reset();
  meHBEnergyRMSThresh_1D->Reset();
  meHEEnergyThresh_1D->Reset();
  meHEEnergyRMSThresh_1D->Reset();
  meHOEnergyThresh_1D->Reset();
  meHOEnergyRMSThresh_1D->Reset();
  meHFEnergyThresh_1D->Reset();
  meHFEnergyRMSThresh_1D->Reset();

  for (int mydepth=0;mydepth<4;++mydepth)
    {
      for (int eta=0;eta<OccupancyByDepth[mydepth]->GetNbinsX();++eta)
    {
      // eta+1=1:  ieta = -42
      // eta+1=13: ieta = -29

      for (int phi=0;phi<72;++phi)
        {
          if (OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1)>0)
        {
          // fill 1D plots
          if (isHB(eta,mydepth+1)) 
            {
              if (validDetId(HcalBarrel, CalcIeta(HcalBarrel, eta, mydepth+1), phi+1, mydepth+1))
            {
              meHBEnergy_1D->Fill(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
              meHBEnergyRMS_1D->Fill(sqrt(pow(SqrtSumEnergy2ByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1)-pow(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)));
            }
            } 
          else if (isHE(eta,mydepth+1)) 
            {
              if (validDetId(HcalEndcap, CalcIeta(HcalEndcap, eta, mydepth+1), phi+1, mydepth+1)) 
            {
              
              meHEEnergy_1D->Fill(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
              meHEEnergyRMS_1D->Fill(sqrt(pow(SqrtSumEnergy2ByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1)-pow(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)));
            }
            }
          else if (isHO(eta,mydepth+1)) 
            {
              if (validDetId(HcalOuter, CalcIeta(HcalOuter, eta, mydepth+1), phi+1, mydepth+1)) 
            {
              meHOEnergy_1D->Fill(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
              meHOEnergyRMS_1D->Fill(sqrt(pow(SqrtSumEnergy2ByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1)-pow(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)));
            }
            } 
          else if (isHF(eta,mydepth+1)) 
            {
              if (validDetId(HcalForward, CalcIeta(HcalForward, eta, mydepth+1), phi+1, mydepth+1)) 
            {
              meHFEnergy_1D->Fill(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
              meHFEnergyRMS_1D->Fill(sqrt(pow(SqrtSumEnergy2ByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1)-pow(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)));
            }
            }
          // normalize 2D plots by number of events
         
          meEnergyByDepth->depth[mydepth]->setBinContent(eta+1, phi+1, SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
          meTimeByDepth->depth[mydepth]->setBinContent(eta+1, phi+1, SumTimeByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
        } // (OccupancyByDepth[mydepth]->GetBinContent(eta+1,phi+1)>0)
                
          if (OccupancyThreshByDepth[mydepth]==0) continue;
          if (OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)>0)
        {
          // fill 1D plots
          if (isHB(eta,mydepth+1)) 
            {
              if (validDetId(HcalBarrel, CalcIeta(HcalBarrel, eta, mydepth+1), phi+1, mydepth+1))
            {
              meHBEnergyThresh_1D->Fill(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
              double RMS=pow(SqrtSumEnergy2ThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1,phi+1)-pow(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2);
              RMS=pow(fabs(RMS),0.5);
              meHBEnergyRMSThresh_1D->Fill(RMS);
            }
            } 
          else if (isHE(eta,mydepth+1)) 
            {
              if (validDetId(HcalEndcap, CalcIeta(HcalEndcap, eta, mydepth+1), phi+1, mydepth+1)) 
            {
              
              meHEEnergyThresh_1D->Fill(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
              double RMS=pow(SqrtSumEnergy2ThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1,phi+1)-pow(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2);
              RMS=pow(fabs(RMS),0.5);
              meHEEnergyRMSThresh_1D->Fill(RMS);
            }
            }
          else if (isHO(eta,mydepth+1)) 
            {
              if (validDetId(HcalOuter, CalcIeta(HcalOuter, eta, mydepth+1), phi+1, mydepth+1)) 
            {
              meHOEnergyThresh_1D->Fill(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
              double RMS=pow(SqrtSumEnergy2ThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1,phi+1)-pow(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2);
              RMS=pow(fabs(RMS),0.5);
              meHOEnergyRMSThresh_1D->Fill(RMS);
            }
            } 
          else if (isHF(eta,mydepth+1)) 
            {
              if (validDetId(HcalForward, CalcIeta(HcalForward, eta, mydepth+1), phi+1, mydepth+1)) 
            {
              meHFEnergyThresh_1D->Fill(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
              double RMS=pow(SqrtSumEnergy2ThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1,phi+1)-pow(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2);
              RMS=pow(fabs(RMS),0.5);
              meHFEnergyRMSThresh_1D->Fill(RMS);          
            }
            }
          // fill 2D plots
          meEnergyThreshByDepth->depth[mydepth]->setBinContent(eta+1, phi+1, SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
          meTimeThreshByDepth->depth[mydepth]->setBinContent(eta+1, phi+1, SumTimeThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
        }
        } // for (int phi=0;phi<72;++phi)
    } // for (int eta=0;eta<OccupancyByDepth->..;++eta)
    } // for (int mydepth=0;...)

  FillUnphysicalHEHFBins(*meEnergyByDepth);
  FillUnphysicalHEHFBins(*meTimeByDepth);
  FillUnphysicalHEHFBins(*meEnergyThreshByDepth);
  FillUnphysicalHEHFBins(*meTimeThreshByDepth);

  calculateProblems(ib,ig);
}

void HcalRecHitClient::beginRun ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 444 of file HcalRecHitClient.cc.

References gather_cfg::cout, and HcalBaseDQClient::debug_.

{
  if (debug_>1)  std::cout <<"<HcalRecHitClient::endRun>"<<std::endl;

}

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

Definition at line 256 of file HcalRecHitClient.cc.

References HcalBaseDQClient::badstatusmap, HcalObjRepresent::CalcIeta(), gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, HcalBaseDQClient::enoughevents_, eta(), MonitorElement::Fill(), HcalObjRepresent::FillUnphysicalHEHFBins(), MonitorElement::getBinContent(), MonitorElement::getTH2F(), HcalBarrel, HcalEmpty, HcalEndcap, HcalEtaPhiHistNames(), HcalForward, HcalOuter, isHB(), isHE(), isHF(), isHO(), min(), HcalBaseDQClient::minerrorrate_, HcalBaseDQClient::minevents_, HcalBaseDQClient::name(), phi, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, MonitorElement::Reset(), MonitorElement::setBinContent(), and ecaldqm::zside().

{
 if (debug_>2) std::cout <<"\t\tHcalRecHitClient::calculateProblems()"<<std::endl;
  double totalevents=0;
  int etabins=0, phibins=0, zside=0;
  double problemvalue=0;
  std::vector<std::string> name = HcalEtaPhiHistNames(); // use this to get EtaPhiHistograms that feed problem calculation, once they exist (see analyze function for example usage of HcalEtaPhiHistNames())

  // Get histograms used in determining rechit problem rate, 
  // and get totalevents from their underflow bins.
  // No such problem histograms are defined so far.
  
  enoughevents_=true;

  if (totalevents==0)
    return;

  // Clear away old problems
  if (ProblemCells!=0)
    {
      ProblemCells->Reset();
      (ProblemCells->getTH2F())->SetMaximum(1.05);
      (ProblemCells->getTH2F())->SetMinimum(0.);
    }
  for  (unsigned int d=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.);
    }
    }
  
  for (unsigned int d=0;d<ProblemCellsByDepth->depth.size();++d)
    {
      if (ProblemCellsByDepth->depth[d]==0) continue;

      // Get total events from some future histogram
      //totalevents=DigiPresentByDepth[d]->GetBinContent(0); // get totalevents from each depth, in case they differ
      if (totalevents==0 || totalevents<minevents_) continue;
      etabins=(ProblemCellsByDepth->depth[d]->getTH2F())->GetNbinsX();
      phibins=(ProblemCellsByDepth->depth[d]->getTH2F())->GetNbinsY();
      problemvalue=0;
      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;
          // Add histograms here when rechit testing decided upon
          /*if (DigiPresentByDepth[d]!=0 && DigiPresentByDepth[d]->GetBinContent(eta+1,phi+1)>0) 
        problemvalue=totalevents; */
          
          if (problemvalue==0) continue;
          problemvalue/=totalevents; // problem value is a rate; should be between 0 and 1
          problemvalue = std::min(1.,problemvalue);
          
          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 999 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 <<"<HcalRecHitClient::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 HcalRecHitClient::cleanup ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 457 of file HcalRecHitClient.cc.

{}

void HcalRecHitClient::endJob ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 363 of file HcalRecHitClient.cc.

{}

bool HcalRecHitClient::hasErrors_Temp ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 459 of file HcalRecHitClient.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 <<"<HcalRecHitClient::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 HcalRecHitClient::hasOther_Temp ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 492 of file HcalRecHitClient.cc.

{return false;}

bool HcalRecHitClient::hasWarnings_Temp ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 491 of file HcalRecHitClient.cc.

{return false;}

void HcalRecHitClient::setup ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 456 of file HcalRecHitClient.cc.

{}

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

private

Definition at line 365 of file HcalRecHitClient.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, doProblemCellSetup_, MonitorElement::getName(), i, meEnergyByDepth, meEnergyThreshByDepth, meHBEnergy_1D, meHBEnergyRMS_1D, meHBEnergyRMSThresh_1D, meHBEnergyThresh_1D, meHEEnergy_1D, meHEEnergyRMS_1D, meHEEnergyRMSThresh_1D, meHEEnergyThresh_1D, meHFEnergy_1D, meHFEnergyRMS_1D, meHFEnergyRMSThresh_1D, meHFEnergyThresh_1D, meHOEnergy_1D, meHOEnergyRMS_1D, meHOEnergyRMSThresh_1D, meHOEnergyThresh_1D, meTimeByDepth, meTimeThreshByDepth, nevts_, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, HcalBaseDQClient::problemnames_, DQMStore::IBooker::setCurrentFolder(), EtaPhiHists::setup(), HcalBaseDQClient::subdir_, x, and detailsBasic3DVector::y.

Referenced by analyze().

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

  nevts_=0;

  ib.setCurrentFolder(subdir_+"Distributions_AllRecHits");
  meEnergyByDepth = new EtaPhiHists();
  meEnergyByDepth->setup(ib,"RecHit Average Energy","GeV");
  meTimeByDepth = new EtaPhiHists();
  meTimeByDepth->setup(ib,"RecHit Average Time","nS");
  // set all average times to -1000 by default (so that they don't show up on plots
  for (unsigned int i=0;i<meTimeByDepth->depth.size();++i)
    {
      (meTimeByDepth->depth[i]->getTH2F())->SetMinimum(-150);
      (meTimeByDepth->depth[i]->getTH2F())->SetMaximum(150);
      int etabins=(meTimeByDepth->depth[i]->getTH2F())->GetNbinsX();
      int phibins=(meTimeByDepth->depth[i]->getTH2F())->GetNbinsY();
      for (int x=1;x<=etabins;++x)
    for (int y=1;y<=phibins;++y)
      meTimeByDepth->depth[i]->setBinContent(x,y,-1000);
    }

  ib.setCurrentFolder(subdir_+"Distributions_PassedMinBias");
  meEnergyThreshByDepth = new EtaPhiHists();
  meEnergyThreshByDepth->setup(ib,"Above Threshold RecHit Average Energy","GeV");
  meTimeThreshByDepth = new EtaPhiHists();
  meTimeThreshByDepth->setup(ib,"Above Threshold RecHit Average Time","nS");
 // set all average times to -1000 by default (so that they don't show up on plots
  for (unsigned int i=0;i<meTimeThreshByDepth->depth.size();++i)
    {
      (meTimeThreshByDepth->depth[i]->getTH2F())->SetMinimum(-150);
      (meTimeThreshByDepth->depth[i]->getTH2F())->SetMaximum(150);
      int etabins=(meTimeThreshByDepth->depth[i]->getTH2F())->GetNbinsX();
      int phibins=(meTimeThreshByDepth->depth[i]->getTH2F())->GetNbinsY();
      for (int x=1;x<=etabins;++x)
    for (int y=1;y<=phibins;++y)
      meTimeThreshByDepth->depth[i]->setBinContent(x,y,-1000);
    }

  ib.setCurrentFolder(subdir_+"Distributions_AllRecHits/rechit_1D_plots/");
  meHBEnergy_1D=ib.book1D("HB_energy_1D","HB Average Energy Per RecHit;Energy (GeV)",200,-5,15);
  meHEEnergy_1D=ib.book1D("HE_energy_1D","HE Average Energy Per RecHit;Energy (GeV)",200,-5,15);
  meHOEnergy_1D=ib.book1D("HO_energy_1D","HO Average Energy Per RecHit;Energy (GeV)",200,-10,20);
  meHFEnergy_1D=ib.book1D("HF_energy_1D","HF Average Energy Per RecHit;Energy (GeV)",200,-5,15);

  meHBEnergyRMS_1D=ib.book1D("HB_energy_RMS_1D","HB Energy RMS Per RecHit;Energy (GeV)",250,0,5);
  meHEEnergyRMS_1D=ib.book1D("HE_energy_RMS_1D","HE Energy RMS Per RecHit;Energy (GeV)",250,0,5);
  meHOEnergyRMS_1D=ib.book1D("HO_energy_RMS_1D","HO Energy RMS Per RecHit;Energy (GeV)",250,0,5);
  meHFEnergyRMS_1D=ib.book1D("HF_energy_RMS_1D","HF Energy RMS Per RecHit;Energy (GeV)",250,0,5);

  ib.setCurrentFolder(subdir_+"Distributions_PassedMinBias/rechit_1D_plots/");
  meHBEnergyThresh_1D=ib.book1D("HB_energyThresh_1D","HB Average Energy Per RecHit Above Threshold;Energy (GeV)",200,-5,35);
  meHEEnergyThresh_1D=ib.book1D("HE_energyThresh_1D","HE Average Energy Per RecHit Above Threshold;Energy (GeV)",200,-5,35);
  meHOEnergyThresh_1D=ib.book1D("HO_energyThresh_1D","HO Average Energy Per RecHit Above Threshold;Energy (GeV)",300,-10,50);
  meHFEnergyThresh_1D=ib.book1D("HF_energyThresh_1D","HF Average Energy Per RecHit Above Threshold;Energy (GeV)",200,-5,95);

  meHBEnergyRMSThresh_1D=ib.book1D("HB_energy_RMSThresh_1D","HB Energy RMS Per RecHit Above Threshold;Energy (GeV)",200,0,10);
  meHEEnergyRMSThresh_1D=ib.book1D("HE_energy_RMSThresh_1D","HE Energy RMS Per RecHit Above Threshold;Energy (GeV)",200,0,10);
  meHOEnergyRMSThresh_1D=ib.book1D("HO_energy_RMSThresh_1D","HO Energy RMS Per RecHit Above Threshold;Energy (GeV)",200,0,10);
  meHFEnergyRMSThresh_1D=ib.book1D("HF_energy_RMSThresh_1D","HF Energy RMS Per RecHit Above Threshold;Energy (GeV)",200,0,20);

  doProblemCellSetup_ = false;

}

bool HcalRecHitClient::test_enabled ( void  )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 493 of file HcalRecHitClient.cc.

{return true;}

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

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 497 of file HcalRecHitClient.cc.

{
  // rechit quality not used to update channel status yet; see dead cell client for example



} //void HcalRecHitClient::updateChannelStatus

Member Data Documentation

bool HcalRecHitClient::doProblemCellSetup_

private

Definition at line 67 of file HcalRecHitClient.h.

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

EtaPhiHists* HcalRecHitClient::meEnergyByDepth

private

Definition at line 39 of file HcalRecHitClient.h.

Referenced by analyze(), HcalRecHitClient(), setupProblemCells(), and ~HcalRecHitClient().

EtaPhiHists* HcalRecHitClient::meEnergyThreshByDepth

private

Definition at line 40 of file HcalRecHitClient.h.

Referenced by analyze(), HcalRecHitClient(), setupProblemCells(), and ~HcalRecHitClient().

MonitorElement* HcalRecHitClient::meHBEnergy_1D

private

Definition at line 46 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHBEnergyRMS_1D

private

Definition at line 51 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHBEnergyRMSThresh_1D

private

Definition at line 61 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHBEnergyThresh_1D

private

Definition at line 56 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHEEnergy_1D

private

Definition at line 47 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHEEnergyRMS_1D

private

Definition at line 52 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHEEnergyRMSThresh_1D

private

Definition at line 62 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHEEnergyThresh_1D

private

Definition at line 57 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHFEnergy_1D

private

Definition at line 49 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHFEnergyRMS_1D

private

Definition at line 54 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHFEnergyRMSThresh_1D

private

Definition at line 64 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHFEnergyThresh_1D

private

Definition at line 59 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHOEnergy_1D

private

Definition at line 48 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHOEnergyRMS_1D

private

Definition at line 53 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHOEnergyRMSThresh_1D

private

Definition at line 63 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

MonitorElement* HcalRecHitClient::meHOEnergyThresh_1D

private

Definition at line 58 of file HcalRecHitClient.h.

Referenced by analyze(), and setupProblemCells().

EtaPhiHists* HcalRecHitClient::meSqrtSumEnergy2ByDepth

private

Definition at line 43 of file HcalRecHitClient.h.

Referenced by HcalRecHitClient(), and ~HcalRecHitClient().

EtaPhiHists* HcalRecHitClient::meSqrtSumEnergy2ThreshByDepth

private

Definition at line 44 of file HcalRecHitClient.h.

Referenced by HcalRecHitClient(), and ~HcalRecHitClient().

EtaPhiHists* HcalRecHitClient::meTimeByDepth

private

Definition at line 41 of file HcalRecHitClient.h.

Referenced by analyze(), HcalRecHitClient(), setupProblemCells(), and ~HcalRecHitClient().

EtaPhiHists* HcalRecHitClient::meTimeThreshByDepth

private

Definition at line 42 of file HcalRecHitClient.h.

Referenced by analyze(), HcalRecHitClient(), setupProblemCells(), and ~HcalRecHitClient().

int HcalRecHitClient::nevts_

private

Definition at line 37 of file HcalRecHitClient.h.

Referenced by setupProblemCells().