1 #ifndef __JFIT__JMATRIXNZ__
2 #define __JFIT__JMATRIXNZ__
16 namespace JPP {
using namespace JFIT; }
66 set(pos, __begin, __end, alpha, sigma);
100 const double ta = alpha *
PI / 180.0;
101 const double ct = cos(ta);
102 const double st =
sin(ta);
107 for (
T i = __begin; i != __end; ++i) {
109 const double dx = i->getX() - pos.
getX();
110 const double dy = i->getY() - pos.
getY();
111 const double dz = i->getZ() - pos.
getZ();
113 const double R = sqrt(dx*dx + dy*dy);
125 x *= (dz*ct - dx*st);
126 y *= (dz*ct - dy*st);
127 v *= -(dx*ct + dz*st);
128 w *= -(dy*ct + dz*st);
133 for (
int i = 0; i !=
N; ++i) {
135 for (
int j = 0;
j != i; ++
j) {
137 (*this)(
j, i) = (*
this)(i,
j);
198 return (first.
x * second.
x +
void set(const JVector3D &pos, T __begin, T __end, const double alpha, const double sigma)
Set co-variance matrix.
double getKappaC()
Get average R-dependence of arrival time of Cherenkov light (a.k.a.
JMatrixNZ()
Default contructor.
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
static double getDot(const variance &first, const variance &second)
Get dot product.
then echo The file $DIR KM3NeT_00000001_00000000 root already please rename or remove it first
JMatrixNZ(const JVector3D &pos, T __begin, T __end, const double alpha, const double sigma)
Constructor.
variance(const double __x, const double __y, const double __v, const double __w)
Constructor.
Data structure for vector in three dimensions.
Determination of the co-variance matrix of hits for a track along z-axis (JFIT::JLine1Z).
do set_variable OUTPUT_DIRECTORY $WORKDIR T
std::vector< variance > buffer
static const double PI
Mathematical constants.
double getY() const
Get y position.
then JCookie sh JDataQuality D $DETECTOR_ID R
then usage $script< input file >[option[primary[working directory]]] nWhere option can be N
const double getInverseSpeedOfLight()
Get inverse speed of light.
double getX() const
Get x position.
variance()
Default constructor.
double getDot(const JVectorND &v) const
Get dot product.
Auxiliary data structure for co-variance calculation.
double getZ() const
Get z position.