Algorithm.h

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

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

      bool Run(hipo::bank& particleBank, hipo::bank const& configBank) const;

      concurrent_key_t PrepareEvent(int const runnum) const;

      bool Filter(double const zvertex, int const pid, int const status, concurrent_key_t const key) const;

      int GetRunNum(concurrent_key_t const key) const;

      std::vector<double> GetElectronZcuts(concurrent_key_t const key) const;

      void SetElectronZcuts(double zcut_lower, double zcut_upper, concurrent_key_t const key);
 
    private:
      hipo::banklist::size_type b_particle, b_config;
 
      // Reload function
      void Reload(int const runnum, concurrent_key_t key) const;

      mutable std::unique_ptr<ConcurrentParam<int>> o_runnum;

      mutable std::unique_ptr<ConcurrentParam<std::vector<double>>> o_electron_vz_cuts;
  };
 
}