|
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
|
#include "CalibMuon/CSCCalibration/interface/CSCFakeGainsConditions.h"
CSCGains *CSCFakeGainsConditions::prefillGains() {
CSCGains *cngains = new CSCGains();
const CSCDetId &detId = CSCDetId();
int max_istrip, id_layer, max_ring, max_cham;
seed = 10000;
srand(seed);
mean = 6.8, min = -10.0, minchi = 1.0, 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;
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<CSCGains::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++) {
// itemvector[istrip].gain_slope = 7.55;
// itemvector[istrip].gain_intercept= -10.00;
// itemvector[istrip].gain_chi2 = 2.00;
itemvector[istrip].gain_slope = ((double)rand() / ((double)(RAND_MAX) + (double)(1))) + mean;
itemvector[istrip].gain_intercept = ((double)rand() / ((double)(RAND_MAX) + (double)(1))) + min;
itemvector[istrip].gain_chi2 = ((double)rand() / ((double)(RAND_MAX) + (double)(1))) + minchi;
cngains->gains[id_layer] = itemvector;
}
}
}
}
}
}
return cngains;
}
CSCFakeGainsConditions::CSCFakeGainsConditions(const edm::ParameterSet &iConfig) {
// the following line is needed to tell the framework what
// data is being produced
setWhatProduced(this, &CSCFakeGainsConditions::produceGains);
findingRecord<CSCGainsRcd>();
// now do what ever other initialization is needed
}
CSCFakeGainsConditions::~CSCFakeGainsConditions() {
// 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 ------------
CSCFakeGainsConditions::ReturnType CSCFakeGainsConditions::produceGains(const CSCGainsRcd &iRecord) {
return CSCFakeGainsConditions::ReturnType(prefillGains());
}
void CSCFakeGainsConditions::setIntervalFor(const edm::eventsetup::EventSetupRecordKey &,
const edm::IOVSyncValue &,
edm::ValidityInterval &oValidity) {
oValidity = edm::ValidityInterval(edm::IOVSyncValue::beginOfTime(), edm::IOVSyncValue::endOfTime());
}
|