ATLAS Offline Software
Loading...
Searching...
No Matches
testRead.cxx
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
3*/
4
5// $Id: testRead.cxx 770805 2016-08-30 14:03:33Z ssnyder $
11
12// STL includes:
13#include <iostream>
14#include <limits>
15
16// Root includes:
17#include <TFile.h>
18#include <THashList.h>
19
20// Gaudi includes:
21#if !defined(RINGER_STANDALONE) && !defined(XAOD_STANDALONE)
22# include "GaudiKernel/SystemOfUnits.h"
23using namespace Gaudi::Units;
24// We just do this ugly thing so that we can have a message stream
25// FIXME How can I create a dummy message stream?
26#define protected public
27#define private public
28# include "RingerSelectorTools/AsgElectronRingerSelector.h"
29#undef protected
30#undef private
31Ringer::AsgElectronRingerSelector __AsgRingerSel__("testRead");
32MsgStream &msg = __AsgRingerSel__.msg();
33#else
34#include "RingerSelectorTools/AsgElectronRingerSelector.h"
35#include "AsgMessaging/MsgStream.h"
36MsgStream msg("testRead");
37#endif
38
39#define BREAKER \
40 "================================================================================"
41
42// Local includes:
43#include "RingerSelectorTools/procedures/RingerProcedureWrapper.h"
44#include "RingerSelectorTools/procedures/Normalizations.h"
45#include "RingerSelectorTools/procedures/NeuralNetwork.h"
46#include "RingerSelectorTools/procedures/Thresholds.h"
47#include "RingerSelectorTools/tools/IOHelperFcns.h"
48#include "RingerSelectorTools/tools/cxx/remove_pointer.h"
49
50using namespace Ringer;
51
55template< class wrapper_t >
56void readWrapperFromFile(wrapper_t *&wrapper, const char* fileName);
57
62template<typename wrapperBase_t>
63void readCollectionFromFile(const char* fileName);
64
65template<>
66void readCollectionFromFile<IThresWrapper>(const char* fileName);
67
68int main(){
69
70 // Change message level and greet user
71 msg.setLevel(MSG::DEBUG);
72
73 msg << MSG::INFO << BREAKER << endmsg;
74 msg << MSG::INFO << "Reading discriminator using interface "
75 "without using normalization." << endmsg;
76 msg << MSG::INFO << BREAKER << endmsg;
77
78 // ----------------------------------------------------------------------------
79 // Read simple discriminator using interface (without using normalization)
80 // ----------------------------------------------------------------------------
81 // Read "basicWrapperWrong" using correct template
82 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
83 EtaIndependent,
84 EtIndependent,
85 NoSegmentation > *basicWrapper(nullptr);
86 readWrapperFromFile<
87 typename Ringer::remove_pointer<decltype(basicWrapper)>::type
88 >( basicWrapper, "basicWrapper.root" );
89 msg << MSG::INFO << BREAKER << endmsg;
90
91 // Read as "basicEtaDepWrapper" using correct template
92 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
93 EtaDependent,
94 EtIndependent,
95 NoSegmentation > *basicEtaDepWrapper(nullptr);
96 readWrapperFromFile<
97 typename Ringer::remove_pointer<decltype(basicEtaDepWrapper)>::type
98 >( basicEtaDepWrapper, "basicEtaDepWrapper.root");
99 msg << MSG::INFO << BREAKER << endmsg;
100
101 // Read "basicEtDepWrapper" using correct template
102 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
103 EtaIndependent,
104 EtDependent,
105 NoSegmentation > *basicEtDepWrapper(nullptr);
106 readWrapperFromFile<
107 typename Ringer::remove_pointer<decltype(basicEtDepWrapper)>::type
108 >( basicEtDepWrapper, "basicEtDepWrapper.root");
109 msg << MSG::INFO << BREAKER << endmsg;
110
111 // Read "basicEtaEtDepWrapper.root" using correct template
112 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
113 EtaDependent,
114 EtDependent,
115 NoSegmentation > *basicEtaEtDepWrapper(nullptr);
116 readWrapperFromFile<
117 typename Ringer::remove_pointer<decltype(basicEtaEtDepWrapper)>::type
118 >( basicEtaEtDepWrapper, "basicEtaEtDepWrapper.root");
119 // -------------------------- END --- END --- END -----------------------------
120
121 msg.setLevel(MSG::ERROR);
122
123 msg << MSG::INFO << BREAKER << endmsg;
124 msg << MSG::INFO << "Reading discriminator with wrong types!" << endmsg;
125 msg << MSG::INFO << BREAKER << endmsg;
126 // ----------------------------------------------------------------------------
127 // Read simple discriminator using interface with WRONG types
128 // ----------------------------------------------------------------------------
129 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
130 EtaDependent,
131 EtDependent,
132 NoSegmentation > *basicWrapperWrong(nullptr);
133 readWrapperFromFile<
134 typename Ringer::remove_pointer<decltype(basicWrapperWrong)>::type
135 >( basicWrapperWrong, "basicWrapper.root" );
136 msg << MSG::INFO << BREAKER << endmsg;
137
138 // Read as "basicEtaDepWrapperWrong" using wrong template
139 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
140 EtaIndependent,
141 EtIndependent,
142 NoSegmentation > *basicEtaDepWrapperWrong(nullptr);
143 readWrapperFromFile<
144 typename Ringer::remove_pointer<decltype(basicEtaDepWrapperWrong)>::type
145 >( basicEtaDepWrapperWrong, "basicEtaDepWrapper.root");
146 msg << MSG::INFO << BREAKER << endmsg;
147
148 // Read "basicEtDepWrapperWrong" using wrong template
149 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
150 EtaIndependent,
151 EtIndependent,
152 NoSegmentation > *basicEtDepWrapperWrong(nullptr);
153 readWrapperFromFile<
154 typename Ringer::remove_pointer<decltype(basicEtDepWrapperWrong)>::type
155 >( basicEtDepWrapperWrong, "basicEtDepWrapper.root");
156 msg << MSG::INFO << BREAKER << endmsg;
157
158 // Read "basicEtaEtDepWrapperWrong.root" using wrong template
159 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
160 EtaDependent,
161 EtDependent,
162 TrackCalJointLayers > *basicEtaEtDepWrapperWrong(nullptr);
163 readWrapperFromFile<
164 typename Ringer::remove_pointer<decltype(basicEtaEtDepWrapperWrong)>::type
165 >( basicEtaEtDepWrapperWrong, "basicEtaEtDepWrapper.root");
166 // -------------------------- END --- END --- END -----------------------------
167 msg.setLevel(MSG::DEBUG);
168
169 msg << MSG::INFO << BREAKER << endmsg;
170 msg << MSG::INFO << "Reading pre-processings!" << endmsg;
171 msg << MSG::INFO << BREAKER << endmsg;
172 // ----------------------------------------------------------------------------
173 // Start testing pre-processings
174 // ----------------------------------------------------------------------------
175 // Read "basicPPWrapper.root"
176 RingerProcedureWrapper<PreProcessing::IPreProcessorVarDep,
177 EtaIndependent,
178 EtIndependent,
179 NoSegmentation > *basicPPWrapper(nullptr);
180 readWrapperFromFile<
181 typename Ringer::remove_pointer<decltype(basicPPWrapper)>::type
182 >( basicPPWrapper, "basicPPWrapper.root");
183 msg << MSG::INFO << BREAKER << endmsg;
184
185 // Read "basicEtaDepPPWrapper.root"
186 RingerProcedureWrapper<PreProcessing::IPreProcessorVarDep,
187 EtaDependent,
188 EtIndependent,
189 NoSegmentation > *basicEtaDepPPWrapper(nullptr);
190 readWrapperFromFile<
191 typename Ringer::remove_pointer<decltype(basicEtaDepPPWrapper)>::type
192 >( basicEtaDepPPWrapper, "basicEtaDepPPWrapper.root");
193 msg << MSG::INFO << BREAKER << endmsg;
194
195 // Read "basicEtDepPPWrapper.root"
196 RingerProcedureWrapper<PreProcessing::IPreProcessorVarDep,
197 EtaIndependent,
198 EtDependent,
199 NoSegmentation > *basicEtDepPPWrapper(nullptr);
200 readWrapperFromFile<
201 typename Ringer::remove_pointer<decltype(basicEtDepPPWrapper)>::type
202 >( basicEtDepPPWrapper, "basicEtDepPPWrapper.root");
203 msg << MSG::INFO << BREAKER << endmsg;
204
205 // Read "basicEtEtaDepPPWrapper.root"
206 RingerProcedureWrapper<PreProcessing::IPreProcessorVarDep,
207 EtaDependent,
208 EtDependent,
209 NoSegmentation > *basicEtEtaDepPPWrapper(nullptr);
210 readWrapperFromFile<
211 typename Ringer::remove_pointer<decltype(basicEtEtaDepPPWrapper)>::type
212 >( basicEtEtaDepPPWrapper, "basicEtEtaDepPPWrapper.root");
213 msg << MSG::INFO << BREAKER << endmsg;
214
215 // Read "basicEtEtaTrackCalDepPPWrapper.root"
216 RingerProcedureWrapper<PreProcessing::IPreProcessorVarDep,
217 EtaDependent,
218 EtDependent,
219 TrackCalSegmentation > *basicEtEtaTrackCalDepPPWrapper(nullptr);
220 readWrapperFromFile<
221 typename Ringer::remove_pointer<decltype(basicEtEtaTrackCalDepPPWrapper)>::type
222 >( basicEtEtaTrackCalDepPPWrapper, "basicEtEtaTrackCalDepPPWrapper.root");
223 msg << MSG::INFO << BREAKER << endmsg;
224
225 // Read "basicEtEtaTrackCalPatTypeDepPPWrapper.root"
226 RingerProcedureWrapper<PreProcessing::IPreProcessorVarDep,
227 EtaDependent,
228 EtDependent,
229 TrackCalPatTypeSegmentation > *basicEtEtaTrackCalPatTypeDepPPWrapper(nullptr);
230 readWrapperFromFile<
231 typename Ringer::remove_pointer<decltype(basicEtEtaTrackCalPatTypeDepPPWrapper)>::type
232 >( basicEtEtaTrackCalPatTypeDepPPWrapper, "basicEtEtaTrackCalPatTypeDepPPWrapper.root");
233 msg << MSG::INFO << BREAKER << endmsg;
234
235 // Read "basicEtaEtSectionDepPPWrapper.root"
236 RingerProcedureWrapper<PreProcessing::IPreProcessorVarDep,
237 EtaDependent,
238 EtDependent,
239 TrackCalJointSections > *basicEtaEtSectionDepPPWrapper(nullptr);
240 readWrapperFromFile<
241 typename Ringer::remove_pointer<decltype(basicEtaEtSectionDepPPWrapper)>::type
242 >( basicEtaEtSectionDepPPWrapper, "basicEtaEtSectionDepPPWrapper.root");
243 msg << MSG::INFO << BREAKER << endmsg;
244
245 // Read "basicEtaEtLayerDepPPWrapper.root"
246 RingerProcedureWrapper<PreProcessing::IPreProcessorVarDep,
247 EtaDependent,
248 EtDependent,
249 TrackCalJointLayers > *basicEtaEtLayerDepPPWrapper(nullptr);
250 readWrapperFromFile<
251 typename Ringer::remove_pointer<decltype(basicEtaEtLayerDepPPWrapper)>::type
252 >( basicEtaEtLayerDepPPWrapper, "basicEtaEtLayerDepPPWrapper.root");
253 msg << MSG::INFO << BREAKER << endmsg;
254
255 // Read "basicNorm1Wrapper.root"
256 RingerProcedureWrapper<PreProcessing::IPreProcessorVarDep,
257 EtaIndependent,
258 EtIndependent,
259 NoSegmentation > *basicNorm1Wrapper(nullptr);
260 readWrapperFromFile<
261 typename Ringer::remove_pointer<decltype(basicNorm1Wrapper)>::type
262 >( basicNorm1Wrapper, "basicNorm1Wrapper.root");
263 msg << MSG::INFO << BREAKER << endmsg;
264
265 // Read "basicNorm1Wrapper.root"
266 RingerProcedureWrapper<PreProcessing::Norm::Norm1VarDep,
267 EtaDependent,
268 EtDependent,
269 TrackCalJointLayers > *basicLayerDepNorm1Wrapper(nullptr);
270 readWrapperFromFile<
271 typename Ringer::remove_pointer<decltype(basicLayerDepNorm1Wrapper)>::type
272 >( basicLayerDepNorm1Wrapper, "basicLayerDepNorm1Wrapper.root");
273 msg << MSG::INFO << BREAKER << endmsg;
274
275 // Read "basicNorm2Wrapper.root"
276 RingerProcedureWrapper<PreProcessing::Norm::Norm2VarDep,
277 EtaIndependent,
278 EtIndependent,
279 TrackCalJointLayers > *basicNorm2Wrapper(nullptr);
280 readWrapperFromFile<
281 typename Ringer::remove_pointer<decltype(basicNorm2Wrapper)>::type
282 >( basicNorm2Wrapper, "basicNorm2Wrapper.root");
283 msg << MSG::INFO << BREAKER << endmsg;
284
285 // Read "basicSqrtWrapper.root"
286 RingerProcedureWrapper<PreProcessing::Norm::SqrtVarDep,
287 EtaIndependent,
288 EtIndependent,
289 NoSegmentation > *basicSqrtWrapper(nullptr);
290 readWrapperFromFile<
291 typename Ringer::remove_pointer<decltype(basicSqrtWrapper)>::type
292 >( basicSqrtWrapper, "basicSqrtWrapper.root");
293 msg << MSG::INFO << BREAKER << endmsg;
294
295 // Read "basicConstantValueWrapper.root"
296 RingerProcedureWrapper<PreProcessing::Norm::ConstantValueVarDep,
297 EtaIndependent,
298 EtIndependent,
299 NoSegmentation > *basicConstantValueWrapper(nullptr);
300 readWrapperFromFile<
301 typename Ringer::remove_pointer<decltype(basicConstantValueWrapper)>::type
302 >( basicConstantValueWrapper, "basicConstantValueWrapper.root");
303 msg << MSG::INFO << BREAKER << endmsg;
304
305 // Read "basicMevToGevWrapper.root"
306 RingerProcedureWrapper<PreProcessing::Norm::MevToGevVarDep,
307 EtaIndependent,
308 EtIndependent,
309 NoSegmentation > *basicMevToGevWrapper(nullptr);
310 readWrapperFromFile<
311 typename Ringer::remove_pointer<decltype(basicMevToGevWrapper)>::type
312 >( basicMevToGevWrapper, "basicMevToGevWrapper.root");
313 msg << MSG::INFO << BREAKER << endmsg;
314
315 // Read "basicSequentialWrapper.root"
316 RingerProcedureWrapper<PreProcessing::Norm::SequentialVarDep,
317 EtaIndependent,
318 EtIndependent,
319 NoSegmentation > *basicSequentialWrapper(nullptr);
320 readWrapperFromFile<
321 typename Ringer::remove_pointer<decltype(basicSequentialWrapper)>::type
322 >( basicSequentialWrapper, "basicSequentialWrapper.root");
323 msg << MSG::INFO << BREAKER << endmsg;
324
325 // Read "basicSpherizationWrapper.root"
326 RingerProcedureWrapper<PreProcessing::Norm::SpherizationVarDep,
327 EtaIndependent,
328 EtIndependent,
329 NoSegmentation > *basicSpherizationWrapper(nullptr);
330 readWrapperFromFile<
331 typename Ringer::remove_pointer<decltype(basicSpherizationWrapper)>::type
332 >( basicSpherizationWrapper, "basicSpherizationWrapper.root");
333 msg << MSG::INFO << BREAKER << endmsg;
334
335 // Read "basicMinMaxWrapper.root"
336 RingerProcedureWrapper<PreProcessing::Norm::MinMaxVarDep,
337 EtaIndependent,
338 EtIndependent,
339 NoSegmentation > *basicMinMaxWrapper(nullptr);
340 readWrapperFromFile<
341 typename Ringer::remove_pointer<decltype(basicMinMaxWrapper)>::type
342 >( basicMinMaxWrapper, "basicMinMaxWrapper.root");
343 msg << MSG::INFO << BREAKER << endmsg;
344 // -------------------------- END --- END --- END -----------------------------
345
346 msg << MSG::INFO << BREAKER << endmsg;
347 msg << MSG::INFO << "Reading discriminators with pre-processings!" << endmsg;
348 msg << MSG::INFO << BREAKER << endmsg;
349
350 // ----------------------------------------------------------------------------
351 // Read discriminators with pre-processings
352 // ----------------------------------------------------------------------------
353 // Read "randomNorm1PPDiscrWrapper.root"
354 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
355 EtaDependent,
356 EtDependent,
357 NoSegmentation > *randomNorm1PPDiscrWrapper(nullptr);
358 readWrapperFromFile<
359 typename Ringer::remove_pointer<decltype(randomNorm1PPDiscrWrapper)>::type
360 >( randomNorm1PPDiscrWrapper, "randomNorm1PPDiscrWrapper.root");
361 msg << MSG::INFO << BREAKER << endmsg;
362
363 // Read "randomNorm1PPDiscrWrapper.root"
364 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
365 EtaDependent,
366 EtDependent,
367 NoSegmentation > *randomDepPPrandomDepPPDiscrWrapper(nullptr);
368 readWrapperFromFile<
369 typename Ringer::remove_pointer<decltype(randomDepPPrandomDepPPDiscrWrapper)>::type
370 >( randomDepPPrandomDepPPDiscrWrapper, "randomDepPPrandomDepPPDiscrWrapper.root");
371 msg << MSG::INFO << BREAKER << endmsg;
372 // ----------------------------------------------------------------------------
373
374 msg << MSG::INFO << BREAKER << endmsg;
375 msg << MSG::INFO << "Reading segmented discriminator PP wrappers:" << endmsg;
376 msg << MSG::INFO << BREAKER << endmsg;
377
378 // ----------------------------------------------------------------------------
379 // Read segmented discriminators with pre-processings
380 // ----------------------------------------------------------------------------
381 // Read "randomDepPPrandomDepPPDiscrSegmentedWrapper.root"
382 RingerProcedureWrapper<Discrimination::IDiscriminatorVarDep,
383 EtaDependent,
384 EtDependent,
385 TrackCalJointLayers > *randomDepPPrandomDepPPDiscrSegmentedWrapper(nullptr);
386 readWrapperFromFile<
387 typename Ringer::remove_pointer<decltype(randomDepPPrandomDepPPDiscrSegmentedWrapper)>::type
388 >( randomDepPPrandomDepPPDiscrSegmentedWrapper, "randomDepPPrandomDepPPDiscrSegmentedWrapper.root");
389 // -------------------------- END --- END --- END -----------------------------
390
391 msg << MSG::INFO << BREAKER << endmsg;
392 msg << MSG::INFO << "Reading specialized discriminator:" << endmsg;
393 msg << MSG::INFO << BREAKER << endmsg;
394
395 // ----------------------------------------------------------------------------
396 // Read specialized discriminator wrapper with no pre-processings
397 // ----------------------------------------------------------------------------
398 // Read "NNWrapper.root"
399 RingerProcedureWrapper<Discrimination::NNFeedForwardVarDep,
400 EtaIndependent,
401 EtIndependent,
402 NoSegmentation > *NNWrapper(nullptr);
403 readWrapperFromFile<
404 typename Ringer::remove_pointer<decltype(NNWrapper)>::type
405 >( NNWrapper, "NNWrapper.root");
406 msg << MSG::INFO << BREAKER << endmsg;
407
408 // Read "NNWrapper.root"
409 RingerProcedureWrapper<Discrimination::NNFeedForwardVarDep,
410 EtaDependent,
411 EtDependent,
412 TrackCalJointLayers > *FullDepNNWrapper(nullptr);
413 readWrapperFromFile<
414 typename Ringer::remove_pointer<decltype(FullDepNNWrapper)>::type
415 >( FullDepNNWrapper, "FullDepNNWrapper.root");
416 // -------------------------- END --- END --- END -----------------------------
417
418 msg << MSG::INFO << BREAKER << endmsg;
419 msg << MSG::INFO << "Reading specialized discriminator with PP wrappers:" << endmsg;
420 msg << MSG::INFO << BREAKER << endmsg;
421
422 // ----------------------------------------------------------------------------
423 // Read specialized discriminator wrapper with specialized pre-processings
424 // ----------------------------------------------------------------------------
425 // Read "Norm1NNWrapper.root"
426 RingerProcedureWrapper<Discrimination::NNFeedForwardVarDep,
427 EtaIndependent,
428 EtIndependent,
429 NoSegmentation > *Norm1NNWrapper(nullptr);
430 readWrapperFromFile<
431 typename Ringer::remove_pointer<decltype(Norm1NNWrapper)>::type
432 >( Norm1NNWrapper, "Norm1NNWrapper.root");
433 msg << MSG::INFO << BREAKER << endmsg;
434
435 // Read "Norm1Norm1FullDepNNWrapper.root"
436 RingerProcedureWrapper<Discrimination::NNFeedForwardVarDep,
437 EtaDependent,
438 EtDependent,
439 TrackCalJointLayers > *Norm1Norm1FullDepNNWrapper(nullptr);
440 readWrapperFromFile<
441 typename Ringer::remove_pointer<decltype(Norm1Norm1FullDepNNWrapper)>::type
442 >( Norm1Norm1FullDepNNWrapper, "Norm1Norm1FullDepNNWrapper.root");
443 // -------------------------- END --- END --- END -----------------------------
444
445 msg << MSG::INFO << BREAKER << endmsg;
446 msg << MSG::INFO << "Reading threshold wrappers:" << endmsg;
447 msg << MSG::INFO << BREAKER << endmsg;
448
449 // ----------------------------------------------------------------------------
450 // Read thresholds
451 // ----------------------------------------------------------------------------
452 // Read "basicThres.root"
453 RingerProcedureWrapper<Discrimination::IThresholdVarDep,
454 EtaIndependent,
455 EtIndependent > *basicThres(nullptr);
456 readWrapperFromFile<
457 typename Ringer::remove_pointer<decltype(basicThres)>::type
458 >( basicThres, "basicThres.root");
459 msg << MSG::INFO << BREAKER << endmsg;
460
461 // Read "basicEtaDepThres.root"
462 RingerProcedureWrapper<Discrimination::IThresholdVarDep,
463 EtaDependent,
464 EtIndependent > *basicEtaDepThres(nullptr);
465 readWrapperFromFile<
466 typename Ringer::remove_pointer<decltype(basicEtaDepThres)>::type
467 >( basicEtaDepThres, "basicEtaDepThres.root");
468 msg << MSG::INFO << BREAKER << endmsg;
469
470 // Read "basicEtDepThres.root"
471 RingerProcedureWrapper<Discrimination::IThresholdVarDep,
472 EtaIndependent,
473 EtDependent > *basicEtDepThres(nullptr);
474 readWrapperFromFile<
475 typename Ringer::remove_pointer<decltype(basicEtDepThres)>::type
476 >( basicEtDepThres, "basicEtDepThres.root");
477 msg << MSG::INFO << BREAKER << endmsg;
478
479 // Read "basicEtaEtDepThres.root"
480 RingerProcedureWrapper<Discrimination::IThresholdVarDep,
481 EtaDependent,
482 EtDependent > *basicEtaEtDepThres(nullptr);
483 readWrapperFromFile<
484 typename Ringer::remove_pointer<decltype(basicEtaEtDepThres)>::type
485 >( basicEtaEtDepThres, "basicEtaEtDepThres.root");
486 // -------------------------- END --- END --- END -----------------------------
487
488 msg << MSG::INFO << BREAKER << endmsg;
489 msg << MSG::INFO << "Reading specialized threshold wrappers:" << endmsg;
490 msg << MSG::INFO << BREAKER << endmsg;
491
492 // ----------------------------------------------------------------------------
493 // Read speaclized thresholds
494 // ----------------------------------------------------------------------------
495 // Read "uniqueThres.root"
496 RingerProcedureWrapper<Discrimination::UniqueThresholdVarDep,
497 EtaIndependent,
498 EtIndependent > *uniqueThres(nullptr);
499 readWrapperFromFile<
500 typename Ringer::remove_pointer<decltype(uniqueThres)>::type
501 >( uniqueThres, "uniqueThres.root");
502 msg << MSG::INFO << BREAKER << endmsg;
503
504 // Read "uniqueEtaDepThres.root"
505 RingerProcedureWrapper<Discrimination::UniqueThresholdVarDep,
506 EtaDependent,
507 EtIndependent > *uniqueEtaDepThres(nullptr);
508 readWrapperFromFile<
509 typename Ringer::remove_pointer<decltype(uniqueEtaDepThres)>::type
510 >( uniqueEtaDepThres, "uniqueEtaDepThres.root");
511 msg << MSG::INFO << BREAKER << endmsg;
512
513 // Read "uniqueEtDepThres.root"
514 RingerProcedureWrapper<Discrimination::UniqueThresholdVarDep,
515 EtaIndependent,
516 EtDependent > *uniqueEtDepThres(nullptr);
517 readWrapperFromFile<
518 typename Ringer::remove_pointer<decltype(uniqueEtDepThres)>::type
519 >( uniqueEtDepThres, "uniqueEtDepThres.root");
520 msg << MSG::INFO << BREAKER << endmsg;
521
522 // Read "uniqueEtaEtDepThres.root"
523 RingerProcedureWrapper<Discrimination::UniqueThresholdVarDep,
524 EtaDependent,
525 EtDependent > *uniqueEtaEtDepThres(nullptr);
526 readWrapperFromFile<
527 typename Ringer::remove_pointer<decltype(uniqueEtaEtDepThres)>::type
528 >( uniqueEtaEtDepThres, "uniqueEtaEtDepThres.root");
529 // -------------------------- END --- END --- END -----------------------------
530
531
532 msg << MSG::INFO << BREAKER << endmsg;
533 msg << MSG::INFO << "Special reading" << endmsg;
534 msg << MSG::INFO << BREAKER << endmsg;
535
536 // ----------------------------------------------------------------------------
537 // Test reading multiple wrappers
538 // ----------------------------------------------------------------------------
539 readCollectionFromFile<IDiscrWrapper>("Norm1NNWrapper.root");
540 msg << MSG::INFO << BREAKER << endmsg;
541 readCollectionFromFile<IDiscrWrapper>("Norm1Norm1FullDepNNWrapper.root");
542 msg << MSG::INFO << BREAKER << endmsg;
543 readCollectionFromFile<IThresWrapper>("uniqueEtaEtDepThres.root");
544 // -------------------------- END --- END --- END -----------------------------
545
546 msg << MSG::INFO << BREAKER << endmsg;
547 msg << MSG::INFO << "Finished!" << endmsg;
548 msg << MSG::INFO << BREAKER << endmsg;
549
550 return 0;
551}
552
553// =============================================================================
554template< class wrapper_t >
555void readWrapperFromFile(wrapper_t *&wrapper, const char* fileName)
556{
557
558 msg << MSG::INFO << "Reading wrapper: " << wrapper_t::staticFullName()
559 << ", at file named: " << fileName << endmsg;
560
561 // Read file and bring everything to memory:
562 TFile wrapperFile(fileName, "READ");
563 wrapperFile.ReadAll();
564 TDirectory *baseDir = wrapperFile.GetDirectory("");
565
566 std::shared_ptr<THashList> wrapperDirList(nullptr);
567 if ( !( wrapperDirList = IOHelperFcns::getDirList(baseDir) ) ) {
568 throw std::runtime_error(std::string("There are no directories available"));
569 }
570
571 msg << MSG::INFO << "Number(s) of folder(s) within this TFile is "
572 << wrapperDirList->GetSize()
573 << endmsg;
574
575 try {
576
577 TIter iter( wrapperDirList.get() );
578 TDirectory* wrapperDir(nullptr);
579
580 bool foundDir = false;
581
582 while ( (wrapperDir = static_cast<TDirectoryFile*>(iter())) ){
583 if ( wrapperDir->GetName() ==
584 ( std::string(wrapper_t::wrapName) +
585 IOHelperFcns::makeIdxStr(0) ) )
586 {
587 msg << MSG::INFO << "Found directory named "
588 << wrapperDir->GetName() << endmsg;
589 foundDir = true;
590 break;
591 }
592 }
593
594 if ( !foundDir ){
595 throw std::runtime_error(std::string("Couldn't find the wrapper directory"));
596 }
597
598 wrapper = wrapper_t::read(wrapperDir,IOHelperFcns::version());
599
600 if (wrapper == nullptr){
601 throw std::runtime_error(std::string("Returned void pointer."));
602 }
603
604 // Use dummy message stream:
605 wrapper->setMsgStream(&msg);
606
607 msg << MSG::INFO << "Printing wrapper " << endmsg;
608
609 // to print wrapper information:
610 wrapper->print();
611
612 } catch ( const std::runtime_error &e){
613 msg << MSG::ERROR << "Couldn't read due to error: "
614 << e.what() << endmsg;
615 }
616
617 // Close file
618 wrapperFile.Close();
619}
620
621// =============================================================================
622template<typename wrapperBase_t>
623void readCollectionFromFile(const char* fileName)
624{
625
626 msg << MSG::INFO << "Reading file: " << fileName << endmsg;
627
628 // Get wrapper collection and read it:
629 typename wrapperBase_t::WrapperCollection col;
630 wrapperBase_t::read( col, fileName );
631
632 unsigned counter(0);
633 for ( auto *wrapper : col ) {
634 if (wrapper){
635 msg << MSG::INFO << "Printing wrapper " << wrapper->fullName()
636 << " at position " << IOHelperFcns::makeIdxStr(counter++) << endmsg;
637 wrapper->setMsgStream(&msg);wrapper->print();
638 } else {
639 msg << MSG::ERROR << "Couldn't print wrapper on file " << fileName << endmsg;
640 }
641 }
642
643 msg << MSG::INFO << "File configuration: " << endmsg;
644 AsgElectronRingerSelector::IOConfStruct ioConf;
645 AsgElectronRingerSelector::retrieveFileConf( fileName, ioConf );
646 AsgElectronRingerSelector::printConf(ioConf, &msg);
647}
648
649// =============================================================================
650template<>
651void readCollectionFromFile<IThresWrapper>(const char* fileName)
652{
653
654 msg << MSG::INFO << "Reading file: " << fileName << endmsg;
655
656 // Get wrapper collection and read it:
657 IThresWrapper *wrapper;
658 IThresWrapper::read( wrapper, fileName );
659 if (wrapper) {wrapper->setMsgStream(&msg); wrapper->print();}
660 else msg << MSG::ERROR << "Couldn't print wrapper on file " << fileName << endmsg;
661}
662
663
664
endmsg
#define endmsg
Definition AnalysisConfig_Ntuple.cxx:63
Ringer::Discrimination::IDiscriminatorVarDep
Discriminator interface to be used by Ringer Wrapper.
Definition IRingerProcedure.h:311
Ringer::Discrimination::IThresholdVarDep
Threshold interface to be used by Ringer Wrapper.
Definition IRingerProcedure.h:349
Ringer::PreProcessing::IPreProcessorVarDep
PreProcessor interface to be used by Ringer Wrapper.
Definition IRingerProcedure.h:269
Ringer::PreProcessing::Norm::ConstantValueVarDep
Use constant value as rings normalization factor.
Definition Normalizations.h:365
Ringer::PreProcessing::Norm::MinMaxVarDep
Normalize data to be bounded within [-1,1] range.
Definition Normalizations.h:461
Ringer::PreProcessing::Norm::Norm2VarDep
Use norm-2 as rings normalization factor.
Definition Normalizations.h:329
Ringer::PreProcessing::Norm::SpherizationVarDep
Normalize data to have zero mean and unitary standard deviation.
Definition Normalizations.h:438
Ringer::PreProcessing::Norm::SqrtVarDep
Use SqrtVarDep of norm-1 as rings normalization factor.
Definition Normalizations.h:347
Ringer::RingerProcedureWrapper
Specialization for Discriminator procedure.
Definition RingerProcedureWrapper.h:171
Ringer
Namespace dedicated for Ringer utilities.
Definition CaloRingsDefs.h:9
type
msg
MsgStream & msg
Definition testRead.cxx:32
readWrapperFromFile
void readWrapperFromFile(wrapper_t *&wrapper, const char *fileName)
Read first wrapper on file with type wrapper_t.
Definition testRead.cxx:555
readCollectionFromFile
void readCollectionFromFile(const char *fileName)
Read wrapper collection at file named fileName with base wrapper type wrapperBase_t.
Definition testRead.cxx:623
__AsgRingerSel__
Ringer::AsgElectronRingerSelector __AsgRingerSel__("testRead")
readCollectionFromFile< IThresWrapper >
void readCollectionFromFile< IThresWrapper >(const char *fileName)
Definition testRead.cxx:651
BREAKER
#define BREAKER
Definition testRead.cxx:39
main
int main()
Definition testRead.cxx:68
Generated on for ATLAS Offline Software by doxygen 1.14.0