{

gROOT->Reset() ;

Int_t threshold=30 , nPix=16 ;

FILE *fin = fopen("myProg.dat","r") ; 
if ( fin==NULL ) 
  { fprintf (stderr, "!!! Error: 'myProg.dat' not readable\n") ; exit(1) ; }

TCanvas *c1 = new TCanvas ("c1","",10,10,800,800) ;
gStyle->SetOptStat(0) ;

TH2I *h_mimo5 = new TH2I ("h_mimo5","",nPix,0,nPix , nPix,0,nPix) ;
h_mimo5->SetMaximum(+50.) ; 
h_mimo5->SetMinimum(-10.) ;

Int_t myPalette[6] ;  // 1=black / 9=blue / 0=white
  myPalette[0] =  1 ; // -10 ->  0
  myPalette[1] =  9 ; //   0 -> 10 
  myPalette[2] =  9 ; //  10 -> 20
  myPalette[3] =  9 ; //  20 -> 30
  myPalette[4] =  0 ; //  30 -> 40
  myPalette[5] =  0 ; //  40 -> 50
gStyle->SetPalette(6,myPalette) ; 

Int_t raw=0,col=0 , data ;
for (Int_t i=0 ; i<nPix**2 ; i++)
{
  fscanf (fin,"%i",&data) ; 
  h_mimo5->SetBinContent(raw+1,col+1,data) ;
  col ++ ;
  if ( col==nPix ) { col=0 ; raw++ ; }
}
h_mimo5->Draw("colz") ; 
c1->SaveAs("myProg.gif") ;

fclose (fin) ;

}