void tablex()
{
   const Int_t nx = 5;
   const Int_t ny = 3;
   char *mult_n[nx]  = {"p_evap","p_cas","total","n_evap","n_cas"};
   char *energia[ny] = {"800MeV","1200MeV","1800MeV"};
   TCanvas *c1 = new TCanvas("c1","multiplicity",10,10,500,500);
   c1->SetGrid();
   c1->SetLeftMargin(0.15);
   c1->SetBottomMargin(0.15);
   TH2F *h = new TH2F("h","Multiplicity table(n,p)-800 MeV-Pb",3,0,3,2,0,2);
   h->SetBit(TH1::kCanRebin);
   h->SetStats(0);
   Int_t vetor[500][3] ;
   Int_t i = 0 ;
   Int_t j = 0 ;
ifstream in;
in.open("tab.dat");
      while (1)
    {
     for(j=0; j<3; j++) in >> vetor[i][j] ;
     if (!in.good()) break;
     i++;
    }
         for (Int_t i=0;i<5;i++)
          {
           for (Int_t j=0;j<3;j++) 
            {
             h->Fill(energia[j],mult_n[i],vetor[4-i][j]);
            }
          }
h->LabelsDeflate("Y");
h->LabelsOption("h");
h->Draw("text");
}