#include <HcalCoarsePedestalClient.h>

Inheritance diagram for HcalCoarsePedestalClient:

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)

	HcalCoarsePedestalClient ()
	Constructors. More...

	HcalCoarsePedestalClient (std::string myname)

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

void	setup (void)

bool	test_enabled (void)

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

	~HcalCoarsePedestalClient ()
	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	setupCoarsePedestal (DQMStore::IBooker &, DQMStore::IGetter &)

Private Attributes
double	ADCDiffThresh_

MonitorElement *	CoarsePedDiff

EtaPhiHists *	CoarsePedestalsByDepth

TH2F *	DatabasePedestalsADCByDepth [4]

bool	doCoarseSetup_

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 HcalCoarsePedestalClient.h.

Constructor & Destructor Documentation

HcalCoarsePedestalClient::HcalCoarsePedestalClient ( )

inline

Constructors.

Definition at line 13 of file HcalCoarsePedestalClient.h.

References HcalBaseDQClient::name_.

13 {name_="";};

HcalBaseDQClient::name_

std::string name_

Definition: HcalBaseDQClient.h:57

HcalCoarsePedestalClient::HcalCoarsePedestalClient ( std::string myname )

Definition at line 18 of file HcalCoarsePedestalClient.cc.

References HcalBaseDQClient::name_.

 {
   name_=myname;
 }

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

Definition at line 23 of file HcalCoarsePedestalClient.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>("CoarsePedestalFolder","CoarsePedestalMonitor_Hcal/"); // CoarsePedestalMonitor_Hcal  
   if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
     subdir_.append("/");
   subdir_=prefixME_+subdir_;
 
   validHtmlOutput_       = ps.getUntrackedParameter<bool>("CoarsePedestal_validHtmlOutput",true);
   cloneME_ = ps.getUntrackedParameter<bool>("cloneME", true);
   badChannelStatusMask_   = ps.getUntrackedParameter<int>("CoarsePedestal_BadChannelStatusMask",
                                                           ps.getUntrackedParameter<int>("BadChannelStatusMask",
                                             ((1<<HcalChannelStatus::HcalCellDead)|(1<<HcalChannelStatus::HcalCellHot)))); // identify channel status values to mask
   
   minerrorrate_ = ps.getUntrackedParameter<double>("CoarsePedestal_minerrorrate",
                            ps.getUntrackedParameter<double>("minerrorrate",0.05));
 
 
   // minevents_ canbe overwritten by monitor task value
   minevents_    = ps.getUntrackedParameter<int>("CoarsePedestal_minevents",
                         ps.getUntrackedParameter<int>("minevents",1));
 
   CoarsePedestalsByDepth=0;
   ProblemCellsByDepth=0;
 
   doCoarseSetup_ = true;
 }

HcalCoarsePedestalClient::~HcalCoarsePedestalClient ( )

Destructor.

Definition at line 342 of file HcalCoarsePedestalClient.cc.

References CoarsePedestalsByDepth, and HcalBaseDQClient::ProblemCellsByDepth.

 {
 
   if ( CoarsePedestalsByDepth ) delete CoarsePedestalsByDepth;
   if ( ProblemCellsByDepth ) delete ProblemCellsByDepth;
 
 }

Member Function Documentation

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

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 56 of file HcalCoarsePedestalClient.cc.

References HcalBaseDQClient::calculateProblems(), gather_cfg::cout, HcalBaseDQClient::debug_, doCoarseSetup_, and setupCoarsePedestal().

 {
   if (debug_>2) std::cout <<"\tHcalCoarsePedestalClient::analyze()"<<std::endl;
   if ( doCoarseSetup_ ) setupCoarsePedestal(ib,ig);
   calculateProblems(ib,ig);
 } // void HcalCoarsePedestalClient::analyze()

void HcalCoarsePedestalClient::beginRun ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 219 of file HcalCoarsePedestalClient.cc.

References HcalBaseDQClient::enoughevents_.

 {
   enoughevents_=false;
 } // void HcalCoarsePedestalClient::beginRun()

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

