#include "TFile.h"
#include "TSpectrum.h"
#include "TH1.h"
void doit()
{
	TFile *f = new TFile("spectrum.root");
	TH1F *pSpectrum = (TH1F *)f->Get("pSpec");
	TH1F *pResponse = (TH1F *)f->Get("pResp");
	TH1F *pDeconv = (TH1F *&)pSpectrum->Clone("pDeconv");
	pDeconv->SetLineColor(2);
	TSpectrum s;
	s.DeconvolutionRL(&pDeconv->fArray[1],&pResponse->fArray[1],pSpectrum->GetNbinsX(),100,100,1.2);
	pDeconv->SetMinimum(1.);
	pDeconv->Draw();
	pSpectrum->Draw("same");
	gPad->GetCanvas()->SetLogy(1);
}