// A reta de calibracao me fornece informacaoes a respeito da energia de uma fonte desconhecida, se eu tenho uma reta (minimo 3 pontos) eu posso medir o canal (que nao tem significado fisico nenhum eh apenas um numero). Assim sabendo o canal pelo multicanal podemos descobrir a energia da fonte.


{

float channel1 [5000], entries1[5000],channel,entries,erro_en1[5000];

FILE *co60_dates = NULL;
co60_dates = fopen("Dados_Fontes/co60/co60.txt","r");

FILE *cs137_dates = NULL;
cs137_dates = fopen("Dados_Fontes/cs137/cs137.txt","r");

//======================================================================

int i =0;
while (fscanf(co60_dates, "%f %f", &channel, &entries)!=EOF){

channel1[i]=channel;
entries1[i]=entries;
erro_en1[i]=sqrt(entries);



i++;

}//fim while


  TCanvas *c1 = new TCanvas("co60","co60",600,600);
  c1 -> cd();
  TGraphErrors *gr1 = new TGraphErrors(i,channel1,entries1,0,erro_en1);
  gr1 -> GetXaxis() -> SetTitle("Canal");
  gr1 -> GetYaxis() -> SetTitle("Entries");
  gr1 -> SetTitle("Co-60");
  gr1->Draw("A*");
  gr1 -> SetMarkerStyle(21);
  c1 -> Print("Espec_co60.pdf");



//===========================================================================

int i =0;
while (fscanf(cs137_dates, "%f %f", &channel, &entries)!=EOF){

channel1[i]=channel;
entries1[i]=entries;
erro_en1[i]=sqrt(entries);


i++;

}//fim while

  printf (" topo_co = %d \n topo_co1 = %d \n topo_cs = %d \n topo_cs1 = %d \n",topo_co60,topo_co60_1,topo_cs137,topo_cs137_1);
  TCanvas *c2 = new TCanvas("cs137","cs137",600,600);
  c2 -> cd();
  TGraphErrors *gr2 = new TGraphErrors(i,channel1,entries1,0,erro_en1);
  gr2 -> GetXaxis() -> SetTitle("Canal");
  gr2 -> GetYaxis() -> SetTitle("Entries");
  gr2 -> SetTitle("Cs-137");
  gr2->Draw("A*");
  gr2 -> SetMarkerStyle(21);
  c2 -> Print("Espec_cs137.pdf");






}//fim main