```
p0.perp()/tan(vec20.theta())
```
==>
```
//define outside vec20R = vec20.perp(); vec20Z = vec20.z();
p0.perp()/vec20R*vec20Z
```

probably the rest outweighs this cost, still, something to keep in mind for the future