Definition at line 63 of file HcalCoarsePedestalClient.cc.

 {
  if (debug_>2) std::cout <<"\t\tHcalCoarsePedestalClient::calculateProblems()"<<std::endl;
   //int totalevents=0; // events checked on a channel-by-channel basis
   int etabins=0, phibins=0, zside=0;
   double problemvalue=0;
 
   if (CoarsePedDiff!=0) CoarsePedDiff->Reset();
 
   // Clear away old problems
   if (ProblemCells!=0)
     {
       ProblemCells->Reset();
       (ProblemCells->getTH2F())->SetMaximum(1.05);
       (ProblemCells->getTH2F())->SetMinimum(0.);
       (ProblemCells->getTH2F())->SetOption("colz");
     }
   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.);
       (ProblemCellsByDepth->depth[d]->getTH2F())->SetOption("colz");
     }
     }
 
   // Get histograms that are used in testing
   TH2F* CoarsePedestalsSumByDepth[4];
   TH2F* CoarsePedestalsOccByDepth[4];
 
   std::vector<std::string> name = HcalEtaPhiHistNames();
   bool gothistos=true;
 
   MonitorElement* me;
   for (int i=0;i<4;++i)
     {
       std::string s=subdir_+"CoarsePedestalSumPlots/"+name[i]+"Coarse Pedestal Summed Map";
       me=ig.get(s.c_str());
       if (me==0) 
     {
       gothistos=false;
       if (debug_>0) std::cout <<"<HcalCoarsePedestalClient::calculateProblems> Could not get histogram with name "<<s<<std::endl;
     }
       CoarsePedestalsSumByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, CoarsePedestalsSumByDepth[i], debug_);
 
       s=subdir_+"CoarsePedestalSumPlots/"+name[i]+"Coarse Pedestal Occupancy Map";
       me=ig.get(s.c_str());
       if (me==0) 
     {
       gothistos=false;
       if (debug_>0) std::cout <<"<HcalCoarsePedestalClient::calculateProblems> Could not get histogram with name "<<s<<std::endl;
     }
       CoarsePedestalsOccByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, CoarsePedestalsOccByDepth[i], debug_);
 
     } // for (int i=0;i<4;++i)
 
   if (gothistos==false)
     {
       if (debug_>0) std::cout <<"<HcalCoarsePedestalClient::calculateProblems> Unable to get all necessary histograms to evaluate problem rate"<<std::endl;
       return;
     }
 
   enoughevents_=true;  // Always set this to true, so that pedestal monitoring doesn't hold up 
 
   int numevents=0;
   for (unsigned int d=0;ProblemCellsByDepth!=0 && d<ProblemCellsByDepth->depth.size();++d)
     {
       if (ProblemCellsByDepth->depth[d]==0) continue;
       if (CoarsePedestalsSumByDepth[d]==0 || 
       CoarsePedestalsOccByDepth[d]==0 ||
       DatabasePedestalsADCByDepth[d]==0) continue;
 
       if (CoarsePedestalsByDepth->depth[d]!=0)
     CoarsePedestalsByDepth->depth[d]->Reset();
 
       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)
         {
           if (abs(ieta)>20 && (phi+1)%2==0)
         continue;
           if (abs(ieta)>39 && (phi+1)%4!=3)
         continue;
           numevents=(int)CoarsePedestalsOccByDepth[d]->GetBinContent(eta+1,phi+1);
           if (numevents==0 || numevents<minevents_)
         {
           if (debug_>1)
             std::cout <<"NOT ENOUGH EVENTS for channel ("<<ieta<<", "<<phi+1<<", "<<d+1<<")  numevents = "<<numevents<<"  minevents = "<<minevents_<<std::endl;
           continue;  // insufficient pedestal information available for this channel; continue on?
         }
 
           problemvalue=1.*CoarsePedestalsSumByDepth[d]->GetBinContent(eta+1,phi+1)/numevents;
           CoarsePedestalsByDepth->depth[d]->setBinContent(eta+1,phi+1,problemvalue);
           problemvalue=(problemvalue-DatabasePedestalsADCByDepth[d]->GetBinContent(eta+1,phi+1));
           CoarsePedDiff->Fill(problemvalue);
           problemvalue=fabs(problemvalue);
           // Pedestal status is cumulative (DetDiag pedestals perform resets after each calibration cycle, and save each output)
           // Either a channels pedestal is 'bad' (outside of allowed range) or 'good' (in range)
           problemvalue>ADCDiffThresh_ ? problemvalue=1 : problemvalue=0;
           if (debug_>0 && problemvalue==1)
         std::cout <<"<HcalCoarsePedestalClient> Problem found for channel ("<<ieta<<", "<<phi+1<<", "<<d+1<<")"<<std::endl;
           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 <<"<HcalCoarsePedestalClient::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);
   FillUnphysicalHEHFBins(*CoarsePedestalsByDepth);
   return;
 }

void HcalCoarsePedestalClient::cleanup ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 296 of file HcalCoarsePedestalClient.cc.

296 {}

void HcalCoarsePedestalClient::endJob ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 217 of file HcalCoarsePedestalClient.cc.

217 {}

