TH1D           *ehi[NGE];
  TH1D           *eside[NGE];
  TH1D           *tge[NGE];
  TH1D           *ebgo[NGE];
  TH1D           *ebgo_seg[NGE][8];
  TH1D           *tbgo[NGE];
  TH1D           *elo[NGE];
  TH1D           *tac1;
  TH1D           *tac2;
  TH1D           *sge;
  TH1D           *sbgo;
  TH1D           *ring[NANG + 1];

  TH1D           *sumehi;
  TH1D           *sumelo;
  TH1D           *sumeBGO;
  TH1D           *sumeside;
  TH1D           *sumtge;
  TH1D           *sumtbgo;
  TH1D           *gehitpat;
  TH1D           *sidehitpat;
  TH1D           *puhitpat;
  TH1D           *orhitpat;
  TH1D           *bgohitpat;
  TH1D           *len;
  TH1D           *len_clean;
  TH1D           *len_dirty;
  TH1D           *len_bgo;
  TH1D           *dtbtev;
  TH1D           *dtbtgm;
  
  /* GS data structures etc */

  struct EVHDR    hdr;
  struct EVENT    ev[NGE];
  unsigned int    ext[5000] = {0};
  unsigned int    nEXT, fera_type[EXTLEN], fera_vsn[EXTLEN], fera_ch[EXTLEN];
  unsigned int    fera_data[EXTLEN], nbadfera = 0;
  unsigned int    tgerf[NGE], tbgorf[NGE];
  int             angno[NGE], fera_st;
  float           ring_angle[NANG + 1];
  int             GSSortError[NERR];
  double          hdr_ttL_last = 0;
  double          dclock1, dclock2;