#ifndef EVENT_H
#define EVENT_H

#include "opt_parameter.h"

class sptr{
public:
  double xdet,ydet,thdet,phdet;
  double xfp,yfp,thfp,phfp;
  double xtg,ytg,thtg,phtg;
  double beamx,beamy;
  double dp;
  double dp_c;
  //
  double zreact(void);
  double dpkin(int peak);
  double ytg_yi(int peak);
  double xsiev(void);
  double ysiev(void);
  //  double delth();
  //  double deldp();
};

class Event {
public:
  sptr sp[3];
  double pe0_loss,pe_loss,px_loss;
  double pmiss;
  double tc_cor;
  int ipeak;
  int bad;

  double emiss(int peak);  
  const char* id() const;
};

class data{
public:
  double ave_beamy,ave_beamx;
  int tot_ev;
  Event ev[NEVENT+1];
  //
  void plot_data(void);
  void read_data(void);
};

#endif