2020-06-10 16:23:33 +02:00
|
|
|
# Landau sample
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
## Sample
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
Sample N = 50'000 random points following $L(x)$
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
$$
|
|
|
|
L(x) = \frac{1}{\pi} \int \limits_{0}^{+ \infty}
|
|
|
|
dt \, e^{-t \ln(t) -xt} \sin (\pi t)
|
|
|
|
$$
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
. . .
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
gsl_ran_Landau(gsl_rng)
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
## Compatibility results:
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
Median:
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
:::: {.columns}
|
|
|
|
::: {.column width=50%}
|
|
|
|
- $t = 0.761$
|
|
|
|
- $p = 0.446$
|
|
|
|
:::
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
::: {.column width=50%}
|
|
|
|
$$
|
|
|
|
\hence \text{Compatible!}
|
|
|
|
$$
|
|
|
|
:::
|
|
|
|
::::
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
\vspace{10pt}
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
. . .
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
Mode:
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
:::: {.columns}
|
|
|
|
::: {.column width=50%}
|
|
|
|
- $t = 1.012$
|
|
|
|
- $p = 0.311$
|
2020-06-09 16:53:16 +02:00
|
|
|
:::
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
::: {.column width=50%}
|
2020-06-09 16:53:16 +02:00
|
|
|
$$
|
2020-06-10 16:23:33 +02:00
|
|
|
\hence \text{Compatible!}
|
2020-06-09 16:53:16 +02:00
|
|
|
$$
|
2020-06-10 16:23:33 +02:00
|
|
|
:::
|
|
|
|
::::
|
|
|
|
|
|
|
|
\vspace{10pt}
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
. . .
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
FWHM:
|
|
|
|
|
|
|
|
:::: {.columns}
|
|
|
|
::: {.column width=50%}
|
|
|
|
- $t=1.338$
|
|
|
|
- $p=0.181$
|
|
|
|
:::
|
|
|
|
|
|
|
|
::: {.column width=50%}
|
2020-06-09 16:53:16 +02:00
|
|
|
$$
|
2020-06-10 16:23:33 +02:00
|
|
|
\hence \text{Compatible!}
|
2020-06-09 16:53:16 +02:00
|
|
|
$$
|
2020-06-10 16:23:33 +02:00
|
|
|
:::
|
|
|
|
::::
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
# Moyal sample
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
## Sample
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
Sample N = 50'000 random points following $M_{\mu \sigma}(x)$
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
$$
|
|
|
|
M_{\mu \sigma}(x) = \frac{1}{\sqrt{2 \pi} \sigma} \exp
|
|
|
|
\left[ - \frac{1}{2} \left(
|
|
|
|
\frac{x - \mu}{\sigma}
|
|
|
|
+ e^{-\frac{x - \mu}{\sigma}} \right) \right]
|
|
|
|
$$
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
. . .
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
reverse sampling
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
- sampling $y$ uniformly in [0, 1] $\hence x = Q_M(y)$
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
## Compatibility results:
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
Median:
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
:::: {.columns}
|
|
|
|
::: {.column width=50%}
|
|
|
|
- $t = 669.940$
|
|
|
|
- $p = 0.000$
|
|
|
|
:::
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
::: {.column width=50%}
|
|
|
|
$$
|
|
|
|
\hence \text{Not compatible!}
|
|
|
|
$$
|
|
|
|
:::
|
|
|
|
::::
|
|
|
|
|
|
|
|
\vspace{10pt}
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
. . .
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
Mode:
|
|
|
|
|
|
|
|
:::: {.columns}
|
|
|
|
::: {.column width=50%}
|
|
|
|
- $t = 0.732$
|
|
|
|
- $p = 0.464$
|
|
|
|
:::
|
|
|
|
|
|
|
|
::: {.column width=50%}
|
|
|
|
$$
|
|
|
|
\hence \text{Compatible!}
|
|
|
|
$$
|
|
|
|
:::
|
|
|
|
::::
|
|
|
|
|
|
|
|
\vspace{10pt}
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
. . .
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
FWHM:
|
|
|
|
|
|
|
|
:::: {.columns}
|
|
|
|
::: {.column width=50%}
|
|
|
|
- $t = 1.329$
|
|
|
|
- $p = 0.184$
|
|
|
|
:::
|
|
|
|
|
|
|
|
::: {.column width=50%}
|
|
|
|
$$
|
|
|
|
\hence \text{Compatible!}
|
|
|
|
$$
|
|
|
|
:::
|
|
|
|
::::
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
# KS results
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
## Samples results
|
|
|
|
|
|
|
|
$N = 50000$ sampled points
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
. . .
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
Landau sample:
|
|
|
|
|
|
|
|
:::: {.columns}
|
|
|
|
::: {.column width=50%}
|
|
|
|
- $D = 0.004$
|
|
|
|
- $p = 0.379$
|
|
|
|
:::
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
::: {.column width=50%}
|
|
|
|
$$
|
|
|
|
\hence \text{Compatible!}
|
|
|
|
$$
|
|
|
|
:::
|
|
|
|
::::
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
\vspace{10pt}
|
2020-06-09 16:53:16 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
. . .
|
2020-06-09 18:28:53 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
Moyal sample:
|
2020-06-09 18:28:53 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
:::: {.columns}
|
|
|
|
::: {.column width=50%}
|
|
|
|
- $D = 0.153$
|
|
|
|
- $p = 0.000$
|
|
|
|
:::
|
2020-06-09 18:28:53 +02:00
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
::: {.column width=50%}
|
|
|
|
$$
|
|
|
|
\hence \text{Not compatible!}
|
|
|
|
$$
|
|
|
|
:::
|
|
|
|
::::
|
2020-06-09 18:28:53 +02:00
|
|
|
|
|
|
|
|
2020-06-10 16:23:33 +02:00
|
|
|
# Trapani results
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
|
|
|
|
## Samples results
|
|
|
|
|
|
|
|
. . .
|
|
|
|
|
|
|
|
Landau sample:
|
|
|
|
|
|
|
|
:::: {.columns}
|
2020-06-09 18:28:53 +02:00
|
|
|
::: {.column width=33%}
|
|
|
|
$$
|
|
|
|
\mu_1
|
|
|
|
\begin{cases}
|
|
|
|
\Theta = 0.255 \\
|
|
|
|
p = 0.614
|
|
|
|
\end{cases}
|
|
|
|
$$
|
2020-06-09 16:53:16 +02:00
|
|
|
:::
|
|
|
|
|
2020-06-09 18:28:53 +02:00
|
|
|
::: {.column width=33% .c}
|
2020-06-09 16:53:16 +02:00
|
|
|
$$
|
2020-06-09 18:28:53 +02:00
|
|
|
\mu_2
|
|
|
|
\begin{cases}
|
|
|
|
\Theta = 0.432 \\
|
|
|
|
p = 0.511
|
|
|
|
\end{cases}
|
2020-06-09 16:53:16 +02:00
|
|
|
$$
|
|
|
|
:::
|
|
|
|
|
2020-06-09 18:28:53 +02:00
|
|
|
::: {.column width=33% .c}
|
|
|
|
$$
|
|
|
|
\hence \text{Infinite!}
|
|
|
|
$$
|
|
|
|
:::
|
|
|
|
::::
|
2020-06-09 16:53:16 +02:00
|
|
|
|
|
|
|
. . .
|
|
|
|
|
2020-06-09 18:28:53 +02:00
|
|
|
\vspace{20pt}
|
|
|
|
|
2020-06-09 16:53:16 +02:00
|
|
|
Moyal sample:
|
|
|
|
|
|
|
|
:::: {.columns}
|
2020-06-09 18:28:53 +02:00
|
|
|
::: {.column width=33%}
|
|
|
|
$$
|
|
|
|
\mu_1
|
|
|
|
\begin{cases}
|
|
|
|
\Theta^2 = 106 \\
|
|
|
|
p = 0.000
|
|
|
|
\end{cases}
|
|
|
|
$$
|
|
|
|
:::
|
|
|
|
|
|
|
|
::: {.column width=33%}
|
|
|
|
$$
|
|
|
|
\mu_2
|
|
|
|
\begin{cases}
|
|
|
|
\Theta^2 = 162 \\
|
|
|
|
p = 0.000
|
|
|
|
\end{cases}
|
|
|
|
$$
|
2020-06-09 16:53:16 +02:00
|
|
|
:::
|
|
|
|
|
2020-06-09 18:28:53 +02:00
|
|
|
::: {.column width=33% .c}
|
2020-06-09 16:53:16 +02:00
|
|
|
$$
|
2020-06-09 18:28:53 +02:00
|
|
|
\hence \text{Finite!}
|
2020-06-09 16:53:16 +02:00
|
|
|
$$
|
|
|
|
:::
|
|
|
|
::::
|