|
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
|
#include <iostream>
#include "DQMOffline/JetMET/plugins/SusyPostProcessor.h"
#include "DQMOffline/JetMET/interface/SusyDQM/Quantile.h"
#include "FWCore/PluginManager/interface/ModuleDef.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
using namespace std;
const char* SusyPostProcessor::messageLoggerCatregory = "SusyDQMPostProcessor";
SusyPostProcessor::SusyPostProcessor(const edm::ParameterSet& pSet) {
iConfig = pSet;
SUSYFolder = iConfig.getParameter<string>("folderName");
_quantile = iConfig.getParameter<double>("quantile");
}
SusyPostProcessor::~SusyPostProcessor() {}
void SusyPostProcessor::QuantilePlots(MonitorElement*& ME, double q_value, DQMStore::IBooker& ibooker_) {
if (ME->getTH1()->GetEntries() > 0.) {
Quantile q(static_cast<const TH1*>(ME->getTH1()));
Float_t mean = q[q_value].first;
Float_t RMS = q[q_value].second;
Float_t xLow = -5.5;
Float_t xUp = 9.5;
Int_t NBin = 15;
if (mean > 0.) {
Float_t DBin = RMS * TMath::Sqrt(12.) / 2.;
xLow = mean - int(mean / DBin + 2) * DBin;
xUp = int(0.2 * mean / DBin) * DBin + mean + 5 * DBin;
NBin = (xUp - xLow) / DBin;
}
ibooker_.setCurrentFolder(ME->getPathname());
TString name = ME->getTH1()->GetName();
name += "_quant";
ME = ibooker_.book1D(name, "", NBin, xLow, xUp);
ME->Fill(mean - RMS);
ME->Fill(mean + RMS);
}
}
void SusyPostProcessor::dqmEndJob(DQMStore::IBooker& ibook_, DQMStore::IGetter& iget_) {
// MET
//----------------------------------------------------------------------------
iget_.setCurrentFolder("JetMET/MET");
Dirs = iget_.getSubdirs();
std::vector<std::string> metFolders;
metFolders.push_back("Uncleaned/");
metFolders.push_back("Cleaned/");
//Need our own copy for thread safety
TF1 mygaus("mygaus", "gaus");
for (int i = 0; i < int(Dirs.size()); i++) {
std::string prefix = "dummy";
if (size_t(Dirs[i].find("met")) != string::npos)
prefix = "met";
if (size_t(Dirs[i].find("pfMet")) != string::npos)
prefix = "pfMET";
for (std::vector<std::string>::const_iterator ic = metFolders.begin(); ic != metFolders.end(); ic++) {
std::string dirName = Dirs[i] + "/" + *ic;
MEx = iget_.get(dirName + "/" + "MEx");
MEy = iget_.get(dirName + "/" + "MEy");
if (MEx && MEx->kind() == MonitorElement::Kind::TH1F) {
if (MEx->getTH1F()->GetEntries() > 50)
MEx->getTH1F()->Fit(&mygaus, "q");
}
if (MEy && MEy->kind() == MonitorElement::Kind::TH1F) {
if (MEy->getTH1F()->GetEntries() > 50)
MEy->getTH1F()->Fit(&mygaus, "q");
}
}
}
// SUSY
//----------------------------------------------------------------------------
iget_.setCurrentFolder(SUSYFolder);
Dirs = iget_.getSubdirs();
for (int i = 0; i < int(Dirs.size()); i++) {
size_t found = Dirs[i].find("Alpha");
if (found != string::npos)
continue;
if (!iget_.dirExists(Dirs[i])) {
edm::LogError(messageLoggerCatregory) << "Directory " << Dirs[i] << " doesn't exist!!";
continue;
}
vector<MonitorElement*> histoVector = iget_.getContents(Dirs[i]);
for (int i = 0; i < int(histoVector.size()); i++) {
QuantilePlots(histoVector[i], _quantile, ibook_);
}
}
}
|