algorithms/reco/InclusiveKinematicsJB.h

0001 // SPDX-License-Identifier: LGPL-3.0-or-later
0002 // Copyright (C) 2022, 2023 Sylvester Joosten, Dmitry Romanov, Wouter Deconinck
0003
0004 #pragma once
0005
0006 #include <algorithms/algorithm.h>
0007 #include <edm4eic/HadronicFinalStateCollection.h>
0008 #include <edm4eic/InclusiveKinematicsCollection.h>
0009 #include <edm4eic/ReconstructedParticleCollection.h>
0010 #include <edm4hep/MCParticleCollection.h>
0011 #include <string>
0012 #include <string_view>
0013
0014 #include "algorithms/interfaces/WithPodConfig.h"
0015 #include "services/particle/ParticleSvc.h"
0016
0017 namespace eicrecon {
0018
0019 using InclusiveKinematicsJBAlgorithm = algorithms::Algorithm<
0020     algorithms::Input<edm4hep::MCParticleCollection, edm4eic::ReconstructedParticleCollection,
0021                       edm4eic::HadronicFinalStateCollection>,
0022     algorithms::Output<edm4eic::InclusiveKinematicsCollection>>;
0023
0024 class InclusiveKinematicsJB : public InclusiveKinematicsJBAlgorithm,
0025                               public WithPodConfig<NoConfig> {
0026
0027 public:
0028   InclusiveKinematicsJB(std::string_view name)
0029       : InclusiveKinematicsJBAlgorithm{
0030             name,
0031             {"MCParticles", "scatteredElectron", "hadronicFinalState"},
0032             {"inclusiveKinematics"},
0033             "Determine inclusive kinematics using Jacquet-Blondel method."} {}
0034
0035   void init() final;
0036   void process(const Input&, const Output&) const final;
0037
0038 private:
0039   const algorithms::ParticleSvc& m_particleSvc = algorithms::ParticleSvc::instance();
0040   double m_crossingAngle{-0.025};
0041 };
0042
0043 } // namespace eicrecon