#ifndef DSelector_cascade_mm_h
#define DSelector_cascade_mm_h

#include <iostream>

#include "DSelector/DSelector.h"
#include "DSelector/DHistogramActions.h"
#include "DSelector/DCutActions.h"

#include "TH1I.h"
#include "TH2I.h"

class DSelector_cascade_mm : public DSelector
{
	public:

		DSelector_cascade_mm(TTree* locTree = NULL) : DSelector(locTree){}
		virtual ~DSelector_cascade_mm(){}

		void Init(TTree *tree);
		Bool_t Process(Long64_t entry);

	private:

		void Get_ComboWrappers(void);
		void Finalize(void);

		// BEAM POLARIZATION INFORMATION
		UInt_t dPreviousRunNumber;
		bool dIsPolarizedFlag; //else is AMO
		bool dIsPARAFlag; //else is PERP or AMO
		double dBeamBunchPeriod;

		//CREATE REACTION-SPECIFIC PARTICLE ARRAYS

		//Step 0
		DParticleComboStep* dStep0Wrapper;
		DBeamParticle* dComboBeamWrapper;
		DChargedTrackHypothesis* dKPlus1Wrapper;
		DChargedTrackHypothesis* dKPlus2Wrapper;

		TH2D* dMissMassOffKPluses_BestCombo_RFSignal;
		TH2D* dMissMassOffKPluses_BestCombo_RFSideband;
		TH2D* dMissMassOffKPluses_BestCombo_RFCut;

	ClassDef(DSelector_cascade_mm, 0);
};

void DSelector_cascade_mm::Get_ComboWrappers(void)
{
	//Step 0
	dStep0Wrapper = dComboWrapper->Get_ParticleComboStep(0);
	dComboBeamWrapper = static_cast<DBeamParticle*>(dStep0Wrapper->Get_InitialParticle());
	dKPlus1Wrapper = static_cast<DChargedTrackHypothesis*>(dStep0Wrapper->Get_FinalParticle(0));
	dKPlus2Wrapper = static_cast<DChargedTrackHypothesis*>(dStep0Wrapper->Get_FinalParticle(1));
}

#endif // DSelector_cascade_mm_h