helper function returning output parameters which an expr depends on

meta/Parameters.m

@@ -77,6 +77,9 @@
 GetDependenceSPhenoRules::usage="Returns list of replacement rules for
  symbols for which a DependenceSPheno rule is defined";
 
+GetOutputParameterDependencies::usage="Returns list of output
+ parameters which appear in the given expression";
+
 CreateLocalConstRefs::usage="creates local const references to model
 parameters / input parameters.";
 
@@ -1240,6 +1243,15 @@
     Cases[SARAH`ParameterDefinitions,
           {parameter_, {___, SARAH`DependenceSPheno -> value:Except[None], ___}} :> RuleDelayed[parameter, value]];
 
+GetSARAHParameters[] :=
+    (#[[1]])& /@ SARAH`SARAHparameters;
+
+GetOutputParameterDependencies[expr_] :=
+    Complement[Select[Join[GetSARAHParameters[],
+                           GetDependenceSPhenoSymbols[]],
+                      (!FreeQ[expr,#])&],
+               GetModelParameters[]];
+
 CreateInputParameterArrayGetter[inputParameters_List] :=
     Module[{get = "", paramCount = 0, name = "", par,
             type, i, assignment = "", nAssignments = 0},