# NOT RUN {
o <- r4vel(5)
vel_to_4mom(o) # five objects, random velocities, all unit mass
vel_to_4mom(o, 1:5) # five objects with identical 4vel, masses 1..5
vel_to_4mom(o[1],1:5) # five objects with random velocities, masses 1..5
p3 <- r3vel(5)
p_to_4mom(p3,E=1)
p_to_4mom(p3,E=10) # slower
p_to_4mom(p3,E=100) # even slower
P <- vel_to_4mom(as.3vel(c(0.8,0,0)) + r3vel(7,0.01))
reflect(P)
reflect(P,c(1,1,1))
sum(P)
# }
Run the code above in your browser using DataCamp Workspace