{
  
 ifstream in;
 in.open("152Eu_END_A_B_1.mca",ios::in);
 gROOT->ProcessLine(".L FitPeaks.C");
 
 
 int x=0;
 int channels = 4096;
 char dummy[30];
 double a,b;
 
 char s2[]="<<DATA>>";
   
   
 
    while(strcmp(dummy,s2))
      {in >>dummy;}
    
     b = (1408-121.8)/(2033 - 173.7); //b = (E2-E1)/(ch2 - ch1)
     a = 1408 - b*2033; //a = E2 - b*ch2
    
       
    
    TH1F *h1 =new TH1F("h1","Spectrum",channels,a,a+b*(channels+1));
         

    for (int i =0; i<channels;i++){
   in>>x;
 
   h1->SetBinContent(i,x);
   
   }
    
       h1->Draw();     
 
}