Pystan log-probability confusing!

See how lp__ is calculated here: How exactly is lp__ being calculated - #4 by Bob_Carpenter

The important piece you are missing here is

So if you would not have constraints then lp__ would probably match exactly. But that will obviously not help you fit the model.