Analysis Software
Documentation for
sPHENIX
simulation software
Home page
Related Pages
Modules
Namespaces
Classes
Files
Examples
External Links
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
RKTools.h
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file RKTools.h
1
/* Copyright 2008-2009, Technische Universitaet Muenchen,
2
Authors: Christian Hoeppner & Sebastian Neubert & Johannes Rauch
3
4
This file is part of GENFIT.
5
6
GENFIT is free software: you can redistribute it and/or modify
7
it under the terms of the GNU Lesser General Public License as published
8
by the Free Software Foundation, either version 3 of the License, or
9
(at your option) any later version.
10
11
GENFIT is distributed in the hope that it will be useful,
12
but WITHOUT ANY WARRANTY; without even the implied warranty of
13
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14
GNU Lesser General Public License for more details.
15
16
You should have received a copy of the GNU Lesser General Public License
17
along with GENFIT. If not, see <http://www.gnu.org/licenses/>.
18
*/
19
24
#ifndef genfit_RKTools_h
25
#define genfit_RKTools_h
26
27
#include <stddef.h>
28
#include <algorithm>
29
#include <initializer_list>
30
31
namespace
genfit {
32
33
template
<
size_t
nRows,
size_t
nCols>
34
struct
RKMatrix
{
35
double
vals
[nRows * nCols];
36
37
RKMatrix
() =
default
;
38
RKMatrix
(
const
RKMatrix
&) =
default
;
39
RKMatrix
(std::initializer_list<double> initList) {
40
std::copy(initList.begin(), initList.end(),
vals
);
41
};
42
43
double
&
operator()
(
size_t
iRow,
size_t
iCol) {
44
return
vals
[nCols*iRow + iCol];
45
}
46
double
&
operator[]
(
size_t
n
) {
47
return
vals
[
n
];
48
}
49
const
double
&
operator[]
(
size_t
n
)
const
{
50
return
vals
[
n
];
51
}
52
53
double
*
begin
() {
return
vals
; }
54
double
*
end
() {
return
vals
+ nRows * nCols; }
55
const
double
*
begin
()
const
{
return
vals
; }
56
const
double
*
end
()
const
{
return
vals
+ nRows * nCols; }
57
58
RKMatrix<nRows, nCols>
&
operator=
(
const
RKMatrix<nRows, nCols>
& o) {
59
std::copy(o.
begin
(), o.
end
(), this->
begin
());
60
return
*
this
;
61
}
62
63
void
print
();
64
};
65
66
typedef
RKMatrix<1, 3>
M1x3
;
67
typedef
RKMatrix<1, 4>
M1x4
;
68
typedef
RKMatrix<1, 7>
M1x7
;
69
typedef
RKMatrix<5, 5>
M5x5
;
70
typedef
RKMatrix<6, 6>
M6x6
;
71
typedef
RKMatrix<7, 7>
M7x7
;
72
typedef
RKMatrix<6, 5>
M6x5
;
73
typedef
RKMatrix<7, 5>
M7x5
;
74
typedef
RKMatrix<5, 6>
M5x6
;
75
typedef
RKMatrix<5, 7>
M5x7
;
76
80
namespace
RKTools {
81
82
void
J_pMTxcov5xJ_pM
(
const
M5x7
& J_pM,
const
M5x5
& cov5,
M7x7
& out7);
83
void
J_pMTxcov5xJ_pM
(
const
M5x6
& J_pM,
const
M5x5
& cov5,
M6x6
& out6);
84
85
void
J_MpTxcov7xJ_Mp
(
const
M7x5
& J_Mp,
const
M7x7
& cov7,
M5x5
& out5);
86
void
J_MpTxcov6xJ_Mp
(
const
M6x5
& J_Mp,
const
M6x6
& cov6,
M5x5
& out5);
87
88
void
J_pMTTxJ_MMTTxJ_MpTT
(
const
M7x5
& J_pMT,
const
M7x7
& J_MMT,
const
M5x7
& J_MpT,
M5x5
& J_pp);
89
90
void
Np_N_NpT
(
const
M7x7
& Np,
M7x7
&
N
);
91
92
void
printDim
(
const
double
* mat,
unsigned
int
dimX,
unsigned
int
dimY);
93
94
}
95
96
template
<
size_t
nRows,
size_t
nCols>
97
inline
void
98
RKMatrix<nRows, nCols>::print
() {
99
RKTools::printDim
(this->vals, nRows, nCols);
100
}
101
102
}
/* End of namespace genfit */
105
#endif // genfit_RKTools_h
106
GenFit
blob
master
trackReps
include
RKTools.h
Built by
Jin Huang
. updated:
Sat Feb 17 2024 22:18:17
using
1.8.2 with
sPHENIX GitHub integration