|
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
|
#ifndef MLClient_H
#define MLClient_H
#include "DQWorkerClient.h"
#include <vector>
#include <valarray>
#include <deque>
namespace ecaldqm {
class MLClient : public DQWorkerClient {
public:
MLClient();
~MLClient() override {}
void producePlots(ProcessType) override;
private:
void setParams(edm::ParameterSet const&) override;
//Occcupancy map is plotted at the tower level
//For the EB: 34 towers in the eta and 72 towers in the phi directions
static const int nEBEtaTowers = 34;
static const int nEBPhiTowers = 72;
//For the EE: 20 towers in the eta and 20 towers in the phi directions
static const int nEEEtaTowers = 20;
static const int nEEPhiTowers = 20;
//After padding with two rows above and below to prevent the edge effect
//For EB: 36 towers in eta direction
//For EE: padding on all four sides, 22 towers in both eta and phi directions.
static const int nEBEtaTowersPad = 36;
static const int nEETowersPad = 22;
float EBThreshold_;
float EEpThreshold_;
float EEmThreshold_;
float EB_PUcorr_slope_;
float EEp_PUcorr_slope_;
float EEm_PUcorr_slope_;
float EB_PUcorr_intercept_;
float EEp_PUcorr_intercept_;
float EEm_PUcorr_intercept_;
size_t nLS = 3; //No.of lumisections to add the occupancy over
size_t nLSloss = 6; //No.of lumisections to multiply the loss over
int nbadtowerEB; //count the no.of bad towers flagged by the ML model.
int nbadtowerEE;
int LScount = 0; //count no.of lumisections over which the MLquality is made.
std::deque<int> NEventQ; //To keep the no.of events in each occupancy plot
//To keep the input occupancy plots to be summed
std::deque<std::valarray<float>> ebOccMap1dQ;
std::deque<std::valarray<float>> eepOccMap1dQ;
std::deque<std::valarray<float>> eemOccMap1dQ;
//To keep the average occupancy to do response correction
std::vector<double> EBavgOcc;
std::vector<double> EEpavgOcc;
std::vector<double> EEmavgOcc;
//To keep the ML losses to be multiplied
std::deque<std::valarray<std::valarray<float>>> EBlossMap2dQ;
std::deque<std::valarray<std::valarray<float>>> EEplossMap2dQ;
std::deque<std::valarray<std::valarray<float>>> EEmlossMap2dQ;
};
} // namespace ecaldqm
#endif
|