I'm trying to create a lognormal mock using a linear power spectrum made using the Eisenstein-Hu fitting functions. When I compute the power spectrum of the matter field, does that yield P_linear(k) or b_1^2 P_linear(k)? (b_1 is the bias parameter for the lognormal field transformation)
I saw the lognormal transformation in the source code, but I'm unable to figure out how the power spectrum transforms to get a factor of b_1^2.