factories/calorimetry/EnergyPositionClusterMerger_factory.h

0001 // SPDX-License-Identifier: LGPL-3.0-or-later
0002 // Copyright (C) 2023 Wouter Deconinck
0003
0004 #pragma once
0005
0006 #include <edm4eic/EDM4eicVersion.h>
0007 #include "algorithms/calorimetry/EnergyPositionClusterMerger.h"
0008 #include "extensions/jana/JOmniFactory.h"
0009 #include "services/algorithms_init/AlgorithmsInit_service.h"
0010
0011 namespace eicrecon {
0012
0013 class EnergyPositionClusterMerger_factory
0014     : public JOmniFactory<EnergyPositionClusterMerger_factory, EnergyPositionClusterMergerConfig> {
0015
0016 public:
0017   using AlgoT = eicrecon::EnergyPositionClusterMerger;
0018
0019 private:
0020   std::unique_ptr<AlgoT> m_algo;
0021
0022   PodioInput<edm4eic::Cluster> m_energy_cluster_input{this};
0023   PodioInput<edm4eic::MCRecoClusterParticleAssociation> m_energy_assoc_input{this};
0024   PodioInput<edm4eic::Cluster> m_position_cluster_input{this};
0025   PodioInput<edm4eic::MCRecoClusterParticleAssociation> m_position_assoc_input{this};
0026
0027   PodioOutput<edm4eic::Cluster> m_cluster_output{this};
0028 #if EDM4EIC_BUILD_VERSION >= EDM4EIC_VERSION(8, 7, 0)
0029   PodioOutput<edm4eic::MCRecoClusterParticleLink> m_links_output{this};
0030 #endif
0031   PodioOutput<edm4eic::MCRecoClusterParticleAssociation> m_assoc_output{this};
0032
0033   ParameterRef<double> m_energyRelTolerance{this, "energyRelTolerance",
0034                                             config().energyRelTolerance};
0035   ParameterRef<double> m_phiTolerance{this, "phiTolerance", config().phiTolerance};
0036   ParameterRef<double> m_etaTolerance{this, "etaTolerance", config().etaTolerance};
0037
0038   Service<AlgorithmsInit_service> m_algorithmsInit{this};
0039
0040 public:
0041   void Configure() {
0042     m_algo = std::make_unique<AlgoT>(GetPrefix());
0043     m_algo->level(static_cast<algorithms::LogLevel>(logger()->level()));
0044     m_algo->applyConfig(config());
0045     m_algo->init();
0046   }
0047
0048   void Process(int32_t /* run_number */, uint64_t /* event_number */) {
0049     m_algo->process({m_energy_cluster_input(), m_energy_assoc_input(), m_position_cluster_input(),
0050                      m_position_assoc_input()},
0051                     {m_cluster_output().get(),
0052 #if EDM4EIC_BUILD_VERSION >= EDM4EIC_VERSION(8, 7, 0)
0053                      m_links_output().get(),
0054 #endif
0055                      m_assoc_output().get()});
0056   }
0057 };
0058
0059 } // namespace eicrecon