|
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
|
#include <fstream>
#include "CalibMuon/CSCCalibration/interface/CSCFakeCrosstalkConditions.h"
CSCcrosstalk *CSCFakeCrosstalkConditions::prefillCrosstalk() {
const CSCDetId &detId = CSCDetId();
CSCcrosstalk *cncrosstalk = new CSCcrosstalk();
int max_istrip, id_layer, max_ring, max_cham;
seed = 10000;
srand(seed);
mean = -0.0009, min = 0.035, minchi = 1.5, M = 1000;
// endcap=1 to 2,station=1 to 4, ring=1 to 4,chamber=1 to 36,layer=1 to 6
for (int iendcap = detId.minEndcapId(); iendcap <= detId.maxEndcapId(); iendcap++) {
for (int istation = detId.minStationId(); istation <= detId.maxStationId(); istation++) {
max_ring = detId.maxRingId();
// station 4 ring 4 not there(36 chambers*2 missing)
// 3 rings max this way of counting (ME1a & b)
if (istation == 1)
max_ring = 3;
if (istation == 2)
max_ring = 2;
if (istation == 3)
max_ring = 2;
if (istation == 4)
max_ring = 1;
for (int iring = detId.minRingId(); iring <= max_ring; iring++) {
max_istrip = 80;
max_cham = detId.maxChamberId();
if (istation == 1 && iring == 1)
max_cham = 36;
if (istation == 1 && iring == 2)
max_cham = 36;
if (istation == 1 && iring == 3)
max_cham = 36;
if (istation == 2 && iring == 1)
max_cham = 18;
if (istation == 2 && iring == 2)
max_cham = 36;
if (istation == 3 && iring == 1)
max_cham = 18;
if (istation == 3 && iring == 2)
max_cham = 36;
if (istation == 4 && iring == 1)
max_cham = 18;
// station 1 ring 3 has 64 strips per layer instead of 80(minus & plus
// side!!!)
for (int ichamber = detId.minChamberId(); ichamber <= max_cham; ichamber++) {
for (int ilayer = detId.minLayerId(); ilayer <= detId.maxLayerId(); ilayer++) {
// station 1 ring 3 has 64 strips per layer instead of 80
if (istation == 1 && iring == 3)
max_istrip = 64;
std::vector<CSCcrosstalk::Item> itemvector;
itemvector.resize(max_istrip);
id_layer = 100000 * iendcap + 10000 * istation + 1000 * iring + 10 * ichamber + ilayer;
for (int istrip = 0; istrip < max_istrip; istrip++) {
// create fake values
itemvector[istrip].xtalk_slope_right =
-((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 10000 + mean;
itemvector[istrip].xtalk_intercept_right =
((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 100 + min;
itemvector[istrip].xtalk_chi2_right = ((double)rand() / ((double)(RAND_MAX) + (double)(1))) + minchi;
itemvector[istrip].xtalk_slope_left =
-((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 10000 + mean;
itemvector[istrip].xtalk_intercept_left =
((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 100 + min;
itemvector[istrip].xtalk_chi2_left = ((double)rand() / ((double)(RAND_MAX) + (double)(1))) + minchi;
cncrosstalk->crosstalk[id_layer] = itemvector;
if (istrip == 0) {
itemvector[istrip].xtalk_slope_right =
-((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 10000 + mean;
itemvector[istrip].xtalk_intercept_right =
((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 100 + min;
itemvector[istrip].xtalk_chi2_right = ((double)rand() / ((double)(RAND_MAX) + (double)(1))) + minchi;
itemvector[istrip].xtalk_slope_left = 0.0;
itemvector[istrip].xtalk_intercept_left = 0.0;
itemvector[istrip].xtalk_chi2_left = 0.0;
cncrosstalk->crosstalk[id_layer] = itemvector;
}
if (istrip == 79) {
itemvector[istrip].xtalk_slope_right = 0.0;
itemvector[istrip].xtalk_intercept_right = 0.0;
itemvector[istrip].xtalk_chi2_right = 0.0;
itemvector[istrip].xtalk_slope_left =
-((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 10000 + mean;
itemvector[istrip].xtalk_intercept_left =
((double)rand() / ((double)(RAND_MAX) + (double)(1))) / 100 + min;
itemvector[istrip].xtalk_chi2_left = ((double)rand() / ((double)(RAND_MAX) + (double)(1))) + minchi;
cncrosstalk->crosstalk[id_layer] = itemvector;
}
}
}
}
}
}
}
return cncrosstalk;
}
CSCFakeCrosstalkConditions::CSCFakeCrosstalkConditions(const edm::ParameterSet &iConfig) {
// the following line is needed to tell the framework what
// data is being produced
setWhatProduced(this, &CSCFakeCrosstalkConditions::produceCrosstalk);
findingRecord<CSCcrosstalkRcd>();
// now do what ever other initialization is needed
}
CSCFakeCrosstalkConditions::~CSCFakeCrosstalkConditions() {
// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)
}
//
// member functions
//
// ------------ method called to produce the data ------------
CSCFakeCrosstalkConditions::ReturnType CSCFakeCrosstalkConditions::produceCrosstalk(const CSCcrosstalkRcd &iRecord) {
return CSCFakeCrosstalkConditions::ReturnType(prefillCrosstalk());
}
void CSCFakeCrosstalkConditions::setIntervalFor(const edm::eventsetup::EventSetupRecordKey &,
const edm::IOVSyncValue &,
edm::ValidityInterval &oValidity) {
oValidity = edm::ValidityInterval(edm::IOVSyncValue::beginOfTime(), edm::IOVSyncValue::endOfTime());
}
|