{
   TH1F *histo = new TH1F("histo", "Spectrum",16385,0.,16385);
   TH1F *Energy_cal = new TH1F("Energy_cal", "Energy_calibrated spectrum",3000,0.,3000.);
   
   ifstream inp; double x;
  int i;
   inp.open("marneli_0.dat");
      for (i=1; i<=16385; i++){
      inp >> x;
      histo->SetBinContent(i,x);
      histo->GetXaxis()->SetTitle("Channel");
      histo->GetYaxis()->SetTitle("Count (sec)");
      
        double energy;
        energy = 0.3724941161*i - 14.5476358221;
        for(energy; energy<=5000; energy++)
        Energy_cal->SetBinContent(energy,x);
     Energy_cal->GetXaxis()->SetTitle("Energy (keV)");
     Energy_cal->GetYaxis()->SetTitle("Count (sec)");
  }
 //  TFile *f = new TFile("mca.root","recreate");
   TFile *f = new TFile("marneli.root","recreate");
 //  TFile *f = new TFile("mca.root","update"); 
   histo->Write();
   Energy_cal->Write();
   
   histo->Draw();
   Energy_cal->Draw();
   
}