From 2df8b068386e6a5216508a82941d61575cb18d38 Mon Sep 17 00:00:00 2001 From: Louis Portales Date: Mon, 13 May 2024 14:42:27 +0200 Subject: [PATCH] WbWb mass scan cards - 340 GeV & 365 GeV --- .../Whizard/v3.0.3/wzp6_ee_WbWb_ecm340.sin | 141 ++++++++++++++++++ .../Whizard/v3.0.3/wzp6_ee_WbWb_ecm365.sin | 141 ++++++++++++++++++ .../v3.0.3/wzp6_ee_WbWb_had_ecm340.sin | 140 +++++++++++++++++ .../v3.0.3/wzp6_ee_WbWb_had_ecm365.sin | 140 +++++++++++++++++ .../v3.0.3/wzp6_ee_WbWb_lep_ecm340.sin | 141 ++++++++++++++++++ .../v3.0.3/wzp6_ee_WbWb_lep_ecm365.sin | 141 ++++++++++++++++++ .../v3.0.3/wzp6_ee_WbWb_semihad_ecm340.sin | 141 ++++++++++++++++++ .../v3.0.3/wzp6_ee_WbWb_semihad_ecm365.sin | 141 ++++++++++++++++++ 8 files changed, 1126 insertions(+) create mode 100644 FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_ecm340.sin create mode 100644 FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_ecm365.sin create mode 100644 FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_had_ecm340.sin create mode 100644 FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_had_ecm365.sin create mode 100644 FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_lep_ecm340.sin create mode 100644 FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_lep_ecm365.sin create mode 100644 FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_semihad_ecm340.sin create mode 100644 FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_semihad_ecm365.sin diff --git a/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_ecm340.sin b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_ecm340.sin new file mode 100644 index 0000000..aeca7c4 --- /dev/null +++ b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_ecm340.sin @@ -0,0 +1,141 @@ +######################################################################## +# +# Copyright (C) 1999-2019 by +# Wolfgang Kilian +# Thorsten Ohl +# Juergen Reuter +# with contributions from +# cf. main AUTHORS file +# +# WHIZARD is free software; you can redistribute it and/or modify it +# under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 2, or (at your option) +# any later version. +# +# WHIZARD is distributed in the hope that it will be useful, but +# WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program; if not, write to the Free Software +# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. +# +######################################################################## +# +# e+ e- -> Z*/gamma* -> t tbar at 365 GeV, t+ -> W+ b, t- -> W- sbar. Polarization kept, CKM enabled +# +# Whizard version is 3.1.0, as in key4hep-stack/2023-04-08 +# +######################################################################## +# +# Card provided by M. Vos - marcel.vos@ific.uv.es +# +######################################################################## + + +model = SM_CKM + +# Center of mass energy +sqrts = 340 GeV + +alphas = 1.e-6 +logical ?is_higgs_process = true + +if ?is_higgs_process then + mH = 125 GeV + wH = 4.1 MeV +else + mH = 2000 GeV + wH = 10 GeV +endif + +logical ?is_higgs_exclusive = false +# logical ?is_higgs_exclusive = true + + +mZ = 91.1876 +mW = 80.385 +mb = 0.0 +me = 0.0 +ms = 0.0 +mc = 0.0 +mmu = 0.0 + +wtop = 1.33 +mtop = 172.5 +!! widths +wZ = 2.443 +wW = 2.049000000000E+00 + +alias quark = u:d:s:c:b:U:D:S:C:B +alias lepton = e1:e2:e3:E1:E2:E3 + +# Processes +process proc = e1, E1 => "W+", "W-", b, B { $omega_flags = "-model:constant_width" } + +?resonance_history = true +resonance_on_shell_limit = 16 +resonance_on_shell_turnoff = 2 +resonance_background_factor = 1 + +beams = e1, E1 => gaussian => isr +?keep_beams = true +?keep_remnants = true + +! beams = e1, E1 => gaussian +gaussian_spread1 = 0.221% +gaussian_spread2 = 0.221% + + +?isr_handler = true +$isr_handler_mode = "recoil" +isr_alpha = 0.0072993 +isr_mass = 0.000511 + +real default_Mqq_cut = 10 GeV + +real default_Mll_cut = 10 GeV + +cuts = all M > default_Mqq_cut [ quark, quark ] and + all M > default_Mll_cut [ lepton, lepton ] and + all M < -default_Mll_cut [ incoming e1, e1 ] and + all M < -default_Mll_cut [ incoming E1, E1 ] + +iterations = 1:500000, 10:500000, 1:1000000 + +relative_error_goal = 0.01 + +$integration_method= "vamp2" +$rng_method= "rng_stream" + + +! Parton shower and hadronization +?ps_fsr_active = true +?ps_isr_active = false +?hadronization_active = true +$shower_method = "PYTHIA6" +!?ps_PYTHIA_verbose = true + +# Process with Higgs +if ?is_higgs_process then + $ps_PYTHIA_PYGIVE = "PMAS(25,1)=125.0;PMAS(25,2)=0.0043;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2;BRAT(210)=0.000052;BRAT(211)=0.000012;BRAT(212)=0.000246;BRAT(213)=0.0289;BRAT(214)=0.582;BRAT(219)=0.000218;BRAT(220)=0.0627;BRAT(222)=0.0819;BRAT(223)=0.00227;BRAT(224)=0.00153;BRAT(225)=0.0262;BRAT(226)=0.214;MDME(210,1)=0;MDME(211,1)=0; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +# +# Processes without Higgs +else + $ps_PYTHIA_PYGIVE = "MSTJ(28)=2;PMAS(25,1)=2000.0;PMAS(25,2)=10.0;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +endif +##$ps_PYTHIA_PYGIVE = "MSTJ(28)=0; PMAS(25,1)=125.; PMAS(25,2)=0.4143E-02; MSTJ(41)=2; MSTU(22)=2000; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200; PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;MSTP(71)=1; MSTP(151)=1; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337; MSTJ(22)=4; PARJ(73)=2250; PARJ(74)=2500" + +integrate (proc) +#{ iterations = 10:100000:"gw", 10:200000:"" } + +# n_events should be passed by the EventProducer +n_events = 100000 + +#sample_format = lhef +sample_format = stdhep +$extension_stdhep = "stdhep" + + +simulate (proc) {checkpoint = 100} diff --git a/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_ecm365.sin b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_ecm365.sin new file mode 100644 index 0000000..3565a72 --- /dev/null +++ b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_ecm365.sin @@ -0,0 +1,141 @@ +######################################################################## +# +# Copyright (C) 1999-2019 by +# Wolfgang Kilian +# Thorsten Ohl +# Juergen Reuter +# with contributions from +# cf. main AUTHORS file +# +# WHIZARD is free software; you can redistribute it and/or modify it +# under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 2, or (at your option) +# any later version. +# +# WHIZARD is distributed in the hope that it will be useful, but +# WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program; if not, write to the Free Software +# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. +# +######################################################################## +# +# e+ e- -> Z*/gamma* -> t tbar at 365 GeV, t+ -> W+ b, t- -> W- sbar. Polarization kept, CKM enabled +# +# Whizard version is 3.1.0, as in key4hep-stack/2023-04-08 +# +######################################################################## +# +# Card provided by M. Vos - marcel.vos@ific.uv.es +# +######################################################################## + + +model = SM_CKM + +# Center of mass energy +sqrts = 365 GeV + +alphas = 1.e-6 +logical ?is_higgs_process = true + +if ?is_higgs_process then + mH = 125 GeV + wH = 4.1 MeV +else + mH = 2000 GeV + wH = 10 GeV +endif + +logical ?is_higgs_exclusive = false +# logical ?is_higgs_exclusive = true + + +mZ = 91.1876 +mW = 80.385 +mb = 0.0 +me = 0.0 +ms = 0.0 +mc = 0.0 +mmu = 0.0 + +wtop = 1.33 +mtop = 172.5 +!! widths +wZ = 2.443 +wW = 2.049000000000E+00 + +alias quark = u:d:s:c:b:U:D:S:C:B +alias lepton = e1:e2:e3:E1:E2:E3 + +# Processes +process proc = e1, E1 => "W+", "W-", b, B { $omega_flags = "-model:constant_width" } + +?resonance_history = true +resonance_on_shell_limit = 16 +resonance_on_shell_turnoff = 2 +resonance_background_factor = 1 + +beams = e1, E1 => gaussian => isr +?keep_beams = true +?keep_remnants = true + +! beams = e1, E1 => gaussian +gaussian_spread1 = 0.221% +gaussian_spread2 = 0.221% + + +?isr_handler = true +$isr_handler_mode = "recoil" +isr_alpha = 0.0072993 +isr_mass = 0.000511 + +real default_Mqq_cut = 10 GeV + +real default_Mll_cut = 10 GeV + +cuts = all M > default_Mqq_cut [ quark, quark ] and + all M > default_Mll_cut [ lepton, lepton ] and + all M < -default_Mll_cut [ incoming e1, e1 ] and + all M < -default_Mll_cut [ incoming E1, E1 ] + +iterations = 1:500000, 10:500000, 1:1000000 + +relative_error_goal = 0.01 + +$integration_method= "vamp2" +$rng_method= "rng_stream" + + +! Parton shower and hadronization +?ps_fsr_active = true +?ps_isr_active = false +?hadronization_active = true +$shower_method = "PYTHIA6" +!?ps_PYTHIA_verbose = true + +# Process with Higgs +if ?is_higgs_process then + $ps_PYTHIA_PYGIVE = "PMAS(25,1)=125.0;PMAS(25,2)=0.0043;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2;BRAT(210)=0.000052;BRAT(211)=0.000012;BRAT(212)=0.000246;BRAT(213)=0.0289;BRAT(214)=0.582;BRAT(219)=0.000218;BRAT(220)=0.0627;BRAT(222)=0.0819;BRAT(223)=0.00227;BRAT(224)=0.00153;BRAT(225)=0.0262;BRAT(226)=0.214;MDME(210,1)=0;MDME(211,1)=0; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +# +# Processes without Higgs +else + $ps_PYTHIA_PYGIVE = "MSTJ(28)=2;PMAS(25,1)=2000.0;PMAS(25,2)=10.0;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +endif +##$ps_PYTHIA_PYGIVE = "MSTJ(28)=0; PMAS(25,1)=125.; PMAS(25,2)=0.4143E-02; MSTJ(41)=2; MSTU(22)=2000; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200; PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;MSTP(71)=1; MSTP(151)=1; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337; MSTJ(22)=4; PARJ(73)=2250; PARJ(74)=2500" + +integrate (proc) +#{ iterations = 10:100000:"gw", 10:200000:"" } + +# n_events should be passed by the EventProducer +n_events = 100000 + +#sample_format = lhef +sample_format = stdhep +$extension_stdhep = "stdhep" + + +simulate (proc) {checkpoint = 100} diff --git a/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_had_ecm340.sin b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_had_ecm340.sin new file mode 100644 index 0000000..742d15c --- /dev/null +++ b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_had_ecm340.sin @@ -0,0 +1,140 @@ +######################################################################## +# +# Copyright (C) 1999-2019 by +# Wolfgang Kilian +# Thorsten Ohl +# Juergen Reuter +# with contributions from +# cf. main AUTHORS file +# +# WHIZARD is free software; you can redistribute it and/or modify it +# under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 2, or (at your option) +# any later version. +# +# WHIZARD is distributed in the hope that it will be useful, but +# WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program; if not, write to the Free Software +# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. +# +######################################################################## +# +# e+ e- -> Z*/gamma* -> t tbar at 365 GeV, t+ -> W+ b, t- -> W- sbar. Polarization kept, CKM enabled +# +# Whizard version is 3.1.0, as in key4hep-stack/2023-04-08 +# +######################################################################## +# +# Card provided by M. Vos - marcel.vos@ific.uv.es +# +######################################################################## + + +model = SM_CKM + +# Center of mass energy +sqrts = 340 GeV + +alphas = 1.e-6 +logical ?is_higgs_process = true + +if ?is_higgs_process then + mH = 125 GeV + wH = 4.1 MeV +else + mH = 2000 GeV + wH = 10 GeV +endif + +logical ?is_higgs_exclusive = false +# logical ?is_higgs_exclusive = true + + +mZ = 91.1876 +mW = 80.385 +mb = 0.0 +me = 0.0 +ms = 0.0 +mc = 0.0 +mmu = 0.0 + +wtop = 1.33 +mtop = 172.5 +!! widths +wZ = 2.443 +wW = 2.049000000000E+00 + +alias quark = u:d:s:c:b:U:D:S:C:B +alias lepton = e1:e2:e3:E1:E2:E3 + +# Processes +process proc = e1, E1 => "W+", "W-", b, B { $omega_flags = "-model:constant_width" } + +?resonance_history = true +resonance_on_shell_limit = 16 +resonance_on_shell_turnoff = 2 +resonance_background_factor = 1 + +beams = e1, E1 => gaussian => isr +?keep_beams = true +?keep_remnants = true + +! beams = e1, E1 => gaussian +gaussian_spread1 = 0.221% +gaussian_spread2 = 0.221% + + +?isr_handler = true +$isr_handler_mode = "recoil" +isr_alpha = 0.0072993 +isr_mass = 0.000511 + +real default_Mqq_cut = 10 GeV + +real default_Mll_cut = 10 GeV + +cuts = all M > default_Mqq_cut [ quark, quark ] and + all M > default_Mll_cut [ lepton, lepton ] and + all M < -default_Mll_cut [ incoming e1, e1 ] and + all M < -default_Mll_cut [ incoming E1, E1 ] + +iterations = 1:500000, 10:500000, 1:1000000 +relative_error_goal = 0.01 + +$integration_method= "vamp2" +$rng_method= "rng_stream" + + +! Parton shower and hadronization +?ps_fsr_active = true +?ps_isr_active = false +?hadronization_active = true +$shower_method = "PYTHIA6" +!?ps_PYTHIA_verbose = true + +# Process with Higgs +if ?is_higgs_process then + $ps_PYTHIA_PYGIVE = "MDME(190,1)=1;MDME(191,1)=1;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=1;MDME(195,1)=1;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=1;MDME(199,1)=1;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=0;MDME(207,1)=0;MDME(208,1)=0;MDME(209,1)=0;PMAS(25,1)=125.0;PMAS(25,2)=0.0043;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2;BRAT(210)=0.000052;BRAT(211)=0.000012;BRAT(212)=0.000246;BRAT(213)=0.0289;BRAT(214)=0.582;BRAT(219)=0.000218;BRAT(220)=0.0627;BRAT(222)=0.0819;BRAT(223)=0.00227;BRAT(224)=0.00153;BRAT(225)=0.0262;BRAT(226)=0.214;MDME(210,1)=0;MDME(211,1)=0; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +# +# Processes without Higgs +else + $ps_PYTHIA_PYGIVE = "MDME(190,1)=1;MDME(191,1)=1;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=1;MDME(195,1)=1;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=1;MDME(199,1)=1;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=0;MDME(207,1)=0;MDME(208,1)=0;MDME(209,1)=0;MSTJ(28)=2;PMAS(25,1)=2000.0;PMAS(25,2)=10.0;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +endif +##$ps_PYTHIA_PYGIVE = "MDME(190,1)=1;MDME(191,1)=1;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=1;MDME(195,1)=1;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=1;MDME(199,1)=1;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=0;MDME(207,1)=0;MDME(208,1)=0;MDME(209,1)=0;MSTJ(28)=0; PMAS(25,1)=125.; PMAS(25,2)=0.4143E-02; MSTJ(41)=2; MSTU(22)=2000; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200; PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;MSTP(71)=1; MSTP(151)=1; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337; MSTJ(22)=4; PARJ(73)=2250; PARJ(74)=2500" + +integrate (proc) +#{ iterations = 10:100000:"gw", 10:200000:"" } + +# n_events should be passed by the EventProducer +n_events = 100000 + +#sample_format = lhef +sample_format = stdhep +$extension_stdhep = "stdhep" + + +simulate (proc) {checkpoint = 100} diff --git a/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_had_ecm365.sin b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_had_ecm365.sin new file mode 100644 index 0000000..be697ce --- /dev/null +++ b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_had_ecm365.sin @@ -0,0 +1,140 @@ +######################################################################## +# +# Copyright (C) 1999-2019 by +# Wolfgang Kilian +# Thorsten Ohl +# Juergen Reuter +# with contributions from +# cf. main AUTHORS file +# +# WHIZARD is free software; you can redistribute it and/or modify it +# under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 2, or (at your option) +# any later version. +# +# WHIZARD is distributed in the hope that it will be useful, but +# WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program; if not, write to the Free Software +# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. +# +######################################################################## +# +# e+ e- -> Z*/gamma* -> t tbar at 365 GeV, t+ -> W+ b, t- -> W- sbar. Polarization kept, CKM enabled +# +# Whizard version is 3.1.0, as in key4hep-stack/2023-04-08 +# +######################################################################## +# +# Card provided by M. Vos - marcel.vos@ific.uv.es +# +######################################################################## + + +model = SM_CKM + +# Center of mass energy +sqrts = 365 GeV + +alphas = 1.e-6 +logical ?is_higgs_process = true + +if ?is_higgs_process then + mH = 125 GeV + wH = 4.1 MeV +else + mH = 2000 GeV + wH = 10 GeV +endif + +logical ?is_higgs_exclusive = false +# logical ?is_higgs_exclusive = true + + +mZ = 91.1876 +mW = 80.385 +mb = 0.0 +me = 0.0 +ms = 0.0 +mc = 0.0 +mmu = 0.0 + +wtop = 1.33 +mtop = 172.5 +!! widths +wZ = 2.443 +wW = 2.049000000000E+00 + +alias quark = u:d:s:c:b:U:D:S:C:B +alias lepton = e1:e2:e3:E1:E2:E3 + +# Processes +process proc = e1, E1 => "W+", "W-", b, B { $omega_flags = "-model:constant_width" } + +?resonance_history = true +resonance_on_shell_limit = 16 +resonance_on_shell_turnoff = 2 +resonance_background_factor = 1 + +beams = e1, E1 => gaussian => isr +?keep_beams = true +?keep_remnants = true + +! beams = e1, E1 => gaussian +gaussian_spread1 = 0.221% +gaussian_spread2 = 0.221% + + +?isr_handler = true +$isr_handler_mode = "recoil" +isr_alpha = 0.0072993 +isr_mass = 0.000511 + +real default_Mqq_cut = 10 GeV + +real default_Mll_cut = 10 GeV + +cuts = all M > default_Mqq_cut [ quark, quark ] and + all M > default_Mll_cut [ lepton, lepton ] and + all M < -default_Mll_cut [ incoming e1, e1 ] and + all M < -default_Mll_cut [ incoming E1, E1 ] + +iterations = 1:500000, 10:500000, 1:1000000 +relative_error_goal = 0.01 + +$integration_method= "vamp2" +$rng_method= "rng_stream" + + +! Parton shower and hadronization +?ps_fsr_active = true +?ps_isr_active = false +?hadronization_active = true +$shower_method = "PYTHIA6" +!?ps_PYTHIA_verbose = true + +# Process with Higgs +if ?is_higgs_process then + $ps_PYTHIA_PYGIVE = "MDME(190,1)=1;MDME(191,1)=1;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=1;MDME(195,1)=1;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=1;MDME(199,1)=1;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=0;MDME(207,1)=0;MDME(208,1)=0;MDME(209,1)=0;PMAS(25,1)=125.0;PMAS(25,2)=0.0043;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2;BRAT(210)=0.000052;BRAT(211)=0.000012;BRAT(212)=0.000246;BRAT(213)=0.0289;BRAT(214)=0.582;BRAT(219)=0.000218;BRAT(220)=0.0627;BRAT(222)=0.0819;BRAT(223)=0.00227;BRAT(224)=0.00153;BRAT(225)=0.0262;BRAT(226)=0.214;MDME(210,1)=0;MDME(211,1)=0; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +# +# Processes without Higgs +else + $ps_PYTHIA_PYGIVE = "MDME(190,1)=1;MDME(191,1)=1;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=1;MDME(195,1)=1;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=1;MDME(199,1)=1;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=0;MDME(207,1)=0;MDME(208,1)=0;MDME(209,1)=0;MSTJ(28)=2;PMAS(25,1)=2000.0;PMAS(25,2)=10.0;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +endif +##$ps_PYTHIA_PYGIVE = "MDME(190,1)=1;MDME(191,1)=1;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=1;MDME(195,1)=1;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=1;MDME(199,1)=1;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=0;MDME(207,1)=0;MDME(208,1)=0;MDME(209,1)=0;MSTJ(28)=0; PMAS(25,1)=125.; PMAS(25,2)=0.4143E-02; MSTJ(41)=2; MSTU(22)=2000; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200; PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;MSTP(71)=1; MSTP(151)=1; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337; MSTJ(22)=4; PARJ(73)=2250; PARJ(74)=2500" + +integrate (proc) +#{ iterations = 10:100000:"gw", 10:200000:"" } + +# n_events should be passed by the EventProducer +n_events = 100000 + +#sample_format = lhef +sample_format = stdhep +$extension_stdhep = "stdhep" + + +simulate (proc) {checkpoint = 100} diff --git a/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_lep_ecm340.sin b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_lep_ecm340.sin new file mode 100644 index 0000000..7ac9484 --- /dev/null +++ b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_lep_ecm340.sin @@ -0,0 +1,141 @@ +######################################################################## +# +# Copyright (C) 1999-2019 by +# Wolfgang Kilian +# Thorsten Ohl +# Juergen Reuter +# with contributions from +# cf. main AUTHORS file +# +# WHIZARD is free software; you can redistribute it and/or modify it +# under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 2, or (at your option) +# any later version. +# +# WHIZARD is distributed in the hope that it will be useful, but +# WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program; if not, write to the Free Software +# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. +# +######################################################################## +# +# e+ e- -> Z*/gamma* -> t tbar at 365 GeV, t+ -> W+ b, t- -> W- sbar. Polarization kept, CKM enabled +# +# Whizard version is 3.1.0, as in key4hep-stack/2023-04-08 +# +######################################################################## +# +# Card provided by M. Vos - marcel.vos@ific.uv.es +# +######################################################################## + + +model = SM_CKM + +# Center of mass energy +sqrts = 340 GeV + +alphas = 1.e-6 +logical ?is_higgs_process = true + +if ?is_higgs_process then + mH = 125 GeV + wH = 4.1 MeV +else + mH = 2000 GeV + wH = 10 GeV +endif + +logical ?is_higgs_exclusive = false +# logical ?is_higgs_exclusive = true + + +mZ = 91.1876 +mW = 80.385 +mb = 0.0 +me = 0.0 +ms = 0.0 +mc = 0.0 +mmu = 0.0 + +wtop = 1.33 +mtop = 172.5 +!! widths +wZ = 2.443 +wW = 2.049000000000E+00 + +alias quark = u:d:s:c:b:U:D:S:C:B +alias lepton = e1:e2:e3:E1:E2:E3 + +# Processes +process proc = e1, E1 => "W+", "W-", b, B { $omega_flags = "-model:constant_width" } + +?resonance_history = true +resonance_on_shell_limit = 16 +resonance_on_shell_turnoff = 2 +resonance_background_factor = 1 + +beams = e1, E1 => gaussian => isr +?keep_beams = true +?keep_remnants = true + +! beams = e1, E1 => gaussian +gaussian_spread1 = 0.221% +gaussian_spread2 = 0.221% + + +?isr_handler = true +$isr_handler_mode = "recoil" +isr_alpha = 0.0072993 +isr_mass = 0.000511 + +real default_Mqq_cut = 10 GeV + +real default_Mll_cut = 10 GeV + +cuts = all M > default_Mqq_cut [ quark, quark ] and + all M > default_Mll_cut [ lepton, lepton ] and + all M < -default_Mll_cut [ incoming e1, e1 ] and + all M < -default_Mll_cut [ incoming E1, E1 ] + +iterations = 1:500000, 10:500000, 1:1000000 + +relative_error_goal = 0.01 + +$integration_method= "vamp2" +$rng_method= "rng_stream" + + +! Parton shower and hadronization +?ps_fsr_active = true +?ps_isr_active = false +?hadronization_active = true +$shower_method = "PYTHIA6" +!?ps_PYTHIA_verbose = true + +# Process with Higgs +if ?is_higgs_process then + $ps_PYTHIA_PYGIVE = "MDME(190,1)=0;MDME(191,1)=0;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=0;MDME(195,1)=0;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=0;MDME(199,1)=0;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=1;MDME(207,1)=1;MDME(208,1)=1;MDME(209,1)=0;PMAS(25,1)=125.0;PMAS(25,2)=0.0043;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2;BRAT(210)=0.000052;BRAT(211)=0.000012;BRAT(212)=0.000246;BRAT(213)=0.0289;BRAT(214)=0.582;BRAT(219)=0.000218;BRAT(220)=0.0627;BRAT(222)=0.0819;BRAT(223)=0.00227;BRAT(224)=0.00153;BRAT(225)=0.0262;BRAT(226)=0.214;MDME(210,1)=0;MDME(211,1)=0; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +# +# Processes without Higgs +else + $ps_PYTHIA_PYGIVE = "MDME(190,1)=0;MDME(191,1)=0;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=0;MDME(195,1)=0;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=0;MDME(199,1)=0;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=1;MDME(207,1)=1;MDME(208,1)=1;MDME(209,1)=0;MSTJ(28)=2;PMAS(25,1)=2000.0;PMAS(25,2)=10.0;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +endif +##$ps_PYTHIA_PYGIVE = "MDME(190,1)=0;MDME(191,1)=0;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=0;MDME(195,1)=0;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=0;MDME(199,1)=0;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=1;MDME(207,1)=1;MDME(208,1)=1;MDME(209,1)=0;MSTJ(28)=0; PMAS(25,1)=125.; PMAS(25,2)=0.4143E-02; MSTJ(41)=2; MSTU(22)=2000; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200; PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;MSTP(71)=1; MSTP(151)=1; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337; MSTJ(22)=4; PARJ(73)=2250; PARJ(74)=2500" + +integrate (proc) +#{ iterations = 10:100000:"gw", 10:200000:"" } + +# n_events should be passed by the EventProducer +n_events = 100000 + +#sample_format = lhef +sample_format = stdhep +$extension_stdhep = "stdhep" + + +simulate (proc) {checkpoint = 100} diff --git a/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_lep_ecm365.sin b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_lep_ecm365.sin new file mode 100644 index 0000000..6ffe6bd --- /dev/null +++ b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_lep_ecm365.sin @@ -0,0 +1,141 @@ +######################################################################## +# +# Copyright (C) 1999-2019 by +# Wolfgang Kilian +# Thorsten Ohl +# Juergen Reuter +# with contributions from +# cf. main AUTHORS file +# +# WHIZARD is free software; you can redistribute it and/or modify it +# under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 2, or (at your option) +# any later version. +# +# WHIZARD is distributed in the hope that it will be useful, but +# WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program; if not, write to the Free Software +# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. +# +######################################################################## +# +# e+ e- -> Z*/gamma* -> t tbar at 365 GeV, t+ -> W+ b, t- -> W- sbar. Polarization kept, CKM enabled +# +# Whizard version is 3.1.0, as in key4hep-stack/2023-04-08 +# +######################################################################## +# +# Card provided by M. Vos - marcel.vos@ific.uv.es +# +######################################################################## + + +model = SM_CKM + +# Center of mass energy +sqrts = 365 GeV + +alphas = 1.e-6 +logical ?is_higgs_process = true + +if ?is_higgs_process then + mH = 125 GeV + wH = 4.1 MeV +else + mH = 2000 GeV + wH = 10 GeV +endif + +logical ?is_higgs_exclusive = false +# logical ?is_higgs_exclusive = true + + +mZ = 91.1876 +mW = 80.385 +mb = 0.0 +me = 0.0 +ms = 0.0 +mc = 0.0 +mmu = 0.0 + +wtop = 1.33 +mtop = 172.5 +!! widths +wZ = 2.443 +wW = 2.049000000000E+00 + +alias quark = u:d:s:c:b:U:D:S:C:B +alias lepton = e1:e2:e3:E1:E2:E3 + +# Processes +process proc = e1, E1 => "W+", "W-", b, B { $omega_flags = "-model:constant_width" } + +?resonance_history = true +resonance_on_shell_limit = 16 +resonance_on_shell_turnoff = 2 +resonance_background_factor = 1 + +beams = e1, E1 => gaussian => isr +?keep_beams = true +?keep_remnants = true + +! beams = e1, E1 => gaussian +gaussian_spread1 = 0.221% +gaussian_spread2 = 0.221% + + +?isr_handler = true +$isr_handler_mode = "recoil" +isr_alpha = 0.0072993 +isr_mass = 0.000511 + +real default_Mqq_cut = 10 GeV + +real default_Mll_cut = 10 GeV + +cuts = all M > default_Mqq_cut [ quark, quark ] and + all M > default_Mll_cut [ lepton, lepton ] and + all M < -default_Mll_cut [ incoming e1, e1 ] and + all M < -default_Mll_cut [ incoming E1, E1 ] + +iterations = 1:500000, 10:500000, 1:1000000 + +relative_error_goal = 0.01 + +$integration_method= "vamp2" +$rng_method= "rng_stream" + + +! Parton shower and hadronization +?ps_fsr_active = true +?ps_isr_active = false +?hadronization_active = true +$shower_method = "PYTHIA6" +!?ps_PYTHIA_verbose = true + +# Process with Higgs +if ?is_higgs_process then + $ps_PYTHIA_PYGIVE = "MDME(190,1)=0;MDME(191,1)=0;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=0;MDME(195,1)=0;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=0;MDME(199,1)=0;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=1;MDME(207,1)=1;MDME(208,1)=1;MDME(209,1)=0;PMAS(25,1)=125.0;PMAS(25,2)=0.0043;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2;BRAT(210)=0.000052;BRAT(211)=0.000012;BRAT(212)=0.000246;BRAT(213)=0.0289;BRAT(214)=0.582;BRAT(219)=0.000218;BRAT(220)=0.0627;BRAT(222)=0.0819;BRAT(223)=0.00227;BRAT(224)=0.00153;BRAT(225)=0.0262;BRAT(226)=0.214;MDME(210,1)=0;MDME(211,1)=0; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +# +# Processes without Higgs +else + $ps_PYTHIA_PYGIVE = "MDME(190,1)=0;MDME(191,1)=0;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=0;MDME(195,1)=0;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=0;MDME(199,1)=0;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=1;MDME(207,1)=1;MDME(208,1)=1;MDME(209,1)=0;MSTJ(28)=2;PMAS(25,1)=2000.0;PMAS(25,2)=10.0;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +endif +##$ps_PYTHIA_PYGIVE = "MDME(190,1)=0;MDME(191,1)=0;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=0;MDME(195,1)=0;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=0;MDME(199,1)=0;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=1;MDME(207,1)=1;MDME(208,1)=1;MDME(209,1)=0;MSTJ(28)=0; PMAS(25,1)=125.; PMAS(25,2)=0.4143E-02; MSTJ(41)=2; MSTU(22)=2000; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200; PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;MSTP(71)=1; MSTP(151)=1; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337; MSTJ(22)=4; PARJ(73)=2250; PARJ(74)=2500" + +integrate (proc) +#{ iterations = 10:100000:"gw", 10:200000:"" } + +# n_events should be passed by the EventProducer +n_events = 100000 + +#sample_format = lhef +sample_format = stdhep +$extension_stdhep = "stdhep" + + +simulate (proc) {checkpoint = 100} diff --git a/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_semihad_ecm340.sin b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_semihad_ecm340.sin new file mode 100644 index 0000000..ccfb771 --- /dev/null +++ b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_semihad_ecm340.sin @@ -0,0 +1,141 @@ +######################################################################## +# +# Copyright (C) 1999-2019 by +# Wolfgang Kilian +# Thorsten Ohl +# Juergen Reuter +# with contributions from +# cf. main AUTHORS file +# +# WHIZARD is free software; you can redistribute it and/or modify it +# under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 2, or (at your option) +# any later version. +# +# WHIZARD is distributed in the hope that it will be useful, but +# WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program; if not, write to the Free Software +# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. +# +######################################################################## +# +# e+ e- -> Z*/gamma* -> t tbar at 365 GeV, t+ -> W+ b, t- -> W- sbar. Polarization kept, CKM enabled +# +# Whizard version is 3.1.0, as in key4hep-stack/2023-04-08 +# +######################################################################## +# +# Card provided by M. Vos - marcel.vos@ific.uv.es +# +######################################################################## + + +model = SM_CKM + +# Center of mass energy +sqrts = 340 GeV + +alphas = 1.e-6 +logical ?is_higgs_process = true + +if ?is_higgs_process then + mH = 125 GeV + wH = 4.1 MeV +else + mH = 2000 GeV + wH = 10 GeV +endif + +logical ?is_higgs_exclusive = false +# logical ?is_higgs_exclusive = true + + +mZ = 91.1876 +mW = 80.385 +mb = 0.0 +me = 0.0 +ms = 0.0 +mc = 0.0 +mmu = 0.0 + +wtop = 1.33 +mtop = 172.5 +!! widths +wZ = 2.443 +wW = 2.049000000000E+00 + +alias quark = u:d:s:c:b:U:D:S:C:B +alias lepton = e1:e2:e3:E1:E2:E3 + +# Processes +process proc = e1, E1 => "W+", "W-", b, B { $omega_flags = "-model:constant_width" } + +?resonance_history = true +resonance_on_shell_limit = 16 +resonance_on_shell_turnoff = 2 +resonance_background_factor = 1 + +beams = e1, E1 => gaussian => isr +?keep_beams = true +?keep_remnants = true + +! beams = e1, E1 => gaussian +gaussian_spread1 = 0.221% +gaussian_spread2 = 0.221% + + +?isr_handler = true +$isr_handler_mode = "recoil" +isr_alpha = 0.0072993 +isr_mass = 0.000511 + +real default_Mqq_cut = 10 GeV + +real default_Mll_cut = 10 GeV + +cuts = all M > default_Mqq_cut [ quark, quark ] and + all M > default_Mll_cut [ lepton, lepton ] and + all M < -default_Mll_cut [ incoming e1, e1 ] and + all M < -default_Mll_cut [ incoming E1, E1 ] + +iterations = 1:500000, 10:500000, 1:1000000 + +relative_error_goal = 0.01 + +$integration_method= "vamp2" +$rng_method= "rng_stream" + + +! Parton shower and hadronization +?ps_fsr_active = true +?ps_isr_active = false +?hadronization_active = true +$shower_method = "PYTHIA6" +!?ps_PYTHIA_verbose = true + +# Process with Higgs +if ?is_higgs_process then + $ps_PYTHIA_PYGIVE = "MDME(190,1)=5;MDME(191,1)=5;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=5;MDME(195,1)=5;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=5;MDME(199,1)=5;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=4;MDME(207,1)=4;MDME(208,1)=4;MDME(209,1)=0;PMAS(25,1)=125.0;PMAS(25,2)=0.0043;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2;BRAT(210)=0.000052;BRAT(211)=0.000012;BRAT(212)=0.000246;BRAT(213)=0.0289;BRAT(214)=0.582;BRAT(219)=0.000218;BRAT(220)=0.0627;BRAT(222)=0.0819;BRAT(223)=0.00227;BRAT(224)=0.00153;BRAT(225)=0.0262;BRAT(226)=0.214;MDME(210,1)=0;MDME(211,1)=0; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +# +# Processes without Higgs +else + $ps_PYTHIA_PYGIVE = "MDME(190,1)=5;MDME(191,1)=5;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=5;MDME(195,1)=5;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=5;MDME(199,1)=5;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=4;MDME(207,1)=4;MDME(208,1)=4;MDME(209,1)=0;MSTJ(28)=2;PMAS(25,1)=2000.0;PMAS(25,2)=10.0;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +endif +##$ps_PYTHIA_PYGIVE = "MDME(190,1)=5;MDME(191,1)=5;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=5;MDME(195,1)=5;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=5;MDME(199,1)=5;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=4;MDME(207,1)=4;MDME(208,1)=4;MDME(209,1)=0;MSTJ(28)=0; PMAS(25,1)=125.; PMAS(25,2)=0.4143E-02; MSTJ(41)=2; MSTU(22)=2000; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200; PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;MSTP(71)=1; MSTP(151)=1; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337; MSTJ(22)=4; PARJ(73)=2250; PARJ(74)=2500" + +integrate (proc) +#{ iterations = 10:100000:"gw", 10:200000:"" } + +# n_events should be passed by the EventProducer +n_events = 100000 + +sample_format = lhef +sample_format = stdhep +$extension_stdhep = "stdhep" + + +simulate (proc) {checkpoint = 100} diff --git a/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_semihad_ecm365.sin b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_semihad_ecm365.sin new file mode 100644 index 0000000..aae1317 --- /dev/null +++ b/FCCee/Generator/Whizard/v3.0.3/wzp6_ee_WbWb_semihad_ecm365.sin @@ -0,0 +1,141 @@ +######################################################################## +# +# Copyright (C) 1999-2019 by +# Wolfgang Kilian +# Thorsten Ohl +# Juergen Reuter +# with contributions from +# cf. main AUTHORS file +# +# WHIZARD is free software; you can redistribute it and/or modify it +# under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 2, or (at your option) +# any later version. +# +# WHIZARD is distributed in the hope that it will be useful, but +# WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program; if not, write to the Free Software +# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. +# +######################################################################## +# +# e+ e- -> Z*/gamma* -> t tbar at 365 GeV, t+ -> W+ b, t- -> W- sbar. Polarization kept, CKM enabled +# +# Whizard version is 3.1.0, as in key4hep-stack/2023-04-08 +# +######################################################################## +# +# Card provided by M. Vos - marcel.vos@ific.uv.es +# +######################################################################## + + +model = SM_CKM + +# Center of mass energy +sqrts = 365 GeV + +alphas = 1.e-6 +logical ?is_higgs_process = true + +if ?is_higgs_process then + mH = 125 GeV + wH = 4.1 MeV +else + mH = 2000 GeV + wH = 10 GeV +endif + +logical ?is_higgs_exclusive = false +# logical ?is_higgs_exclusive = true + + +mZ = 91.1876 +mW = 80.385 +mb = 0.0 +me = 0.0 +ms = 0.0 +mc = 0.0 +mmu = 0.0 + +wtop = 1.33 +mtop = 172.5 +!! widths +wZ = 2.443 +wW = 2.049000000000E+00 + +alias quark = u:d:s:c:b:U:D:S:C:B +alias lepton = e1:e2:e3:E1:E2:E3 + +# Processes +process proc = e1, E1 => "W+", "W-", b, B { $omega_flags = "-model:constant_width" } + +?resonance_history = true +resonance_on_shell_limit = 16 +resonance_on_shell_turnoff = 2 +resonance_background_factor = 1 + +beams = e1, E1 => gaussian => isr +?keep_beams = true +?keep_remnants = true + +! beams = e1, E1 => gaussian +gaussian_spread1 = 0.221% +gaussian_spread2 = 0.221% + + +?isr_handler = true +$isr_handler_mode = "recoil" +isr_alpha = 0.0072993 +isr_mass = 0.000511 + +real default_Mqq_cut = 10 GeV + +real default_Mll_cut = 10 GeV + +cuts = all M > default_Mqq_cut [ quark, quark ] and + all M > default_Mll_cut [ lepton, lepton ] and + all M < -default_Mll_cut [ incoming e1, e1 ] and + all M < -default_Mll_cut [ incoming E1, E1 ] + +iterations = 1:500000, 10:500000, 1:1000000 + +relative_error_goal = 0.01 + +$integration_method= "vamp2" +$rng_method= "rng_stream" + + +! Parton shower and hadronization +?ps_fsr_active = true +?ps_isr_active = false +?hadronization_active = true +$shower_method = "PYTHIA6" +!?ps_PYTHIA_verbose = true + +# Process with Higgs +if ?is_higgs_process then + $ps_PYTHIA_PYGIVE = "MDME(190,1)=5;MDME(191,1)=5;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=5;MDME(195,1)=5;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=5;MDME(199,1)=5;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=4;MDME(207,1)=4;MDME(208,1)=4;MDME(209,1)=0;PMAS(25,1)=125.0;PMAS(25,2)=0.0043;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2;BRAT(210)=0.000052;BRAT(211)=0.000012;BRAT(212)=0.000246;BRAT(213)=0.0289;BRAT(214)=0.582;BRAT(219)=0.000218;BRAT(220)=0.0627;BRAT(222)=0.0819;BRAT(223)=0.00227;BRAT(224)=0.00153;BRAT(225)=0.0262;BRAT(226)=0.214;MDME(210,1)=0;MDME(211,1)=0; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +# +# Processes without Higgs +else + $ps_PYTHIA_PYGIVE = "MDME(190,1)=5;MDME(191,1)=5;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=5;MDME(195,1)=5;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=5;MDME(199,1)=5;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=4;MDME(207,1)=4;MDME(208,1)=4;MDME(209,1)=0;MSTJ(28)=2;PMAS(25,1)=2000.0;PMAS(25,2)=10.0;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1; MSTP(125)=2; MWID(25)=2; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337" +endif +##$ps_PYTHIA_PYGIVE = "MDME(190,1)=5;MDME(191,1)=5;MDME(192,1)=0;MDME(193,1)=0;MDME(194,1)=5;MDME(195,1)=5;MDME(196,1)=0;MDME(197,1)=0;MDME(198,1)=5;MDME(199,1)=5;MDME(200,1)=0;MDME(201,1)=0;MDME(202,1)=0;MDME(203,1)=0;MDME(204,1)=0;MDME(205,1)=0;MDME(206,1)=4;MDME(207,1)=4;MDME(208,1)=4;MDME(209,1)=0;MSTJ(28)=0; PMAS(25,1)=125.; PMAS(25,2)=0.4143E-02; MSTJ(41)=2; MSTU(22)=2000; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000; PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200; PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000; PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000; PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;MSTP(71)=1; MSTP(151)=1; PARP(151)=0.0273; PARP(152)=4.88e-5; PARP(153)=1.33; PARP(154)=1.337; MSTJ(22)=4; PARJ(73)=2250; PARJ(74)=2500" + +integrate (proc) +#{ iterations = 10:100000:"gw", 10:200000:"" } + +# n_events should be passed by the EventProducer +n_events = 100000 + +sample_format = lhef +sample_format = stdhep +$extension_stdhep = "stdhep" + + +simulate (proc) {checkpoint = 100}