Fit.
1213 using namespace std;
1214 using namespace JPP;
1227 for (data_type::const_iterator ix =
data.begin(); ix !=
data.end(); ++ix) {
1234 ndf += ix->second.size();
1243 double precessor = numeric_limits<double>::max();
1276 return { precessor /
estimator->getRho(1.0), ndf };
1303 for (
size_t i = 0; i != N; ++i) {
1309 h[i] = 1.0 / sqrt(
V(i,i));
1314 for (
size_t i = 0; i != N; ++i) {
1315 for (
size_t j = 0;
j != i; ++
j) {
1316 V(
j,i) *=
h[i] *
h[
j];
1321 for (
size_t i = 0; i != N; ++i) {
1327 for (
size_t col = 0; col != N; ++col) {
1334 catch (
const exception& error) {
1336 ERROR(
"JGandalf: " << error.what() << endl <<
V << endl);
1343 const double factor = 2.0;
1362 return { precessor /
estimator->getRho(1.0), ndf };
#define DEBUG(A)
Message macros.
static constexpr double LAMBDA_MIN
minimal value control parameter
static constexpr double LAMBDA_DOWN
multiplication factor control parameter
static constexpr double LAMBDA_MAX
maximal value control parameter
static constexpr double LAMBDA_UP
multiplication factor control parameter
static constexpr double EPSILON
maximal distance to minimum.
void evaluate(const data_type &data)
Evaluation of fit.
static constexpr int MAXIMUM_ITERATIONS
maximal number of iterations.
static constexpr double PIVOT
minimal value diagonal element of matrix
bool isFree() const
Check if parameter is free.
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Auxiliary data structure for floating point format specification.
JParameter_t R
maximal coincidence rate [Hz]
JParameter_t p1
1st order angle dependence coincidence rate
JParameter_t p2
2nd order angle dependence coincidence rate
JParameter_t p3
3rd order angle dependence coincidence rate
JParameter_t p4
4th order angle dependence coincidence rate
JParameter_t cc
fraction of signal correlated background
size_t getN() const
Get number of fit parameters.
JPMTParameters_t parameters[NUMBER_OF_PMTS]
void setIndex()
Set index of PMT used for fixed time offset.
JParameter_t t0
time offset [ns]
JParameter_t TTS
transition-time spread [ns]
JParameter_t bg
background [Hz/ns]
JParameter_t QE
relative quantum efficiency [unit]
void resize(const size_t size)
Resize matrix.
void solve(JVectorND_t &u)
Get solution of equation A x = b.