I think P1 should just be p in your subtitle?

I assume theta_1 and theta_2 are a function of the material? Or of the normal? Or are they essentially just parameters to the BDRF that are used to compute radiance?