|
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
|
#include "DQM/SiStripCommissioningClients/interface/SamplingHistograms.h"
#include "CondFormats/SiStripObjects/interface/SamplingAnalysis.h"
#include "DQM/SiStripCommissioningAnalysis/interface/SamplingAlgorithm.h"
#include "DQM/SiStripCommissioningSummary/interface/SummaryGenerator.h"
#include "DQM/SiStripCommon/interface/ExtractTObject.h"
#include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
#include "DataFormats/SiStripCommon/interface/SiStripEnumsAndStrings.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iostream>
#include <memory>
#include <sstream>
#include <iomanip>
#include "TProfile.h"
using namespace std;
using namespace sistrip;
// -----------------------------------------------------------------------------
/** */
SamplingHistograms::SamplingHistograms(const edm::ParameterSet& pset, DQMStore* bei, const sistrip::RunType& task)
: CommissioningHistograms(pset.getParameter<edm::ParameterSet>("SamplingParameters"), bei, task), sOnCut_(3) {
LogTrace(mlDqmClient_) << "[SamplingHistograms::" << __func__ << "]"
<< " Constructing object...";
factory_ = std::make_unique<SamplingSummaryFactory>();
// retreive the latency code from the root file
std::string dataPath = std::string(sistrip::collate_) + "/" + sistrip::root_ + "/latencyCode";
MonitorElement* codeElement = bei->get(dataPath);
if (codeElement)
latencyCode_ = codeElement->getIntValue();
else
latencyCode_ = 0;
}
// -----------------------------------------------------------------------------
/** */
SamplingHistograms::~SamplingHistograms() {
LogTrace(mlDqmClient_) << "[SamplingHistograms::" << __func__ << "]"
<< " Deleting object...";
}
// -----------------------------------------------------------------------------
/** */
void SamplingHistograms::histoAnalysis(bool debug) {
// Clear map holding analysis objects
Analyses::iterator ianal;
for (ianal = data().begin(); ianal != data().end(); ianal++) {
if (ianal->second) {
delete ianal->second;
}
}
data().clear();
// Iterate through map containing vectors of profile histograms
HistosMap::const_iterator iter = histos().begin();
for (; iter != histos().end(); iter++) {
// Check vector of histos is not empty (should be 1 histo)
if (iter->second.empty()) {
edm::LogWarning(mlDqmClient_) << "[SamplingHistograms::" << __func__ << "]"
<< " Zero collation histograms found!";
continue;
}
// Retrieve pointers to profile histos for this FED channel
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);
}
}
// Perform histo analysis
SamplingAnalysis* anal = new SamplingAnalysis(iter->first);
anal->setSoNcut(sOnCut_);
SamplingAlgorithm algo(this->pset(), anal, latencyCode_);
algo.analysis(profs);
data()[iter->first] = anal;
}
}
void SamplingHistograms::configure(const edm::ParameterSet& pset, const edm::EventSetup&) {
//TODO: should use the parameter set. Why is this crashing ???
// sOnCut_ = pset.getParameter<double>("SignalToNoiseCut");
sOnCut_ = 3.;
}
|