{
  TFile *f = new TFile("/home/muzalevsky/AculUti/exp1803/cali1mm/calDL/SQX_L_full_calibrated_spectra.root");
  const Int_t nhists = 32;
  TH1F *h[nhists];

  TString hname,cName;
  TTree *t = (TTree*)f->Get("AnalysisxTree");
  for(Int_t i=0;i<nhists;i++) {
    hname.Form("HistSQX_L[%d]Efull",i);
    h[i] = (TH1F*)f->Get(hname.Data()); 
  }
  
  TF1* g1 = new TF1("g1", "gaus", 4.4, 4.6);
  g1->SetParLimits(0,1.,150.);
  g1->SetParLimits(1,4.,5.);
  g1->SetParLimits(2,0.1,3.);

  TF1* g2 = new TF1("g2", "gaus", 5.12, 5.35);
  g2->SetParLimits(0,1.,150.);
  g2->SetParLimits(1,5.12,5.35);
  g2->SetParLimits(2,0.1,3.);

  TF1* g3 = new TF1("g3", "gaus", 5.65, 5.95);
  g3->SetParLimits(0,1.,150.);
  g3->SetParLimits(1,5.65,5.95);
  g3->SetParLimits(2,0.1,3.);

  TF1* g4 = new TF1("g4", "gaus", 7.25, 7.65);
  g4->SetParLimits(0,1.,150.);
  g4->SetParLimits(1,7.25,7.65);
  g4->SetParLimits(2,0.1,3.);
  
  TCanvas *c[8];
  for(Int_t i=0;i<8;i++){
    cName.Form("c%d",i+1);
    TCanvas *c[i] = new TCanvas(cName.Data(),"calibrated spectra",1000,1000);
    c[i]->Divide(2,2);
  }

  /*c1->cd(1);
  h[0]->Draw();
  h[0]->Rebin(8);
  h[0]->GetXaxis()->SetRangeUser(4.,8.);
  h[0]->Fit("g1","R+");
  h[0]->Fit("g2","R+");
  h[0]->Fit("g3","R+");
  h[0]->Fit("g4","R+");*/

	ofstream outcalfile;
	outcalfile.open("/home/muzalevsky/AculUti/exp1803/cali1mm/calDL/check_cali.txt");
	if (!outcalfile.is_open()) {
		cout <<"Output file was not opened" << endl;
		return;
	}

  Int_t count=-1;
  Int_t nPad;
  for(Int_t i=0;i<32;i++) {
    count = i/4;
    nPad = (i%4)+1;
    c[count]->cd(nPad);

    h[i]->Draw();
    h[i]->Rebin(8);
    h[i]->Fit("g1","R+");
    h[i]->Fit("g2","R+");
    h[i]->Fit("g3","R+");
    h[i]->Fit("g4","R+");
    h[i]->GetXaxis()->SetRangeUser(4.,8.);
    c[count]->Update();
    outcalfile << g1->GetParameter(1) << "\t" << g2->GetParameter(1) << "\t" << g3->GetParameter(1) << "\t" << g4->GetParameter(1) << endl; 
  }

  return;
}