|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
#ifndef GenericMinL3Algorithm_H
#define GenericMinL3Algorithm_H
/** \class GenericMinL3Algorithm
* Implementation of the L3 Collaboration algorithm to solve a system Ax = B
* by minimization of |Ax-B| using an iterative linear approach
*
* \author R.Ofierzynski, CERN
*/
#include <vector>
#include <iostream>
class GenericMinL3Algorithm {
public:
/// Default constructor
/// CAVEAT: use normalise = true only if you know what you're doing!
GenericMinL3Algorithm(bool normalise = false);
/// Destructor
~GenericMinL3Algorithm();
/// run the Minimization L3 Algorithm "nIter" number of times, recalibrating the event matrix after each iteration with the new solution
/// returns the vector of calibration coefficients built from all iteration solutions
std::vector<float> iterate(const std::vector<std::vector<float> >& eventMatrix,
const std::vector<float>& energyVector,
int nIter);
/// perform one iteration using the Minimization L3 Algorithm
/// returns the vector of calibration coefficients
std::vector<float> iterate(const std::vector<std::vector<float> >& eventMatrix,
const std::vector<float>& energyVector);
protected:
private:
bool normaliseFlag;
};
#endif // GenericMinL3Algorithm_H
|