Skip to content

Samples for Higgs Zgamma effective coupling re-analysis (with Tagger) #13

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
kjvbrt wants to merge 2 commits into
base: master
Choose a base branch
from
Open

Samples for Higgs Zgamma effective coupling re-analysis (with Tagger) #13

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
99eb86e
Samples for Higgs Zgamma effective coupling re-analysis (with Tagger …
kjvbrt May 5, 2025
bcd7676
Enabling the actual tagger
kjvbrt May 5, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
44 changes: 44 additions & 0 deletions ...oductions/spring2024ConfigFiles/ee_ZH/HiggsZgammaEffCouplingTagger/whizardddsimgaudi.conf
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,44 @@
[whizard2]
Version = 2.8.3

[ddsim]
SteeringFile = cld_steer.py
ExtraCLIArguments = --vertexSigma 0.0098 2.54e-5 0.646 0.0063 --physics.rejectPDGs 1,2,3,4,5,6,12,14,16,21,23,24,25

[gaudiapp]
ExecutableName = k4run
SteeringFile = CLDReconstruction_jet_tagger_impl.py
InputFileFlag = --inputFiles
OutputFileFlag = --IOSvc.Output
Arguments = --enableLCFIJet --enableMLJetTagger --cms 240

[Production Parameters]
machine = ee
prodGroup = fcc_ee_spring2024_HiggsZgammaEffCouplingTagger1

softwareVersion = key4hep_nightly
generatorApplication = whizard2
simulationApplication = ddsim
reconstructionApplication = gaudiapp
generatorSteeringFile = wzp6_ee_ZH_ecm240.sin

configVersion = key4hep-devel-6
configPackage = fccConfig
eventsPerJobs = 1000

numberOfTasks = 1

campaign = test_spring2024
energies = 240
processes = ZH
ProdIDs = 16880
detectorModel = CLD_o2_v05
additionalName = HiggsZgammaEffCouplingTagger1
productionLogLevel = VERBOSE
outputSE = CERN-DST-EOS

finalOutputSE = CERN-SRM
MoveStatus = Stopped
MoveGroupSize = 10

ProdTypes = Rec
104 changes: 104 additions & 0 deletions ...roductions/spring2024ConfigFiles/ee_ZH/HiggsZgammaEffCouplingTagger/wzp6_ee_ZH_ecm240.sin
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,104 @@
########################################################################
#
# Copyright (C) 1999-2019 by
# Wolfgang Kilian <kilian@physik.uni-siegen.de>
# Thorsten Ohl <ohl@physik.uni-wuerzburg.de>
# Juergen Reuter <juergen.reuter@desy.de>
# 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.
#
########################################################################

model = SM

# Center of mass energy
sqrts = 240 GeV

# Processes

?vis_diags = true

# ZH
# Z decays are handled by Pythia

process proc = e1, E1 => Z, H

beams = e1, E1 => gaussian => isr
?keep_beams = true
?keep_remnants = true

gaussian_spread1 = 0.221%
gaussian_spread2 = 0.221%


?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

# I 210 d + dbar I 2.142D-09 I 7.849D-07 I off I 0.000D+00 I
# I 211 u + ubar I 6.855D-10 I 2.511D-07 I off I 0.000D+00 I
# I 212 s + sbar I 8.656D-07 I 3.171D-04 I off I 0.000D+00 I
# I 213 c + cbar I 1.070D-04 I 3.919D-02 I off I 0.000D+00 I
# I 214 b + bbar I 1.952D-03 I 7.150D-01 I on I 1.000D+00 I
# I 215 t + tbar I 0.000D+00 I 0.000D+00 I off I 0.000D+00 I
# I 218 e- + e+ I 1.938D-11 I 7.101D-09 I off I 0.000D+00 I
# I 219 mu- + mu+ I 8.319D-07 I 3.048D-04 I off I 0.000D+00 I
# I 220 tau- + tau+ I 2.350D-04 I 8.608D-02 I off I 0.000D+00 I
# I 222 g + g I 1.490D-04 I 5.460D-02 I off I 0.000D+00 I
# I 223 gamma + gamma I 7.496D-06 I 2.746D-03 I off I 0.000D+00 I
# I 224 gamma + Z0 I 2.255D-06 I 8.263D-04 I off I 0.000D+00 I
# I 225 Z0 + Z0 I 2.561D-05 I 9.383D-03 I off I 0.000D+00 I
# I 226 W+ + W- I 2.498D-04 I 9.153D-02 I off I 0.000D+00 I

# Z decays:
# MDME(174,1) !Z decay into d dbar
# MDME(175,1) !Z decay into u ubar
# MDME(176,1) !Z decay into s sbar
# MDME(177,1) !Z decay into c cbar
# MDME(178,1) !Z decay into b bbar
# MDME(179,1) !Z decay into t tbar
# MDME(180,1) !Z decay into b4 b4bar
# MDME(181,1) !Z decay into t4 t4bar
# MDME(182,1) !Z decay into e- e+
# MDME(183,1) !Z decay into nu_e nu_ebar
# MDME(184,1) !Z decay into mu- mu+
# MDME(185,1) !Z decay into nu_mu nu_mubar
# MDME(186,1) !Z decay into tau- tau+
# MDME(187,1) !Z decay into nu_tau nu_taubar
# MDME(188,1) !Z decay into tau4 tau4bar
# MDME(189,1) !Z decay into nu_tau4 nu_tau4bar

# Z -> tau
#$ps_PYTHIA_PYGIVE = "MDME(174,1)=0; MDME(175,1)=0; MDME(176,1)=0; MDME(177,1)=0; MDME(178,1)=0; MDME(179,1)=0; MDME(180,1)=0; MDME(181,1)=0; MDME(182,1)=0; MDME(183,1)=0; MDME(184,1)=0; MDME(185,1)=0; MDME(186,1)=1; MDME(187,1)=0; MDME(188,1)=0; MDME(189,1)=0; MDME(210,1)=0; MDME(211,1)=0; MDME(212,1)=0; MDME(213,1)=0; MDME(214,1)=0; MDME(215,1)=0; MDME(218,1)=0; MDME(219,1)=0; MDME(220,1)=0; MDME(222,1)=0; MDME(223,1)=0; MDME(224,1)=0; MDME(225,1)=1; MDME(226,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"
$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; MSTJ(22)=4; PARJ(73)=2250; PARJ(74)=2500"


integrate (proc) { iterations = 10:100000:"gw", 5:200000:"" }

sample_format = stdhep

# Number of events and simulation step is passed in by the DIRAC
44 changes: 44 additions & 0 deletions ...ions/spring2024ConfigFiles/ee_xxgamma/HiggsZgammaEffCouplingTagger/whizardddsimgaudi.conf
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,44 @@
[whizard2]
Version = 2.8.3

[ddsim]
SteeringFile = cld_steer.py
ExtraCLIArguments = --vertexSigma 0.0098 2.54e-5 0.646 0.0063 --physics.rejectPDGs 1,2,3,4,5,6,12,14,16,21,23,24,25

[gaudiapp]
ExecutableName = k4run
SteeringFile = CLDReconstruction_jet_tagger_impl.py
InputFileFlag = --inputFiles
OutputFileFlag = --IOSvc.Output
Arguments = --enableLCFIJet --enableMLJetTagger --cms 240

[Production Parameters]
machine = ee
prodGroup = fcc_ee_spring2024_HiggsZgammaEffCouplingTagger1

softwareVersion = key4hep_nightly
generatorApplication = whizard2
simulationApplication = ddsim
reconstructionApplication = gaudiapp
generatorSteeringFile = wzp6_ee_eegamma_ecm240.sin, wzp6_ee_mumugamma_ecm240.sin, wzp6_ee_tautaugamma_ecm240.sin, wzp6_ee_qqgamma_ecm240.sin, wzp6_ee_ccgamma_ecm240.sin, wzp6_ee_bbgamma_ecm240.sin

configVersion = key4hep-devel-6
configPackage = fccConfig
eventsPerJobs = 1000, 1000, 1000, 1000, 1000, 1000

numberOfTasks = 1, 1, 1, 1, 1, 1

campaign = test_spring2024
energies = 240, 240, 240, 240, 240, 240
processes = eegamma, mumugamma, tautaugamma, qqgamma, ccgamma, bbgamma
ProdIDs = 16883, 16886, 16889, 16892, 16895, 16898
detectorModel = CLD_o2_v05
additionalName = HiggsZgammaEffCouplingTagger1
productionLogLevel = VERBOSE
outputSE = CERN-DST-EOS

finalOutputSE = CERN-SRM
MoveStatus = Stopped
MoveGroupSize = 10

ProdTypes = Rec
121 changes: 121 additions & 0 deletions .../spring2024ConfigFiles/ee_xxgamma/HiggsZgammaEffCouplingTagger/wzp6_ee_bbgamma_ecm240.sin
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,121 @@
########################################################################
#
# Copyright (C) 1999-2019 by
# Wolfgang Kilian <kilian@physik.uni-siegen.de>
# Thorsten Ohl <ohl@physik.uni-wuerzburg.de>
# Juergen Reuter <juergen.reuter@desy.de>
# 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.
#
########################################################################

model = SM

# Center of mass energy
sqrts = 240 GeV

# Processes

?vis_diags = true

# Z(->bb)+gamma
# Z decays are handled by Pythia


process proc = e1, E1 => b, B, A


cuts = let subevt @meA = select if Index > 2 [A] in
all E > 50 GeV [@meA]
and all 0.01 < Dist [@meA, b:B]
and all -0.998750 < cos (Theta) < 0.998750 [@meA]
and all -0.998750 < cos (Theta) < 0.998750 [b:B]
and all M > 5 GeV [b, B]

# cuts = let subevt @meA = select if Index > 2 [A]
# in all E > 50 GeV [@meA]
# and all -0.998750 < cos (Theta) < 0.998750 [@meA]
# and all M > 3GeV [@meA, e2:E2]
# and all M > 5 GeV [e2, E2]
# and all -0.998750 < cos (Theta) < 0.998750 [e2:E2]


beams = e1, E1 => gaussian => isr
?keep_beams = true
?keep_remnants = true

gaussian_spread1 = 0.221%
gaussian_spread2 = 0.221%


?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

# I 210 d + dbar I 2.142D-09 I 7.849D-07 I off I 0.000D+00 I
# I 211 u + ubar I 6.855D-10 I 2.511D-07 I off I 0.000D+00 I
# I 212 s + sbar I 8.656D-07 I 3.171D-04 I off I 0.000D+00 I
# I 213 c + cbar I 1.070D-04 I 3.919D-02 I off I 0.000D+00 I
# I 214 b + bbar I 1.952D-03 I 7.150D-01 I on I 1.000D+00 I
# I 215 t + tbar I 0.000D+00 I 0.000D+00 I off I 0.000D+00 I
# I 218 e- + e+ I 1.938D-11 I 7.101D-09 I off I 0.000D+00 I
# I 219 mu- + mu+ I 8.319D-07 I 3.048D-04 I off I 0.000D+00 I
# I 220 tau- + tau+ I 2.350D-04 I 8.608D-02 I off I 0.000D+00 I
# I 222 g + g I 1.490D-04 I 5.460D-02 I off I 0.000D+00 I
# I 223 gamma + gamma I 7.496D-06 I 2.746D-03 I off I 0.000D+00 I
# I 224 gamma + Z0 I 2.255D-06 I 8.263D-04 I off I 0.000D+00 I
# I 225 Z0 + Z0 I 2.561D-05 I 9.383D-03 I off I 0.000D+00 I
# I 226 W+ + W- I 2.498D-04 I 9.153D-02 I off I 0.000D+00 I

# Z decays:
# MDME(174,1) !Z decay into d dbar
# MDME(175,1) !Z decay into u ubar
# MDME(176,1) !Z decay into s sbar
# MDME(177,1) !Z decay into c cbar
# MDME(178,1) !Z decay into b bbar
# MDME(179,1) !Z decay into t tbar
# MDME(180,1) !Z decay into b4 b4bar
# MDME(181,1) !Z decay into t4 t4bar
# MDME(182,1) !Z decay into e- e+
# MDME(183,1) !Z decay into nu_e nu_ebar
# MDME(184,1) !Z decay into mu- mu+
# MDME(185,1) !Z decay into nu_mu nu_mubar
# MDME(186,1) !Z decay into tau- tau+
# MDME(187,1) !Z decay into nu_tau nu_taubar
# MDME(188,1) !Z decay into tau4 tau4bar
# MDME(189,1) !Z decay into nu_tau4 nu_tau4bar

# Z -> tau
#$ps_PYTHIA_PYGIVE = "MDME(174,1)=0; MDME(175,1)=0; MDME(176,1)=0; MDME(177,1)=0; MDME(178,1)=0; MDME(179,1)=0; MDME(180,1)=0; MDME(181,1)=0; MDME(182,1)=0; MDME(183,1)=0; MDME(184,1)=0; MDME(185,1)=0; MDME(186,1)=1; MDME(187,1)=0; MDME(188,1)=0; MDME(189,1)=0; MDME(210,1)=0; MDME(211,1)=0; MDME(212,1)=0; MDME(213,1)=0; MDME(214,1)=0; MDME(215,1)=0; MDME(218,1)=0; MDME(219,1)=0; MDME(220,1)=0; MDME(222,1)=0; MDME(223,1)=0; MDME(224,1)=0; MDME(225,1)=1; MDME(226,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"
$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; MSTJ(22)=4; PARJ(73)=2250; PARJ(74)=2500"


integrate (proc) { iterations = 10:100000:"gw", 5:200000:"" }

sample_format = stdhep

# Number of events and simulation step is passed in by the DIRAC
Loading