TAxis *axis = h1->GetXaxis();
double bmin = axis->FindBin(123);
double bmax = axis->FindBin(456);
double sum1 = h1->Integral(bmin,bmax);