plPed (int i=4, int j=0, char *name="outP.hist.root"){
  TFile f=TFile(name);
  int ic=0;
  int nh=0;
  // int i,j,k;
  int k;
 
  //  for(i=4;i<5;i++)  // sectors
  {
  //  for(j=0;j<5;j++) // subsectors
    {
      c=new TCanvas("bb","ccc",400,500);
      c->Divide(3,4);
      for(k=0;k<12;k++) { // pseudorapidity
	char tt1[100];
	sprintf(tt1,"%2.2dT%c%2.2dc%d",i+1,'A'+j,k+1,ic);
	h= (TH1F*) f->Get(tt1); // find histrogram in file with this name
	nh++;
	c->cd(k+1);
	h->Draw();
	h->SetAxisRange(5.,30.);
	printf("try %s ...%p\n",tt1,h);
	
      }
      char text[100];
      sprintf(text,"%2.2dT%c.ps",i+1,'A'+j);
      c->Print(text);
      delete c;
      //      break;
    }
    //  break; 

  }
  printf("nh=%d\n",nh);
}

//--------------------------- Plot Ped summary -----
plPedSum(char *name="ped5"){
//eemcDb -D Calibrations_eemc -p VerD/sector07/eemcPMTped -s -g

FILE *f=fopen(name,"r");

char text[100];
TString tit;
int i=fscanf(f,"%s",text);
assert(i==1);
c=new TCanvas(name,text);

tit="Mean value " ; tit+=text;
hx=new TH1F("ped",tit,50,0,50);

tit="Sigma " ; tit+=text;
hsx=new TH1F("sig",tit,60,0,3);

int k;
for(k=0;k<60;k++) {
  float x, sx;
  i=fscanf(f,"%s%f%f",text,&x,&sx);
  // printf("%d %s%f%f",i,text,x,sx);
  assert(i==3);
  hx->Fill(x);
  hsx->Fill(sx);
}

c->Divide(1,2);
c->cd(1);
hx->Fit("gaus");
c->cd(2);
hsx->Fit("gaus");;
}