OK.
Found the problem in the code. Very silly (put embraced face here:)
I have accidentally fed the LL function with the shock vector instead of the stim and vise verse. Once this was corrected, 30 trials are enough to recover alpha relatively accurate.
Now, back to how to implement this in pymc