A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [@Chatrchyan:2008zzk].

Simulated event samples {#sec:mcsamples}
=======================

Signal and background processes are simulated using various event generators, while the CMS detector response is modelled with  [@GEANT4]. The quark-induced [$\PZ\PH$]{}and [$\PW\PH$]{}signal processes are generated at next-to-leading order (NLO) accuracy in QCD using the v2 [@Nason:2004rx; @POWHEG; @Alioli:2010xd] event generator extended with the Multi-scale improved NLO (MiNLO) procedure [@Hamilton2012; @Luisoni:2013kna], while the gluon-induced [$\PZ\PH$]{}process is generated at leading order (LO) accuracy with v2. The boson mass is set to 125for all signal samples. The production cross sections of the signal processes [@deFlorian:2016spz] are corrected as a function of [$\pt(\PV)$]{}to next-to-next-to-leading order (NNLO) QCD + NLO electroweak (EW) accuracy combining the <span style="font-variant:small-caps;">vhnnlo</span> [@Ferrera:2013yga; @Ferrera:2014lca; @Ferrera:2011bk; @Ferrera:2017zex], <span style="font-variant:small-caps;">vh@nnlo</span> [@Brein:2012ne; @Harlander:2013mla], and <span style="font-variant:small-caps;">hawk</span> v2.0 [@Denner:2014cla] generators as described in Ref. [@deFlorian:2016spz].

The +jets events are generated with v2.4.2 [@Alwall:2014hca] at NLO with up to two additional partons, and at LO accuracy with up to four additional partons. The production cross sections for the +jets samples are scaled to the NNLO cross sections obtained using 3.1 [@Li:2012wna]. 