//
// TAGMboundaries.h - class to compute the boundaries in focal plane position
//                    and energy of the tagger microscope fiber columns.
//
// author: richard.t.jones at uconn.edu
// version: june 27, 2014

#ifndef TAGMBOUNDARIES_H
#define TAGMBOUNDAIRES_H 1

#include <TFile.h>
#include <vector>

class TAGMboundaries {
 public:
   TAGMboundaries(const char *rootfile=0);
   ~TAGMboundaries();

   int Write(const char *textfile=0);
   double fXorigin_cm;
   double fYorigin_cm;

   double PhotonEnergy_GeV(double x_cm, double y_cm);
   double fEendpoint_GeV;

   int fTAGM_sums_ncolumns;
   int fTAGM_singles_ncolumns;
   int *fTAGM_singles_column;
   double fTAGM_fiber_pitch_cm;
   double fTAGM_fiber_width_cm;
   double *fTAGM_fiber_beta_deg;

   ClassDef(TAGMboundaries,0);

 private:
   struct FPtable_t {
      double FPx_cm;
      double Ephoton_GeV;
      double beta_deg;
   };
   static const char *fFPtable_struct;
   std::vector<FPtable_t> fFPtable;
   FPtable_t fFPtable_entry;

   int loadFPtable_from_textfile(const char *textfile);
   int loadFPtable_from_rootfile(const char *rootfile=0);
   int saveFPtable_to_rootfile(const char *rootfile=0);

   TFile *fRootfile;
};

#endif