With[{T=10},
sol=ParametricNDSolve[{
(x^\[Prime])[t]==a[t],(y^\[Prime])[t]==b[t],
x[0]==-5,y[0]==0,
a[0]==Cos[\[Theta]],b[0]==Sin[\[Theta]],
WhenEvent[-2-x[t]+x[t]^3+y[t]^2==0,
{vx=(4 b[t] (1-3 x[t]^2) y[t]+a[t] (-1+6 x[t]^2-9 x[t]^4+4 y[t]^2))/((1-3 x[t]^2)^2+4 y[t]^2),
vy=(4 a[t] (1-3 x[t]^2) y[t]+b[t] ((1-3 x[t]^2)^2-4 y[t]^2))/((1-3 x[t]^2)^2+4 y[t]^2),
a[t]->vx,b[t]->vy}]
},{x,y},{t,0,T},{\[Theta]},DiscreteVariables->{a,b}];
Show[{
ParametricPlot[Evaluate[Table[{x[\[Theta]][t],y[\[Theta]][t]}/.sol,{\[Theta],0\[Degree],30\[Degree],.3\[Degree]}]],{t,0,T},

With[{T = 10},
 sol = ParametricNDSolve[{
    Derivative[1][x][t] == a[t], Derivative[1][y][t] == b[t],
    x[0] == -5, y[0] == 0,
    a[0] == Cos[\[Theta]], b[0] == Sin[\[Theta]],
    WhenEvent[-2 - x[t] + x[t]^3 + y[t]^2 == 0,
     {vx = (
       4 b[t] (1 - 3 x[t]^2) y[t] +
        a[t] (-1 + 6 x[t]^2 - 9 x[t]^4 + 4 y[t]^2))/((1 -
          3 x[t]^2)^2 + 4 y[t]^2),
      vy = (
       4 a[t] (1 - 3 x[t]^2) y[t] +
        b[t] ((1 - 3 x[t]^2)^2 - 4 y[t]^2))/((1 - 3 x[t]^2)^2 +
        4 y[t]^2),