#ifndef DATABASE_H
#define DATABASE_H

#include "opt_parameter.h"

class database {
public:
  const char* infil;    // job input file name
  char db_in[100];      // input database file name
  char db_out[20];      // output database file name
  char dat_file[20];    // fort.51 event input file name
  class spectrometer {
  public:
    double B0[MAXPEAK];
    float Gamma[4];
    float th;
    float th_0;
    bool vary[150];
    float dsiev;
    int lines;
  }sp[3];
  //  char line[1650][190];
  int spectro;
  int hit_arm[3];
  float zoffset;
  float z_off;
  float z_step;
  float target_angle;
  float sieve_y_step;
  int fit_indx;
  float spect_h_off,spect_v_off,tg_off;
  float sieve_x_off,sieve_y_off;
  float vert_D;
  float D;
  float peak_0[MAXPEAK];
  float peak_1[MAXPEAK],start_1[MAXPEAK];
  float peak_2[MAXPEAK],start_2[MAXPEAK];
  float peak_3[MAXPEAK];
  double m_A[MAXPEAK];
  float m_B[MAXPEAK];
  float E0;
  float eloss1[MAXPEAK],eloss2[MAXPEAK],eloss3[MAXPEAK];
  int npeak;
  int nzpeak;
  int use_this[MAXPEAK];
  int pk_y[MAXPEAK],pk_z[MAXPEAK],pk_x[MAXPEAK];
  int nysievpeak,nxsievpeak;
  float xsievpeak[10],ysievpeak[10],zpeak[10];
  int init();
  int read_infile();
  int read_dp();
  int read_emiss();
  int read_database(int arm);
  int dpk_peak();
  const char* id() const;
};

#endif