Commit 7f1e890f authored by Breno Rilho Lemos's avatar Breno Rilho Lemos 💬
Browse files

Fixed Lorentz Factor issue

parent 2c7812ba
......@@ -11,7 +11,6 @@
// number of tracks, vx, vy, vz, magnetic field => esd-summary
// mass, charge, px, py, pz, Vx, Vy, Vz => esd-detail
//
// For now, we are considering classical mechanics, so p = mv (without relativistic considerations)
//
......@@ -88,7 +87,7 @@ int main(int argc, char **argv) {
fscanf(detail,"%lf %lf %lf %lf %lf %lf %lf %lf \n",&m[i],&q[i],&px[i],&py[i],&pz[i],&Vx[i],&Vy[i],&Vz[i]);
lorentz[i] = 1 / sqrt( 1 - Vx[i]*Vx[i]+Vy[i]*Vy[i]+Vz[i]*Vz[i] ); // Lorentz factor [SOMETHING FUNNY HERE...]
lorentz[i] = 1 / sqrt( 1 - (Vx[i]*Vx[i]+Vy[i]*Vy[i]+Vz[i]*Vz[i]) ); // Lorentz factor [SOMETHING FUNNY HERE...]
}
......@@ -161,7 +160,7 @@ int main(int argc, char **argv) {
j++;
} while( radius < TPCRADIUS && radius >= prevRadius ); // if the instant of time equals half the period, stop writing! (w = 2*PI / period) OR stop writing if it has taken half a lap (previous radius is, then, greater than current radius)
} while( radius < TPCRADIUS && radius >= prevRadius ); // if the instant of time equals half the period, stop writing! (w = 2*PI / period) OR stop writing if it has taken half a lap (previous vertex distance is, then, greater than current vertex distance)
}
......
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