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

Author: maxgalli

src/utils.cc

@@ -765,6 +765,11 @@ 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 +796,14 @@ 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] 
+          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());