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;
      void Run(hipo::banklist& banks) const override;
      void Stop() override;
 
      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;
 
  };
 
}