// $Id$
//
//    File: DCustomAction_jpsi_hists.h
// Created: Wed Jan 21 16:53:41 EST 2015
// Creator: jrsteven (on Linux ifarm1401 2.6.32-431.el6.x86_64 x86_64)
//

#ifndef _DCustomAction_jpsi_hists_
#define _DCustomAction_jpsi_hists_

#include <string>
#include <iostream>

#include "TH1.h"

#include "JANA/JEventLoop.h"
#include "JANA/JApplication.h"

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

using namespace std;
using namespace jana;

class DCustomAction_jpsi_hists : public DAnalysisAction
{
	public:

DCustomAction_jpsi_hists(const DReaction* locReaction, bool locUseKinFitResultsFlag, string locActionUniqueString = "", double in_EbeamMin=0., double in_EbeamMax=12., bool in_VertexReject=false) : 
	DAnalysisAction(locReaction, "Custom_jpsi_hists", locUseKinFitResultsFlag, locActionUniqueString),dEbeamMin(in_EbeamMin),dEbeamMax(in_EbeamMax),dVertexReject(in_VertexReject) {}

		void Initialize(JEventLoop* locEventLoop);

	private:

		bool Perform_Action(JEventLoop* locEventLoop, const DParticleCombo* locParticleCombo);

		// Optional: Useful utility functions.
		const DAnalysisUtilities* dAnalysisUtilities;

		Particle_t dMissingPID;
		
		double dEbeamMin, dEbeamMax;
		bool dVertexReject;

		//Store any histograms as member variables here
		TH1I *dEgamma;
        TH1I *dMeep;
        TH1I *dMee;


        TH1D *dBCAL_EP_pos;
        TH2D *dBCAL_EP_p_pos;
        TH1D *dBCAL_EP_neg;
        TH2D *dBCAL_EP_p_neg;
        TH1D *dBCAL_EP_proton;
        TH2D *dBCAL_EP_p_proton;

        TH1D *dFCAL_EP_pos;
        TH2D *dFCAL_EP_p_pos;
        TH1D *dFCAL_EP_neg;
        TH2D *dFCAL_EP_p_neg;
        TH1D *dFCAL_EP_proton;
        TH2D *dFCAL_EP_p_proton;

	//TH1D *h1d_MCWeight;

};

#endif // _DCustomAction_jpsi_hists_