EnergyPositionClusterMerger_factory.h

Go to the documentation of this file.

// SPDX-License-Identifier: LGPL-3.0-or-later
// Copyright (C) 2023 Wouter Deconinck
 
#pragma once
 
#include <edm4eic/EDM4eicVersion.h>
#include "algorithms/calorimetry/EnergyPositionClusterMerger.h"
#include "extensions/jana/JOmniFactory.h"
#include "services/algorithms_init/AlgorithmsInit_service.h"
 
namespace eicrecon {
 
class EnergyPositionClusterMerger_factory
    : public JOmniFactory<EnergyPositionClusterMerger_factory, EnergyPositionClusterMergerConfig> {
 
public:
  using AlgoT = eicrecon::EnergyPositionClusterMerger;
 
private:
  std::unique_ptr<AlgoT> m_algo;
 
  PodioInput<edm4eic::Cluster> m_energy_cluster_input{this};
  PodioInput<edm4eic::MCRecoClusterParticleAssociation> m_energy_assoc_input{this};
  PodioInput<edm4eic::Cluster> m_position_cluster_input{this};
  PodioInput<edm4eic::MCRecoClusterParticleAssociation> m_position_assoc_input{this};
 
  PodioOutput<edm4eic::Cluster> m_cluster_output{this};
#if EDM4EIC_BUILD_VERSION >= EDM4EIC_VERSION(8, 7, 0)
  PodioOutput<edm4eic::MCRecoClusterParticleLink> m_links_output{this};
#endif
  PodioOutput<edm4eic::MCRecoClusterParticleAssociation> m_assoc_output{this};
 
  ParameterRef<double> m_energyRelTolerance{this, "energyRelTolerance",
                                            config().energyRelTolerance};
  ParameterRef<double> m_phiTolerance{this, "phiTolerance", config().phiTolerance};
  ParameterRef<double> m_etaTolerance{this, "etaTolerance", config().etaTolerance};
 
  Service<AlgorithmsInit_service> m_algorithmsInit{this};
 
public:
  void Configure() {
    m_algo = std::make_unique<AlgoT>(GetPrefix());
    m_algo->level(static_cast<algorithms::LogLevel>(logger()->level()));
    m_algo->applyConfig(config());
    m_algo->init();
  }
 
  void Process(int32_t /* run_number */, uint64_t /* event_number */) {
    m_algo->process({m_energy_cluster_input(), m_energy_assoc_input(), m_position_cluster_input(),
                     m_position_assoc_input()},
                    {m_cluster_output().get(),
#if EDM4EIC_BUILD_VERSION >= EDM4EIC_VERSION(8, 7, 0)
                     m_links_output().get(),
#endif
                     m_assoc_output().get()});
  }
};
 
} // namespace eicrecon