#ifndef _DAnalysisActionEtaU_
#define _DAnalysisActionEtaU_

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

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

    //define any histograms here
    TH1D *uboson_mass,*ep_em_gamma_mass,*ep_em_gamma_p_MM2,*ep_em_gamma_p_MM2_cut,*uboson_mass_cut,*ep_em_gamma_mass_cut,*mandelstam_t,*mandelstam_t_kin,*ep_em_gamma_mass_kin;
    TH2D *vertex_kin,*ep_em_vertex_kin;
    TH2D *electron_E_over_p,*positron_E_over_p,*ep_em_vertex_hist,*vertex_hist;
    TH1D *p_missing_mass,*ep_em_gamma_p_missing_E,*kinfit_cl,*uboson_mass_kin;

};

#endif // _DAnalysisActionEtaU_