128 {
129
130
131 m_encoded_streamLength =
new TH1F(
"m_encoded_streamLength",
"Encoded Stream Length in bits", 2000, 0., 2000.);
132 if ((
m_thistSvc->regHist(
m_path + m_encoded_streamLength->GetName(), m_encoded_streamLength)).isFailure())
133 return StatusCode::FAILURE;
134
135 ATH_MSG_DEBUG(
"Histogram " << m_encoded_streamLength->GetName() <<
" successfully registered.");
136
137
138
139 std::array<std::vector<std::vector<double>>,
N_REGIONS>
bins;
141
142
144 if (not barrel_z[layer].
empty()) {
145 for (
unsigned int z_bin = 0; z_bin<(barrel_z[
layer].size()-1); z_bin++) {
146
147 float interModulePoint = 0.5 * ( barrel_z[
layer].at(z_bin) + barrel_z[layer].
at(z_bin+1) );
148 if (z_bin==0)
151 }
152
153 double last_value = 2.*barrel_z[
layer].back() - barrel_z[
layer].at( barrel_z[layer].
size() - 2 );
155 }
156
157 if (not endcap_z[layer].
empty()) {
158 for (
unsigned int z_bin = 0; z_bin<(endcap_z[
layer].size()-1); z_bin++) {
159
160 float interModulePoint = 0.5 * ( endcap_z[
layer].at(z_bin) + endcap_z[layer].
at(z_bin+1) );
161 if (z_bin==0) {
162 float initialValue = 2.*endcap_z[
layer].at(z_bin)-interModulePoint;
164 }
167 }
168
169 double last_value = 2.*endcap_z[
layer].back() - endcap_z[
layer].at( endcap_z[layer].
size() - 2);
171 }
172 }
173
174 for (
unsigned int region=0; region<
N_REGIONS; region++) {
176
177
179 continue;
180
182
185
186
187 m_p_streamLength[
layer][region][
side] =
new TProfile((
"m_p_streamLength_" + name_base).c_str(),
188 (title_base + " Stream Length; z[mm]; <stream length> [bits]").c_str(),
189 int(
bins[region][layer].
size()-1), &
bins[region][layer][0]);
190 if ((
m_thistSvc->regHist(
m_path + m_p_streamLength[layer][region][side]->GetName(), m_p_streamLength[layer][region][side])).isFailure())
191 return StatusCode::FAILURE;
192 ATH_MSG_DEBUG(
"Histogram " << m_p_streamLength[layer][region][side]->GetName() <<
" successfully registered.");
193
194 m_h2_streamLength[
layer][region][
side] =
new TH2F((
"m_h2_streamLength_" + name_base).c_str(),
195 (title_base + " Stream Length; z[mm]; <stream length> [bits]").c_str(),
196 int(
bins[region][layer].
size()-1), &
bins[region][layer][0],
197 1000, 0., 10000.);
198 if ((
m_thistSvc->regHist(
m_path + m_h2_streamLength[layer][region][side]->GetName(), m_h2_streamLength[layer][region][side])).isFailure())
199 return StatusCode::FAILURE;
200 ATH_MSG_DEBUG(
"Histogram " << m_h2_streamLength[layer][region][side]->GetName() <<
" successfully registered.");
201
202
203
204
205 m_p_dataRate[
layer][region][
side] =
new TProfile((
"m_p_dataRate_" + name_base).c_str(),
206 (title_base + " Data rate per chip; z[mm]; <data rate per chip> [b/s]").c_str(),
207 int(
bins[region][layer].
size()-1), &
bins[region][layer][0]);
208 if ((
m_thistSvc->regHist(
m_path + m_p_dataRate[layer][region][side]->GetName(), m_p_dataRate[layer][region][side])).isFailure())
209 return StatusCode::FAILURE;
210 ATH_MSG_DEBUG(
"Histogram " << m_p_dataRate[layer][region][side]->GetName() <<
" successfully registered.");
211
212 m_h2_dataRate[
layer][region][
side] =
new TH2F((
"m_h2_dataRate_" + name_base).c_str(),
213 (title_base + " Data rate per chip; z[mm]; <data rate per chip> [b/s]").c_str(),
214 int(
bins[region][layer].
size()-1), &
bins[region][layer][0],
215 1000, 0., 10000.);
216 if ((
m_thistSvc->regHist(
m_path + m_h2_dataRate[layer][region][side]->GetName(), m_h2_dataRate[layer][region][side])).isFailure())
217 return StatusCode::FAILURE;
218 ATH_MSG_DEBUG(
"Histogram " << m_h2_dataRate[layer][region][side]->GetName() <<
" successfully registered.");
219
220
222 (title_base + " Hits per chip; z[mm]; <hits/chip>").c_str(),
223 int(
bins[region][layer].
size()-1), &
bins[region][layer][0]);
224 if ((
m_thistSvc->regHist(
m_path + m_p_hits[layer][region][side]->GetName(), m_p_hits[layer][region][side])).isFailure())
225 return StatusCode::FAILURE;
226 ATH_MSG_DEBUG(
"Histogram " << m_p_hits[layer][region][side]->GetName() <<
" successfully registered.");
227
228 m_h2_hits[
layer][region][
side] =
new TH2F((
"m_h2_hits_" + name_base).c_str(),
229 (title_base + " Hits per chip; z[mm]; <hits/chip>").c_str(),
230 int(
bins[region][layer].
size()-1), &
bins[region][layer][0],
231 1000, 0., 10000.);
232 if ((
m_thistSvc->regHist(
m_path + m_h2_hits[layer][region][side]->GetName(), m_h2_hits[layer][region][side])).isFailure())
233 return StatusCode::FAILURE;
234 ATH_MSG_DEBUG(
"Histogram " << m_h2_hits[layer][region][side]->GetName() <<
" successfully registered.");
235 }
236 }
237 }
238
239 return StatusCode::SUCCESS;
240
241}
static const std::vector< std::string > bins
size_t size() const
Number of registered mappings.
static const Attributes_t empty
return m_collEvts back().back().max_entries
mapped_type at(key_type key) const
Look up an element in the map.
TH2F(name, title, nxbins, bins_par2, bins_par3, bins_par4, bins_par5=None, bins_par6=None, path='', **kwargs)
TProfile(*args, **kwargs)
TH1F(name, title, nxbins, bins_par2, bins_par3=None, path='', **kwargs)