Algorithm.h

#pragma once
 
#include "iguana/algorithms/Algorithm.h"
#include "iguana/algorithms/TypeDefs.h"
#include "iguana/services/ConcurrentParam.h"
#include "iguana/services/RCDBReader.h"
 
namespace iguana::physics {

  class InclusiveKinematics : public Algorithm
  {
 
      DEFINE_IGUANA_ALGORITHM(InclusiveKinematics, physics::InclusiveKinematics)
 
    public:
 
      void Start(hipo::banklist& banks) override;
      bool Run(hipo::banklist& banks) const override;
      void Stop() override;

      bool Run(
          hipo::bank const& particle_bank,
          hipo::bank const& config_bank,
          hipo::bank& result_bank) const;

      concurrent_key_t PrepareEvent(int const runnum, double const beam_energy = -1) const;

      InclusiveKinematicsVars ComputeFromLepton(
          vector_element_t const lepton_px,
          vector_element_t const lepton_py,
          vector_element_t const lepton_pz,
          concurrent_key_t const key) const;
 
    private:

      int FindScatteredLepton(hipo::bank const& particle_bank, concurrent_key_t const key) const;
 
      void Reload(int const runnum, double const user_beam_energy, concurrent_key_t key) const;
 
      // banklist indices
      hipo::banklist::size_type b_particle;
      hipo::banklist::size_type b_config;
      hipo::banklist::size_type b_result;
 
      // `b_result` bank item indices
      int i_pindex;
      int i_Q2;
      int i_x;
      int i_y;
      int i_W;
      int i_nu;
      int i_qx;
      int i_qy;
      int i_qz;
      int i_qE;
      int i_beamPz;
      int i_targetM;
 
      // config options
      mutable std::unique_ptr<ConcurrentParam<int>> o_runnum;
      mutable std::unique_ptr<ConcurrentParam<std::vector<double>>> o_target_PxPyPzM;
      mutable std::unique_ptr<ConcurrentParam<std::vector<double>>> o_beam_PxPyPzM;
      double o_beam_mass; // unlikely to change
      int o_beam_pdg; // unlikely to change
 
      enum method_reconstruction { scattered_lepton };
      enum method_lepton_finder { highest_energy_FD_trigger };
      method_reconstruction o_method_reconstruction;
      method_lepton_finder o_method_lepton_finder;
 
      std::unique_ptr<RCDBReader> m_rcdb;
  };
 
}