#include <string>
#include <TTree.h>
#include <TMatrixD.h>
#include <TVectorD.h>
#include "run.h"
#include "stat.h"

#ifndef BIN
#define BIN

using std::string;
using namespace std;

class Bin {

public:
  Bin();
  ~Bin();

  void add(double, Run*);
  void add(Run*);
//  void out();

  int get_n();
  int get_np();
  int get_nm();

  double get_x();

  double get_Px(); double get_dPx();
  double get_Py(); double get_dPy();
  double get_Pz(); double get_dPz();
  double get_R(); double get_dR();
  double get_Ac(); double get_dAc();
  double get_Ac(double, double); double get_dAc(double, double);
  double get_Ac_calc();
  double get_px_prime_calc();
  double get_pz_prime_calc();
  double get_r_calc();

  double get_a1(); double get_da1();
  double get_b1(); double get_db1();

private:
  int no_of_events;
  int no_of_events_plus;
  int no_of_events_minus;

  TMatrixD A_Ac;
  TMatrixD b_Ac;
  TMatrixD A_noAc;
  TMatrixD b_noAc;

  Stat x;
  Stat Ac_calc;
  Stat kin_Ebeam;
  Stat kin_Ee;
  Stat kin_theta_e;
  Stat kin_tau;
  Stat kin_epsilon;
  Stat px_prime;
  Stat pz_prime;

  Double_t sum_sin;
  Double_t sum_sin2;
  Double_t sum_cos;
  Double_t sum_cos2;

  bool is_calc;
  void do_calc();

  // results

  double Px, dPx;
  double Py, dPy;
  double Pz, dPz;
  double R, dR;
  double PxAc, dPxAc;
  double PzAc, dPzAc;
  double Ac, dAc;

  double a1, da1;
  double b1, db1;

  double R_calc;
};

#endif