---

mMuiFastRoadFinder.h

Go to the documentation of this file.

#ifndef __MMUIFASTROADFINDER_H__
#define __MMUIFASTROADFINDER_H__

#include "PHModule.h"
#include "PHTimeServer.h"
#include <MUIOO.h>

#include <mMuiFastRoadFinderPar.h>
#include <TMuiRoadMapO.h>
#include <TMuiHitMapO.h>

#include <vector>

struct mui_symset
{
  UShort_t idx;    // Symset index corresponding to gap0 twopack  
  Float_t cos;     // Slope of 1D projection
  UShort_t depth;  // Tracked to gap 1-4
  UShort_t nhit;   // How many hits?
  UShort_t dof;    // Degrees of freedom (nhit - 2)
  UShort_t gapbit;      // which gaps were hit (bit coded)
  // start coordinates
  Float_t x;
  Float_t y;
  Float_t z;
  Float_t chisq;
  
  mui_symset() : idx(0), cos(0.), depth(0), nhit(0), dof(1), gapbit(0), 
    x(0.0), y(0.0), z(0.0), chisq(1.0) {}
};
typedef std::vector<mui_symset> symset_list;
typedef symset_list::iterator symset_iterator;

class mMuiFastRoadFinder : public PHModule
{
  enum{kTwoPacksPerPlaneMax=326};
 public:
  mMuiFastRoadFinder();
  ~mMuiFastRoadFinder();
  
  PHBoolean event(PHCompositeNode* top_node);
  
 private:
  
  void set_interface_ptrs(PHCompositeNode* top_node);
  
  PHBoolean init();
  
  void read_idx_XoZ_map(const char* filename = "L2MuiIdxXoZ.adb");
  
  void read_twopack_idx_map(const char* filename = "muioomap.adb");
  
  void find_tracks();
  
  void hit_to_logical(UShort_t arm, UShort_t orient);
  
  void calc_symset(UShort_t arm, UShort_t orient);

  void hit_to_logical_per_panel(UShort_t arm, UShort_t orient,
                                UShort_t virt_panel);
  
  void mui_trigger_edge_finder(int nsymsets, 
                               const float idx_XoZ[kTwoPacksPerPlaneMax],
                               bool *gap0bit,
                               bool *gap1bit, bool *gap2bit, 
                               bool *gap3bit, bool *gap4bit,
                               UShort_t arm, UShort_t orient);
  
  void simple_trigger_road_finder(bool **p,
                                  bool& gap1bit, bool& gap2bit, 
                                  bool& gap3bit, bool& gap4bit);
  void orig_trigger_road_finder(bool **p,
                                  bool& gap1bit, bool& gap2bit, 
                                  bool& gap3bit, bool& gap4bit);
  
  /* Methods for orig road finding (copied from LVL2) */
  void gap0(bool **p, short* showermax);
  
  void gap1(bool **p, short swi, short* showermax, bool* roadfound);
                   
  void gap2(bool **p, short swi, short* showermax, bool* roadfound);
                   
  void gap3(bool **p, short swi, short* showermax, bool* roadfound);

  void gap4(bool **p, short swi, short* showermax, bool* roadfound);

  short bsum (bool *st, long n);

  bool _signal[kTwoPacksPerPlaneMax][MUIOO::MAX_PLANE];
  bool* _trigger_signal[kTwoPacksPerPlaneMax][mMuiFastRoadFinderPar::MAX_SYMSET];
  
  Float_t _idx_XoZ[MUIOO::MAX_ARM][MUIOO::MAX_ORIENTATION][kTwoPacksPerPlaneMax];
  UShort_t _twopack_idx[MUIOO::MAX_ARM][MUIOO::MAX_PLANE][MUIOO::MAX_PANEL][MUIOO::MAX_ORIENTATION][kTwoPacksPerPlaneMax];
  
  symset_list _symsets[MUIOO::MAX_ARM][MUIOO::MAX_ORIENTATION];
  
  const mMuiFastRoadFinderPar* _mod_par;   
  TMuiRoadMapO* _roadmap;
  TMuiHitMapO* _hitmap;
  PHTimeServer::timer _timer;  
};


#endif /*__MMUIFASTROADFINDER_H__*/

---