Likelihood ratio test for bi dimensional data

ROOT Version: 6.16
_Platform:_MacOs Mojave
Compiler: gcc

Hello Everyone. This might not be the best place to post this question but I canβt find a way out of this problem.
1- bosons have an angular distribution in ΞΈ and π

πΉ(π,π)=34π[0.5(1βπΌ)+(0.5)(3πΌβ1)πππ (π)2βπ½π ππ(π)2πππ (2π)β2βΎβπΎπ ππ(2π)πππ (π)]

With πΌ=0.65,π½=0.06,πΎ=β0.18

I sample some Monte Carlo events from this distribution (50000). Now 0- bosons are uniformly distributed π,π,
that means that the pdf is a constant:
π(π,π)=12π2πππβ[0,π],πβ[0,2π]

Iβm asked to test the hypothesis that the Montecarlo data comes from 1- or 0- distribution. This sounds really trivial of course. I opted for a simple likelihood ratio. Note now that πΉ(π,π) reduces to uniform when πΌ=πΌ0=1./3,π½=π½0=0,πΎ=πΎ0=0:

π=βπ(π,π|πΌ0,π½0,πΎ0) / βπΉ(π,π|πΌ,π½,πΎ)

So basically what I do now is simple taking every π and π in my 50000 events and plug it into each pdf and get the value or the likelihood for that instance to be sampled given the uniform or alternative hypothesis.

Since both values are, by definition of pdf minor than 1, then this multiplication computation exceeds the numerical precision of a long double in c++. What I did then is turn into a logarithm:
πππ(π)=βπππ(π(ππ,ππ|πΌ0,π½0,πΎ0))ββπππ(πΉ(ππ,ππ|πΌ,π½,πΎ))

This value turns out to be:
πππ(π)=β40051βπ=πxp(β40051)

with the method I used before.

I expected a really small value for the statistics since a sample from one of the hypothesis Iβm testing on but can it be this small?
Furthermore how do I compute the rejection region according to the Neyman Pearson Lemma?
Is this the best way do do a statistical hypothesis test in 2 dimensions?

Iβll post some snapshot of the distributions and their fits. Thank you everyone who can help, Iβm kind of desperate ahah.

result_ML_unbinned_2.pdf (104.0 KB)
result_ML_unbinned.pdf (76.3 KB)

@moneta perhaps you can help here please?