void test() {
   // === Provokes error
   // root -l -q test.C
   // or
   // lldb -- root.exe -l -q test.C
   // break set --name TH1::KolmogorovTest

   const size_t NBINS = 200;

   TH1F * h0 = new TH1F("h0", "h0", NBINS, 0., 1.);
   TH1F * h1 = new TH1F("h1", "h1", NBINS, 0., 1.);

   for (size_t i = 0; i<NBINS-1; ++i) {
      h0->AddBinContent(i, i);
      h1->AddBinContent(i, i);
   }
   h0->AddBinContent(1, -20.);
   h1->AddBinContent(1, -20.);

   // std::cout << "Integral: " << h0->ComputeIntegral(true) << std::endl;
   std::cout << "Before : h0[1]: " << h0->GetBinContent(1) << std::endl;
   std::cout << "Res    : " << h0->KolmogorovTest(h1, "X") << std::endl;
   std::cout << "After  : h0[1]: " << h0->GetBinContent(1) << std::endl;
}