#ifndef DStatUtilities_h
#define DStatUtilities_h

#include <vector>
#include <utility>
#include <string>
#include <sstream>
#include <iostream>

#include "TObject.h"
#include "TClass.h"
#include "TFile.h"
#include "TF1.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TObjArray.h"
#include "TDirectory.h"
#include "TKey.h"
#include "TMap.h"

using namespace std;

class DStatUtilities : public TObject
{
	//Pairs: value, uncertainty
	public:

		//CONTROL VARIABLES
		static double dMinNumTrialsForAcceptance;

		//POINT: BINOMIAL
		static double Calc_BinomialSuccessUncertainty(double locNumSuccesses, double locNumTrials);
		static pair<double, double> Calc_BinomialRate(double locNumSuccesses, double locNumTrials); //e.g. efficiency, acceptance

		//POINTS: SUM
		static pair<double, double> Calc_Sum(pair<double, double> locTerm1, pair<double, double> locTerm2, double locFactor1 = 1.0, double locFactor2 = 1.0);
		static pair<double, double> Calc_Sum(const vector<pair<double, double> >& locTerms);

		//POINTS: PRODUCT
		static pair<double, double> Calc_Product(pair<double, double> locTerm1, pair<double, double> locTerm2, double locFactor = 1.0);
		static pair<double, double> Calc_Product(pair<double, double> locTerm, double locFactor = 1.0);
		static pair<double, double> Calc_Product(const vector<pair<double, double> >& locTerms, double locFactor = 1.0);

		//POINTS: RATIO
		static pair<double, double> Calc_Ratio(pair<double, double> locNumerator, pair<double, double> locDenominator, double locFactor = 1.0);
		static pair<double, double> Calc_Ratio(pair<double, double> locDenominator, double locFactor = 1.0);
		static pair<double, double> Calc_Ratio(pair<double, double> locNumeratorTerm, const vector<pair<double, double> >& locDenominatorTerms, double locFactor = 1.0);
		static pair<double, double> Calc_Ratio(const vector<pair<double, double> >& locNumeratorTerms, pair<double, double> locDenominatorTerm, double locFactor = 1.0);
		static pair<double, double> Calc_Ratio(const vector<pair<double, double> >& locNumeratorTerms, const vector<pair<double, double> >& locDenominatorTerms, double locFactor = 1.0);

		//POINT COMPARISON: UNWEIGHTED
		static double Calc_Average(const vector<double>& locData);
		static double Calc_SampleSigma(const vector<double>& locData);
		static double Calc_SampleSigma(double locAverage, const vector<double>& locData);

		//POINT COMPARISON: WEIGHTED
		static pair<double, double> Calc_WeightedAverage(const vector<pair<double, double> >& locData);
		static double Calc_WeightedSampleSigma(const vector<pair<double, double> >& locData);
		static double Calc_WeightedSampleSigma(double locWeightedAverage, const vector<pair<double, double> >& locData);
		static pair<double, double> Calc_SystematicUncertainty(pair<double, double> locTerm1, pair<double, double> locTerm2);

		//HISTOGRAMS //Works for 1D, 2D, and 3D
		static TH1* Calc_Acceptance(TH1* locSuccessesHist, TH1* locTrialsHist, pair<string, string> locStringReplacement = pair<string, string>("", ""));
		static TH1* Calc_Acceptance_SeparateTotals(TH1* locSuccessesHist, TH1* locFailuresHist, pair<string, string> locStringReplacement = pair<string, string>("", ""));
		static TH1* Calc_Ratio(TH1* locNumeratorHist, TH1* locDenominatorHist, pair<string, string> locStringReplacement = pair<string, string>("", ""));
		static TH1* Compare(TH1* locCheckHist, TH1* locCompareToHist, pair<string, string> locStringReplacement = pair<string, string>("", ""));
		using TObject::Compare; //this is necessary because the above declaration hides the base class function

		static TPair* Calc_SystematicUncertainty(TH1* locResultHist, TH1* locAlternativeHist);

		//GROUP OF HISTOGRAMS //Works for 1D, 2D, and 3D
		static TH1* Calc_TotalSystematicUncertainty(TObjArray* locUncertaintyHists);
		static TH1* Calc_WeightedSampleSigmaFraction(TObjArray* locHistArray);
		static TH1* Calc_WeightedAverage(TObjArray* locHistArray);

		//THESE BELONG SOMEWHERE ELSE
		static TH1* CloneAndReset_Hist(TH1* locReferenceHist, string locNameSuffix, bool locCreateDoubleFlag = false, pair<string, string> locStringReplacement = pair<string, string>("", ""));
		static void Replace_TitleSubString(TH1* locHist, pair<string, string> locStringReplacement);

		//Recognized actions: "Calc_Acceptance," "Calc_Ratio," "Compare," "Calc_SystematicUncertainty"
			//In the future, try: http://stackoverflow.com/questions/4013391/is-it-possible-to-have-a-function-name-as-a-template-parameter-in-c
		static TObjArray* Perform_Action(const TObjArray* locHistArray1, const TObjArray* locHistArray2, string locAction, pair<string, string> locStringReplacement = pair<string, string>("", ""));
		static void Perform_Action(TDirectory* locDir1, TDirectory* locDir2, TDirectory* locOutputDir, string locAction, pair<string, string> locStringReplacement = pair<string, string>("", ""));

	ClassDef(DStatUtilities, 1)
};

#endif