| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | #pragma once | ||
| 2 | |||
| 3 | #include "iguana/algorithms/Algorithm.h" | ||
| 4 | |||
| 5 | namespace iguana::clas12 { | ||
| 6 | |||
| 7 | /// @brief_algo Find the sector for all rows in `REC::Particle` | ||
| 8 | /// | ||
| 9 | /// @begin_doc_algo{clas12::SectorFinder | Creator} | ||
| 10 | /// @input_banks{REC::Particle} | ||
| 11 | /// @output_banks{%REC::Particle::Sector} | ||
| 12 | /// @end_doc | ||
| 13 | /// | ||
| 14 | /// @begin_doc_config{clas12/SectorFinder} | ||
| 15 | /// @config_param{bank_charged | string | if not `default`, use this bank for sector finding of charged particles} | ||
| 16 | /// @config_param{bank_uncharged | string | if not `default`, use this bank for sector finding of neutral particles} | ||
| 17 | /// @end_doc | ||
| 18 | /// | ||
| 19 | /// @creator_note | ||
| 20 | class SectorFinder : public Algorithm | ||
| 21 | { | ||
| 22 | |||
| 23 |
6/14✓ Branch 0 (2→3) taken 5 times.
✗ Branch 1 (2→6) not taken.
✓ Branch 2 (7→8) taken 5 times.
✗ Branch 3 (7→19) not taken.
✓ Branch 4 (8→9) taken 5 times.
✗ Branch 5 (8→11) not taken.
✓ Branch 6 (11→12) taken 5 times.
✗ Branch 7 (11→25) not taken.
✓ Branch 8 (12→13) taken 5 times.
✗ Branch 9 (12→23) not taken.
✓ Branch 10 (15→16) taken 5 times.
✗ Branch 11 (15→23) not taken.
✗ Branch 12 (19→20) not taken.
✗ Branch 13 (19→22) not taken.
|
40 | DEFINE_IGUANA_ALGORITHM(SectorFinder, clas12::SectorFinder) |
| 24 | |||
| 25 | public: | ||
| 26 | |||
| 27 | void Start(hipo::banklist& banks) override; | ||
| 28 | void Run(hipo::banklist& banks) const override; | ||
| 29 | void Stop() override; | ||
| 30 | |||
| 31 | /// @action_function{vector creator} for a given particle with index `pindex`, get its sector from a detector bank's list of `sectors` and `pindices` (both must be ordered in the same way); | ||
| 32 | /// _nb_: this is done instead of finding the `pindex` in the bank directly, to have an action function | ||
| 33 | /// | ||
| 34 | /// **Example:** | ||
| 35 | /// ```cpp | ||
| 36 | /// | ||
| 37 | /// //... Initialise algorithms & banks ... | ||
| 38 | /// | ||
| 39 | /// //For each event, do: | ||
| 40 | /// | ||
| 41 | /// std::vector<int> sectors; | ||
| 42 | /// std::vector<int> pindices | ||
| 43 | /// | ||
| 44 | /// //bank is a hipo::bank object from which we want to get the sectors | ||
| 45 | /// //for example the bank object related to REC::Calorimeter | ||
| 46 | /// for(auto const& row : bank.getRowList()) { | ||
| 47 | /// | ||
| 48 | /// int det=bank.getInt("detector",row); | ||
| 49 | /// | ||
| 50 | /// //NB: you should check you read from an FD detector | ||
| 51 | /// // eg det 7 is the ECAL | ||
| 52 | /// if(det==7){ | ||
| 53 | /// sectors.push_back(bank.getInt("sector", row)); | ||
| 54 | /// pindices.push_back(bank.getShort("pindex", row)); | ||
| 55 | /// } | ||
| 56 | /// } | ||
| 57 | /// | ||
| 58 | /// //partbank is a hipo::bank object related to REC::Particle | ||
| 59 | /// //algo_sector_finder is the iguana::clas12::SectorFinder object | ||
| 60 | /// for(auto const& row : partbank.getRowList()) { | ||
| 61 | /// int sector = algo_sector_finder.GetSector(sectors, pindices, row); | ||
| 62 | /// } | ||
| 63 | /// ``` | ||
| 64 | /// | ||
| 65 | /// @param sectors list of sectors in a detector bank | ||
| 66 | /// @param pindices list of pindices in a detector bank | ||
| 67 | /// @param pindex index in `REC::Particle` bank for which to get sector | ||
| 68 | /// @returns sector for `pindex` in list, -1 if `pindex` not in inputted list | ||
| 69 | int GetSector(std::vector<int> const& sectors, std::vector<int> const& pindices, int const pindex) const; | ||
| 70 | |||
| 71 | /// fill lists of sectors and pindices present in the input bank | ||
| 72 | /// @param bank bank from which to get lists of sectors and pindices | ||
| 73 | /// @param sectors list to fill with sectors in the bank | ||
| 74 | /// @param pindices list to fill with pindices in the bank | ||
| 75 | void GetListsSectorPindex(hipo::bank const& bank, std::vector<int>& sectors, std::vector<int>& pindices) const; | ||
| 76 | |||
| 77 | private: | ||
| 78 | |||
| 79 | /// `hipo::banklist` index for the particle bank | ||
| 80 | hipo::banklist::size_type b_particle; | ||
| 81 | hipo::banklist::size_type b_calorimeter; | ||
| 82 | hipo::banklist::size_type b_track; | ||
| 83 | hipo::banklist::size_type b_scint; | ||
| 84 | hipo::banklist::size_type b_user_charged; | ||
| 85 | hipo::banklist::size_type b_user_uncharged; | ||
| 86 | hipo::banklist::size_type b_result; | ||
| 87 | bool userSpecifiedBank_charged{false}; | ||
| 88 | bool userSpecifiedBank_uncharged{true}; | ||
| 89 | |||
| 90 | // `b_result` bank item indices | ||
| 91 | int i_sector; | ||
| 92 | int i_pindex; | ||
| 93 | |||
| 94 | /// Configuration options | ||
| 95 | std::string o_bankname_charged; | ||
| 96 | std::string o_bankname_uncharged; | ||
| 97 | |||
| 98 | //only want sectors from FD detectors | ||
| 99 | std::set<int> const listFDDets{6,7,12,15,16,18}; | ||
| 100 | }; | ||
| 101 | |||
| 102 | } | ||
| 103 |