#include <TObjArray.h>
#include <TFractionFitter.h>
#include <TFile.h>
#include <TH1D.h>

void tffit ()
{

  TFile *_file0=new TFile("testmc.root","read");

  TObjArray *mc = new TObjArray(4);
  mc->Add(_file0->Get("hee_energytot_ext-exbg11-bi214")->Clone());
  mc->Add(_file0->Get("hee_energytot_ext-exbg11-tl208")->Clone());
  mc->Add(_file0->Get("hee_energytot_ext-exbg4-bi214")->Clone());
  mc->Add(_file0->Get("fitted_tot")->Clone());
  TH1D* h=(TH1D*)(_file0->Get("hee_energytot_ext-betabeta")->Clone());
  h->Smooth(100);
  TFractionFitter* tf=new TFractionFitter (h,mc);
  tf->SetRangeX(10,80);
  tf->GetFitter()->SetParameter( 0, "0", 0.7, 0.001, 0.001, 1.5);
  tf->GetFitter()->SetParameter( 1, "1", 0.7, 0.001, 0.001, 1.5);
  tf->GetFitter()->SetParameter( 2, "2", 0.7, 0.001, 0.001, 1.5);
  tf->GetFitter()->SetParameter( 3, "0", 0.7, 0.01, 0.8, 1.1);

  tf->Fit();
}