bool HcalCoarsePedestalClient::hasErrors_Temp ( void )

virtual

Reimplemented from HcalBaseDQClient.

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

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 332 of file HcalCoarsePedestalClient.cc.

332 {return false;}

bool HcalCoarsePedestalClient::hasWarnings_Temp ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 331 of file HcalCoarsePedestalClient.cc.

331 {return false;}

void HcalCoarsePedestalClient::setup ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 295 of file HcalCoarsePedestalClient.cc.

295 {}

void HcalCoarsePedestalClient::setupCoarsePedestal	(	DQMStore::IBooker &	ib,
		DQMStore::IGetter &	ig
	)

private

Definition at line 224 of file HcalCoarsePedestalClient.cc.

References ADCDiffThresh_, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), HcalBaseDQClient::cloneME_, CoarsePedDiff, CoarsePedestalsByDepth, gather_cfg::cout, DatabasePedestalsADCByDepth, HcalBaseDQClient::debug_, EtaPhiHists::depth, doCoarseSetup_, DQMStore::IGetter::get(), MonitorElement::getFloatValue(), MonitorElement::getIntValue(), MonitorElement::getName(), HcalEtaPhiHistNames(), i, HcalBaseDQClient::minevents_, HcalBaseDQClient::name(), nevts_, HcalBaseDQClient::prefixME_, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, HcalBaseDQClient::problemnames_, alignCSCRings::s, DQMStore::IBooker::setCurrentFolder(), EtaPhiHists::setup(), AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQClient::subdir_, and groupFilesInBlocks::temp.

Referenced by analyze().

 {
 
   ib.setCurrentFolder(subdir_);
   CoarsePedestalsByDepth = new EtaPhiHists();
   CoarsePedestalsByDepth->setup(ib," Coarse Pedestal Map");
 
   CoarsePedDiff=ib.book1D("PedRefDiff","(Pedestal-Reference)",200,-10,10);
 
   problemnames_.clear();
   ProblemCells=ib.book2D(" ProblemCoarsePedestals",
                  " Problem Coarse Pedestal 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_coarsepedestals");
   ProblemCellsByDepth = new EtaPhiHists();
   ProblemCellsByDepth->setup(ib," Problem Coarse Pedestal Rate");
   for (unsigned int i=0; i<ProblemCellsByDepth->depth.size();++i)
     problemnames_.push_back(ProblemCellsByDepth->depth[i]->getName());
   nevts_=0;
   
   std::vector<std::string> name = HcalEtaPhiHistNames();
   MonitorElement* me;
   for (int i=0;i<4;++i)
     {
       std::string s=prefixME_+"HcalInfo/PedestalsFromCondDB/"+name[i]+"ADC Pedestals From Conditions DB";
       me=ig.get(s.c_str());
       if (me==0) 
     {
       if (debug_>0) std::cout <<"<HcalCoarsePedestalClient::beginRun> Could not get histogram with name "<<s<<std::endl;
     }
       DatabasePedestalsADCByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, DatabasePedestalsADCByDepth[i], debug_);
     }
   std::string s=subdir_+"CoarsePedestal_parameters/ADCdiff_Problem_Threshold";
   me=ig.get(s.c_str());
   if (me==0)
     {
       if (debug_>0) 
     std::cout <<"<HcalCoarsePedestalClient::beginRun> Could not get value with name "<<s<<std::endl;
     }
   else 
     ADCDiffThresh_=me->getFloatValue();
   s=subdir_+"CoarsePedestal_parameters/minEventsNeededForPedestalCalculation";
 
   me=ig.get(s.c_str());
   int temp = 0;
   if (me==0) 
     {
       if (debug_>0)
     {
       std::cout <<"<HcalCoarsePedestalClient::beginRun> Could not get value with name "<<s<<"\n\t  Continuing on using default 'minevents' value of "<<minevents_<<std::endl;
     }
     }
   else 
     temp=me->getIntValue();
   if (temp>minevents_)
     {
       if (debug_>0)
     std::cout <<"<HcalCoarsePedestalClient::beginRun>  Specified client 'minevents' value of "<<minevents_<<"  is less than minimum task 'minevents' value of "<<temp<<"\n\t  Setting client 'minevents' to "<<temp<<std::endl;
       minevents_=temp;
     }
 
   doCoarseSetup_ = false;
 
 }

bool HcalCoarsePedestalClient::test_enabled ( void )

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 333 of file HcalCoarsePedestalClient.cc.

333 {return true;}

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

virtual

Reimplemented from HcalBaseDQClient.

Definition at line 336 of file HcalCoarsePedestalClient.cc.

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