{

Double_t E1_l, E1_r, E2_l, E2_r, centr1, centr2, centr_diff, centr1_err, centr2_err, centr_diff_err ;
double fit_limit1_l, fit_limit1_r, fit_limit2_l, fit_limit2_r;
centr1=0;
centr2=0;
centr_diff=0;



E1_l =500;
E1_r =520;
E2_l =770;
E2_r =800;


fit_limit1_l = -2000;
fit_limit1_r = 2000;
fit_limit2_l = -1000;
fit_limit2_r = 1000;

TCanvas * canv1 = new TCanvas("Time_distr","",1);
TCanvas * canv2 = new TCanvas("start_stop_matrix","",1);

TFile *f=new TFile("tspec_91_cube.root");
TH3D *cube_1 = (TH3D*)ggdt_cube->Clone("cube_1");
TH3D *cube_2 = (TH3D*)ggdt_cube->Clone("cube_2");
TH3D *cube_3 = (TH3D*)ggdt_cube->Clone("cube_3");
TH3D *cube_4 = (TH3D*)ggdt_cube->Clone("cube_4");

TF1* gaus1 = new TF1("gauss1","[0]*exp(-(pow(x-[1],2))/(2*pow([2],2))) +[3]");
gaus1->SetParameter(0,500);
gaus1->SetParameter(1,0);
gaus1->SetParameter(2,300);
gaus1->SetParameter(3,10);
TF1* gaus2 = new TF1("gauss2","[0]*exp(-(pow(x-[1],2))/(2*pow([2],2))) +[3]");
gaus2->SetParameter(0,100);
gaus2->SetParameter(1,0);
gaus2->SetParameter(2,300);
gaus2->SetParameter(3,10); 




cube_1->GetXaxis()->SetRangeUser(E1_l,E1_r);
TH2D* matrix1 = (TH2D*)cube_1->Project3D("yz");
matrix1->GetYaxis()->SetRangeUser(E2_l,E2_r);
TH1D* delta1 = matrix1 -> ProjectionX();

canv1->cd();
delta1->Draw("");
gaus1->SetLineColor(kBlue);
delta1->Fit("gauss1","L M E", "R", fit_limit1_l,fit_limit1_r);
centr1=gaus1->GetParameter(1);
centr1_err=gaus1->GetParError(1);


cube_2->GetYaxis()->SetRangeUser(E1_l,E1_r);
TH2D* matrix2 = (TH2D*) cube_2->Project3D("xz");
matrix2->GetYaxis()->SetRangeUser(E2_l,E2_r);
TH1D* delta2 = matrix2 -> ProjectionX();
delta2->Draw("same");
delta2->SetLineColor(kRed);
delta2->Fit("gauss2","L M E", "R", fit_limit2_l,fit_limit2_r);
gaus2->SetLineColor(kRed);
centr2=gaus2->GetParameter(1);
centr2_err=gaus2->GetParError(1);


centr_diff=centr2-centr1;
centr_diff_err = sqrt(pow(centr1_err,2) + pow(centr2_err,2));

canv2->cd();
cube_3->GetZaxis()->SetRangeUser(-3000,3000);
TH2D* matrix3 = (TH2D*) cube_3->Project3D("xy");
matrix3->Draw("colz");

TCanvas * canv3 = new TCanvas("deltat","",1);
canv3->cd();
cube_4->GetXaxis()->SetRangeUser(0,1500);
cube_4->GetYaxis()->SetRangeUser(0,1500);
TH2D* matrix4 = (TH2D*) cube_4->Project3D("yz");
matrix4->GetYaxis()->SetRangeUser(0,1500);
TH1D* deltasum = matrix4->ProjectionX();
deltasum->Draw("");
cout<<centr_diff<<"("<<centr_diff_err<<") ps"<<endl;
}