maxRapidity = 5;
deltaRapidity = .2;
bigBangEvent = {0 Sinh[rapidity], 0 Cosh[rapidity]};
e0 = Table[
bigBangEvent, {rapidity, -maxRapidity, maxRapidity, deltaRapidity}];
(*The e0 values here will all be {0,0}*)
e1 = Table[{1 Sinh[rapidity],
1 Cosh[rapidity]}, {rapidity, -maxRapidity, maxRapidity,
deltaRapidity}];
firstHalfLife = 0.4; fHL = firstHalfLife;
secondHalfLife = 1 - firstHalfLife; sHL = secondHalfLife;
initialRapidity = -3; iR = initialRapidity;
deltaRapidity = .1; dR = deltaRapidity;
resultingParticles = 12; rP = resultingParticles;
decayEvent = {fHL Sinh[iR], fHL Cosh[iR]};
nextWorldLinesBegin = Table[decayEvent, {n, -rP, +rP}];
nextWorldLinesEnd =
Table[decayEvent + {sHL*Sinh[iR + n*dR],
sHL*Cosh[iR + n*dR]}, {n, -rP, rP}];
e0 = Join[e0, nextWorldLinesBegin];
e1 = Join[e1, nextWorldLinesEnd];
(*Apply the Lorentz Transformation around the decayEvent.*)
decayEventList = Table[decayEvent, {n, 1, Length[e0]}];
e0 = e0 - decayEventList;
e1 = e1 - decayEventList;
LT[theta_] := {{Cosh[theta], -Sinh[theta]}, {-Sinh[theta],
Cosh[theta]}};
Manipulate[
ePrime0 = Transpose[LT[theta].Transpose[e0]];
ePrime1 = Transpose[LT[theta].Transpose[e1]];
milneWorldLines = Transpose[{ePrime0, ePrime1}];
ListLinePlot[milneWorldLines,
PlotRange -> {{-2, 2}, {-.5, 2}}], {{theta, iR}, iR - (rP*dR)/2,
iR + (rP*dR)/2}]