factories/tracking/CKFTracking_factory.h

0001 // Created by Dmitry Romanov
0002 // Subject to the terms in the LICENSE file found in the top-level directory.
0003 //
0004
0005 #pragma once
0006
0007 #include <ActsExamples/EventData/Track.hpp>
0008 #include <JANA/JEvent.h>
0009 #include <edm4eic/TrackParametersCollection.h>
0010 #include <edm4eic/TrajectoryCollection.h>
0011 #include <memory>
0012 #include <string>
0013 #include <utility>
0014 #include <vector>
0015
0016 #include "ActsExamples/EventData/Trajectories.hpp"
0017 #include "algorithms/tracking/CKFTracking.h"
0018 #include "algorithms/tracking/CKFTrackingConfig.h"
0019 #include "extensions/jana/JOmniFactory.h"
0020 #include "services/geometry/acts/ACTSGeo_service.h"
0021
0022 namespace eicrecon {
0023
0024 class CKFTracking_factory : public JOmniFactory<CKFTracking_factory, CKFTrackingConfig> {
0025
0026 private:
0027   using AlgoT = eicrecon::CKFTracking;
0028   std::unique_ptr<AlgoT> m_algo;
0029
0030   PodioInput<edm4eic::TrackParameters> m_parameters_input{this};
0031   PodioInput<edm4eic::Measurement2D> m_measurements_input{this};
0032   Output<ActsExamples::Trajectories> m_acts_trajectories_output{this};
0033   Output<ActsExamples::ConstTrackContainer> m_acts_tracks_output{this};
0034
0035   ParameterRef<std::vector<double>> m_etaBins{this, "EtaBins", config().etaBins,
0036                                               "Eta Bins for ACTS CKF tracking reco"};
0037   ParameterRef<std::vector<double>> m_chi2CutOff{this, "Chi2CutOff", config().chi2CutOff,
0038                                                  "Chi2 Cut Off for ACTS CKF tracking"};
0039   ParameterRef<std::vector<std::size_t>> m_numMeasurementsCutOff{
0040       this, "NumMeasurementsCutOff", config().numMeasurementsCutOff,
0041       "Number of measurements Cut Off for ACTS CKF tracking"};
0042
0043   Service<ACTSGeo_service> m_ACTSGeoSvc{this};
0044
0045 public:
0046   void Configure() {
0047     m_algo = std::make_unique<AlgoT>();
0048     m_algo->applyConfig(config());
0049     m_algo->init(m_ACTSGeoSvc().actsGeoProvider(), logger());
0050   }
0051
0052   void ChangeRun(int32_t /* run_number */) {}
0053
0054   void Process(int32_t /* run_number */, uint64_t /* event_number */) {
0055     std::tie(m_acts_trajectories_output(), m_acts_tracks_output()) =
0056         m_algo->process(*m_parameters_input(), *m_measurements_input());
0057   }
0058 };
0059
0060 } // namespace eicrecon