void RecoH2TB::initFADCtable()
{
  //   Updated,  2002.9.  S.Kunori
  float rangefactor[4]={1.,5.,25.,125.};
  float binwidth[5]={1.0,2.0,3.0,4.0,5.0};
  int numbins[5]={15,7,4,3,3};
  float charge=-1.0;
  int m=0;
  int nbins;

  // loop over 4 ranges
  for(int i=0; i<4;i++){
    // loop over 5 slopes;
    for(int j=0; j<5; j++) {
      nbins=numbins[j];
      // loop over bins
      for(int k=0; k<nbins; k++) {
        charge=charge+binwidth[j]*rangefactor[i];
        fadc2charge[m]=charge-binwidth[j]*rangefactor[i]/2.0;
        m=m+1;
      }   //  (end loop over bins)
    }    //  (end loop over 5 slopes)
    // shift next start value by two bins to take into account
    // overlap at the range boundaries.
    charge=charge-binwidth[4]*rangefactor[i]*2.0;
  }   // (end loop over 4 ranges)

  //    cout<<" FADC code "<<endl;
  //    int j1,j2,j3,j4;
  //   for(int i=0; i<32; i++) {
  //      j1=i;
  //      j2=j1+32;
  //     j3=j2+32;
  //     j4=j3+32;
  //     cout<<"   "<< j1 <<"  "<<fadc2charge[j1]
  //        <<"   "<< j2 <<"  "<<fadc2charge[j2]
  //        <<"   "<< j3 <<"  "<<fadc2charge[j3]
  //         <<"   "<< j4 <<"  "<<fadc2charge[j4]<<endl;
  //    }
}