Algorithm.h

#pragma once
 
#include "iguana/algorithms/Algorithm.h"
 
namespace iguana::clas12 {
 
  class SectorFinder : public Algorithm
  {
 
      DEFINE_IGUANA_ALGORITHM(SectorFinder, clas12::SectorFinder)
 
    public:
 
      void Start(hipo::banklist& banks) override;
      void Run(hipo::banklist& banks) const override;
      void Stop() override;
 
      int GetSector(std::vector<int> const& sectors, std::vector<int> const& pindices, int const pindex) const;
 
      void GetListsSectorPindex(hipo::bank const& bank, std::vector<int>& sectors, std::vector<int>& pindices) const;
 
    private:
 
      hipo::banklist::size_type b_particle;
      hipo::banklist::size_type b_calorimeter;
      hipo::banklist::size_type b_track;
      hipo::banklist::size_type b_scint;
      hipo::banklist::size_type b_user_charged;
      hipo::banklist::size_type b_user_uncharged;
      hipo::banklist::size_type b_result;
      bool userSpecifiedBank_charged{false};
      bool userSpecifiedBank_uncharged{true};
 
      // `b_result` bank item indices
      int i_sector;
      int i_pindex;
 
      std::string o_bankname_charged;
      std::string o_bankname_uncharged;
 
      //only want sectors from FD detectors
      std::set<int> const listFDDets{6,7,12,15,16,18};
  };
 
}