
A high statistics data sample of the $K^{+}\to\mu^{+}\nu_{\mu}$ decay was
accumulated by the OKA experiment in 2012. The missing mass analysis was
performed to search for the decay channel $K^{+}\to\mu^{+}\nu_{H}$ with a
hypothetic stable heavy neutrino in the final state. The obtained missing mass
spectrum does not show peaks which could be attributed to existence of stable
heavy neutrinos in the mass range (220 $< m_{\nu_{H}} <$ 375) MeV$/c^{2}$. As a
result, we obtain upper limits on the branching ratio and on the value of the
mixing element $U_{\mu H}^{2}$.

Recent results from OKA setup concerning form factor studies in $K_{e3}$
decay are presented. About 5.25M events are selected for the analysis. The
linear and quadratic slopes for the decay formfactor $f_{+}(t)$ are measured:
$\lambda'_{+}=(26.1 \pm 0.35 \pm 0.28 )\times 10^{3}$, $\lambda"_{+}=(1.91 \pm
0.19 \pm 0.14)\times 10^{3}$. The scalar and tensor contributions are
compatible with zero. Several alternative parametrizations are tried: the Pole
fit parameter is found to be $M_V = 891 \pm 2.0$ MeV ; the parameter of the
Dispersive parametrization is measured to be $\Lambda_+ =(24.58 \pm 0.18)
\times 10^{3}$. The presented results are considered as preliminary.

Results of study of the radiative kaon decay K>pi0 e nu gamma at ISTRA+
setup are presented. 4476 events of this decay have been observed. The
branching ratio is measured.

The branching ratio in the chargedpion kinetic energy region of 55 to 90 MeV
for the direct photon emission in the K > pi pi0 gamma decay has been
measured using inflight decays detected with the ISTRA+ setup operating in the
25 GeV/c negative secondary beam of the U70 PS. The value
Br(DE)=[0.37+0.39(stat)+0.10(syst)]*10^(5) obtained from the analysis of 930
completely reconstructed events is consistent with the average value of two
stoppedkaon experiments, but it differs by 2.5 standard deviations from the
average value of three inflightkaon experiments. The result is also compared
with recent theoretical predictions.

The decay K > pi0 mu nu has been studied using inflight decays detected
with the "ISTRA+" spectrometer. About 540K events were collected for the
analysis. The lambda+ and lambda0 slope parameters of the decay formfactors
f+(t), f0(t) have been measured : lambda+ = 0.0277+0.0013 (stat)+0.0009
(syst), lambda0 = 0.0183+0.0011(stat)+0.0006(syst), and
d(lambda0)/d(lambda+)=0.348. The limits on the possible tensor and scalar
couplings have been derived: fT/f+(0)=0.0007 + 0.0071, fS/f+(0)=0.0017 +
0.0014. No visible nonlinearity in the formfactors have been observed.

The decay K > pi0 e nu has been studied using inflight decays detected
with the "ISTRA+" spectrometer working at the 25 GeV negative secondary beam of
the U70 PS. About 550K events were used for the analysis. The lambda+
parameter of the vector formfactor has been measured: lambda+ = 0.0286 +
0.0008 (stat) + 0.0006(syst). The limits on the possible tensor and scalar
couplings have been obtained: f(T)/f+(0)=0.021 +0.064 0.075 (stat) +
0.026(syst) ; f(S)/f+(0)=0.002 +0.020 0.022 (stat) + 0.003(syst)

The Dalitz plot slope parameters g, h and k for the K > pi0 pi0 pi decay
have been measured using inflight decays detected with the ISTRA+ setup
operating in the 25 GeV negative secondary beam of the U70 PS. About 252 K
events with fourmomenta measured for the pi and four involved photons were
used for the analysis. The values obtained g=0.627+/0.004(stat)+/0.010(syst),
h=0.046+/0.004(stat)+/0.012(syst), k=0.001+/0.001(stat)+/0.002(syst) are
consistent with the world averages dominated by K+ data, but have significantly
smaller errors.

An explicit solution of timehomogeneous pure birth branching processes is
described. It gives alternative extensions for the negative binomial
distribution (branching processes with immigration) and for the FurryYule
distribution (branching processes without immigration).

An explicit solution of noncritical timehomogeneous branching processes is
described.

Recursive solution for a general homogeneous in time pure birth branching
process with simultaneous production of any number of particles and with
continuous evolution parameter is given. Calculational algorithm based on the
use of Koenigs function and functional Schroder equation is described. It is
shown that multiplicity distributions in e+e annihilation into hadrons for
c.m. energies up to 189 GeV are well described by the modified negative
binomial distribution, explained by simple pure birth branching process without
multiple simultaneous particle production. The energy dependence of the
evolution parameter is also discussed.

It is shown that at collider energies experimental multiplicity distributions
are well parameterized by a sum of GuptaSarma distributions. This extends
earlier description of the lower energy data by the two parameter sum of
Poissonians. Implications of the proposed parametrization for LHC are
discussed.

It is shown that simple extension of the modified negative binomial
distribution describes negatively charged particle multiplicity distributions
in e+e annihilation, measured in the whole phase space, as well as the
modified negative binomial.