| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | /// @begin_doc_example{cpp} | ||
| 2 | /// @file iguana_ex_cpp_00_run_functions.cc | ||
| 3 | /// @brief Example using **full HIPO banks** with Iguana algorithms' `Run` functions. This example requires the | ||
| 4 | /// user to have the C++ `hipo::bank` objects; see other examples if you do not have banks in this format. | ||
| 5 | /// @par Usage | ||
| 6 | /// ```bash | ||
| 7 | /// iguana_ex_cpp_00_run_functions [HIPO_FILE] [NUM_EVENTS] | ||
| 8 | /// | ||
| 9 | /// HIPO_FILE the HIPO file to analyze | ||
| 10 | /// | ||
| 11 | /// NUM_EVENTS the number of events to analyze; | ||
| 12 | /// set to zero to analyze all events | ||
| 13 | /// ``` | ||
| 14 | /// @end_doc_example | ||
| 15 | |||
| 16 | #include <hipo4/reader.h> | ||
| 17 | #include <iguana/algorithms/AlgorithmSequence.h> | ||
| 18 | |||
| 19 | /// @brief show a bank along with a header | ||
| 20 | /// @param header the header to print above the bank | ||
| 21 | /// @param bank the bank to show | ||
| 22 | 200 | void prettyPrint(std::string header, hipo::bank& bank) | |
| 23 | { | ||
| 24 | 400 | fmt::print("{:=^70}\n", " " + header + " "); | |
| 25 | 200 | bank.show(); | |
| 26 | 200 | } | |
| 27 | |||
| 28 | /// main function | ||
| 29 | 1 | int main(int argc, char** argv) | |
| 30 | { | ||
| 31 | |||
| 32 | // parse arguments | ||
| 33 | int argi = 1; | ||
| 34 |
1/2✓ Branch 0 (2→3) taken 1 times.
✗ Branch 1 (2→8) not taken.
|
1 | char const* inFileName = argc > argi ? argv[argi++] : "data.hipo"; |
| 35 |
1/2✓ Branch 0 (3→4) taken 1 times.
✗ Branch 1 (3→8) not taken.
|
1 | int const numEvents = argc > argi ? std::stoi(argv[argi++]) : 1; |
| 36 | |||
| 37 | // read input file | ||
| 38 |
1/2✓ Branch 0 (9→10) taken 1 times.
✗ Branch 1 (9→81) not taken.
|
1 | hipo::reader reader(inFileName,{0}); |
| 39 | |||
| 40 | // set list of banks to be read | ||
| 41 |
1/2✓ Branch 0 (13→14) taken 1 times.
✗ Branch 1 (13→125) not taken.
|
2 | hipo::banklist banks = reader.getBanks({"RUN::config", |
| 42 | "REC::Particle", | ||
| 43 | "REC::Calorimeter", | ||
| 44 | "REC::Track", | ||
| 45 |
1/2✓ Branch 0 (14→15) taken 1 times.
✗ Branch 1 (14→85) not taken.
|
1 | "REC::Scintillator"}); |
| 46 | |||
| 47 | // get bank index, for each bank we want to use after Iguana algorithms run | ||
| 48 |
3/6✓ Branch 0 (16→17) taken 1 times.
✗ Branch 1 (16→123) not taken.
✓ Branch 2 (17→18) taken 1 times.
✗ Branch 3 (17→87) not taken.
✓ Branch 4 (19→20) taken 1 times.
✗ Branch 5 (19→123) not taken.
|
2 | auto b_particle = hipo::getBanklistIndex(banks, "REC::Particle"); |
| 49 | |||
| 50 | // iguana algorithm sequence | ||
| 51 |
1/2✓ Branch 0 (19→20) taken 1 times.
✗ Branch 1 (19→123) not taken.
|
1 | iguana::AlgorithmSequence seq; |
| 52 |
3/6✓ Branch 0 (20→21) taken 1 times.
✗ Branch 1 (20→121) not taken.
✓ Branch 2 (21→22) taken 1 times.
✗ Branch 3 (21→91) not taken.
✓ Branch 4 (22→23) taken 1 times.
✗ Branch 5 (22→89) not taken.
|
2 | seq.Add("clas12::EventBuilderFilter"); // filter by Event Builder PID |
| 53 |
3/6✓ Branch 0 (25→26) taken 1 times.
✗ Branch 1 (25→121) not taken.
✓ Branch 2 (26→27) taken 1 times.
✗ Branch 3 (26→95) not taken.
✓ Branch 4 (27→28) taken 1 times.
✗ Branch 5 (27→93) not taken.
|
2 | seq.Add("clas12::SectorFinder"); // get the sector for each particle |
| 54 |
3/6✓ Branch 0 (30→31) taken 1 times.
✗ Branch 1 (30→121) not taken.
✓ Branch 2 (31→32) taken 1 times.
✗ Branch 3 (31→99) not taken.
✓ Branch 4 (32→33) taken 1 times.
✗ Branch 5 (32→97) not taken.
|
2 | seq.Add("clas12::MomentumCorrection"); // momentum corrections |
| 55 | // TODO: | ||
| 56 | // - getRowList for filter | ||
| 57 | // - add a creator algo | ||
| 58 | |||
| 59 | // set log levels | ||
| 60 |
2/4✓ Branch 0 (35→36) taken 1 times.
✗ Branch 1 (35→121) not taken.
✓ Branch 2 (36→37) taken 1 times.
✗ Branch 3 (36→103) not taken.
|
2 | seq.SetOption("clas12::EventBuilderFilter", "log", "debug"); |
| 61 |
2/4✓ Branch 0 (42→43) taken 1 times.
✗ Branch 1 (42→121) not taken.
✓ Branch 2 (43→44) taken 1 times.
✗ Branch 3 (43→107) not taken.
|
2 | seq.SetOption("clas12::MomentumCorrection", "log", "debug"); |
| 62 | |||
| 63 | // set algorithm options (overrides configuration files) | ||
| 64 |
5/12✓ Branch 0 (49→50) taken 1 times.
✗ Branch 1 (49→121) not taken.
✓ Branch 2 (50→51) taken 1 times.
✗ Branch 3 (50→113) not taken.
✓ Branch 4 (51→52) taken 1 times.
✗ Branch 5 (51→111) not taken.
✓ Branch 6 (52→53) taken 1 times.
✗ Branch 7 (52→109) not taken.
✓ Branch 8 (55→56) taken 1 times.
✗ Branch 9 (55→58) not taken.
✗ Branch 10 (113→114) not taken.
✗ Branch 11 (113→116) not taken.
|
2 | seq.SetOption<std::vector<int>>("clas12::EventBuilderFilter", "pids", {11, 211, -211}); |
| 65 | |||
| 66 | // start the algorithms | ||
| 67 |
1/2✓ Branch 0 (58→59) taken 1 times.
✗ Branch 1 (58→121) not taken.
|
1 | seq.Start(banks); |
| 68 | |||
| 69 | // run the algorithm sequence on each event | ||
| 70 | int iEvent = 0; | ||
| 71 |
5/8✓ Branch 0 (69→70) taken 101 times.
✗ Branch 1 (69→121) not taken.
✓ Branch 2 (70→71) taken 101 times.
✗ Branch 3 (70→74) not taken.
✓ Branch 4 (71→72) taken 101 times.
✗ Branch 5 (71→73) not taken.
✓ Branch 6 (72→73) taken 100 times.
✓ Branch 7 (72→74) taken 1 times.
|
101 | while(reader.next(banks) && (numEvents == 0 || iEvent++ < numEvents)) { |
| 72 |
3/6✓ Branch 0 (60→61) taken 100 times.
✗ Branch 1 (60→121) not taken.
✓ Branch 2 (61→62) taken 100 times.
✗ Branch 3 (61→117) not taken.
✓ Branch 4 (73→60) taken 100 times.
✗ Branch 5 (73→121) not taken.
|
100 | prettyPrint("BEFORE", banks.at(b_particle)); |
| 73 |
1/2✓ Branch 0 (63→64) taken 100 times.
✗ Branch 1 (63→121) not taken.
|
100 | seq.Run(banks); |
| 74 |
3/6✓ Branch 0 (64→65) taken 100 times.
✗ Branch 1 (64→121) not taken.
✓ Branch 2 (65→66) taken 100 times.
✗ Branch 3 (65→121) not taken.
✓ Branch 4 (66→67) taken 100 times.
✗ Branch 5 (66→119) not taken.
|
200 | prettyPrint("AFTER", banks.at(b_particle)); |
| 75 | } | ||
| 76 | |||
| 77 | // stop algorithms | ||
| 78 |
1/2✓ Branch 0 (74→75) taken 1 times.
✗ Branch 1 (74→121) not taken.
|
1 | seq.Stop(); |
| 79 | return 0; | ||
| 80 | 1 | } | |
| 81 |