Commit 702edd51 authored by kraus's avatar kraus

fix position of reference particle in temporal monitor

parent 3f0bdc0a
......@@ -103,7 +103,7 @@ bool Monitor::applyToReferenceParticle(const Vector_t &R,
dt * (R(2) + P(2) * recpgamma) > dt * middle) {
double frac = (middle - R(2)) / (P(2) * recpgamma);
double time = t + frac * dt;
Vector_t dR = (0.5 + frac) * P * recpgamma;
Vector_t dR = frac * P * recpgamma;
double ds = euclidean_norm(dR);
lossDs_m->addReferenceParticle(csTrafoGlobal2Local_m.transformFrom(R + dR),
csTrafoGlobal2Local_m.rotateFrom(P),
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment