Allow robustHesse without a fit

cms-analysis · Aug 30, 2023 · 5b061b9 · 5b061b9
1 parent 9e40e3c
commit 5b061b9
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/src/RobustHesse.cc b/src/RobustHesse.cc
@@ -101,6 +101,13 @@ int RobustHesse::setParameterStencil(unsigned i) {
   double boundLo = rrv->getMin();
   double boundHi = rrv->getMax();
 
+  // If we skip initial fit (e.g. to compute from asimov point)
+  // need a guess for the initial step size
+  if ( rrv->getError() == 0.0 ) {
+    valLo = x + 1e-3*(boundLo - x);
+    valHi = x + 1e-3*(boundHi - x);
+  }
+
   bool closeToLo = valLo < boundLo;
   bool closeToHi = valHi > boundHi;