Algorithm.h

#pragma once
 
#include "iguana/algorithms/Algorithm.h"
 
namespace iguana::clas12 {
 
  class FiducialFilter : public Algorithm
  {
 
      DEFINE_IGUANA_ALGORITHM(FiducialFilter, clas12::FiducialFilter)
 
    public:
 
      void Start(hipo::banklist& banks) override;
      void Run(hipo::banklist& banks) const override;
      void Stop() override;
 
      struct traj_row_data {
        double x1 = -999;
        double x2 = -999;
        double x3 = -999;
        double y1 = -999;
        double y2 = -999;
        double y3 = -999;
        double z1 = -999;
        double z2 = -999;
        double z3 = -999;
        int sector= -1;
      };
 
      static std::map<int, FiducialFilter::traj_row_data> GetTrajMap(hipo::bank const &bank);
 
      static int determineSectorDC(float x, float y, float z);
 
    private:
 
      bool Filter(FiducialFilter::traj_row_data const traj_row, float const torus, int const pid) const;
 
      
      bool DC_fiducial_cut_XY_pass1(FiducialFilter::traj_row_data const traj_row, float const torus, int const pid) const;
      
      
      bool DC_fiducial_cut_theta_phi_pass1(FiducialFilter::traj_row_data const traj_row, float const torus, int const pid) const;
      
      hipo::banklist::size_type b_particle;
      hipo::banklist::size_type b_traj;
      hipo::banklist::size_type b_config;
      
      int o_pass = 1;
      
  };
 
}