|
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
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
|
#include "CondFormats/EcalObjects/interface/EcalSRSettings.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <sstream>
#include <iostream>
#include <cstdlib>
#include <cassert>
#include <algorithm>
using namespace std;
EcalSRSettings::EcalSRSettings()
: ebDccAdcToGeV_(0.), eeDccAdcToGeV_(0.), bxGlobalOffset_(0), automaticMasks_(0), automaticSrpSelect_(0) {}
#define SR_PRINT(a) o << #a ": " << val.a##_ << "\n";
#define SR_VPRINT(a) \
o << #a; \
if (val.a##_.empty()) \
o << "[0.." << (val.a##_.size() - 1) << "]:"; \
else \
o << "[]: <empty>"; \
for (size_t i = 0; i < val.a##_.size(); ++i) \
o << "\t" << val.a##_[i]; \
o << "\n";
#define SR_VVPRINT(a) \
if (val.a##_.empty()) \
o << #a "[][]: <empty>\n"; \
for (size_t i = 0; i < val.a##_.size(); ++i) { \
o << #a "[" << i << "]"; \
if (val.a##_.empty()) \
o << "[0.." << (val.a##_[i].size() - 1) << "]:"; \
else \
o << "[]: <empty>"; \
for (size_t j = 0; j < val.a##_[i].size(); ++j) \
o << "\t" << val.a##_[i][j]; \
o << "\n"; \
}
std::ostream& operator<<(std::ostream& o, const EcalSRSettings& val) {
o << "# Neighbour eta range, neighborhood: (2*deltaEta+1)*(2*deltaPhi+1)\n"
"# In the vector contains:\n"
"# - 1 element, then value applies to whole ECAL\n"
"# - 2 elements, then element 0 applies to EB, element 1 to EE\n"
"# - 12 elements, then element i applied to SRP (i+1)\n"
"# SRP emulation (see SimCalorimetry/EcalSelectiveReadoutProcuders) supports\n"
"# only 1 element mode.\n";
SR_VPRINT(deltaEta);
o << "\n# Neighbouring eta range, neighborhood: (2*deltaEta+1)*(2*deltaPhi+1)\n"
"# If the vector contains...\n"
"# ... 1 element, then value applies to whole ECAL\n"
"# ... 2 elements, then element 0 applies to EB, element 1 to EE\n"
"# ... 12 elements, then element i applied to SRP (i+1)\n"
"# If the vector contains...\n"
"# ... 1 element, then value applies to whole ECAL\n"
"# ... 2 elements, then element 0 applies to EB, element 1 to EE\n"
"# ... 12 elements, then element i applied to SRP (i+1)\n"
"# SRP emulation (see SimCalorimetry/EcalSelectiveReadoutProcuders) supports\n"
"# only the single-element mode.\n";
SR_VPRINT(deltaPhi);
o << "\n# Index of time sample (staring from 1) the first DCC weights is implied\n"
"# If the vector contains:\n"
"# ... 1 element, then value applies to whole ECAL\n"
"# ... 2 elements, then element 0 applies to EB, element 1 to EE\n"
"# ... 54 elements, then element i applied to DCC (i+1) (FED ID 651+i)\n"
"# SRP emulation (see SimCalorimetry/EcalSelectiveReadoutProcuders) supports\n"
"# only the single-element mode.\n";
SR_VPRINT(ecalDccZs1stSample);
o << "\n# ADC to GeV conversion factor used in ZS filter for EB\n";
SR_PRINT(ebDccAdcToGeV);
o << "\n# ADC to GeV conversion factor used in ZS filter for EE\n";
SR_PRINT(eeDccAdcToGeV);
o << "\n# DCC ZS FIR weights: weights are rounded in such way that in Hw\n"
"# representation (weigth*1024 rounded to nearest integer) the sum is null:\n"
"# Each element is a vector of 6 values, the 6 weights\n"
"# If the vector contains...\n"
"# ... 1 element, then the weight set applies to whole ECAL\n"
"# ... 2 elements, then element 0 applies to EB, element 1 to EE\n"
"# ... 54 elements, then element i applied to DCC (i+1) (FED ID 651+i)\n";
SR_VVPRINT(dccNormalizedWeights);
o << "\n# Switch to use a symetric zero suppression (cut on absolute value). For\n"
"# studies only, for time being it is not supported by the hardware.\n"
"# having troubles for vector<bool> with coral (3.8.0pre1), using vector<int> instead,\n"
"# 0 means false, a value different than 0 means true.\n"
"# If the vector contains...\n"
"# ... 1 element, then the weight set applies to whole ECAL\n"
"# ... 2 elements, then element 0 applies to EB, element 1 to EE\n"
"# ... 54 elements, then element i applied to DCC (i+1) (FED ID 651+i)\n"
"# ... 75848 elements, then:\n"
"# for i < 61200, element i applies to EB crystal with denseIndex i\n"
"# (see EBDetId::denseIndex())\n"
"# for i >= 61200, element i applies to EE crystal with denseIndex (i+61200)\n"
"# (see EBDetId::denseIndex())\n"
"# SRP emulation supports only 1 element mode. Hardware does not support\n"
"# the symetric ZS, so symetricZS = 0 for real data.\n";
SR_VPRINT(symetricZS);
o << "\n# ZS energy threshold in GeV to apply to low interest channels of barrel\n"
"# If the vector contains...\n"
"# ... 1 element, then the weight set applies to whole ECAL\n"
"# ... 2 elements, then element 0 applies to EB, element 1 to EE\n"
"# ... 54 elements, then element i applied to DCC (i+1) (FED ID 651+i)\n"
"# SRP emulation supports only the 2-element mode.\n"
"# Corresponds to srpBarrelLowInterestChannelZS and srpEndcapLowInterestChannelZS\n"
"# of python configuration file parameters\n";
SR_VPRINT(srpLowInterestChannelZS);
o << "\n# ZS energy threshold in GeV to apply to high interest channels of endcap\n"
"# If the vector contains...\n"
"# ... 1 element, then the weight set applies to whole ECAL\n"
"# ... 2 elements, then element 0 applies to EB, element 1 to EE\n"
"# ... 54 elements, then element i applied to DCC (i+1) (FED ID 651+i)\n"
"# SRP emulation supports only the 2-element mode.\n"
"# Corresponds to srpBarrelLowInterestChannelZS and srpEndcapLowInterestChannelZS\n"
"# of python configuration file parameters\n";
SR_VPRINT(srpHighInterestChannelZS);
// o << "\n# Switch to run w/o trigger primitive. For debug use only\n"
// "# having troubles for vector<bool> with coral (3.8.0pre1), using vector<int> instead\n"
// "# Parameter only relevant for emulation. For real data, must be contains 1 element with\n"
// "# value 0.\n"
// "# ... 1 element, then the weight set applies to whole ECAL\n"
// "# ... 2 elements, then element 0 applies to EB, element 1 to EE\n"
// "# ... 54 elements, then element i applied to DCC (i+1) (FED ID 651+i)\n"
// "# SRP emulation supports only the single-element mode.\n";
// SR_VPRINT(trigPrimBypass);\n"
//
// o << "\n# Mode selection for "# Trig bypass" mode\n"
// "# 0: TT thresholds applied on sum of crystal Et's\n"
// "# 1: TT thresholds applies on compressed Et from Trigger primitive\n"
// "# @see trigPrimByPass switch\n"
// "# Parameter only relevant for \n";
// SR_VPRINT(trigPrimBypassMode);
//
// o << "\n# for debug mode only:\n";
// SR_VPRINT( trigPrimBypassLTH);
//
// o << "\n# for debug mode only:\n";
// SR_VPRINT(trigPrimBypassHTH);
//
// o << "\n# for debug mode only\n"
// "# having troubles for vector<bool> with coral (3.8.0pre1), using vector<int> instead\n";
// SR_VPRINT( trigPrimBypassWithPeakFinder);
//
// o << "\n# Trigger Tower Flag to use when a flag is not found from the input\n"
// "# Trigger Primitive collection. Must be one of the following values:\n"
// "# 0: low interest, 1: mid interest, 3: high interest\n"
// "# 4: forced low interest, 5: forced mid interest, 7: forced high interest\n";
// SR_VPRINT(defaultTtf);\n"
o << "\n# SR->action flag map. 4 elements\n"
"# action_[i]: action for flag value i\n";
SR_VPRINT(actions);
o << "\n# Masks for TTC inputs of SRP cards\n"
"# One element per TCC, that is 108 elements: element i applies to TCC (i+1)\n";
SR_VPRINT(tccMasksFromConfig);
o << "\n# Masks for SRP-SRP inputs of SRP cards\n"
"# One element per SRP, that is 12 elements: element i applies to SRP (i+1)\n"
"# indices: [iSrp][iCh]\n";
SR_VVPRINT(srpMasksFromConfig);
o << "\n# Masks for DCC output of SRP cards\n"
"# One element per DCC, that is 54 elements: element i applies to DCC (i+1)\n";
SR_VPRINT(dccMasks);
o << "\n# Mask to enable pattern test. Typical value: 0.\n"
"# One element per SRP, that is 12 elements: element i applies to SRP (i+1)\n";
SR_VPRINT(srfMasks);
o << "\n# Substitution flags used in patterm mode\n"
"# indices: [iSrp][iFlag]\n";
SR_VVPRINT(substitutionSrfs);
o << "\n# Tester mode configuration\n";
SR_VPRINT(testerTccEmuSrpIds);
SR_VPRINT(testerSrpEmuSrpIds);
SR_VPRINT(testerDccTestSrpIds);
SR_VPRINT(testerSrpTestSrpIds);
//@}
o << "\n# Per SRP card bunch crossing counter offset.\n"
"# This offset is added to the bxGlobalOffset\n";
SR_VPRINT(bxOffsets);
o << "\n# SRP system bunch crossing counter offset.\n"
"# For each card the bxOffset[i]\n"
"# is added to this one.\n";
SR_PRINT(bxGlobalOffset);
o << "\n# Switch for automatic channel masking. 0: disabled; 1: enabled. Standard configuration: 1.\n"
"# When enabled, if a FED is excluded from the run, the corresponding TCC inputs is automatically\n"
"# masked (overwrites the tccInputMasks).\n";
SR_PRINT(automaticMasks);
o << "\n# Switch for automatic SRP card selection. 0: disabled; 1 : enabled..\n"
"# When enabled, if all the FEDs corresponding to a given SRP is excluded from the run,\n"
"# Then the corresponding SRP card is automatically excluded.\n";
SR_PRINT(automaticSrpSelect);
return o;
}
|