{


TFile *f=new TFile("gate244_Eu.root");
TH1D *hist=(TH1D*)f->Get("et_1");
TF1  *fit = new TF1("fit","([0]/(2*[3]))*exp((pow([2],2))/(2*pow([3],2))-(x-[1])/[3])*TMath::Erfc((1/sqrt(2))*(([2]/[3])-((x-[1])/[2])))",-30000,30000);

fit->SetParameters(1.868E04,350,0,0);
fit->SetParLimits(0,1.487E04,1.868E04);
fit->SetParLimits(1,250,350);
 
et_1->Draw();

et_1->Fit("fit");


//e_t2_1_py__1 ->Draw("AP");
}