Algorithm.h

#pragma once
 
#include "iguana/algorithms/Algorithm.h"
 
namespace iguana::physics {

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

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

      DepolarizationVars Compute(double const Q2, double const x, double const y, double const targetM) const;
 
    private:
 
      // banklist indices
      hipo::banklist::size_type b_inc_kin;
      hipo::banklist::size_type b_result;
 
      // `b_result` bank item indices
      int i_epsilon;
      int i_A;
      int i_B;
      int i_C;
      int i_V;
      int i_W;
  };
 
}