void cofit()
{
#include <fstream>


   TH1F *hist1 = new TH1F("hist1","title1", 100, 0., 2000);
   TH1F *hist2 = new TH1F("hist2","title2", 100, 0., 2000);

   std::ifstream infile("Co60.dat");
   int x1,x2;
   while (infile >> x1 >> x2) {
      hist1->Fill(x1);
      hist2->Fill(x2);
   }
   
   
   
   /*
   double p0,p1;

 TF1 *f3 = new TF1("f3","[0]*exp(-0.5*pow((x-[1])/[2],2))+[3]*exp(-0.5*pow((x-[4])/[5],2))+[6]*exp(-[7]*x)",1000,1600);
 
   f3->SetParLimits(1,1220,1260);
   f3->SetParLimits(4,1380,1430);
   f3->SetParLimits(2,20,60);
   f3->SetParLimits(5,20,60);
 
   

  hist1->Draw();

  hist1->Fit("f3","","",1000,1600);
  f3->Draw("Same"); 
  
 
  p0=f3->GetParameter(6);
  p1=f3->GetParameter(7);
  
  TF1 *fbg = new TF1("fbg","[0]*exp(-[1]*x)",1000,1600);
  

  fbg->SetParameter(0,f3->GetParameter(6));
  fbg->SetParameter(1,f3->GetParameter(7));
  
  fbg->SetLineColor(5);
  fbg->Draw("Same");
  */
  
  
  double p2,p3;
  TF1 *f4 = new TF1("f4","[8]*exp(-0.5*pow((x-[9])/[10],2))+[11]*exp(-0.5*pow((x-[12])/[13],2))+[14]*exp(-[15]*x)",1000,1600);
 // TF1 *f4 = new TF1("f4","[8]*exp(-0.5*pow((x-[9])/[10],2))+[11]*exp(-0.5*pow((x-[12])/[13],2))+[14]+[15]*x+[16]*x*x",1000,1600);
   f4->SetParLimits(9,1060,1180);
   f4->SetParLimits(11,1205,1430);
   f4->SetParLimits(10,10,80);
   f4->SetParLimits(13,10,80);
   hist2->Draw();
   
   hist2->Fit("f4","","",1000,1600);
  f4->Draw("Same"); 
  p2=f4->GetParameter(14);
  p3=f4->GetParameter(15);
//  p4=f4->GetParameter(16);
 // p5=f4->GetParameter(17);
  
 // TF1 *fbg2 = new TF1("fbg2","[0]+[1]*x+[2]*x*x+[3]*x*x*x",1000,1600);
//  fbg2->SetParameter(0,f4->GetParameter(14));
//  fbg2->SetParameter(1,f4->GetParameter(15));
  //fbg2->SetParameter(2,f4->GetParameter(16));
 // fbg2->SetParameter(3,f4->GetParameter(17));
 // fbg2->SetLineColor(6);
 // fbg2->Draw("Same");
  
  
 // auto c = new TCanvas("canvas", "title", 800, 600);
//   c->Divide(1,2);
 //  c->cd(1);
  
 //  hist1->Draw();
 //  c->cd(2);
 //  hist2->Draw();
   
   
  
  
   
  
   
   }