| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | #pragma once | ||
| 2 | |||
| 3 | #include "iguana/algorithms/Algorithm.h" | ||
| 4 | #include "iguana/algorithms/physics/Tools.h" | ||
| 5 | #include <Math/Vector3D.h> | ||
| 6 | #include <Math/Vector4D.h> | ||
| 7 | |||
| 8 | namespace iguana::physics { | ||
| 9 | |||
| 10 | /// @brief_algo Calculate semi-inclusive dihadron kinematic quantities defined in `iguana::physics::DihadronKinematicsVars` | ||
| 11 | /// | ||
| 12 | /// @begin_doc_algo{physics::DihadronKinematics | Creator} | ||
| 13 | /// @input_banks{REC::Particle, %physics::InclusiveKinematics} | ||
| 14 | /// @output_banks{%physics::DihadronKinematics} | ||
| 15 | /// @end_doc | ||
| 16 | /// | ||
| 17 | /// @begin_doc_config{physics/DihadronKinematics} | ||
| 18 | /// @config_param{hadron_a_list | list[int] | list of "hadron A" PDGs} | ||
| 19 | /// @config_param{hadron_b_list | list[int] | list of "hadron B" PDGs} | ||
| 20 | /// @config_param{phi_r_method | string | method used to calculate @latex{\phi_R} (see section "phiR calculation methods" below)} | ||
| 21 | /// @config_param{theta_method | string | method used to calculate @latex{\theta} (see section "theta calculation methods" below)} | ||
| 22 | /// @end_doc | ||
| 23 | /// | ||
| 24 | /// Dihadron PDGs will be formed from pairs from `hadron_a_list` and `hadron_b_list`. For example, | ||
| 25 | /// if you define: | ||
| 26 | /// ```yaml | ||
| 27 | /// hadron_a_list: [ 211 ] | ||
| 28 | /// hadron_b_list: [ -211, 2212 ] | ||
| 29 | /// ``` | ||
| 30 | /// then the algorithm will calculate kinematics for @latex{\pi^+\pi^-} and @latex{\pi^+p} dihadrons; hadron A | ||
| 31 | /// is the @latex{\pi^+} for both of these, whereas hadron B is the @latex{\pi^-} for the former and the proton | ||
| 32 | /// for the latter. | ||
| 33 | /// | ||
| 34 | /// @par phiR calculation methods | ||
| 35 | /// - `"RT_via_covariant_kT"`: use @latex{R_T} computed via covariant @latex{k_T} formula | ||
| 36 | /// | ||
| 37 | /// @par theta calculation methods | ||
| 38 | /// - `"hadron_a"`: use hadron A's "decay angle" in the dihadron rest frame | ||
| 39 | /// | ||
| 40 | /// @creator_note | ||
| 41 | class DihadronKinematics : public Algorithm | ||
| 42 | { | ||
| 43 | |||
| 44 |
5/12✓ Branch 0 (2→3) taken 2002 times.
✗ Branch 1 (2→6) not taken.
✓ Branch 2 (7→8) taken 4 times.
✗ Branch 3 (7→18) not taken.
✓ Branch 4 (8→9) taken 2 times.
✗ Branch 5 (8→11) not taken.
✓ Branch 6 (11→12) taken 2 times.
✗ Branch 7 (11→22) not taken.
✓ Branch 8 (14→15) taken 2 times.
✗ Branch 9 (14→22) not taken.
✗ Branch 10 (18→19) not taken.
✗ Branch 11 (18→21) not taken.
|
2020 | DEFINE_IGUANA_ALGORITHM(DihadronKinematics, physics::DihadronKinematics) |
| 45 | |||
| 46 | public: | ||
| 47 | |||
| 48 | void Start(hipo::banklist& banks) override; | ||
| 49 | void Run(hipo::banklist& banks) const override; | ||
| 50 | void Stop() override; | ||
| 51 | |||
| 52 | /// @brief form dihadrons by pairing hadrons | ||
| 53 | /// @param particle_bank the particle bank (`REC::Particle`) | ||
| 54 | /// @returns a list of pairs of hadron rows | ||
| 55 | std::vector<std::pair<int,int>> PairHadrons(hipo::bank const& particle_bank) const; | ||
| 56 | |||
| 57 | private: | ||
| 58 | |||
| 59 | // banklist indices | ||
| 60 | hipo::banklist::size_type b_particle; | ||
| 61 | hipo::banklist::size_type b_inc_kin; | ||
| 62 | hipo::banklist::size_type b_result; | ||
| 63 | |||
| 64 | // `b_result` bank item indices | ||
| 65 | int i_pindex_a; | ||
| 66 | int i_pindex_b; | ||
| 67 | int i_pdg_a; | ||
| 68 | int i_pdg_b; | ||
| 69 | int i_Mh; | ||
| 70 | int i_z; | ||
| 71 | int i_PhPerp; | ||
| 72 | int i_MX2; | ||
| 73 | int i_xF; | ||
| 74 | int i_yB; | ||
| 75 | int i_phiH; | ||
| 76 | int i_phiR; | ||
| 77 | int i_theta; | ||
| 78 | |||
| 79 | // config options | ||
| 80 | std::set<int> o_hadron_a_pdgs; | ||
| 81 | std::set<int> o_hadron_b_pdgs; | ||
| 82 | std::string o_phi_r_method; | ||
| 83 | std::string o_theta_method; | ||
| 84 | enum {e_RT_via_covariant_kT} m_phi_r_method; | ||
| 85 | enum {e_hadron_a} m_theta_method; | ||
| 86 | |||
| 87 | // storage for a single hadron | ||
| 88 | struct Hadron { | ||
| 89 | int row; | ||
| 90 | int pdg; | ||
| 91 | ROOT::Math::PxPyPzMVector p; | ||
| 92 | double z; | ||
| 93 | std::optional<ROOT::Math::XYZVector> p_perp; | ||
| 94 | }; | ||
| 95 | |||
| 96 | }; | ||
| 97 | |||
| 98 | } | ||
| 99 |