Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
ParameterSettings.h
Go to the documentation of this file.
1// @(#)root/mathcore:$Id$
2// Author: L. Moneta Thu Sep 21 16:21:48 2006
3
4/**********************************************************************
5 * *
6 * Copyright (c) 2006 LCG ROOT Math Team, CERN/PH-SFT *
7 * *
8 * *
9 **********************************************************************/
10
11// Header file for class ParameterSettings
12
13#ifndef ROOT_Fit_ParameterSettings
14#define ROOT_Fit_ParameterSettings
15
16#include <string>
17
18#include "Math/Error.h"
19
20
21namespace ROOT {
22
23 namespace Fit {
24
25
26//___________________________________________________________________________________
27/**
28 Class, describing value, limits and step size of the parameters
29 Provides functionality also to set/retrieve values, step sizes, limits and fix the
30 parameters.
31
32 To be done: add constraints (equality and inequality) as functions of the parameters
33
34 @ingroup FitMain
35*/
37
38public:
39
40 /**
41 Default constructor
42 */
44 fValue(0.), fStepSize(0.1), fFix(false),
46 fName("")
47 {}
48
49
50 ///constructor for unlimited named Parameter
51 ParameterSettings(const std::string & name, double val, double err) :
52 fValue(val), fStepSize(err), fFix(false),
54 fName(name)
55 {}
56
57 ///constructor for double limited Parameter. The given value should be within the given limits [min,max]
58 ParameterSettings(const std::string & name, double val, double err,
59 double min, double max) :
60 fValue(val), fStepSize(err), fFix(false),
62 fName(name)
63 {
64 SetLimits(min,max);
65 }
66
67 ///constructor for fixed Parameter
68 ParameterSettings(const std::string & name, double val) :
69 fValue(val), fStepSize(0), fFix(true),
71 fName(name)
72 {}
73
74
75
76
77 /// set value and name (unlimited parameter)
78 void Set(const std::string & name, double value, double step) {
81 SetStepSize(step);
82 }
83
84 /// set a limited parameter. The given value should be within the given limits [min,max]
85 void Set(const std::string & name, double value, double step, double lower, double upper ) {
88 SetStepSize(step);
89 SetLimits(lower,upper);
90 }
91
92 /// set a fixed parameter
93 void Set(const std::string & name, double value) {
96 Fix();
97 }
98
99
100 /**
101 Destructor (no operations)
102 */
104
105 /// copy constructor and assignment operators (leave them to the compiler)
106
107public:
108
109 /// return parameter value
110 double Value() const { return fValue; }
111 /// return step size
112 double StepSize() const { return fStepSize; }
113 /// return lower limit value
114 double LowerLimit() const {return fLowerLimit;}
115 /// return upper limit value
116 double UpperLimit() const {return fUpperLimit;}
117 /// check if is fixed
118 bool IsFixed() const { return fFix; }
119 /// check if parameter has lower limit
120 bool HasLowerLimit() const {return fHasLowerLimit; }
121 /// check if parameter has upper limit
122 bool HasUpperLimit() const {return fHasUpperLimit; }
123 /// check if is bound
124 bool IsBound() const { return fHasLowerLimit || fHasUpperLimit; }
125 /// check if is double bound (upper AND lower limit)
126 bool IsDoubleBound() const { return fHasLowerLimit && fHasUpperLimit; }
127 /// return name
128 const std::string & Name() const { return fName; }
129
130 /** interaction **/
131
132 /// set name
133 void SetName(const std::string & name ) { fName = name; }
134
135 /// fix the parameter
136 void Fix() {fFix = true;}
137 /// release the parameter
138 void Release() {fFix = false;}
139 /// set the value
140 void SetValue(double val) {fValue = val;}
141 /// set the step size
142 void SetStepSize(double err) {fStepSize = err;}
143 /// set a double side limit,
144 /// if low == up the parameter is fixed if low > up the limits are removed
145 /// The current parameter value should be within the given limits [low,up].
146 /// If the value is outside the limits, then a new parameter value is set to = (up+low)/2
147 void SetLimits(double low, double up) {
148
149 if ( low > up ) {
150 RemoveLimits();
151 return;
152 }
153 if (low == up && low == fValue) {
154 Fix();
155 return;
156 }
157 if (low > fValue || up < fValue) {
158 MATH_INFO_MSG("ParameterSettings","lower/upper bounds outside current parameter value. The value will be set to (low+up)/2 ");
159 fValue = 0.5 * (up+low);
160 }
161 fLowerLimit = low;
162 fUpperLimit = up;
163 fHasLowerLimit = true;
164 fHasUpperLimit = true;
165 }
166 /// set a single upper limit
167 void SetUpperLimit(double up) {
168 fLowerLimit = 0.;
169 fUpperLimit = up;
170 fHasLowerLimit = false;
171 fHasUpperLimit = true;
172 }
173 /// set a single lower limit
174 void SetLowerLimit(double low) {
175 fLowerLimit = low;
176 fUpperLimit = 0.;
177 fHasLowerLimit = true;
178 fHasUpperLimit = false;
179 }
180
181 /// remove all limit
183 fLowerLimit = 0.;
184 fUpperLimit = 0.;
185 fHasLowerLimit = false;
186 fHasUpperLimit = false;
187 }
188
189
190
191protected:
192
193
194private:
195
196 double fValue; ///< parameter value
197 double fStepSize; ///< parameter step size (used by minimizer)
198 bool fFix; ///< flag to control if parameter is fixed
199 double fLowerLimit; ///< lower parameter limit
200 double fUpperLimit; ///< upper parameter limit
201 bool fHasLowerLimit; ///< flag to control lower parameter limit
202 bool fHasUpperLimit; ///< flag to control upper parameter limit
203
204 std::string fName; ///< parameter name
205
206};
207
208 } // end namespace Fit
209
210} // end namespace ROOT
211
212
213#endif /* ROOT_Fit_ParameterSettings */
#define MATH_INFO_MSG(loc, str)
Pre-processor macro to report messages which can be configured to use ROOT error or simply an std::io...
Definition Error.h:77
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void value
char name[80]
Definition TGX11.cxx:110
Class, describing value, limits and step size of the parameters Provides functionality also to set/re...
void Release()
release the parameter
bool IsFixed() const
check if is fixed
bool HasUpperLimit() const
check if parameter has upper limit
bool fFix
flag to control if parameter is fixed
ParameterSettings(const std::string &name, double val, double err)
constructor for unlimited named Parameter
ParameterSettings()
Default constructor.
void RemoveLimits()
remove all limit
void SetValue(double val)
set the value
double LowerLimit() const
return lower limit value
void SetName(const std::string &name)
interaction
const std::string & Name() const
return name
bool HasLowerLimit() const
check if parameter has lower limit
double fLowerLimit
lower parameter limit
double fUpperLimit
upper parameter limit
bool fHasUpperLimit
flag to control upper parameter limit
void SetStepSize(double err)
set the step size
void Set(const std::string &name, double value, double step)
set value and name (unlimited parameter)
double fStepSize
parameter step size (used by minimizer)
bool fHasLowerLimit
flag to control lower parameter limit
void SetLimits(double low, double up)
set a double side limit, if low == up the parameter is fixed if low > up the limits are removed The c...
ParameterSettings(const std::string &name, double val)
constructor for fixed Parameter
double Value() const
copy constructor and assignment operators (leave them to the compiler)
void Set(const std::string &name, double value)
set a fixed parameter
void Set(const std::string &name, double value, double step, double lower, double upper)
set a limited parameter. The given value should be within the given limits [min,max]
ParameterSettings(const std::string &name, double val, double err, double min, double max)
constructor for double limited Parameter. The given value should be within the given limits [min,...
~ParameterSettings()
Destructor (no operations)
double StepSize() const
return step size
void SetUpperLimit(double up)
set a single upper limit
std::string fName
parameter name
void Fix()
fix the parameter
double UpperLimit() const
return upper limit value
bool IsDoubleBound() const
check if is double bound (upper AND lower limit)
bool IsBound() const
check if is bound
void SetLowerLimit(double low)
set a single lower limit
TFitResultPtr Fit(FitObject *h1, TF1 *f1, Foption_t &option, const ROOT::Math::MinimizerOptions &moption, const char *goption, ROOT::Fit::DataRange &range)
Definition HFitImpl.cxx:133
This file contains a specialised ROOT message handler to test for diagnostic in unit tests.