Throw runtime error when parameter is set via command line in a range outside its value

Author: maxgalli

src/utils.cc

@@ -765,6 +765,15 @@ void utils::setModelParameters( const std::string & setPhysicsModelParameterExpr
                     cout << "Set Default Value of Parameter " << target
                      << " To : " << PhysicsParameterValue << "\n";
                     tmpParameter->setVal(PhysicsParameterValue);
+                    if (tmpParameter->getVal() != PhysicsParameterValue) {
+                      throw std::runtime_error(
+                          Form("Parameter %s value %g is outside its range [%g, %g]. "
+                               "Use --rMin/--rMax or --setParameterRanges to adjust the range.",
+                               target.c_str(),
+                               PhysicsParameterValue,
+                               tmpParameter->getMin(),
+                               tmpParameter->getMax()));
+                    }
                 }
             } else {
                 RooCategory *tmpCategory  = dynamic_cast<RooCategory*>(tmp);
@@ -791,9 +800,18 @@ void utils::setModelParameters( const std::string & setPhysicsModelParameterExpr
         if (isrvar) {
           RooRealVar *tmpParameter = dynamic_cast<RooRealVar*>(tmp);
           double PhysicsParameterValue = atof(SetParameterExpression[1].c_str());
-          cout << "Set Default Value of Parameter " << SetParameterExpression[0] 
-               << " To : " << PhysicsParameterValue << "\n";
+          cout << "Set Default Value of Parameter " << SetParameterExpression[0] << " To : " << PhysicsParameterValue
+               << "\n";
           tmpParameter->setVal(PhysicsParameterValue);
+          if (tmpParameter->getVal() != PhysicsParameterValue) {
+            throw std::runtime_error(
+                Form("Parameter %s value %g is outside its range [%g, %g]. "
+                     "Use --rMin/--rMax or --setParameterRanges to adjust the range.",
+                     SetParameterExpression[0].c_str(),
+                     PhysicsParameterValue,
+                     tmpParameter->getMin(),
+                     tmpParameter->getMax()));
+          }
         } else {
           RooCategory *tmpCategory  = dynamic_cast<RooCategory*>(tmp);
           int PhysicsParameterValue = atoi(SetParameterExpression[1].c_str());