|
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
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
|
#include <iostream>
#include <fstream>
#include <iomanip>
#include <string>
#include <map>
#include <cmath>
int main(int argc, char **argv) {
//std::string file = "TimingInfo.txt";
std::string file;
file.assign(argv[1]);
std::map<std::string,double> timingPerModule, timingPerLabel;
std::map<unsigned,double> timingPerEvent;
std::ifstream myTimingFile(file.c_str(),std::ifstream::in);
std::string dummy1, label, module;
double timing;
unsigned idummy1,evt;
int nbofevts = 0;
// If the machine is busy, the factor is not 100%.
double factor = 0.995;
if ( myTimingFile ) {
while ( !myTimingFile.eof() ) {
myTimingFile >> dummy1 >> evt >> idummy1 >> label >> module >> timing ;
// std::cout << evt << " " << module << " " << timing << std::endl;
timingPerEvent[evt] += timing * factor * 1000.;
timingPerModule[module] += timing * factor * 1000.;
timingPerLabel[module+":"+label] += timing * factor * 1000.;
}
nbofevts = (int) timingPerEvent.size();
} else {
std::cout << "File " << file << " does not exist!" << std::endl;
}
std::map<std::string,double>::const_iterator modIt = timingPerModule.begin();
std::map<std::string,double>::const_iterator labIt = timingPerLabel.begin();
std::map<std::string,double>::const_iterator modEnd = timingPerModule.end();
std::map<std::string,double>::const_iterator labEnd = timingPerLabel.end();
std::map<double,std::string> modulePerTiming;
std::map<double,std::string> labelPerTiming;
for ( ; modIt != modEnd; ++modIt ) {
double time = modIt->second/((double)nbofevts-1.);
std::string name = modIt->first;
modulePerTiming[time] = name;
}
for ( ; labIt != labEnd; ++labIt ) {
double time = labIt->second/((double)nbofevts-1.);
std::string name = labIt->first;
labelPerTiming[time] = name;
}
std::map<double,std::string>::const_reverse_iterator timeIt = modulePerTiming.rbegin();
std::map<double,std::string>::const_reverse_iterator timeEnd = modulePerTiming.rend();
std::map<double,std::string>::const_reverse_iterator timeIt2 = labelPerTiming.rbegin();
std::map<double,std::string>::const_reverse_iterator timeEnd2 = labelPerTiming.rend();
std::cout << "Timing per module " << std::endl;
std::cout << "================= " << std::endl;
double totalTime = 0.;
unsigned i=1;
for ( ; timeIt != timeEnd; ++timeIt ) {
totalTime += timeIt->first;
std::cout << std::setw(3) << i++
<< std::setw(50) << timeIt->second << " : "
<< std::setw(7) << std::setprecision(3) << timeIt-> first << " ms/event"
<< std::endl;
}
std::cout << "Total time = " << totalTime << " ms/event " << std::endl;
/*
std::cout << "================= " << std::endl;
std::cout << "Timing per label " << std::endl;
std::cout << "================= " << std::endl;
totalTime = 0.;
i = 1;
for ( ; timeIt2 != timeEnd2; ++timeIt2 ) {
totalTime += timeIt2->first;
std::cout << std::setw(3) << i++
<< std::setw(100) << timeIt2->second << " : "
<< std::setw(7) << std::setprecision(3) << timeIt2-> first << " ms/event"
<< std::endl;
}
*/
double subtotaltimepermodule = 0;
double cumulativetimepermodule = 0;
std::cout << "================= " << std::endl;
std::cout << " DQM for collision : Timing per step " << std::endl;
std::cout << "================= " << std::endl;
std::cout << "1. Reconstruction " << std::endl;
std::cout << " 1.1 Raw2Digi+LocalReco : " << std::endl;
std::cout << " 1.1.1 : Raw2Digi " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "SiPixelRawToDigi" << std::setw(30) << "" << std::setw(8) << timingPerModule["SiPixelRawToDigi"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SiPixelRawToDigi"];
cumulativetimepermodule += timingPerModule["SiPixelRawToDigi"];
std::cout << " - " << std::setw(40) << "SiStripRawToDigi" << std::setw(30) << "" << std::setw(8) << timingPerModule["SiStripRawToDigiModule"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SiStripRawToDigiModule"];
cumulativetimepermodule += timingPerModule["SiStripRawToDigiModule"];
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << " 1.1.2 : LocalReco" << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "SiPixelClusterProducer" << std::setw(30) << "" << std::setw(8) << timingPerModule["SiPixelClusterProducer"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SiPixelClusterProducer"];
cumulativetimepermodule += timingPerModule["SiPixelClusterProducer"];
std::cout << " - " << std::setw(40) << "SiPixelRecHitConverter" << std::setw(30) << "" << std::setw(8) << timingPerLabel["SiPixelRecHitConverter:siPixelRecHits"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerLabel["SiPixelRecHitConverter:siPixelRecHits"];
cumulativetimepermodule += timingPerLabel["SiPixelRecHitConverter:siPixelRecHits"];
std::cout << " - " << std::setw(40) << "SiStripZeroSuppression" << std::setw(30) << "" << std::setw(8) << timingPerModule["SiStripZeroSuppression"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SiStripZeroSuppression"];
cumulativetimepermodule += timingPerModule["SiStripZeroSuppression"];
std::cout << " - " << std::setw(40) << "SiStripClusterizer" << std::setw(30) << "" << std::setw(8) << timingPerModule["SiStripClusterizer"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SiStripClusterizer"];
cumulativetimepermodule += timingPerModule["SiStripClusterizer"];
std::cout << " - " << std::setw(40) << "SiStripRecHitConverter" << std::setw(30) << "" << std::setw(8) << timingPerLabel["SiStripRecHitConverter:siStripMatchedRecHits"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerLabel["SiStripRecHitConverter:siStripMatchedRecHits"];
cumulativetimepermodule += timingPerLabel["SiStripRecHitConverter:siStripMatchedRecHits"];
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << " 1.2 BeamSpot+RecoPixelVertexing+CkfTracks :" << std::endl;
std::cout << " 1.2.1 : BeamSpot " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "BeamSpotProducer" << std::setw(30) << "" << std::setw(8) << timingPerModule["BeamSpotProducer"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["BeamSpotProducer"];
cumulativetimepermodule += timingPerModule["BeamSpotProducer"];
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << " 1.2.2 : RecoPixelVertexing " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "PixelTrackProducer" << std::setw(30) << "" << std::setw(8) << timingPerModule["PixelTrackProducer"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["PixelTrackProducer"];
cumulativetimepermodule += timingPerModule["PixelTrackProducer"];
std::cout << " - " << std::setw(40) << "PixelVertexProducer" << std::setw(30) << "" << std::setw(8) << timingPerModule["PixelVertexProducer"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["PixelVertexProducer"];
cumulativetimepermodule += timingPerModule["PixelVertexProducer"];
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << " 1.2.3 : CkfTracks " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "SeedGeneratorFromRegionHitsEDProducer" << std::setw(30) << "" << std::setw(8) << timingPerModule["SeedGeneratorFromRegionHitsEDProducer"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SeedGeneratorFromRegionHitsEDProducer"];
cumulativetimepermodule += timingPerModule["SeedGeneratorFromRegionHitsEDProducer"];
/*
- newSeedFromTriplets
- newSeedFromPairs
- secTriplets
- thPLSeeds
- fourthPLSeeds
- fifthSeeds
*/
std::cout << " - " << std::setw(40) << "CkfTrackCandidateMaker" << std::setw(30) << "" << std::setw(8) << timingPerModule["CkfTrackCandidateMaker"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["CkfTrackCandidateMaker"];
cumulativetimepermodule += timingPerModule["CkfTrackCandidateMaker"];
/*
- newTrackCandidateMaker
- secTrackCandidates
- thTrackCandidates
- fourthTrackCandidates
- fifthTrackCandidates
*/
std::cout << " - " << std::setw(40) << "TrackProducer" << std::setw(30) << "" << std::setw(8) << timingPerModule["TrackProducer"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["TrackProducer"];
cumulativetimepermodule += timingPerModule["TrackProducer"];
/*
- preFilterZeroStepTracks
- preFilterStepOneTracks
- secWithMaterialTracks
- thWithMaterialTracks
- fourthWithMaterialTracks
- fifthWithMaterialTracks
*/
std::cout << " - " << std::setw(40) << "AnalyticalTrackSelector" << std::setw(30) << "" << std::setw(8) << timingPerModule["AnalyticalTrackSelector"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["AnalyticalTrackSelector"];
cumulativetimepermodule += timingPerModule["AnalyticalTrackSelector"];
/*
- zeroStepWithLooseQuality
- zeroStepWithTightQuality
- zeroStepTracksWithQuality
- firstStepWithLooseQuality
- firstStepWithTightQuality
- preMergingFirstStepTracksWithQuality
- secStepVtxLoose
- secStepTrkLoose
- secStepVtxTight
- secStepTrkTight
- secStepVtx
- secStepTrk
- thStepVtxLoose
- thStepTrkLoose
- thStepVtxTight
- thStepTrkTight
- thStepVtx
- thStepTrk
- pixellessStepLoose
- pixellessStepTight
- pixellessStep
- tobtecStepLoose
- tobtecStepTight
- tobtecStep
*/
std::cout << " - " << std::setw(40) << "QualityFilter" << std::setw(30) << "" << std::setw(8) << timingPerModule["QualityFilter"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["QualityFilter"];
cumulativetimepermodule += timingPerModule["QualityFilter"];
/*
- zeroStepFilter
- firstfilter
- secfilter
- thfilter
- fourthfilter
*/
std::cout << " - " << std::setw(40) << "TrackClusterRemover" << std::setw(30) << "" << std::setw(8) << timingPerModule["TrackClusterRemover"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["TrackClusterRemover"];
cumulativetimepermodule += timingPerModule["TrackClusterRemover"];
/*
- newClusters
- secClusters
- thClusters
- fourthClusters
- fifthClusters
*/
std::cout << " - " << std::setw(40) << "SiPixelRecHitConverter" << std::setw(30) << "" << std::setw(8) << timingPerLabel["SiPixelRecHitConverter:newPixelRecHits"]/((double)nbofevts-1.)
+timingPerLabel["SiPixelRecHitConverter:secPixelRecHits"]/((double)nbofevts-1.)
+timingPerLabel["SiPixelRecHitConverter:thPixelRecHits"]/((double)nbofevts-1.)
+timingPerLabel["SiPixelRecHitConverter:fourthPixelRecHits"]/((double)nbofevts-1.)
+timingPerLabel["SiPixelRecHitConverter:fifthPixelRecHits"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerLabel["SiPixelRecHitConverter:newPixelRecHits"];
subtotaltimepermodule += timingPerLabel["SiPixelRecHitConverter:secPixelRecHits"];
subtotaltimepermodule += timingPerLabel["SiPixelRecHitConverter:thPixelRecHits"];
subtotaltimepermodule += timingPerLabel["SiPixelRecHitConverter:fourthPixelRecHits"];
subtotaltimepermodule += timingPerLabel["SiPixelRecHitConverter:fifthPixelRecHits"];
cumulativetimepermodule += timingPerLabel["SiPixelRecHitConverter:newPixelRecHits"];
cumulativetimepermodule += timingPerLabel["SiPixelRecHitConverter:secPixelRecHits"];
cumulativetimepermodule += timingPerLabel["SiPixelRecHitConverter:thPixelRecHits"];
cumulativetimepermodule += timingPerLabel["SiPixelRecHitConverter:fourthPixelRecHits"];
cumulativetimepermodule += timingPerLabel["SiPixelRecHitConverter:fifthPixelRecHits"];
/*
- newPixelRecHits
- secPixelRecHits
- thPixelRecHits
- fourthPixelRecHits
- fifthPixelRecHits
*/
std::cout << " - " << std::setw(40) << "SiStripRecHitConverter" << std::setw(30) << "" << std::setw(8) << timingPerLabel["SiStripRecHitConverter:newStripRecHits"]/((double)nbofevts-1.)
+timingPerLabel["SiStripRecHitConverter:secStripRecHits"]/((double)nbofevts-1.)
+timingPerLabel["SiStripRecHitConverter:thStripRecHits"]/((double)nbofevts-1.)
+timingPerLabel["SiStripRecHitConverter:fourthStripRecHits"]/((double)nbofevts-1.)
+timingPerLabel["SiStripRecHitConverter:fifthStripRecHits"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerLabel["SiStripRecHitConverter:newStripRecHits"];
subtotaltimepermodule += timingPerLabel["SiStripRecHitConverter:secStripRecHits"];
subtotaltimepermodule += timingPerLabel["SiStripRecHitConverter:thStripRecHits"];
subtotaltimepermodule += timingPerLabel["SiStripRecHitConverter:fourthStripRecHits"];
subtotaltimepermodule += timingPerLabel["SiStripRecHitConverter:fifthStripRecHits"];
cumulativetimepermodule += timingPerLabel["SiStripRecHitConverter:newStripRecHits"];
cumulativetimepermodule += timingPerLabel["SiStripRecHitConverter:secStripRecHits"];
cumulativetimepermodule += timingPerLabel["SiStripRecHitConverter:thStripRecHits"];
cumulativetimepermodule += timingPerLabel["SiStripRecHitConverter:fourthStripRecHits"];
cumulativetimepermodule += timingPerLabel["SiStripRecHitConverter:fifthStripRecHits"];
/*
- newStripRecHits
- secStripRecHits
- thStripRecHits
- fourthStripRecHits
- fifthStripRecHits
*/
std::cout << " - " << std::setw(40) << "SimpleTrackListMerger" << std::setw(30) << "" << std::setw(8) << timingPerModule["SimpleTrackListMerger"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SimpleTrackListMerger"];
cumulativetimepermodule += timingPerModule["SimpleTrackListMerger"];
/*
- merge2nd3rdTracks
- merge4th5thTracks
- iterTracks
- generalTracks
*/
std::cout << " - " << std::setw(40) << "SeedCombiner" << std::setw(30) << "" << std::setw(8) << timingPerModule["SeedCombiner"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SeedCombiner"];
cumulativetimepermodule += timingPerModule["SeedCombiner"];
/*
- newCombinedSeeds
*/
std::cout << " - " << std::setw(40) << "DeDxEstimatorProducer" << std::setw(30) << "" << std::setw(8) << timingPerModule["DeDxEstimatorProducer"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["DeDxEstimatorProducer"];
cumulativetimepermodule += timingPerModule["DeDxEstimatorProducer"];
/*
- dedxMedian
- dedxHarmonic2
*/
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << "2. Data quality monitoring : " << std::endl;
std::cout << " 2.1 DQM common modules " << std::endl;
std::cout << " 2.1.1 : Quality tests " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "" << std::setw(30) << "" << std::setw(8) << timingPerModule["QualityTester"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["QualityTester"];
cumulativetimepermodule += timingPerModule["QualityTester"];
/*
- qTester
*/
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << " 2.1.2 : DQM playback environment " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "DQMEventInfo" << std::setw(30) << "" << std::setw(8) << timingPerLabel["DQMEventInfo:dqmEnv"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerLabel["DQMEventInfo:dqmEnv"];
cumulativetimepermodule += timingPerLabel["DQMEventInfo:dqmEnv"];
/*
- dqmEnv
*/
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << " 2.1.3 : DQM playback for Tracking info " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "DQMEventInfo" << std::setw(30) << "" << std::setw(8) << timingPerLabel["DQMEventInfo:dqmEnvTr"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerLabel["DQMEventInfo:dqmEnvTr"];
cumulativetimepermodule += timingPerLabel["DQMEventInfo:dqmEnvTr"];
/*
- dqmEnvTr
*/
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << " 2.1.4 : DQM file saver " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "DQMFileSaver" << std::setw(30) << "" << std::setw(8) << timingPerModule["DQMFileSaver"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["DQMFileSaver"];
cumulativetimepermodule += timingPerModule["DQMFileSaver"];
/*
- dqmSaver
*/
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << " 2.2 DQM monitor " << std::endl;
std::cout << " 2.2.1 : SiStripMonitor " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "SiStripFEDMonitorPlugin" << std::setw(30) << "" << std::setw(8) << timingPerModule["SiStripFEDMonitorPlugin"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SiStripFEDMonitorPlugin"];
cumulativetimepermodule += timingPerModule["SiStripFEDMonitorPlugin"];
/*
- siStripFEDMonitor
*/
std::cout << " - " << std::setw(40) << "SiStripMonitorDigi" << std::setw(30) << "" << std::setw(8) << timingPerModule["SiStripMonitorDigi"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SiStripMonitorDigi"];
cumulativetimepermodule += timingPerModule["SiStripMonitorDigi"];
/*
- SiStripMonitorDigi
*/
std::cout << " - " << std::setw(40) << "SiStripMonitorCluster" << std::setw(30) << "" << std::setw(8) << timingPerModule["SiStripMonitorCluster"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SiStripMonitorCluster"];
cumulativetimepermodule += timingPerModule["SiStripMonitorCluster"];
/*
- SiStripMonitorClusterReal
*/
std::cout << " - " << std::setw(40) << "SiStripMonitorTrack" << std::setw(30) << "" << std::setw(8) << timingPerModule["SiStripMonitorTrack"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SiStripMonitorTrack"];
cumulativetimepermodule += timingPerModule["SiStripMonitorTrack"];
/*
- SiStripMonitorTrack_gentk
*/
std::cout << " - " << std::setw(40) << "MonitorTrackResiduals" << std::setw(30) << "" << std::setw(8) << timingPerModule["MonitorTrackResiduals"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["MonitorTrackResiduals"];
cumulativetimepermodule += timingPerModule["MonitorTrackResiduals"];
/*
- MonitorTrackResiduals_gentk
*/
std::cout << " - " << std::setw(40) << "TrackingMonitor" << std::setw(30) << "" << std::setw(8) << timingPerModule["TrackingMonitor"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["TrackingMonitor"];
cumulativetimepermodule += timingPerModule["TrackingMonitor"];
/*
- TrackMon_gentk
*/
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << " 2.2.2 : SiStripAnalyser " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "SiStripAnalyser" << std::setw(30) << "" << std::setw(8) << timingPerModule["SiStripAnalyser"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["SiStripAnalyser"];
cumulativetimepermodule += timingPerModule["SiStripAnalyser"];
/*
- SiStripAnalyser
*/
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << " 2.2.3 : Miscellaneous " << std::endl;
subtotaltimepermodule = 0;
std::cout << " - " << std::setw(40) << "TriggerResultInserter" << std::setw(30) << "" << std::setw(8) << timingPerModule["TriggerResultInserter"]/((double)nbofevts-1.) << " ms/event" << std::endl;
subtotaltimepermodule += timingPerModule["TriggerResultInserter"];
cumulativetimepermodule += timingPerModule["TriggerResultInserter"];
/*
- TriggerResults
*/
std::cout << " " << std::setw(70) << "" << std::setw(8) << "--------" << std::endl;
std::cout << " " << std::setw(70) << "subtotal : " << std::setw(8) << subtotaltimepermodule/((double)nbofevts-1.) << " ms/event" << " / " << std::setw(8) << cumulativetimepermodule/((double)nbofevts-1.) << " ms/event" << std::endl;
std::cout << "Total nb of events read = " << nbofevts << std::endl;
std::cout << "Total time = " << totalTime << " ms/event " << std::endl;
std::map<unsigned,double>::const_iterator eventIt = timingPerEvent.begin();
std::map<unsigned,double>::const_iterator eventEnd = timingPerEvent.end();
double minEv = 99999999.;
double maxEv = 0.;
double rms = 0.;
double mean = 0.;
double timeEv = 0;
for ( ; eventIt != eventEnd; ++eventIt ) {
if ( eventIt->first == 1 ) continue;
timeEv = eventIt->second;
//std::cout << "Evt nr : " << eventIt->first << " / " << timeEv << " ms" << std::endl;
if ( timeEv > maxEv ) maxEv = timeEv;
if ( timeEv < minEv ) minEv = timeEv;
mean += timeEv;
rms += timeEv*timeEv;
}
mean /= ((double)nbofevts-1.);
rms /= ((double)nbofevts-1.);
rms = std::sqrt(rms-mean*mean);
std::cout << "Total time = " << mean << " +/- " << rms << " ms/event" << std::endl;
std::cout << "Min. time = " << minEv << " ms/event" << std::endl;
std::cout << "Max. time = " << maxEv << " ms/event" << std::endl;
}
|