#ifndef _DAnalysisActionEtaDalitz_
#define _DAnalysisActionEtaDalitz_

#include "TROOT.h"
#include "TDirectoryFile.h"
#include "TH1D.h"
#include "TH2D.h"

#include "JANA/JEventLoop.h"
#include "DANA/DApplication.h"

#include "ANALYSIS/DAnalysisAction.h"
#include "ANALYSIS/DParticleCombo.h"

using namespace std;
using namespace jana;

class DAnalysisActionEtaDalitz : public DAnalysisAction
{
  public:
    DAnalysisActionEtaDalitz(const DReaction* locReaction, string locActionUniqueString = "Dalitz") :
    DAnalysisAction(locReaction, "Eta", false, locActionUniqueString){}

   
    class DalitzInfo{
    public:
      DalitzInfo(double cl,double x_dalitz,double y_dalitz,double z_dalitz,
		 double p, double theta, double E, double InvMass):cl(cl),x_dalitz(x_dalitz),
	y_dalitz(y_dalitz),z_dalitz(z_dalitz),p(p),theta(theta),E(E),InvMass(InvMass){}
      double cl;
      double x_dalitz,y_dalitz,z_dalitz;
      double p,theta,E,InvMass;
    };

  private:
    bool Perform_Action(JEventLoop* locEventLoop, const DParticleCombo* locParticleCombo);
    void Initialize(JEventLoop* locEventLoop);

  
    vector<DalitzInfo>dalitzCombos;

    unsigned int myevt;

    //define any histograms here
    TH1D *dalitzZ;
    TH2D *dalitzXY;
  

};

#endif // _DAnalysisActionEtaDalitz_