Skip to content

Commit

Permalink
Merge pull request #163 from portalesHEP/winter23_BES365
Browse files Browse the repository at this point in the history
adding llH ecm365 cards w/ 10% BES variation; fixing llH ecm365 cards w/ 1% BES variation
  • Loading branch information
portalesHEP authored May 24, 2024
2 parents f81bed7 + a6cb8b8 commit f07287a
Show file tree
Hide file tree
Showing 8 changed files with 244 additions and 8 deletions.
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
########################################################################
#
# Whizard 3.0.3
#
# e e -> mu mu H (only one diagram, e e -> Z H and Z -> mumu )
#
# Hadronisation by Pythia6. Inclusive H decays (by Pythia)
#
########################################################################
#

model = SM

mH = 125 GeV

# Center of mass energy
sqrts = 365 GeV

# Processes

#?vis_diags = true


process proc = e1, E1 => H, e1, E1

beams = e1, E1 => gaussian => isr
?keep_beams = true # do not use this option, makes Pythia crash
?keep_remnants = true

gaussian_spread1 = 0.2431%
gaussian_spread2 = 0.2431%


?isr_handler = true
$isr_handler_mode = "recoil"
isr_alpha = 0.0072993
isr_mass = 0.000511



! Parton shower and hadronization
?ps_fsr_active = true
?ps_isr_active = false
?hadronization_active = true
$shower_method = "PYTHIA6"
!?ps_PYTHIA_verbose = true


$ps_PYTHIA_PYGIVE = "MSTJ(28)=0; PMAS(25,1)=125.1; 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", 5:200000:"" }

# n_events should be passed by the EventProducer
# n_events = 100000

sample_format = stdhep
$extension_stdhep = "stdhep"
simulate (proc) {checkpoint = 100}

Original file line number Diff line number Diff line change
Expand Up @@ -27,8 +27,8 @@ beams = e1, E1 => gaussian => isr
?keep_beams = true # do not use this option, makes Pythia crash
?keep_remnants = true

gaussian_spread1 = 0.18685%
gaussian_spread2 = 0.18685%
gaussian_spread1 = 0.2232%
gaussian_spread2 = 0.2232%


?isr_handler = true
Expand Down
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
########################################################################
#
# Whizard 3.0.3
#
# e e -> mu mu H (only one diagram, e e -> Z H and Z -> mumu )
#
# Hadronisation by Pythia6. Inclusive H decays (by Pythia)
#
########################################################################
#

model = SM

mH = 125 GeV

# Center of mass energy
sqrts = 365 GeV

# Processes

#?vis_diags = true


process proc = e1, E1 => H, e1, E1

beams = e1, E1 => gaussian => isr
?keep_beams = true # do not use this option, makes Pythia crash
?keep_remnants = true

gaussian_spread1 = 0.1989%
gaussian_spread2 = 0.1989%


?isr_handler = true
$isr_handler_mode = "recoil"
isr_alpha = 0.0072993
isr_mass = 0.000511



! Parton shower and hadronization
?ps_fsr_active = true
?ps_isr_active = false
?hadronization_active = true
$shower_method = "PYTHIA6"
!?ps_PYTHIA_verbose = true


$ps_PYTHIA_PYGIVE = "MSTJ(28)=0; PMAS(25,1)=125.1; 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", 5:200000:"" }

# n_events should be passed by the EventProducer
# n_events = 100000

sample_format = stdhep
$extension_stdhep = "stdhep"
simulate (proc) {checkpoint = 100}

Original file line number Diff line number Diff line change
Expand Up @@ -27,8 +27,8 @@ beams = e1, E1 => gaussian => isr
?keep_beams = true # do not use this option, makes Pythia crash
?keep_remnants = true

gaussian_spread1 = 0.18315%
gaussian_spread2 = 0.18315%
gaussian_spread1 = 0.2188%
gaussian_spread2 = 0.2188%


?isr_handler = true
Expand Down
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
########################################################################
#
# Whizard 3.0.3
#
# e e -> mu mu H (only one diagram, e e -> Z H and Z -> mumu )
#
# Hadronisation by Pythia6. Inclusive H decays (by Pythia)
#
########################################################################
#

model = SM

mH = 125 GeV

# Center of mass energy
sqrts = 365 GeV

# Processes

#?vis_diags = true


process proc = e1, E1 => H, e2, E2

beams = e1, E1 => gaussian => isr
?keep_beams = true # do not use this option, makes Pythia crash
?keep_remnants = true

gaussian_spread1 = 0.2431%
gaussian_spread2 = 0.2431%


?isr_handler = true
$isr_handler_mode = "recoil"
isr_alpha = 0.0072993
isr_mass = 0.000511



! Parton shower and hadronization
?ps_fsr_active = true
?ps_isr_active = false
?hadronization_active = true
$shower_method = "PYTHIA6"
!?ps_PYTHIA_verbose = true


$ps_PYTHIA_PYGIVE = "MSTJ(28)=0; PMAS(25,1)=125.1; 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", 5:200000:"" }

# n_events should be passed by the EventProducer
# n_events = 100000

sample_format = stdhep
$extension_stdhep = "stdhep"
simulate (proc) {checkpoint = 100}

Original file line number Diff line number Diff line change
Expand Up @@ -27,8 +27,8 @@ beams = e1, E1 => gaussian => isr
?keep_beams = true # do not use this option, makes Pythia crash
?keep_remnants = true

gaussian_spread1 = 0.18685%
gaussian_spread2 = 0.18685%
gaussian_spread1 = 0.2232%
gaussian_spread2 = 0.2232%


?isr_handler = true
Expand Down
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
########################################################################
#
# Whizard 3.0.3
#
# e e -> mu mu H (only one diagram, e e -> Z H and Z -> mumu )
#
# Hadronisation by Pythia6. Inclusive H decays (by Pythia)
#
########################################################################
#

model = SM

mH = 125 GeV

# Center of mass energy
sqrts = 365 GeV

# Processes

#?vis_diags = true


process proc = e1, E1 => H, e2, E2

beams = e1, E1 => gaussian => isr
?keep_beams = true # do not use this option, makes Pythia crash
?keep_remnants = true

gaussian_spread1 = 0.1989%
gaussian_spread2 = 0.1989%


?isr_handler = true
$isr_handler_mode = "recoil"
isr_alpha = 0.0072993
isr_mass = 0.000511



! Parton shower and hadronization
?ps_fsr_active = true
?ps_isr_active = false
?hadronization_active = true
$shower_method = "PYTHIA6"
!?ps_PYTHIA_verbose = true


$ps_PYTHIA_PYGIVE = "MSTJ(28)=0; PMAS(25,1)=125.1; 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", 5:200000:"" }

# n_events should be passed by the EventProducer
# n_events = 100000

sample_format = stdhep
$extension_stdhep = "stdhep"
simulate (proc) {checkpoint = 100}

Original file line number Diff line number Diff line change
Expand Up @@ -27,8 +27,8 @@ beams = e1, E1 => gaussian => isr
?keep_beams = true # do not use this option, makes Pythia crash
?keep_remnants = true

gaussian_spread1 = 0.18315%
gaussian_spread2 = 0.18315%
gaussian_spread1 = 0.2188%
gaussian_spread2 = 0.2188%


?isr_handler = true
Expand Down

0 comments on commit f07287a

Please sign in to comment.