|
|
|
#include <HcalTrigPrimClient.h>
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) |
| HcalTrigPrimClient (std::string myname) | |
| HcalTrigPrimClient () | |
| Constructors. | |
| HcalTrigPrimClient (std::string myname, const edm::ParameterSet &ps) | |
| void | setup (void) |
| bool | test_enabled (void) |
| void | updateChannelStatus (std::map< HcalDetId, unsigned int > &myqual) |
| ~HcalTrigPrimClient () | |
| Destructor. | |
Private Attributes | |
| int | nevts_ |
| EtaPhiHists * | ProblemsByDepthNZS_ |
| EtaPhiHists * | ProblemsByDepthZS_ |
Definition at line 8 of file HcalTrigPrimClient.h.
| HcalTrigPrimClient::HcalTrigPrimClient | ( | ) | [inline] |
Constructors.
Definition at line 13 of file HcalTrigPrimClient.h.
References HcalBaseDQClient::name_.
{name_="";};
| HcalTrigPrimClient::HcalTrigPrimClient | ( | std::string | myname | ) |
Definition at line 20 of file HcalTrigPrimClient.cc.
References HcalBaseDQClient::name_.
{
name_=myname;
}
| HcalTrigPrimClient::HcalTrigPrimClient | ( | std::string | myname, |
| const edm::ParameterSet & | ps | ||
| ) |
Definition at line 25 of file HcalTrigPrimClient.cc.
References HcalBaseDQClient::badChannelStatusMask_, HcalBaseDQClient::cloneME_, HcalBaseDQClient::debug_, HcalBaseDQClient::enableCleanup_, edm::ParameterSet::getUntrackedParameter(), HcalBaseDQClient::minerrorrate_, HcalBaseDQClient::minevents_, HcalBaseDQClient::name_, HcalBaseDQClient::Online_, HcalBaseDQClient::prefixME_, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, 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>("TrigPrimFolder","TrigPrimMonitor_Hcal/"); // TrigPrimMonitor
if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
subdir_.append("/");
subdir_=prefixME_+subdir_;
validHtmlOutput_ = ps.getUntrackedParameter<bool>("TrigPrim_validHtmlOutput",true);
cloneME_ = ps.getUntrackedParameter<bool>("cloneME", true);
badChannelStatusMask_ = ps.getUntrackedParameter<int>("TrigPrim_BadChannelStatusMask",
ps.getUntrackedParameter<int>("BadChannelStatusMask",0));
// Set error rate to 1%, not (earlier) value of 0.1% -- Jeff, 11 May 2010
minerrorrate_ = ps.getUntrackedParameter<double>("TrigPrim_minerrorrate",
ps.getUntrackedParameter<double>("minerrorrate",0.01));
minevents_ = ps.getUntrackedParameter<int>("TrigPrim_minevents",
ps.getUntrackedParameter<int>("minevents",1));
Online_ = ps.getUntrackedParameter<bool>("online",false);
ProblemCells=0;
ProblemCellsByDepth=0;
}
| HcalTrigPrimClient::~HcalTrigPrimClient | ( | ) |
| void HcalTrigPrimClient::analyze | ( | void | ) | [virtual] |
Reimplemented from HcalBaseDQClient.
Definition at line 53 of file HcalTrigPrimClient.cc.
References calculateProblems(), gather_cfg::cout, and HcalBaseDQClient::debug_.
Referenced by endRun().
{
if (debug_>2) std::cout <<"\tHcalTrigPrimClient::analyze()"<<std::endl;
calculateProblems();
}
| void HcalTrigPrimClient::beginJob | ( | void | ) | [virtual] |
Reimplemented from HcalBaseDQClient.
Definition at line 311 of file HcalTrigPrimClient.cc.
References gather_cfg::cout, HcalBaseDQClient::debug_, HcalBaseDQClient::dqmStore_, cppFunctionSkipper::operator, and DQMStore::showDirStructure().
{
dqmStore_ = edm::Service<DQMStore>().operator->();
if (debug_>0)
{
std::cout <<"<HcalTrigPrimClient::beginJob()> Displaying dqmStore directory structure:"<<std::endl;
dqmStore_->showDirStructure();
}
}
| void HcalTrigPrimClient::beginRun | ( | void | ) | [virtual] |
Reimplemented from HcalBaseDQClient.
Definition at line 322 of file HcalTrigPrimClient.cc.
References DQMStore::book2D(), gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, HcalBaseDQClient::dqmStore_, HcalBaseDQClient::enoughevents_, MonitorElement::getName(), i, nevts_, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, HcalBaseDQClient::problemnames_, ProblemsByDepthNZS_, ProblemsByDepthZS_, DQMStore::setCurrentFolder(), EtaPhiHists::setup(), and HcalBaseDQClient::subdir_.
{
enoughevents_=false;
if (!dqmStore_)
{
if (debug_>0) std::cout <<"<HcalTrigPrimClient::beginRun> dqmStore does not exist!"<<std::endl;
return;
}
dqmStore_->setCurrentFolder(subdir_);
problemnames_.clear();
// Put the appropriate name of your problem summary here
ProblemCells=dqmStore_->book2D(" ProblemTriggerPrimitives",
" Problem Trigger Primitive 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_triggerprimitives");
ProblemCellsByDepth = new EtaPhiHists();
ProblemCellsByDepth->setup(dqmStore_," Problem Trigger Primitive Rate");
for (unsigned int i=0; i<ProblemCellsByDepth->depth.size();++i)
problemnames_.push_back(ProblemCellsByDepth->depth[i]->getName());
nevts_=0;
dqmStore_->setCurrentFolder(subdir_+"problem_ZS");
ProblemsByDepthZS_ = new EtaPhiHists();
ProblemsByDepthZS_->setup(dqmStore_,"ZS Problem Trigger Primitive Rate");
dqmStore_->setCurrentFolder(subdir_+"problem_NZS");
ProblemsByDepthNZS_ = new EtaPhiHists();
ProblemsByDepthNZS_->setup(dqmStore_,"NZS Problem Trigger Primitive Rate");
}
| void HcalTrigPrimClient::calculateProblems | ( | void | ) | [virtual] |
Reimplemented from HcalBaseDQClient.
Definition at line 59 of file HcalTrigPrimClient.cc.
References abs, HcalBaseDQClient::cloneME_, gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, HcalBaseDQClient::dqmStore_, HcalBaseDQClient::enoughevents_, eta(), MonitorElement::Fill(), HcalObjRepresent::FillUnphysicalHEHFBins(), DQMStore::get(), MonitorElement::getBinContent(), MonitorElement::getTH2F(), HcalEtaPhiHistNames(), i, HcalBaseDQClient::minevents_, HcalBaseDQClient::name(), phi, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, ProblemsByDepthNZS_, ProblemsByDepthZS_, MonitorElement::Reset(), MonitorElement::setBinContent(), and HcalBaseDQClient::subdir_.
Referenced by analyze().
{
if (debug_>2) std::cout <<"\t\tHcalTrigPrimClient::calculateProblems()"<<std::endl;
if(!dqmStore_) return;
double totalevents=0;
int etabins=0, phibins=0;
double problemvalue=0;
enoughevents_=false; // assume we lack sufficient events until proven otherwise
// 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<ProblemsByDepthZS_->depth.size();++d)
{
if (ProblemsByDepthZS_->depth[d]!=0)
{
ProblemsByDepthZS_->depth[d]->Reset();
(ProblemsByDepthZS_->depth[d]->getTH2F())->SetMaximum(1.05);
(ProblemsByDepthZS_->depth[d]->getTH2F())->SetMinimum(0.);
}
}
for (unsigned int d=0;d<ProblemsByDepthNZS_->depth.size();++d)
{
if (ProblemsByDepthNZS_->depth[d]!=0)
{
ProblemsByDepthNZS_->depth[d]->Reset();
(ProblemsByDepthNZS_->depth[d]->getTH2F())->SetMaximum(1.05);
(ProblemsByDepthNZS_->depth[d]->getTH2F())->SetMinimum(0.);
}
}
// Get histograms that are used in testing
// currently none used,
std::vector<std::string> name = HcalEtaPhiHistNames();
/*
// This is a sample of how to get a histogram from the task that can then be used for evaluation purposes
*/
MonitorElement* me;
TH2F *goodZS=0;
TH2F *badZS=0;
TH2F* goodNZS=0;
TH2F* badNZS=0;
me=dqmStore_->get(subdir_+"Good TPs_ZS");
if (!me && debug_>0)
std::cout <<"<HcalTrigPrimClient::calculateProblems> Could not get histogram named '"<<subdir_<<"Good TPs_ZS'"<<std::endl;
else goodZS = HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, goodZS, debug_);
me=dqmStore_->get(subdir_+"Bad TPs_ZS");
if (!me && debug_>0)
std::cout <<"<HcalTrigPrimClient::calculateProblems> Could not get histogram named '"<<subdir_<<"Bad TPs_ZS'"<<std::endl;
else badZS = HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, badZS, debug_);
me=dqmStore_->get(subdir_+"noZS/Good TPs_noZS");
if (!me && debug_>0)
std::cout <<"<HcalTrigPrimClient::calculateProblems> Could not get histogram named '"<<subdir_<<"noZS/Good TPs_noZS'"<<std::endl;
else goodNZS = HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, goodNZS, debug_);
me=dqmStore_->get(subdir_+"noZS/Bad TPs_noZS");
if (!me && debug_>0)
std::cout <<"<HcalTrigPrimClient::calculateProblems> Could not get histogram named '"<<subdir_<<"noZS/Bad TPs_noZS'"<<std::endl;
else badNZS = HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, badNZS, debug_);
// get bin info from good histograms
if (goodZS!=0)
{
etabins=goodZS->GetNbinsX();
phibins=goodZS->GetNbinsY();
totalevents=goodNZS->GetBinContent(0);
}
else if (goodNZS!=0)
{
etabins=goodNZS->GetNbinsX();
phibins=goodNZS->GetNbinsY();
totalevents=goodNZS->GetBinContent(0);
}
if (totalevents<minevents_)
{
enoughevents_=false;
if (debug_>2) std::cout <<"<HcalTrigPrimClient::calculateProblems()> Not enough events! events = "<<totalevents<<" minimum required = "<<minevents_<<std::endl;
return;
}
enoughevents_=true;
// got good and bad histograms; now let's loop over them
int ieta=-99, iphi=-99;
int badvalZS=0, goodvalZS=0;
int badvalNZS=0, goodvalNZS=0;
for (int eta=1;eta<=etabins;++eta)
{
ieta=eta-33; // Patrick's eta-phi maps starts at ieta=-32
for (int phi=1;phi<=phibins;++phi)
{
badvalZS=0, goodvalZS=0;
badvalNZS=0, goodvalNZS=0;
iphi=phi;
if (badZS!=0) badvalZS=(int)badZS->GetBinContent(eta,phi);
if (badNZS!=0) badvalNZS=(int)badNZS->GetBinContent(eta,phi);
if (badvalZS+badvalNZS==0) continue;
if (goodZS!=0) goodvalZS=(int)goodZS->GetBinContent(eta,phi);
if (goodNZS!=0) goodvalNZS=(int)goodNZS->GetBinContent(eta,phi);
if (badvalNZS>0)
{
problemvalue=badvalNZS*1./(badvalNZS+goodvalNZS);
if (abs(ieta)<29)
{
ProblemsByDepthNZS_->depth[0]->Fill(ieta,iphi,problemvalue);
if (abs(ieta)==28) // TP 28 spans towers 28 and 29
ProblemsByDepthNZS_->depth[0]->Fill(ieta+abs(ieta)/ieta,iphi,problemvalue);
}
else // HF
{
/* HF TPs:
TP29 = ieta 29-31
TP30 = ieta 32-34
TP31 = ieta 35-37
TP38 = ieta 38-41
iphi = 1, 5, ..., with 1 covering iphi=1 and iphi=71, etc.
*/
int newieta=-99;
for (int i=0;i<3;++i)
{
newieta=i+29+3*(abs(ieta)-29)+1; // shift values by 1 for HF in EtaPhiHistsplot
if (ieta<0) newieta*=-1;
ProblemsByDepthNZS_->depth[0]->Fill(newieta,iphi,problemvalue);
ProblemsByDepthNZS_->depth[0]->Fill(newieta,(iphi-2+72)%72,problemvalue);
}
if (abs(ieta)==32)
{
ProblemsByDepthNZS_->depth[0]->Fill(42*abs(ieta)/ieta,iphi,problemvalue);
ProblemsByDepthNZS_->depth[0]->Fill(newieta,(iphi-2+72)%72,problemvalue);
}
}
} // errors found in NZS;
if (badvalZS>0)
{
problemvalue=badvalZS*1./(badvalZS+goodvalZS);
if (abs(ieta)<29) // Make special case for ieta=16 (HB/HE overlap?)
{
ProblemsByDepthZS_->depth[0]->Fill(ieta,iphi,problemvalue);
if (abs(ieta)==28) // TP 28 spans towers 28 and 29
ProblemsByDepthZS_->depth[0]->Fill(ieta+abs(ieta)/ieta,iphi,problemvalue);
}
else
{
int newieta=-99;
for (int i=0;i<3;++i)
{
newieta=i+29+3*(abs(ieta)-29)+1; // shift values by 1 for HF in EtaPhiHistsplot
if (ieta<0) newieta*=-1;
ProblemsByDepthZS_->depth[0]->Fill(newieta,iphi,problemvalue);
ProblemsByDepthZS_->depth[0]->Fill(newieta,(iphi-2+72)%72,problemvalue);
}
if (abs(ieta)==32)
{
ProblemsByDepthZS_->depth[0]->Fill(42*abs(ieta)/ieta,iphi,problemvalue);
ProblemsByDepthZS_->depth[0]->Fill(42*abs(ieta)/ieta,(iphi-2+72)%72,problemvalue);
}
}
} // errors found in ZS
if (badvalZS>0 || badvalNZS>0)
{
// Fill overall problem histograms with sum from both ZS & NZS, or ZS only?
//problemvalue=(badvalZS+badvalNZS)*1./(badvalZS+badvalNZS+goodvalZS+goodvalNZS);
// Update on 8 March -- NZS shows lots of errors; let's not include that in problem cells just yet
if (badvalZS==0) continue;
problemvalue=(badvalZS*1.)/(badvalZS+goodvalZS);
if (abs(ieta)<29) // Make special case for ieta=16 (HB/HE overlap?)
{
ProblemCellsByDepth->depth[0]->Fill(ieta,iphi,problemvalue);
ProblemCells->Fill(ieta,iphi,problemvalue);
if (abs(ieta)==28) // TP 28 spans towers 28 and 29
{
ProblemCellsByDepth->depth[0]->Fill(ieta+abs(ieta)/ieta,iphi,problemvalue);
ProblemCells->Fill(ieta+abs(ieta)/ieta,iphi,problemvalue);
}
}
else
{
int newieta=-99;
int newiphi=(iphi-2+72)%72; // forward triggers combine two HF cells; *subtract* 2 to the iphi, and allow wraparound
// FIXME:
// iphi seems to start at 1 in Patrick's plots, continues mod 4;
// adjust in far-forward region, where cells start at iphi=3? Check with Patrick.
for (int i=0;i<3;++i)
{
newieta=i+29+3*(abs(ieta)-29)+1; // shift values by 1 for HF in EtaPhiHistsplot
if (ieta<0) newieta*=-1;
ProblemCellsByDepth->depth[0]->Fill(newieta,iphi,problemvalue);
ProblemCells->Fill(newieta,iphi,problemvalue);
ProblemCellsByDepth->depth[0]->Fill(newieta,newiphi,problemvalue);
ProblemCells->Fill(newieta,newiphi,problemvalue);
}
if (abs(ieta)==32)
{
ProblemCellsByDepth->depth[0]->Fill(42*abs(ieta)/ieta,iphi,problemvalue);
ProblemCells->Fill(42*abs(ieta)/ieta,iphi,problemvalue);
ProblemCellsByDepth->depth[0]->Fill(42*abs(ieta)/ieta,newiphi,problemvalue);
ProblemCells->Fill(42*abs(ieta)/ieta,newiphi,problemvalue);
}
}
}
}
} // for (int eta=1;eta<etabins;++eta)
if (ProblemCells==0)
{
if (debug_>0) std::cout <<"<HcalTrigPrimClient::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(*ProblemsByDepthZS_);
FillUnphysicalHEHFBins(*ProblemsByDepthNZS_);
FillUnphysicalHEHFBins(ProblemCells);
return;
}
| void HcalTrigPrimClient::cleanup | ( | void | ) | [virtual] |
| void HcalTrigPrimClient::endJob | ( | void | ) | [virtual] |
| void HcalTrigPrimClient::endRun | ( | void | ) | [virtual] |
Reimplemented from HcalBaseDQClient.
Definition at line 356 of file HcalTrigPrimClient.cc.
References analyze().
{analyze();}
| bool HcalTrigPrimClient::hasErrors_Temp | ( | void | ) | [virtual] |
Reimplemented from HcalBaseDQClient.
Definition at line 361 of file HcalTrigPrimClient.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 <<"<HcalTrigPrimClient::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 HcalTrigPrimClient::hasOther_Temp | ( | void | ) | [virtual] |
Reimplemented from HcalBaseDQClient.
Definition at line 395 of file HcalTrigPrimClient.cc.
{return false;}
| bool HcalTrigPrimClient::hasWarnings_Temp | ( | void | ) | [virtual] |
Reimplemented from HcalBaseDQClient.
Definition at line 394 of file HcalTrigPrimClient.cc.
{return false;}
| void HcalTrigPrimClient::setup | ( | void | ) | [virtual] |
| bool HcalTrigPrimClient::test_enabled | ( | void | ) | [virtual] |
Reimplemented from HcalBaseDQClient.
Definition at line 396 of file HcalTrigPrimClient.cc.
{return true;}
| void HcalTrigPrimClient::updateChannelStatus | ( | std::map< HcalDetId, unsigned int > & | myqual | ) | [virtual] |
Reimplemented from HcalBaseDQClient.
Definition at line 399 of file HcalTrigPrimClient.cc.
{
// This gets called by HcalMonitorClient
// trigger primitives don't yet contribute to channel status (though they could...)
// see dead or hot cell code for an example
} //void HcalTrigPrimClient::updateChannelStatus
int HcalTrigPrimClient::nevts_ [private] |
Definition at line 36 of file HcalTrigPrimClient.h.
Referenced by beginRun().
Definition at line 39 of file HcalTrigPrimClient.h.
Referenced by beginRun(), and calculateProblems().
Definition at line 38 of file HcalTrigPrimClient.h.
Referenced by beginRun(), and calculateProblems().
1.7.3