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/ParticleSvc.h"
0015
0016 namespace eicrecon {
0017
0018 using InclusiveKinematicsJBAlgorithm = algorithms::Algorithm<
0019     algorithms::Input<edm4hep::MCParticleCollection, edm4eic::ReconstructedParticleCollection,
0020                       edm4eic::HadronicFinalStateCollection>,
0021     algorithms::Output<edm4eic::InclusiveKinematicsCollection>>;
0022
0023 class InclusiveKinematicsJB : public InclusiveKinematicsJBAlgorithm {
0024
0025 public:
0026   InclusiveKinematicsJB(std::string_view name)
0027       : InclusiveKinematicsJBAlgorithm{
0028             name,
0029             {"MCParticles", "scatteredElectron", "hadronicFinalState"},
0030             {"inclusiveKinematics"},
0031             "Determine inclusive kinematics using Jacquet-Blondel method."} {}
0032
0033   void init() final;
0034   void process(const Input&, const Output&) const final;
0035
0036 private:
0037   const algorithms::ParticleSvc& m_particleSvc = algorithms::ParticleSvc::instance();
0038   double m_crossingAngle{-0.025};
0039 };
0040
0041 } // namespace eicrecon