|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
|
#include "DQM/SiStripCommissioningClients/interface/NoiseHistograms.h"
#include "CondFormats/SiStripObjects/interface/NoiseAnalysis.h"
#include "DQM/SiStripCommissioningAnalysis/interface/NoiseAlgorithm.h"
#include "DQM/SiStripCommissioningSummary/interface/NoiseSummaryFactory.h"
#include "DQM/SiStripCommon/interface/ExtractTObject.h"
#include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iostream>
#include <memory>
#include <sstream>
#include <iomanip>
#include "TProfile.h"
using namespace std;
using namespace sistrip;
// -----------------------------------------------------------------------------
/** */
NoiseHistograms::NoiseHistograms(const edm::ParameterSet& pset, DQMStore* bei)
: CommissioningHistograms(pset.getParameter<edm::ParameterSet>("NoiseParameters"), bei, sistrip::NOISE) {
factory_ = std::make_unique<NoiseSummaryFactory>();
LogTrace(mlDqmClient_) << "[NoiseHistograms::" << __func__ << "]"
<< " Constructing object...";
}
// -----------------------------------------------------------------------------
/** */
NoiseHistograms::~NoiseHistograms() {
LogTrace(mlDqmClient_) << "[NoiseHistograms::" << __func__ << "]"
<< " Destructing object...";
}
// -----------------------------------------------------------------------------
/** */
void NoiseHistograms::histoAnalysis(bool debug) {
LogTrace(mlDqmClient_) << "[NoiseHistograms::" << __func__ << "]";
// Some initialisation
uint16_t valid = 0;
HistosMap::const_iterator iter;
Analyses::iterator ianal;
std::map<std::string, uint16_t> errors;
// Clear map holding analysis objects
for (ianal = data().begin(); ianal != data().end(); ianal++) {
if (ianal->second) {
delete ianal->second;
}
}
data().clear();
// Iterate through map containing histograms
for (iter = histos().begin(); iter != histos().end(); iter++) {
// Check vector of histos is not empty
if (iter->second.empty()) {
edm::LogWarning(mlDqmClient_) << "[NoiseHistograms::" << __func__ << "]"
<< " Zero histograms found!";
continue;
}
// Retrieve pointers to profile histos
std::vector<TH1*> profs;
Histos::const_iterator ihis = iter->second.begin();
for (; ihis != iter->second.end(); ihis++) {
TProfile* prof = ExtractTObject<TProfile>().extract((*ihis)->me_);
if (prof) {
profs.push_back(prof);
}
//@@ Common mode histos?...
//TH1F* his = ExtractTObject<TH1F>().extract( (*ihis)->me_ );
//if ( his ) { profs.push_back(his); }
}
// Perform histo analysis
NoiseAnalysis* anal = new NoiseAnalysis(iter->first);
NoiseAlgorithm algo(this->pset(), anal);
algo.analysis(profs);
data()[iter->first] = anal;
if (anal->isValid()) {
valid++;
}
if (!anal->getErrorCodes().empty()) {
errors[anal->getErrorCodes()[0]]++;
}
}
if (!histos().empty()) {
edm::LogVerbatim(mlDqmClient_) << "[NoiseHistograms::" << __func__ << "]"
<< " Analyzed histograms for " << histos().size() << " FED channels, of which "
<< valid << " (" << 100 * valid / histos().size() << "%) are valid.";
if (!errors.empty()) {
uint16_t count = 0;
std::stringstream ss;
ss << std::endl;
std::map<std::string, uint16_t>::const_iterator ii;
for (ii = errors.begin(); ii != errors.end(); ++ii) {
ss << " " << ii->first << ": " << ii->second << std::endl;
count += ii->second;
}
edm::LogWarning(mlDqmClient_) << "[NoiseHistograms::" << __func__ << "]"
<< " Found " << count << " errors (" << 100 * count / histos().size()
<< "%): " << ss.str();
}
} else {
edm::LogWarning(mlDqmClient_) << "[NoiseHistograms::" << __func__ << "]"
<< " No histograms to analyze!";
}
}
// -----------------------------------------------------------------------------
/** */
void NoiseHistograms::printAnalyses() {
Analyses::iterator ianal = data().begin();
Analyses::iterator janal = data().end();
for (; ianal != janal; ++ianal) {
if (ianal->second) {
std::stringstream ss;
ianal->second->print(ss, 1);
ianal->second->print(ss, 2);
if (ianal->second->isValid()) {
LogTrace(mlDqmClient_) << ss.str();
} else {
edm::LogWarning(mlDqmClient_) << ss.str();
}
}
}
}
|