Simulation of the process e^+e^- rightarrow W^+W^- with the heavy right-handed neutrino exchange at 1 TeV future lepton colliders

We study potential contribution of the heavy right-handed neutrino exchange in the process $e^{+}e^{-} \rightarrow W^{+}W^{-}$. This process is sensitive to heavy neutrinos with masses larger than $\sqrt{s}$. The Monte Carlo simulation of the studied process is performed assuming the Seesaw type-I model, where heavy right-handed neutrinos (heavy neutral leptons, HNL) are introduced in the leptonic sector. Within the Standard Model (SM), the process has a large cross section described by diagrams with $s$-channel $Z / \gamma$ exchange and $t$-channel active neutrino exchange. Respectively, the $t$-channel right-handed neutrino exchange amplitude will interfere with these SM amplitudes. However, the angular distributions of the $W$ boson production and decay are different for the right-handed neutrino and SM amplitudes. That can be used to evaluate potential HNL contribution using the extended likelihood method. The simulation of the $e^{+}e^{-} \rightarrow W^{+}W^{-}$ process is performed at the 1 TeV center-of-mass energy and polarisation $P_{e^+e^-}$ of (20 %, -80 %), which is a standard option for the future linear $e^+e^-$ collider ILC. Both $W$ bosons are reconstructed from two hadronic jets. Simulation of the SM background processes is also done. The beam-induced backgrounds and the initial state radiation (ISR) effects are taken into account. The majority of background processes are effectively suppressed by the cuts on the invariant masses of two and four jets. Finally, we obtain upper limits on the mixing parameter $|V_{eN}|^2$ as a function of $M(N)$.