factories/reco/FarForwardNeutralsReconstruction_factory.h

0001 // SPDX-License-Identifier: LGPL-3.0-or-later
0002 // Copyright (C) 2025 Sebouh Paul
0003
0004 #pragma once
0005
0006 #include "algorithms/reco/FarForwardNeutralsReconstruction.h"
0007 #include "algorithms/reco/FarForwardNeutralsReconstructionConfig.h"
0008 #include "services/algorithms_init/AlgorithmsInit_service.h"
0009 #include "extensions/jana/JOmniFactory.h"
0010
0011 namespace eicrecon {
0012
0013 class FarForwardNeutralsReconstruction_factory
0014     : public JOmniFactory<FarForwardNeutralsReconstruction_factory,
0015                           FarForwardNeutralsReconstructionConfig> {
0016
0017 public:
0018   using AlgoT = eicrecon::FarForwardNeutralsReconstruction;
0019
0020 private:
0021   std::unique_ptr<AlgoT> m_algo;
0022   PodioInput<edm4eic::Cluster> m_clusters_hcal_input{this};
0023   PodioOutput<edm4eic::ReconstructedParticle> m_neutrals_output{this};
0024
0025   ParameterRef<std::string> m_offset_position_name{this, "offsetPositionName",
0026                                                    config().offsetPositionName};
0027   ParameterRef<std::vector<double>> m_n_scale_corr_coeff_hcal{this, "neutronScaleCorrCoeffHcal",
0028                                                               config().neutronScaleCorrCoeffHcal};
0029   ParameterRef<std::vector<double>> m_gamma_scale_corr_coeff_hcal{this, "gammaScaleCorrCoeffHcal",
0030                                                                   config().gammaScaleCorrCoeffHcal};
0031   ParameterRef<double> m_global_to_proton_rotation{this, "globalToProtonRotation",
0032                                                    config().globalToProtonRotation};
0033   ParameterRef<double> m_gamma_zmax_offset{this, "gammaZMaxOffset", config().gammaZMaxOffset};
0034
0035   ParameterRef<double> m_gamma_max_length{this, "gammaMaxLength", config().gammaMaxLength};
0036   ParameterRef<double> m_gamma_max_width{this, "gammaMaxWidth", config().gammaMaxWidth};
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((algorithms::LogLevel)logger()->level());
0044
0045     m_algo->applyConfig(config());
0046     m_algo->init();
0047   }
0048
0049   void Process(int32_t /* run_number */, uint64_t /* event_number */) {
0050     m_algo->process(
0051         {
0052             m_clusters_hcal_input(),
0053         },
0054         {
0055             m_neutrals_output().get(),
0056         });
0057   }
0058 };
0059
0060 } // namespace eicrecon