Evaluation of fit.
801 for (
T hit = __begin; hit != __end; ++hit) {
807 const double D = hit->getDistance(position);
808 const double Vi = velocity.getInverseVelocity(D, hit->getZ(), position.
getZ());
809 const double toa_s =
value.
emission[hit->getEKey()].t1 + D * Vi;
811 const double u = (toa_s - hit->getValue()) / hit->getSigma();
812 const double W = sqrt(hit->getWeight());
814 successor += (W*W) * estimator->getRho(u);
816 H_t
H(1.0,
string.getGradient(parameters, hit->getPosition(), hit->getFloor()) * Vi);
818 H *= W * estimator->getPsi(u) / hit->getSigma();
829 V(i.t1, i.t1) +=
H.t1 *
H.t1;
833 if (hit->getFloor() != 0) {
835 switch (this->option) {
838 V(i.t1, i.vs) +=
H.t1 *
H.vs;
V(i.tx, i.vs) +=
H.tx *
H.vs;
V(i.ty, i.vs) +=
H.ty *
H.vs;
V(i.tx2, i.vs) +=
H.tx2 *
H.vs;
V(i.ty2, i.vs) +=
H.ty2 *
H.vs;
840 V(i.vs, i.t1) +=
H.vs *
H.t1;
841 V(i.vs, i.tx) +=
H.vs *
H.tx;
842 V(i.vs, i.ty) +=
H.vs *
H.ty;
843 V(i.vs, i.tx2) +=
H.vs *
H.tx2;
844 V(i.vs, i.ty2) +=
H.vs *
H.ty2;
846 V(i.vs, i.vs) +=
H.vs *
H.vs;
851 V(i.t1, i.tx2) +=
H.t1 *
H.tx2;
V(i.tx, i.tx2) +=
H.tx *
H.tx2;
V(i.ty, i.tx2) +=
H.ty *
H.tx2;
853 V(i.tx2, i.t1) +=
H.tx2 *
H.t1;
854 V(i.tx2, i.tx) +=
H.tx2 *
H.tx;
855 V(i.tx2, i.ty) +=
H.tx2 *
H.ty;
857 V(i.t1, i.ty2) +=
H.t1 *
H.ty2;
V(i.tx, i.ty2) +=
H.tx *
H.ty2;
V(i.ty, i.ty2) +=
H.ty *
H.ty2;
859 V(i.ty2, i.t1) +=
H.ty2 *
H.t1;
860 V(i.ty2, i.tx) +=
H.ty2 *
H.tx;
861 V(i.ty2, i.ty) +=
H.ty2 *
H.ty;
863 V(i.tx2, i.tx2) +=
H.tx2 *
H.tx2;
V(i.tx2, i.ty2) +=
H.tx2 *
H.ty2;
864 V(i.ty2, i.tx2) +=
H.ty2 *
H.tx2;
V(i.ty2, i.ty2) +=
H.ty2 *
H.ty2;
866 Y[i.tx2] += W *
H.tx2;
867 Y[i.ty2] += W *
H.ty2;
870 V(i.t1, i.tx) +=
H.t1 *
H.tx;
V(i.t1, i.ty) +=
H.t1 *
H.ty;
871 V(i.tx, i.t1) +=
H.tx *
H.t1;
V(i.ty, i.t1) +=
H.ty *
H.t1;
873 V(i.tx, i.tx) +=
H.tx *
H.tx;
V(i.tx, i.ty) +=
H.tx *
H.ty;
874 V(i.ty, i.tx) +=
H.ty *
H.tx;
V(i.ty, i.ty) +=
H.ty *
H.ty;
fit times of emission of emitters and tilt angles of strings with second order correction and stretch...
static const double H
Planck constant [eV s].
*fatal Wrong number of arguments esac JCookie sh typeset Z DETECTOR typeset Z SOURCE_RUN typeset Z TARGET_RUN set_variable PARAMETERS_FILE $WORKDIR parameters
JMatrixND & reset()
Set matrix to the null matrix.
fit times of emission of emitters and tilt angles of strings
do set_variable OUTPUT_DIRECTORY $WORKDIR T
JACOUSTICS::JModel::emission_type emission
const JPosition3D & getPosition() const
Get position.
fit times of emission of emitters and tilt angles of strings with second order correction ...
size_t getIndex(int id, double JString::*p) const
Get index of fit parameter for given string.
JACOUSTICS::JModel::string_type string
Data structure for position in three dimensions.
do echo Generating $dir eval D
double getZ() const
Get z position.