#ifndef SEPARATIONPOWER_HH
#define SEPARATIONPOWER_HH

class TGraph;

//Separation Power for a specific Particle-Pair
class SeparationPower
{
public:
	SeparationPower(string Name1,string Name2, Color_t color=kRed) : ParticleName1(Name1), ParticleName2(Name2)	{}
	~SeparationPower() {}
	void FillSeparationPower(map<Double_t, map<string, TF1*> > &FitMap);
	void FillSeparationPower(TF1 *FitFunction1, TF1 *FitFunction2, const Double_t X);
	void FillSeparationPower(const Double_t Mean1, const Double_t Mean2,  const Double_t Sigma1, const Double_t Sigma2, const Double_t X);
	void FillSeparationPower(const Double_t SepPo,const Double_t X);
	inline bool IsFor(string ParticleName)	{ return( ParticleName==ParticleName1 || ParticleName==ParticleName2 ); }
	inline bool IsFor(string ParticleName1, string ParticleName2)	{ return( IsFor(ParticleName1) && IsFor(ParticleName2) ); }
	string GetTheOtherParticle(string ParticleName);
	void FillHistogram(TH2F  *pHisto);
	TGraph *CreateGraph();
	inline bool CheckSanity()	{	return(vSepPo.size() == vP.size()); }
	inline int Size()	{	return(vSepPo.size());	}	
private:
	string ParticleName1;
	string ParticleName2;
	vector<Double_t> vSepPo;
	vector<Double_t> vP;
};

#endif