// $Id$
//
//    File: DCustomAction_sdobbs_pi0_dalitz.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_sdobbs_pi0_dalitz_
#define _DCustomAction_sdobbs_pi0_dalitz_

#include <string>
#include <iostream>

#include "TH1.h"
#include "TMath.h"
#include "Math/Point3D.h"
#include "Math/Vector3D.h"
#include "Math/Plane3D.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"
#include "PID/DDetectorMatches.h"
#include "BCAL/DBCALShower.h"
#include "FCAL/DFCALShower.h"
#include "FCAL/DFCALCluster.h"

using namespace std;
using namespace jana;
using namespace ROOT::Math;

class DCustomAction_sdobbs_pi0_dalitz : public DAnalysisAction
{
	public:

		DCustomAction_sdobbs_pi0_dalitz(const DReaction* locReaction, bool locUseKinFitResultsFlag, string locActionUniqueString = "") : 
		DAnalysisAction(locReaction, "Custom_sdobbs_pi0_dalitz", locUseKinFitResultsFlag, locActionUniqueString) {}

		void Initialize(JEventLoop* locEventLoop);

	private:

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

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

		//Store any histograms as member variables here
		TH1I *dEgamma, *dmass_ee, *dmassPi0_ee_cut, *ebyp, *ebyp_BCAL_neg, *ebyp_BCAL_plus, *ebyp_FCAL_neg, *ebyp_FCAL_plus;
		TH2I *dMM2_MPi0, *dProton_dEdx_P, *dProton_P_Theta;
		TH2I *dDeltaE_MPi0;
		TH2I *dMM2_M3pi_ProtonTag, *dDeltaE_M3pi_ProtonTag, *dMM2_DeltaE_ProtonTag;
                TH2I *dMM2_M3pi_CoplanarTag, *dDeltaE_M3pi_CoplanarTag, *dMM2_DeltaE_CoplanarTag;
                TH2I *dDeltaPhi_MPi0, *dPhiPi0_PhiP;
		TH2I *dEgamma_MPi0_ProtonTag;
		TH2I *dMppizero_M2pi, *dMppiplus_M2pi, *dMppiminus_M2pi;
		TH2I *dmass_ee_v_masseeg, *ebyplogic_mee_v_meeg, *ebyplogic_both_mee_v_meeg, *meeg_v_pmag, *meeg_v_shwrE;
		TH2I *BCAL_Ee_minus_pe_vs_pe, *BCAL_Ee_plus_pe_vs_pe, *FCAL_Ee_minus_pe_vs_pe, *FCAL_Ee_plus_pe_vs_pe;

};

#endif // _DCustomAction_sdobbs_pi0_dalitz_