#ifndef _DAnalysisActionEtaPrime_
#define _DAnalysisActionEtaPrime_

#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 DAnalysisActionEtaPrime : public DAnalysisAction
{
  public:
    DAnalysisActionEtaPrime(const DReaction* locReaction, string locActionUniqueString = "EtaPrime") :
    DAnalysisAction(locReaction, "EtaPrime", false, locActionUniqueString){}

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

    unsigned int myevt;

    //define any histograms here
    TH1D *h_p_missing_mass_eta_cut,*h_2gamma_dphi;
    TH1D *h_t_eta,*h_t_pi0,*h_p_eta_mass;
    TH1D *h_2gamma_measured_invariant_mass;
    TH1D *h_p_missing_mass,*h_2gamma_invariant_mass,*h_p_missing_mass_cut;
    TH2D *dalitzXY,*h_theta_vs_p,*h_p_mm_vs_E,*h_2gamma_2d;
    TH1D *h_pip_pim_mass,*h_ppippim_missing_mass_sq;
  
};

#endif // _DAnalysisActionEtaPrime_