|
|
|
#include <HcalRecHitsClient.h>
Public Member Functions | |
| virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
| virtual void | beginJob (void) |
| virtual void | beginRun (const edm::Run &run, const edm::EventSetup &c) |
| virtual void | endJob () |
| virtual void | endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c) |
| virtual void | endRun (const edm::Run &run, const edm::EventSetup &c) |
| HcalRecHitsClient (const edm::ParameterSet &) | |
| int | HcalRecHitsEndjob (const std::vector< MonitorElement * > &hcalMEs) |
| virtual void | runClient_ () |
| virtual | ~HcalRecHitsClient () |
Private Attributes | |
| edm::ParameterSet | conf_ |
| DQMStore * | dbe_ |
| bool | debug_ |
| std::string | dirName_ |
| std::string | dirNameJet_ |
| std::string | dirNameMET_ |
| std::string | outputFile_ |
| bool | verbose_ |
Definition at line 36 of file HcalRecHitsClient.h.
| HcalRecHitsClient::HcalRecHitsClient | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 12 of file HcalRecHitsClient.cc.
References dbe_, debug_, dirName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), cmsCodeRules::cppFunctionSkipper::operator, outputFile_, and verbose_.
:conf_(iConfig) { outputFile_ = iConfig.getUntrackedParameter<std::string>("outputFile", "myfile.root"); dbe_ = edm::Service<DQMStore>().operator->(); if (!dbe_) { edm::LogError("HcalRecHitsClient") << "unable to get DQMStore service, upshot is no client histograms will be made"; } if(iConfig.getUntrackedParameter<bool>("DQMStore", false)) { if(dbe_) dbe_->setVerbose(0); } debug_ = false; verbose_ = false; dirName_=iConfig.getParameter<std::string>("DQMDirName"); if(dbe_) dbe_->setCurrentFolder(dirName_); }
| HcalRecHitsClient::~HcalRecHitsClient | ( | ) | [virtual] |
Definition at line 34 of file HcalRecHitsClient.cc.
{
}
| void HcalRecHitsClient::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
| void HcalRecHitsClient::beginJob | ( | void | ) | [virtual] |
| void HcalRecHitsClient::beginRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | c | ||
| ) | [virtual] |
| void HcalRecHitsClient::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 45 of file HcalRecHitsClient.cc.
References dbe_, outputFile_, and DQMStore::save().
{
if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
}
| void HcalRecHitsClient::endLuminosityBlock | ( | const edm::LuminosityBlock & | lumiSeg, |
| const edm::EventSetup & | c | ||
| ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 67 of file HcalRecHitsClient.cc.
{
// runClient_();
}
| void HcalRecHitsClient::endRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | c | ||
| ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 56 of file HcalRecHitsClient.cc.
References runClient_().
{
runClient_();
}
| int HcalRecHitsClient::HcalRecHitsEndjob | ( | const std::vector< MonitorElement * > & | hcalMEs | ) |
Definition at line 109 of file HcalRecHitsClient.cc.
References gather_cfg::cout, MonitorElement::Fill(), MonitorElement::getBinContent(), MonitorElement::getBinError(), MonitorElement::getEntries(), edm::getName(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), i, getHLTprescales::index, j, MonitorElement::setBinContent(), and verbose_.
Referenced by runClient_().
{
int useAllHistos = 0, subdet =5;
// for ZS ...
// MonitorElement* emap_min_ME =0;
MonitorElement* ZS_HO=0, *ZS_seqHO=0;
MonitorElement* ZS_HB1=0, *ZS_seqHB1=0, *ZS_HB2=0, *ZS_seqHB2=0;
MonitorElement* ZS_HF1=0, *ZS_seqHF1=0, *ZS_HF2=0, *ZS_seqHF2=0;
MonitorElement* ZS_HE1=0, *ZS_seqHE1=0, *ZS_HE2=0, *ZS_seqHE2=0, *ZS_HE3=0, *ZS_seqHE3=0;
MonitorElement* map_depth1 =0, *map_depth2 =0, *map_depth3 =0, *map_depth4 =0;
// others
MonitorElement* Nhf=0;
MonitorElement* emap_depth1 =0, *emap_depth2 =0, *emap_depth3 =0, *emap_depth4 =0;
MonitorElement* occupancy_seqHB1 =0, *occupancy_seqHB2 =0;
MonitorElement* occupancy_seqHE1 =0, *occupancy_seqHE2 =0, *occupancy_seqHE3 =0;
MonitorElement* occupancy_seqHF1 =0, *occupancy_seqHF2 =0;
MonitorElement* occupancy_seqHO =0;
MonitorElement* emean_seqHB1 =0, *emean_seqHB2 =0;
MonitorElement* emean_seqHE1 =0, *emean_seqHE2 =0, *emean_seqHE3 =0;
MonitorElement* emean_seqHF1 =0, *emean_seqHF2 =0;
MonitorElement* emean_seqHO =0;
MonitorElement* RMS_seq_HB1 =0, *RMS_seq_HB2 =0;
MonitorElement* RMS_seq_HE1 =0, *RMS_seq_HE2 =0, *RMS_seq_HE3 =0;
MonitorElement* RMS_seq_HF1 =0, *RMS_seq_HF2 =0;
MonitorElement* RMS_seq_HO =0;
MonitorElement *occupancy_map_HO =0;
MonitorElement* occupancy_map_HB1 =0, *occupancy_map_HB2 =0;
MonitorElement* occupancy_map_HF1 =0, *occupancy_map_HF2 =0;
MonitorElement* occupancy_map_HE1 =0, *occupancy_map_HE2 =0, *occupancy_map_HE3 =0;
MonitorElement* emean_vs_ieta_HB1 =0, *emean_vs_ieta_HB2 =0;
MonitorElement* emean_vs_ieta_HE1 =0, *emean_vs_ieta_HE2 =0, *emean_vs_ieta_HE3 =0;
MonitorElement* emean_vs_ieta_HF1 =0, *emean_vs_ieta_HF2 =0;
MonitorElement* emean_vs_ieta_HO =0;
MonitorElement* RMS_vs_ieta_HB1 =0, *RMS_vs_ieta_HB2 =0;
MonitorElement* RMS_vs_ieta_HE1 =0, *RMS_vs_ieta_HE2 =0, *RMS_vs_ieta_HE3 =0;
MonitorElement* RMS_vs_ieta_HF1 =0, *RMS_vs_ieta_HF2 =0;
MonitorElement* RMS_vs_ieta_HO =0;
MonitorElement* occupancy_vs_ieta_HB1 =0, *occupancy_vs_ieta_HB2 =0;
MonitorElement* occupancy_vs_ieta_HE1 =0, *occupancy_vs_ieta_HE2 =0, *occupancy_vs_ieta_HE3 =0;
MonitorElement* occupancy_vs_ieta_HF1 =0, *occupancy_vs_ieta_HF2 =0;
MonitorElement* occupancy_vs_ieta_HO =0;
MonitorElement* RecHit_StatusWord_HB =0, *RecHit_StatusWord_HE=0, *RecHit_StatusWord_HO =0, *RecHit_StatusWord_HF =0, *RecHit_StatusWord_HF67 =0;
MonitorElement* RecHit_Aux_StatusWord_HB =0, *RecHit_Aux_StatusWord_HE=0, *RecHit_Aux_StatusWord_HO =0, *RecHit_Aux_StatusWord_HF =0;
for(unsigned int ih=0; ih<hcalMEs.size(); ih++){
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_map_depth1") ==0 ){
useAllHistos =1; subdet =6;
}
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_vs_ieta_HB1") ==0 ){
useAllHistos =1;
}
// if( strcmp(hcalMEs[ih]->getName().c_str(), "emap_min_ME") ==0 ){ emap_min_ME = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_simple1D_HO") ==0 ){ ZS_HO = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_sequential1D_HO") ==0 ){ ZS_seqHO = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_simple1D_HB1") ==0 ){ ZS_HB1 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_sequential1D_HB1") ==0 ){ ZS_seqHB1 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_simple1D_HB2") ==0 ){ ZS_HB2 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_sequential1D_HB2") ==0 ){ ZS_seqHB2 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_simple1D_HF1") ==0 ){ ZS_HF1 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_sequential1D_HF1") ==0 ){ ZS_seqHF1 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_simple1D_HF2") ==0 ){ ZS_HF2 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_sequential1D_HF2") ==0 ){ ZS_seqHF2 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_simple1D_HE1") ==0 ){ ZS_HE1 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_sequential1D_HE1") ==0 ){ ZS_seqHE1 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_simple1D_HE2") ==0 ){ ZS_HE2 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_sequential1D_HE2") ==0 ){ ZS_seqHE2 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_simple1D_HE3") ==0 ){ ZS_HE3 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_sequential1D_HE3") ==0 ){ ZS_seqHE3 = hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_map_depth1") ==0 ){ map_depth1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_map_depth2") ==0 ){ map_depth2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_map_depth3") ==0 ){ map_depth3= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "ZSmin_map_depth4") ==0 ){ map_depth4= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "N_HF") ==0 ){ Nhf= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emap_depth1") ==0 ){ emap_depth1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emap_depth2") ==0 ){ emap_depth2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emap_depth3") ==0 ){ emap_depth3= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emap_depth4") ==0 ){ emap_depth4= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occ_sequential1D_HB1") ==0 ){ occupancy_seqHB1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occ_sequential1D_HB2") ==0 ){ occupancy_seqHB2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occ_sequential1D_HE1") ==0 ){ occupancy_seqHE1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occ_sequential1D_HE2") ==0 ){ occupancy_seqHE2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occ_sequential1D_HE3") ==0 ){ occupancy_seqHE3= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occ_sequential1D_HF1") ==0 ){ occupancy_seqHF1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occ_sequential1D_HF2") ==0 ){ occupancy_seqHF2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occ_sequential1D_HO") ==0 ){ occupancy_seqHO= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_seq_HB1") ==0 ){ emean_seqHB1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_seq_HB2") ==0 ){ emean_seqHB2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_seq_HE1") ==0 ){ emean_seqHE1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_seq_HE2") ==0 ){ emean_seqHE2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_seq_HE3") ==0 ){ emean_seqHE3= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_seq_HF1") ==0 ){ emean_seqHF1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_seq_HF2") ==0 ){ emean_seqHF2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_seq_HO") ==0 ){ emean_seqHO= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_seq_HB1") ==0 ){ RMS_seq_HB1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_seq_HB2") ==0 ){ RMS_seq_HB2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_seq_HE1") ==0 ){ RMS_seq_HE1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_seq_HE2") ==0 ){ RMS_seq_HE2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_seq_HE3") ==0 ){ RMS_seq_HE3= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_seq_HF1") ==0 ){ RMS_seq_HF1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_seq_HF2") ==0 ){ RMS_seq_HF2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_seq_HO") ==0 ){ RMS_seq_HO= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_map_HB1") ==0 ){ occupancy_map_HB1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_map_HB2") ==0 ){ occupancy_map_HB2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_map_HE1") ==0 ){ occupancy_map_HE1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_map_HE2") ==0 ){ occupancy_map_HE2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_map_HE3") ==0 ){ occupancy_map_HE3= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_map_HF1") ==0 ){ occupancy_map_HF1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_map_HF2") ==0 ){ occupancy_map_HF2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_map_HO") ==0 ){ occupancy_map_HO= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_vs_ieta_HB1") ==0 ){ emean_vs_ieta_HB1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_vs_ieta_HB2") ==0 ){ emean_vs_ieta_HB2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_vs_ieta_HE1") ==0 ){ emean_vs_ieta_HE1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_vs_ieta_HE2") ==0 ){ emean_vs_ieta_HE2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_vs_ieta_HE3") ==0 ){ emean_vs_ieta_HE3= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_vs_ieta_HF1") ==0 ){ emean_vs_ieta_HF1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_vs_ieta_HF2") ==0 ){ emean_vs_ieta_HF2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "emean_vs_ieta_HO") ==0 ){ emean_vs_ieta_HO= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_vs_ieta_HB1") ==0 ){ RMS_vs_ieta_HB1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_vs_ieta_HB2") ==0 ){ RMS_vs_ieta_HB2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_vs_ieta_HE1") ==0 ){ RMS_vs_ieta_HE1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_vs_ieta_HE2") ==0 ){ RMS_vs_ieta_HE2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_vs_ieta_HE3") ==0 ){ RMS_vs_ieta_HE3= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_vs_ieta_HF1") ==0 ){ RMS_vs_ieta_HF1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_vs_ieta_HF2") ==0 ){ RMS_vs_ieta_HF2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "RMS_vs_ieta_HO") ==0 ){ RMS_vs_ieta_HO= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_vs_ieta_HB1") ==0 ){ occupancy_vs_ieta_HB1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_vs_ieta_HB2") ==0 ){ occupancy_vs_ieta_HB2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_vs_ieta_HE1") ==0 ){ occupancy_vs_ieta_HE1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_vs_ieta_HE2") ==0 ){ occupancy_vs_ieta_HE2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_vs_ieta_HE3") ==0 ){ occupancy_vs_ieta_HE3= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_vs_ieta_HF1") ==0 ){ occupancy_vs_ieta_HF1= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_vs_ieta_HF2") ==0 ){ occupancy_vs_ieta_HF2= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "occupancy_vs_ieta_HO") ==0 ){ occupancy_vs_ieta_HO= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "HcalRecHitTask_RecHit_StatusWord_HB") ==0 ){ RecHit_StatusWord_HB= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "HcalRecHitTask_RecHit_StatusWord_HE") ==0 ){ RecHit_StatusWord_HE= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "HcalRecHitTask_RecHit_StatusWord_HO") ==0 ){ RecHit_StatusWord_HO= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "HcalRecHitTask_RecHit_StatusWord_HF") ==0 ){ RecHit_StatusWord_HF= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "HcalRecHitTask_RecHit_StatusWord_HF67") ==0 ){ RecHit_StatusWord_HF67= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "HcalRecHitTask_RecHit_Aux_StatusWord_HB") ==0 ){ RecHit_Aux_StatusWord_HB= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "HcalRecHitTask_RecHit_Aux_StatusWord_HE") ==0 ){ RecHit_Aux_StatusWord_HE= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "HcalRecHitTask_RecHit_Aux_StatusWord_HO") ==0 ){ RecHit_Aux_StatusWord_HO= hcalMEs[ih]; }
if( strcmp(hcalMEs[ih]->getName().c_str(), "HcalRecHitTask_RecHit_Aux_StatusWord_HF") ==0 ){ RecHit_Aux_StatusWord_HF= hcalMEs[ih]; }
}
if( useAllHistos !=0 && useAllHistos !=1 ) return 0;
//None of the ZS stuff necessary for drawn histograms
if (subdet==6 && useAllHistos) {
// FIXME: a dummy emap_min! No solutions yet found.
// Since useAllHistos is set to disable it in our ususal validation, it's left here to be fixed later...
double emap_min[82][72][4][4];
/* NOT a valid solution
for(unsigned int i1=0; i1 <82; i1++){
for(unsigned int i2=0; i2<72; i2++){
for(unsigned int i3=0; i3<4; i3++){
for(unsigned int i4=0; i4<4; i4++){
int idx = i1 + i2*82 + i3*(82*72) + i4*(82*72*4);
emap_min[i1][i2][i3][i4] = emap_min_ME->GetBinContent(idx+1);
}
}
}
}
*/
for (unsigned int i1 = 0; i1 < 82; i1++) {
for (unsigned int i2 = 0; i2 < 72; i2++) {
int index = (i1-41) * 72 + i2;
double e = emap_min[i1][i2][0][0];
if( e < 10000.) {
ZS_HB1->Fill(e);
ZS_seqHB1->Fill(double(index),e);
}
e = emap_min[i1][i2][1][0];
if( e < 10000.) {
ZS_HB2->Fill(e);
ZS_seqHB2->Fill(double(index),e);
}
e = emap_min[i1][i2][0][1];
if( e < 10000.) {
ZS_HE1->Fill(e);
ZS_seqHE1->Fill(double(index),e);
}
e = emap_min[i1][i2][1][1];
if( e < 10000.) {
ZS_HE2->Fill(e);
ZS_seqHE2->Fill(double(index),e);
}
e = emap_min[i1][i2][2][1];
if( e < 10000.) {
ZS_HE3->Fill(e);
ZS_seqHE3->Fill(double(index),e);
}
e = emap_min[i1][i2][3][2];
if( e < 10000.) {
ZS_HO->Fill(e);
ZS_seqHO->Fill(double(index),e);
}
e = emap_min[i1][i2][0][3];
if( e < 10000.) {
ZS_HF1->Fill(e);
ZS_seqHF1->Fill(double(index),e);
}
e = emap_min[i1][i2][1][3];
if( e < 10000.) {
ZS_HF2->Fill(e);
ZS_seqHF2->Fill(double(index),e);
}
for (unsigned int i3 = 0; i3 < 4; i3++) { // depth
double emin = 100000.;
for (unsigned int i4 = 0; i4 < 4; i4++) { // subdet
/*
std::cout << "* ieta, iphi, depth, sub = "
<< i1 << ", " << i2 << ", " << i3 << ", " << i4
<< " emap_min = " << emap_min [i1][i2][i3][i4]
<< std::endl;
*/
if ( emin > emap_min [i1][i2][i3][i4])
emin = emap_min [i1][i2][i3][i4];
}
int ieta = i1-41;
if( i3 == 0 && emin < 10000.) {
map_depth1->Fill(double(ieta),double(i2),emin);
/*
std::cout << "* Fill map_depth1 " << double(ieta) << " "
<< double(i2) << " with " << emin << std::endl;
*/
}
if( i3 == 1 && emin < 10000.)
map_depth2->Fill(double(ieta),double(i2),emin);
if( i3 == 2 && emin < 10000.)
map_depth3->Fill(double(ieta),double(i2),emin);
if( i3 == 3 && emin < 10000.)
map_depth4->Fill(double(ieta),double(i2),emin);
}
}
}
}
// mean energies and occupancies evaluation
else {
double nevtot = Nhf->getEntries();
if(verbose_) std::cout<<"nevtot : "<<nevtot<<std::endl;
int nx = occupancy_map_HB1->getNbinsX();
int ny = occupancy_map_HB1->getNbinsY();
float cnorm;
float fev = float (nevtot);
// std::cout << "*** nevtot " << nevtot << std::endl;
float sumphi_hb1, sumphi_hb2, sumphi_he1, sumphi_he2, sumphi_he3,
sumphi_ho, sumphi_hf1, sumphi_hf2;
/*
if(nx != 82 || ny != 72)
std::cout << "*** problem with binning " << std::endl;
*/
float phi_factor;
// First - special <E> maps
int nx1 = emap_depth1->getNbinsX();
int ny1 = emap_depth1->getNbinsY();
for (int i = 1; i <= nx1; i++) {
for (int j = 1; j <= ny1; j++) {
cnorm = emap_depth1->getBinContent(i,j) / fev;
emap_depth1->setBinContent(i,j,cnorm);
cnorm = emap_depth2->getBinContent(i,j) / fev;
emap_depth2->setBinContent(i,j,cnorm);
cnorm = emap_depth3->getBinContent(i,j) / fev;
emap_depth3->setBinContent(i,j,cnorm);
cnorm = emap_depth4->getBinContent(i,j) / fev;
emap_depth4->setBinContent(i,j,cnorm);
}
}
// Second: all others regular maps
for (int i = 1; i <= nx; i++) {
sumphi_hb1 = 0.;
sumphi_hb2 = 0.;
sumphi_he1 = 0.;
sumphi_he2 = 0.;
sumphi_he3 = 0.;
sumphi_ho = 0.;
sumphi_hf1 = 0.;
sumphi_hf2 = 0.;
for (int j = 1; j <= ny; j++) {
int index = (i-42) * ny + j-1;
//Occupancies (needed for occ vs ieta histos)
cnorm = occupancy_map_HB1->getBinContent(i,j) / fev;
occupancy_map_HB1->setBinContent(i,j,cnorm);
cnorm = occupancy_map_HB2->getBinContent(i,j) / fev;
occupancy_map_HB2->setBinContent(i,j,cnorm);
cnorm = occupancy_map_HE1->getBinContent(i,j) / fev;
occupancy_map_HE1->setBinContent(i,j,cnorm);
cnorm = occupancy_map_HE2->getBinContent(i,j) / fev;
occupancy_map_HE2->setBinContent(i,j,cnorm);
cnorm = occupancy_map_HE3->getBinContent(i,j) / fev;
occupancy_map_HE3->setBinContent(i,j,cnorm);
cnorm = occupancy_map_HO->getBinContent(i,j) / fev;
occupancy_map_HO->setBinContent(i,j,cnorm);
cnorm = occupancy_map_HF1->getBinContent(i,j) / fev;
occupancy_map_HF1->setBinContent(i,j,cnorm);
cnorm = occupancy_map_HF2->getBinContent(i,j) / fev;
occupancy_map_HF2->setBinContent(i,j,cnorm);
sumphi_hb1 += occupancy_map_HB1->getBinContent(i,j);
sumphi_hb2 += occupancy_map_HB2->getBinContent(i,j);
sumphi_he1 += occupancy_map_HE1->getBinContent(i,j);
sumphi_he2 += occupancy_map_HE2->getBinContent(i,j);
sumphi_he3 += occupancy_map_HE3->getBinContent(i,j);
sumphi_ho += occupancy_map_HO->getBinContent(i,j);
sumphi_hf1 += occupancy_map_HF1->getBinContent(i,j);
sumphi_hf2 += occupancy_map_HF2->getBinContent(i,j);
// Occupancies - not in main drawn set of histos
if(useAllHistos){
occupancy_seqHB1->Fill(double(index),cnorm);
occupancy_seqHB2->Fill(double(index),cnorm);
occupancy_seqHE1->Fill(double(index),cnorm);
occupancy_seqHE2->Fill(double(index),cnorm);
occupancy_seqHE3->Fill(double(index),cnorm);
occupancy_seqHO->Fill(double(index),cnorm);
occupancy_seqHF1->Fill(double(index),cnorm);
occupancy_seqHF2->Fill(double(index),cnorm);
}
}
int ieta = i - 42; // -41 -1, 0 40
if(ieta >=0 ) ieta +=1; // -41 -1, 1 41 - to make it detector-like
if(ieta >= -20 && ieta <= 20 )
{phi_factor = 72.;}
else {
if(ieta >= 40 || ieta <= -40 ) {phi_factor = 18.;}
else
phi_factor = 36.;
}
if(ieta >= 0) ieta -= 1; // -41 -1, 0 40 - to bring back to histo num
/*
std::cout << "*** ieta = " << ieta << " sumphi_hb1, sumphi_hb2, sumphi_he1, sumphi_he2, simphi_he3, sumphi_ho, simphi_hf1, sumphi_hf2" << std::endl
<< sumphi_hb1 << " " << sumphi_hb2 << " " << sumphi_he1 << " "
<< sumphi_he2 << " " << simphi_he3 << " " << sumphi_ho << " "
<< simphi_hf1 << " " << sumphi_hf2 << std::endl << std::endl;
*/
//Occupancy vs. ieta histos are drawn, RMS is not
cnorm = sumphi_hb1 / phi_factor;
occupancy_vs_ieta_HB1->Fill(float(ieta), cnorm);
cnorm = sumphi_hb2 / phi_factor;
occupancy_vs_ieta_HB2->Fill(float(ieta), cnorm);
cnorm = sumphi_he1 / phi_factor;
occupancy_vs_ieta_HE1->Fill(float(ieta), cnorm);
cnorm = sumphi_he2 / phi_factor;
occupancy_vs_ieta_HE2->Fill(float(ieta), cnorm);
cnorm = sumphi_he3 / phi_factor;
occupancy_vs_ieta_HE3->Fill(float(ieta), cnorm);
cnorm = sumphi_ho / phi_factor;
occupancy_vs_ieta_HO->Fill(float(ieta), cnorm);
cnorm = sumphi_hf1 / phi_factor;
occupancy_vs_ieta_HF1->Fill(float(ieta), cnorm);
cnorm = sumphi_hf2 / phi_factor;
occupancy_vs_ieta_HF2->Fill(float(ieta), cnorm);
if (useAllHistos){
// RMS vs ieta (Emean's one)
cnorm = emean_vs_ieta_HB1->getBinError(i);
RMS_vs_ieta_HB1->Fill(ieta,cnorm);
cnorm = emean_vs_ieta_HB2->getBinError(i);
RMS_vs_ieta_HB2->Fill(ieta,cnorm);
cnorm = emean_vs_ieta_HE1->getBinError(i);
RMS_vs_ieta_HE1->Fill(ieta,cnorm);
cnorm = emean_vs_ieta_HE2->getBinError(i);
RMS_vs_ieta_HE2->Fill(ieta,cnorm);
cnorm = emean_vs_ieta_HE1->getBinError(i);
RMS_vs_ieta_HE3->Fill(ieta,cnorm);
cnorm = emean_vs_ieta_HB1->getBinError(i);
RMS_vs_ieta_HO->Fill(ieta,cnorm);
cnorm = emean_vs_ieta_HB1->getBinError(i);
RMS_vs_ieta_HF1->Fill(ieta,cnorm);
cnorm = emean_vs_ieta_HB1->getBinError(i);
RMS_vs_ieta_HF2->Fill(ieta,cnorm);
}
} // end of i-loop
// RMS seq (not drawn)
if(useAllHistos){
nx = emean_seqHB1->getNbinsX();
for(int ibin = 1; ibin <= nx; ibin++ ){
cnorm = emean_seqHB1->getBinError(ibin);
RMS_seq_HB1->setBinContent(ibin, cnorm);
cnorm = emean_seqHB2->getBinError(ibin);
RMS_seq_HB2->setBinContent(ibin, cnorm);
cnorm = emean_seqHO->getBinError(ibin);
RMS_seq_HO->setBinContent(ibin, cnorm);
}
nx = emean_seqHE1->getNbinsX();
for(int ibin = 1; ibin <= nx; ibin++ ){
cnorm = emean_seqHE1->getBinError(ibin);
RMS_seq_HE1->setBinContent(ibin, cnorm);
cnorm = emean_seqHE2->getBinError(ibin);
RMS_seq_HE2->setBinContent(ibin, cnorm);
cnorm = emean_seqHE3->getBinError(ibin);
RMS_seq_HE3->setBinContent(ibin, cnorm);
}
nx = emean_seqHF1->getNbinsX();
for(int ibin = 1; ibin <= nx; ibin++ ){
cnorm = emean_seqHF1->getBinError(ibin);
RMS_seq_HF1->setBinContent(ibin, cnorm);
cnorm = emean_seqHF2->getBinError(ibin);
RMS_seq_HF2->setBinContent(ibin, cnorm);
}
}
//Status Word (drawn)
nx = RecHit_StatusWord_HB->getNbinsX();
for (int ibin = 1; ibin <= nx; ibin++) {
cnorm = RecHit_StatusWord_HB->getBinContent(ibin) / (fev * 2592.);
RecHit_StatusWord_HB->setBinContent(ibin,cnorm);
cnorm = RecHit_StatusWord_HE->getBinContent(ibin) / (fev * 2592.);
RecHit_StatusWord_HE->setBinContent(ibin,cnorm);
cnorm = RecHit_StatusWord_HO->getBinContent(ibin) / (fev * 2160.);
RecHit_StatusWord_HO->setBinContent(ibin,cnorm);
cnorm = RecHit_StatusWord_HF->getBinContent(ibin) / (fev * 1728.);
RecHit_StatusWord_HF->setBinContent(ibin,cnorm);
cnorm = RecHit_Aux_StatusWord_HB->getBinContent(ibin) / (fev * 2592.);
RecHit_Aux_StatusWord_HB->setBinContent(ibin,cnorm);
cnorm = RecHit_Aux_StatusWord_HE->getBinContent(ibin) / (fev * 2592.);
RecHit_Aux_StatusWord_HE->setBinContent(ibin,cnorm);
cnorm = RecHit_Aux_StatusWord_HO->getBinContent(ibin) / (fev * 2160.);
RecHit_Aux_StatusWord_HO->setBinContent(ibin,cnorm);
cnorm = RecHit_Aux_StatusWord_HF->getBinContent(ibin) / (fev * 1728.);
RecHit_Aux_StatusWord_HF->setBinContent(ibin,cnorm);
}
//HF 2-bit status word (not drawn)
if(useAllHistos){
nx = RecHit_StatusWord_HF67->getNbinsX();
for (int ibin = 1; ibin <= nx; ibin++) {
cnorm = RecHit_StatusWord_HF67->getBinContent(ibin) / (fev * 1728.);
RecHit_StatusWord_HF67->setBinContent(ibin,cnorm);
}
}
}
return 1;
}
| void HcalRecHitsClient::runClient_ | ( | ) | [virtual] |
Definition at line 72 of file HcalRecHitsClient.cc.
References gather_cfg::cout, dbe_, dirName_, DQMStore::getContents(), DQMStore::getSubdirs(), HcalRecHitsEndjob(), i, j, DQMStore::setCurrentFolder(), and verbose_.
Referenced by endRun().
{
if(!dbe_) return; //we dont have the DQMStore so we cant do anything
dbe_->setCurrentFolder(dirName_);
if (verbose_) std::cout << "\nrunClient" << std::endl;
std::vector<MonitorElement*> hcalMEs;
// Since out folders are fixed to three, we can just go over these three folders
// i.e., CaloTowersV/CaloTowersTask, HcalRecHitsV/HcalRecHitTask, NoiseRatesV/NoiseRatesTask.
std::vector<std::string> fullPathHLTFolders = dbe_->getSubdirs();
for(unsigned int i=0;i<fullPathHLTFolders.size();i++) {
if (verbose_) std::cout <<"\nfullPath: "<< fullPathHLTFolders[i] << std::endl;
dbe_->setCurrentFolder(fullPathHLTFolders[i]);
std::vector<std::string> fullSubPathHLTFolders = dbe_->getSubdirs();
for(unsigned int j=0;j<fullSubPathHLTFolders.size();j++) {
if (verbose_) std::cout <<"fullSub: "<<fullSubPathHLTFolders[j] << std::endl;
if( strcmp(fullSubPathHLTFolders[j].c_str(), "HcalRecHitsV/HcalRecHitTask") ==0 ){
hcalMEs = dbe_->getContents(fullSubPathHLTFolders[j]);
if (verbose_) std::cout <<"hltMES size : "<<hcalMEs.size()<<std::endl;
if( !HcalRecHitsEndjob(hcalMEs) ) std::cout<<"\nError in HcalRecHitsEndjob!"<<std::endl<<std::endl;
}
}
}
}
edm::ParameterSet HcalRecHitsClient::conf_ [private] |
Definition at line 42 of file HcalRecHitsClient.h.
DQMStore* HcalRecHitsClient::dbe_ [private] |
Definition at line 39 of file HcalRecHitsClient.h.
Referenced by endJob(), HcalRecHitsClient(), and runClient_().
bool HcalRecHitsClient::debug_ [private] |
Definition at line 45 of file HcalRecHitsClient.h.
Referenced by HcalRecHitsClient().
std::string HcalRecHitsClient::dirName_ [private] |
Definition at line 47 of file HcalRecHitsClient.h.
Referenced by HcalRecHitsClient(), and runClient_().
std::string HcalRecHitsClient::dirNameJet_ [private] |
Definition at line 48 of file HcalRecHitsClient.h.
std::string HcalRecHitsClient::dirNameMET_ [private] |
Definition at line 49 of file HcalRecHitsClient.h.
std::string HcalRecHitsClient::outputFile_ [private] |
Definition at line 40 of file HcalRecHitsClient.h.
Referenced by endJob(), and HcalRecHitsClient().
bool HcalRecHitsClient::verbose_ [private] |
Definition at line 44 of file HcalRecHitsClient.h.
Referenced by HcalRecHitsClient(), HcalRecHitsEndjob(), and runClient_().
1.7.3