-
Notifications
You must be signed in to change notification settings - Fork 29
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge branch 'main' into nopayloadclient
- Loading branch information
Showing
51 changed files
with
895 additions
and
355 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,14 @@ | ||
| set(PLUGIN_NAME "algorithms_pid_lut") | ||
|
|
||
| # Function creates ${PLUGIN_NAME}_plugin and ${PLUGIN_NAME}_library targets | ||
| # Setting default includes, libraries and installation paths | ||
| plugin_add(${PLUGIN_NAME} WITH_SHARED_LIBRARY WITHOUT_PLUGIN) | ||
|
|
||
| # The macro grabs sources as *.cc *.cpp *.c and headers as *.h *.hh *.hpp Then | ||
| # correctly sets sources for ${_name}_plugin and ${_name}_library targets Adds | ||
| # headers to the correct installation directory | ||
| plugin_glob_all(${PLUGIN_NAME}) | ||
|
|
||
| # Find dependencies | ||
| plugin_add_algorithms(${PLUGIN_NAME}) | ||
| plugin_add_event_model(${PLUGIN_NAME}) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,142 @@ | ||
| // SPDX-License-Identifier: LGPL-3.0-or-later | ||
| // Copyright (C) 2024, Nathan Brei, Dmitry Kalinkin | ||
|
|
||
| #include <algorithms/service.h> | ||
| #include <edm4eic/MCRecoParticleAssociationCollection.h> | ||
| #include <edm4eic/ReconstructedParticleCollection.h> | ||
| #include <edm4hep/MCParticleCollection.h> | ||
| #include <edm4hep/utils/vector_utils.h> | ||
| #include <fmt/core.h> | ||
| #include <podio/ObjectID.h> | ||
| #include <cmath> | ||
| #include <gsl/pointers> | ||
| #include <stdexcept> | ||
|
|
||
| #include "algorithms/pid_lut/PIDLookup.h" | ||
| #include "algorithms/pid_lut/PIDLookupConfig.h" | ||
| #include "services/pid_lut/PIDLookupTableSvc.h" | ||
|
|
||
| namespace eicrecon { | ||
|
|
||
| void PIDLookup::init() { | ||
| auto& serviceSvc = algorithms::ServiceSvc::instance(); | ||
| auto lut_svc = serviceSvc.service<PIDLookupTableSvc>("PIDLookupTableSvc"); | ||
|
|
||
| m_lut = lut_svc->load(m_cfg.filename, { | ||
| .pdg_values=m_cfg.pdg_values, | ||
| .charge_values=m_cfg.charge_values, | ||
| .momentum_edges=m_cfg.momentum_edges, | ||
| .polar_edges=m_cfg.polar_edges, | ||
| .azimuthal_binning=m_cfg.azimuthal_binning, | ||
| .azimuthal_bin_centers_in_lut=m_cfg.azimuthal_bin_centers_in_lut, | ||
| .momentum_bin_centers_in_lut=m_cfg.momentum_bin_centers_in_lut, | ||
| .polar_bin_centers_in_lut=m_cfg.polar_bin_centers_in_lut, | ||
| .skip_legacy_header=m_cfg.skip_legacy_header, | ||
| .use_radians=m_cfg.use_radians, | ||
| .missing_electron_prob=m_cfg.missing_electron_prob, | ||
| }); | ||
| if (m_lut == nullptr) { | ||
| throw std::runtime_error("LUT not available"); | ||
| } | ||
| } | ||
|
|
||
| void PIDLookup::process(const Input& input, const Output& output) const { | ||
| const auto [recoparts_in, partassocs_in] = input; | ||
| auto [recoparts_out, partassocs_out, partids_out] = output; | ||
|
|
||
| for (const auto& recopart_without_pid : *recoparts_in) { | ||
| edm4hep::MCParticle mcpart; | ||
| auto recopart = recopart_without_pid.clone(); | ||
|
|
||
| // Find MCParticle from associations and propagate the relevant ones further | ||
| bool assoc_found = false; | ||
| for (auto assoc_in : *partassocs_in) { | ||
| if (assoc_in.getRec() == recopart_without_pid) { | ||
| if (assoc_found) { | ||
| warning("Found a duplicate association for ReconstructedParticle at index {}", recopart_without_pid.getObjectID().index); | ||
| warning("The previous MCParticle was at {} and the duplicate is at {}", mcpart.getObjectID().index, assoc_in.getSim().getObjectID().index); | ||
| } | ||
| assoc_found = true; | ||
| mcpart = assoc_in.getSim(); | ||
| auto assoc_out = assoc_in.clone(); | ||
| assoc_out.setRec(recopart); | ||
| partassocs_out->push_back(assoc_out); | ||
| } | ||
| } | ||
| if (not assoc_found) { | ||
| recoparts_out->push_back(recopart); | ||
| continue; | ||
| } | ||
|
|
||
| int true_pdg = mcpart.getPDG(); | ||
| int true_charge = mcpart.getCharge(); | ||
| int charge = recopart.getCharge(); | ||
| double momentum = edm4hep::utils::magnitude(recopart.getMomentum()); | ||
|
|
||
| double theta = edm4hep::utils::anglePolar(recopart.getMomentum()) / M_PI * 180.; | ||
| double phi = edm4hep::utils::angleAzimuthal(recopart.getMomentum()) / M_PI * 180.; | ||
|
|
||
| trace("lookup for true_pdg={}, true_charge={}, momentum={:.2f} GeV, polar={:.2f}, aziumthal={:.2f}", | ||
| true_pdg, true_charge, momentum, theta, phi); | ||
| auto entry = m_lut->Lookup(true_pdg, true_charge, momentum, theta, phi); | ||
|
|
||
| int identified_pdg = 0; // unknown | ||
|
|
||
| if ((entry != nullptr) && ((entry->prob_electron != 0.) || (entry->prob_pion != 0.) || (entry->prob_kaon != 0.) || (entry->prob_electron != 0.))) { | ||
| double random_unit_interval = m_dist(m_gen); | ||
|
|
||
| trace("entry with e:pi:K:P={}:{}:{}:{}", entry->prob_electron, entry->prob_pion, entry->prob_kaon, entry->prob_proton); | ||
|
|
||
| recopart.addToParticleIDs(partids_out->create( | ||
| m_cfg.system, // std::int32_t type | ||
| std::copysign(11, charge), // std::int32_t PDG | ||
| 0, // std::int32_t algorithmType | ||
| static_cast<float>(entry->prob_electron) // float likelihood | ||
| )); | ||
| recopart.addToParticleIDs(partids_out->create( | ||
| m_cfg.system, // std::int32_t type | ||
| std::copysign(211, charge), // std::int32_t PDG | ||
| 0, // std::int32_t algorithmType | ||
| static_cast<float>(entry->prob_pion) // float likelihood | ||
| )); | ||
| recopart.addToParticleIDs(partids_out->create( | ||
| m_cfg.system, // std::int32_t type | ||
| std::copysign(321, charge), // std::int32_t PDG | ||
| 0, // std::int32_t algorithmType | ||
| static_cast<float>(entry->prob_kaon) // float likelihood | ||
| )); | ||
| recopart.addToParticleIDs(partids_out->create( | ||
| m_cfg.system, // std::int32_t type | ||
| std::copysign(2212, charge), // std::int32_t PDG | ||
| 0, // std::int32_t algorithmType | ||
| static_cast<float>(entry->prob_proton) // float likelihood | ||
| )); | ||
|
|
||
| if (random_unit_interval < entry->prob_electron) { | ||
| identified_pdg = 11; // electron | ||
| recopart.setParticleIDUsed((*partids_out)[partids_out->size() - 4]); | ||
| } else if (random_unit_interval < (entry->prob_electron + entry->prob_pion)) { | ||
| identified_pdg = 211; // pion | ||
| recopart.setParticleIDUsed((*partids_out)[partids_out->size() - 3]); | ||
| } else if (random_unit_interval < | ||
| (entry->prob_electron + entry->prob_pion + entry->prob_kaon)) { | ||
| identified_pdg = 321; // kaon | ||
| recopart.setParticleIDUsed((*partids_out)[partids_out->size() - 2]); | ||
| } else if (random_unit_interval < (entry->prob_electron + entry->prob_pion + | ||
| entry->prob_kaon + entry->prob_electron)) { | ||
| identified_pdg = 2212; // proton | ||
| recopart.setParticleIDUsed((*partids_out)[partids_out->size() - 1]); | ||
| } | ||
| } | ||
|
|
||
| recopart.setPDG(std::copysign(identified_pdg, charge)); | ||
|
|
||
| if (identified_pdg != 0) { | ||
| trace("randomized PDG is {}", recopart.getPDG()); | ||
| } | ||
|
|
||
| recoparts_out->push_back(recopart); | ||
| } | ||
| } | ||
|
|
||
| } // namespace eicrecon |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,44 @@ | ||
| // SPDX-License-Identifier: LGPL-3.0-or-later | ||
| // Copyright (C) 2024, Nathan Brei, Dmitry Kalinkin | ||
|
|
||
| #include <algorithms/algorithm.h> | ||
| #include <edm4eic/MCRecoParticleAssociationCollection.h> | ||
| #include <edm4eic/ReconstructedParticleCollection.h> | ||
| #include <edm4hep/ParticleIDCollection.h> | ||
| #include <random> | ||
| #include <string> | ||
| #include <string_view> | ||
|
|
||
| #include "PIDLookupConfig.h" | ||
| #include "algorithms/interfaces/WithPodConfig.h" | ||
| #include "services/pid_lut/PIDLookupTable.h" | ||
|
|
||
| namespace eicrecon { | ||
|
|
||
| using PIDLookupAlgorithm = | ||
| algorithms::Algorithm<algorithms::Input<edm4eic::ReconstructedParticleCollection, | ||
| edm4eic::MCRecoParticleAssociationCollection>, | ||
| algorithms::Output<edm4eic::ReconstructedParticleCollection, | ||
| edm4eic::MCRecoParticleAssociationCollection, | ||
| edm4hep::ParticleIDCollection>>; | ||
|
|
||
| class PIDLookup : public PIDLookupAlgorithm, public WithPodConfig<PIDLookupConfig> { | ||
|
|
||
| public: | ||
| PIDLookup(std::string_view name) | ||
| : PIDLookupAlgorithm{name, | ||
| {"inputParticlesCollection", "inputParticleAssociationsCollection"}, | ||
| {"outputParticlesCollection", "outputParticleAssociationsCollection", | ||
| "outputParticleIDCollection"}, | ||
| ""} {} | ||
|
|
||
| void init() final; | ||
| void process(const Input&, const Output&) const final; | ||
|
|
||
| private: | ||
| mutable std::mt19937 m_gen{}; | ||
| mutable std::uniform_real_distribution<double> m_dist{0, 1}; | ||
| const PIDLookupTable* m_lut; | ||
| }; | ||
|
|
||
| } // namespace eicrecon |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,26 @@ | ||
| // SPDX-License-Identifier: LGPL-3.0-or-later | ||
| // Copyright (C) 2024, Nathan Brei, Dmitry Kalinkin | ||
|
|
||
| #pragma once | ||
|
|
||
| #include <vector> | ||
|
|
||
| namespace eicrecon { | ||
|
|
||
| struct PIDLookupConfig { | ||
| std::string filename; | ||
| int system; | ||
| std::vector<int> pdg_values; | ||
| std::vector<int> charge_values; | ||
| std::vector<double> momentum_edges; | ||
| std::vector<double> polar_edges; | ||
| std::vector<double> azimuthal_binning; | ||
| bool azimuthal_bin_centers_in_lut {false}; | ||
| bool momentum_bin_centers_in_lut {false}; | ||
| bool polar_bin_centers_in_lut {false}; | ||
| bool skip_legacy_header {false}; | ||
| bool use_radians {false}; | ||
| bool missing_electron_prob {false}; | ||
| }; | ||
|
|
||
| } // namespace eicrecon |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Oops, something went wrong.