#include "TH1.h"
#include "TF1.h"
#include "iostream.h"
#include "fstream.h"
#include <TMath.h>
#include <TSpectrum.h>
#include <TRandom.h>
#include <TVirtualFitter.h>

using namespace std;
using std::ofstream;

const int nbin=8192;

void main()
{
	TH1F *h = new TH1F("h","Spectrum",nbin,0,nbin);
	for (int i=1; i<=nbin; i++) {h->SetBinContent(i,0);}
	ifstream if1;
	if1.open("8020cS.txt");
	float content; int i;
	while (if1 >> i >> content) h->SetBinContent(i,content);
	if1.close();
	
	//Search peak;
	TSpectrum *s = new TSpectrum(500);
	Int_t nfound = s -> Search(h,2.5,"",0.05);
	
	Double_t *xpeaks = s->GetPositionX();
	for (Int_t i=0;i<nfound;i++)
	{cout<<xpeaks[i]<<endl;}
	
	//Estimate background;
	TH1 *hb = s->Background(h,20,"same");
	hb->SetFillColor(kYellow);
	hb->Draw("same");
	
	//Reject peaks at background level and print result;
	
}