start writing AlgorithmProcess

src/algorithms/pid/IrtParticleID.cc

@@ -85,6 +85,21 @@ std::vector<edm4hep::ParticleID*> eicrecon::IrtParticleID::AlgorithmProcess(
   std::vector<edm4hep::ParticleID*> particle_id;
   if(m_init_failed) return particle_id;
 
+  // loop over charged particles ********************************************
+  for(const auto& charged_particle : charged_particles) {
+    int charged_particle_pdg = 0; // FIXME: not known, maybe not needed
+    auto irt_particle = std::make_unique<ChargedParticle>(charged_particle_pdg);
+
+    // start radiator history
+    for(auto [name,rad] : m_irt_det->Radiators())
+      irt_particle->StartRadiatorHistory({ rad, new RadiatorHistory() });
+
+    // loop over `charged_particle`'s projected track points
+
+    // loop over raw hits ***************************************************
+    // for(const auto& raw_hit : raw_hits) {
+  }
+
   m_log->trace("{:=^70}"," call IrtParticleID::AlgorithmProcess ");
   m_log->trace("number of raw hits: {}", raw_hits.size());
   m_log->trace("number of propagated tracks: {}", charged_particles.size());