#ifndef TDiPhoMat_hh
#define TDiPhoMat_hh


#include "TMatrix.h"
#include "TVector.h"
#include <iostream>

class TDiPhoMat {
public:
  TDiPhoMat();
  ~TDiPhoMat();
  void Clear();
  void Print();
  void AddNonFiducialEvent() { nonfiducial++; }
  // pass means pass in the sense of the efficiencies
  double AddEvent(double es1, double eb1, double es2, double eb2,
		  bool pass1, bool pass2);
  double AddEvent(TMatrix& eMat,bool pass1, bool pass2);
  int    GetObs(int i) { return int(m4(i));}
  double GetEbe(int i) { return ebew(i);}
  double GetEbeTotal() { return (ebew(0)+ebew(1)+ebew(2)+ebew(3));}
  double GetEbeSig(int i) { return sqrt(ebeww(i));}
  int GetFiducial() { return fiducial;}
  int GetNonFiducial() { return nonfiducial;}
  int FillMat(TMatrix& eMat, double es1, double eb1, double es2, double eb2);

private:
  // holds the sum of the observed events
  TVector m4;

  TVector ebew;  // sum of weights
  TVector ebeww; // sum of sq of weights
  int nonfiducial;
  int fiducial;
};

#endif