|
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
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
|
#ifndef _CSCFAKEDBCROSSTALK_H
#define _CSCFAKEDBCROSSTALK_H
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/SourceFactory.h"
#include <memory>
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "CondFormats/CSCObjects/interface/CSCDBCrosstalk.h"
#include "CondFormats/DataRecord/interface/CSCDBCrosstalkRcd.h"
#include <DataFormats/MuonDetId/interface/CSCDetId.h>
class CSCFakeDBCrosstalk : public edm::ESProducer, public edm::EventSetupRecordIntervalFinder {
public:
CSCFakeDBCrosstalk(const edm::ParameterSet &);
~CSCFakeDBCrosstalk() override;
inline static CSCDBCrosstalk *prefillDBCrosstalk();
typedef std::unique_ptr<CSCDBCrosstalk> Pointer;
Pointer produceDBCrosstalk(const CSCDBCrosstalkRcd &);
private:
// ----------member data ---------------------------
void setIntervalFor(const edm::eventsetup::EventSetupRecordKey &,
const edm::IOVSyncValue &,
edm::ValidityInterval &) override;
Pointer cndbCrosstalk;
};
#include <fstream>
#include <iostream>
#include <vector>
// to workaround plugin library
inline CSCDBCrosstalk *CSCFakeDBCrosstalk::prefillDBCrosstalk() {
int seed;
float mean, min;
int ii, jj, iii, jjj;
const int MAX_SIZE = 217728; // or 252288 for ME4/2 chambers
const int SLOPE_FACTOR = 10000000;
const int INTERCEPT_FACTOR = 100000;
CSCDBCrosstalk *cndbcrosstalk = new CSCDBCrosstalk();
cndbcrosstalk->crosstalk.resize(MAX_SIZE);
seed = 10000;
srand(seed);
mean = -0.0009, min = 0.035;
ii = 0, jj = 0, iii = 0, jjj = 0;
cndbcrosstalk->factor_slope = int(SLOPE_FACTOR);
cndbcrosstalk->factor_intercept = int(INTERCEPT_FACTOR);
for (int i = 0; i < MAX_SIZE; i++) {
cndbcrosstalk->crosstalk[i].xtalk_slope_right =
(short int)((-((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 10000 + mean) * SLOPE_FACTOR + 0.5);
cndbcrosstalk->crosstalk[i].xtalk_intercept_right =
(short int)((((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 100 + min) * INTERCEPT_FACTOR + 0.5);
cndbcrosstalk->crosstalk[i].xtalk_slope_left =
(short int)((-((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 10000 + mean) * SLOPE_FACTOR + 0.5);
cndbcrosstalk->crosstalk[i].xtalk_intercept_left =
(short int)((((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 100 + min) * INTERCEPT_FACTOR + 0.5);
// 80 strips per chamber
if (i < 34561 && i % 80 == 0) {
cndbcrosstalk->crosstalk[i].xtalk_slope_left = 0;
cndbcrosstalk->crosstalk[i].xtalk_intercept_left = 0;
}
if (i != 0 && i < 34561 && (i + 1) % 80 == 0) {
cndbcrosstalk->crosstalk[i].xtalk_slope_right = 0;
cndbcrosstalk->crosstalk[i].xtalk_intercept_right = 0;
}
// 64 strips per chamber
if (i > 34560 && i < 48385 && i % 64 == 0) {
cndbcrosstalk->crosstalk[i].xtalk_slope_left = 0;
cndbcrosstalk->crosstalk[i].xtalk_intercept_left = 0;
}
if (i > 34560 && i < 48385 && (i + 1) % 64 == 0) {
cndbcrosstalk->crosstalk[i].xtalk_slope_right = 0;
cndbcrosstalk->crosstalk[i].xtalk_intercept_right = 0;
}
// 80 strips per chamber again
if (i > 48384 && i < 143425) {
ii++;
if (i > 48384 && i < 143425 && ii % 80 == 0) {
cndbcrosstalk->crosstalk[i].xtalk_slope_left = 0;
cndbcrosstalk->crosstalk[i].xtalk_intercept_left = 0;
}
}
if (i > 48384 && i < 143425) {
jj++;
if (i > 48384 && i < 143425 && (jj + 1) % 80 == 0) {
cndbcrosstalk->crosstalk[i].xtalk_slope_right = 0;
cndbcrosstalk->crosstalk[i].xtalk_intercept_right = 0;
}
}
// 64 strips per chamber again
if (i > 143424 && i < 157249 && i % 64 == 0) {
cndbcrosstalk->crosstalk[i].xtalk_slope_left = 0;
cndbcrosstalk->crosstalk[i].xtalk_intercept_left = 0;
}
if (i > 143424 && i < 157249 && (i + 1) % 64 == 0) {
cndbcrosstalk->crosstalk[i].xtalk_slope_right = 0;
cndbcrosstalk->crosstalk[i].xtalk_intercept_right = 0;
}
// 80 strips per chamber last time
if (i > 157248) {
iii++;
if (i > 157248 && iii % 80 == 0) {
cndbcrosstalk->crosstalk[i].xtalk_slope_left = 0;
cndbcrosstalk->crosstalk[i].xtalk_intercept_left = 0;
}
}
if (i > 157248) {
jjj++;
if (i > 157248 && (jjj + 1) % 80 == 0) {
cndbcrosstalk->crosstalk[i].xtalk_slope_right = 0;
cndbcrosstalk->crosstalk[i].xtalk_intercept_right = 0;
}
}
}
return cndbcrosstalk;
}
#endif
|