#ifndef SimPipeline_h
#define SimPipeline_h

#include "Vector.h"
#include "Sample.cpp"

const Double_t SL_SYS = 0.5;	//Percentage of background events for systematic error

/*************************************************
This object contains backgrounds and various 
signals, which are added by the user.  It can
compute signal limits and determine which signals
are excluded.
*************************************************/
class SimPipeline
{
	public:
		TMatrixT<Double_t> total_background[4];
		TMatrixT<Double_t> signal_limits[4];

		Double_t luminosity;

		SimPipeline();
		~SimPipeline();

		void AddBackground(Sample* bg);
		void AddSignal(Sample* sig);
		void ComputeSignalLimits();
		void GetSignalExclusions();

	private:
		Vector<Sample*> backgrounds;	//Vector containing all your backgrounds
		Vector<Sample*> signals;		//Vector containing all your signals

		

		Double_t CalcSLStat(Double_t);
		Double_t CalcMu(Int_t nm);
};

#endif