Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Sign in
Toggle navigation
Menu
Open sidebar
atlas-physics
pmg
MC Job Options
Commits
bebd61ff
Commit
bebd61ff
authored
Sep 10, 2021
by
Yanlin Liu
Committed by
Jan Kretzschmar
Sep 10, 2021
Browse files
JOs for SM non-resonant HH->bb+ZZ and HH->yy+leptons VBF pythia8 signal
parent
a57083cf
Pipeline
#3006057
passed with stages
in 1 minute and 19 seconds
Changes
5
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
508xxx/508674/mc.MGPy8EG_NNPDF30NLO_hhml_bbZZ2l2q_l1cvv1cv1.py
0 → 100644
View file @
bebd61ff
import
MadGraphControl.MadGraph_NNPDF30NLO_Base_Fragment
from
MadGraphControl.MadGraphUtils
import
*
beamEnergy
=-
999
if
hasattr
(
runArgs
,
'ecmEnergy'
):
beamEnergy
=
runArgs
.
ecmEnergy
/
2.
else
:
raise
RuntimeError
(
"No center of mass energy found."
)
#---------------------------------------------------------------------------------------------------
# Due to the low filter efficiency, the number of generated events are set to safefactor times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor
=
30
nevents
=
5000
*
safefactor
if
runArgs
.
maxEvents
>
0
:
nevents
=
runArgs
.
maxEvents
*
safefactor
mode
=
0
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
=
{}
parameters
[
'NEW'
]
=
{
'CV'
:
'1.0000000'
,
# CV
'C2V'
:
'1.0000000'
,
# C2V
'C3'
:
'1.0000000'
}
# C3
#---------------------------------------------------------------------------------------------------
# Setting higgs mass to 125 GeV for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
[
'MASS'
]
=
{
'25'
:
'1.250000e+02'
}
#MH
#---------------------------------------------------------------------------------------------------
# Setting some parameters for run_card.dat
#---------------------------------------------------------------------------------------------------
extras
=
{
'lhe_version'
:
'2.0'
,
'cut_decays'
:
'F'
,
'scale'
:
'125'
,
'dsqrt_q2fact1'
:
'125'
,
'dsqrt_q2fact2'
:
'125'
,
'nevents'
:
int
(
nevents
)}
#---------------------------------------------------------------------------------------------------
# Generating non-resonant VBF-Only HH process with MadGraph
# Parameters are set above
#---------------------------------------------------------------------------------------------------
process
=
"""
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
import model /cvmfs/atlas.cern.ch/repo/sw/Generators/madgraph/models/latest/HHVBF_UFO
generate p p > h h j j $$ z w+ w- / a j QED=4
output -f"""
process_dir
=
new_process
(
process
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the run_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new run_card.dat from an existing one
# Using the values given in "extras" above for the selected parameters when setting up the run_card
# If not set in "extras", default values are used
#---------------------------------------------------------------------------------------------------
modify_run_card
(
process_dir
=
process_dir
,
runArgs
=
runArgs
,
settings
=
extras
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the param_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new param_card.dat from an existing one
# Used values given in "parameters" for CV, C2V, C3
# Higgs mass is set to 125 GeV by "higgsMass"
#---------------------------------------------------------------------------------------------------
modify_param_card
(
process_dir
=
process_dir
,
params
=
parameters
)
#---------------------------------------------------------------------------------------------------
# Printing cards
#---------------------------------------------------------------------------------------------------
print_cards
()
#---------------------------------------------------------------------------------------------------
# Generate events
#---------------------------------------------------------------------------------------------------
generate
(
process_dir
=
process_dir
,
runArgs
=
runArgs
)
#---------------------------------------------------------------------------------------------------
# Move output files into the appropriate place, with the appropriate name
#---------------------------------------------------------------------------------------------------
arrange_output
(
process_dir
=
process_dir
,
runArgs
=
runArgs
,
lhe_version
=
2
,
saveProcDir
=
True
)
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig
.
generators
=
[
"MadGraph"
,
"Pythia8"
]
evgenConfig
.
description
=
"SM diHiggs production, decay to multi-lepton, with Powheg-Box-V2, at NLO + full top mass."
evgenConfig
.
keywords
=
[
"hh"
,
"SM"
,
"SMHiggs"
,
"nonResonant"
,
"VBF"
,
"bbbar"
,
"ZZ"
]
evgenConfig
.
contact
=
[
'Varsha Senthilkumar <varsha.senthilkumar@cern.ch>'
]
evgenConfig
.
nEventsPerJob
=
10000
#evgenConfig.tune = "MMHT2014"
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include
(
"Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py"
)
include
(
"Pythia8_i/Pythia8_MadGraph.py"
)
# configure Pythia8
genSeq
.
Pythia8
.
Commands
+=
[
"25:oneChannel = on 0.5 100 23 23 "
,
# ZZ decay
"25:addChannel = on 0.5 100 5 -5 "
,
# bb decay
"23:onMode = off"
,
"23:onIfAny= 1 2 3 4 5 11 13 15"
,
"24:mMin = 0"
,
# W minimum mass
"24:mMax = 99999"
,
# W maximum mass
"23:mMin = 0"
,
# Z minimum mass
"23:mMax = 99999"
,
# Z maximum mass "TimeShower:mMaxGamma = 0" ] # Z/gamma* combination scale
"TimeShower:mMaxGamma = 0"
]
# Z/gamma* combination scale
#--------------------------------------------------------------
# Dipole option Pythia8
#--------------------------------------------------------------
genSeq
.
Pythia8
.
Commands
+=
[
"SpaceShower:dipoleRecoil = on"
]
# #---------------------------------------------------------------------------------------------------
# # Generator Filters
# #---------------------------------------------------------------------------------------------------
from
GeneratorFilters.GeneratorFiltersConf
import
ParentChildFilter
filtSeq
+=
ParentChildFilter
(
"hbbFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
5
])
filtSeq
+=
ParentChildFilter
(
"hZZFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
23
])
filtSeq
+=
ParentChildFilter
(
"ZllFilter"
,
PDGParent
=
[
23
],
PDGChild
=
[
11
,
13
,
15
])
filtSeq
+=
ParentChildFilter
(
"ZqqFilter"
,
PDGParent
=
[
23
],
PDGChild
=
[
1
,
2
,
3
,
4
,
5
])
# -- require HH->bbZZ
filtSeq
.
Expression
=
"hbbFilter and hZZFilter and ZllFilter and ZqqFilter"
\ No newline at end of file
508xxx/508675/mc.MGPy8EG_NNPDF30NLO_hhml_bbZZ4l_l1cvv1cv1.py
0 → 100644
View file @
bebd61ff
import
MadGraphControl.MadGraph_NNPDF30NLO_Base_Fragment
from
MadGraphControl.MadGraphUtils
import
*
beamEnergy
=-
999
if
hasattr
(
runArgs
,
'ecmEnergy'
):
beamEnergy
=
runArgs
.
ecmEnergy
/
2.
else
:
raise
RuntimeError
(
"No center of mass energy found."
)
#---------------------------------------------------------------------------------------------------
# Due to the low filter efficiency, the number of generated events are set to safefactor times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor
=
30
nevents
=
5000
*
safefactor
if
runArgs
.
maxEvents
>
0
:
nevents
=
runArgs
.
maxEvents
*
safefactor
mode
=
0
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
=
{}
parameters
[
'NEW'
]
=
{
'CV'
:
'1.0000000'
,
# CV
'C2V'
:
'1.0000000'
,
# C2V
'C3'
:
'1.0000000'
}
# C3
#---------------------------------------------------------------------------------------------------
# Setting higgs mass to 125 GeV for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
[
'MASS'
]
=
{
'25'
:
'1.250000e+02'
}
#MH
#---------------------------------------------------------------------------------------------------
# Setting some parameters for run_card.dat
#---------------------------------------------------------------------------------------------------
extras
=
{
'lhe_version'
:
'2.0'
,
'cut_decays'
:
'F'
,
'scale'
:
'125'
,
'dsqrt_q2fact1'
:
'125'
,
'dsqrt_q2fact2'
:
'125'
,
'nevents'
:
int
(
nevents
)}
#---------------------------------------------------------------------------------------------------
# Generating non-resonant VBF-Only HH process with MadGraph
# Parameters are set above
#---------------------------------------------------------------------------------------------------
process
=
"""
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
import model /cvmfs/atlas.cern.ch/repo/sw/Generators/madgraph/models/latest/HHVBF_UFO
generate p p > h h j j $$ z w+ w- / a j QED=4
output -f"""
process_dir
=
new_process
(
process
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the run_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new run_card.dat from an existing one
# Using the values given in "extras" above for the selected parameters when setting up the run_card
# If not set in "extras", default values are used
#---------------------------------------------------------------------------------------------------
modify_run_card
(
process_dir
=
process_dir
,
runArgs
=
runArgs
,
settings
=
extras
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the param_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new param_card.dat from an existing one
# Used values given in "parameters" for CV, C2V, C3
# Higgs mass is set to 125 GeV by "higgsMass"
#---------------------------------------------------------------------------------------------------
modify_param_card
(
process_dir
=
process_dir
,
params
=
parameters
)
#---------------------------------------------------------------------------------------------------
# Printing cards
#---------------------------------------------------------------------------------------------------
print_cards
()
#---------------------------------------------------------------------------------------------------
# Generate events
#---------------------------------------------------------------------------------------------------
generate
(
process_dir
=
process_dir
,
runArgs
=
runArgs
)
#---------------------------------------------------------------------------------------------------
# Move output files into the appropriate place, with the appropriate name
#---------------------------------------------------------------------------------------------------
arrange_output
(
process_dir
=
process_dir
,
runArgs
=
runArgs
,
lhe_version
=
2
,
saveProcDir
=
True
)
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig
.
generators
=
[
"MadGraph"
,
"Pythia8"
]
evgenConfig
.
description
=
"SM diHiggs production, decay to multi-lepton, with Powheg-Box-V2, at NLO + full top mass."
evgenConfig
.
keywords
=
[
"hh"
,
"SM"
,
"SMHiggs"
,
"nonResonant"
,
"VBF"
,
"bbbar"
,
"ZZ"
]
evgenConfig
.
contact
=
[
'Varsha Senthilkumar <varsha.senthilkumar@cern.ch>'
]
evgenConfig
.
nEventsPerJob
=
10000
#evgenConfig.tune = "MMHT2014"
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include
(
"Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py"
)
include
(
"Pythia8_i/Pythia8_MadGraph.py"
)
# configure Pythia8
genSeq
.
Pythia8
.
Commands
+=
[
"25:oneChannel = on 0.99 100 23 23 "
,
# ZZ decay
"25:addChannel = on 0.01 100 5 -5 "
,
# bb decay
"23:onMode = off"
,
"23:onIfAny= 11 13 15"
,
"24:mMin = 0"
,
# W minimum mass
"24:mMax = 99999"
,
# W maximum mass
"23:mMin = 0"
,
# Z minimum mass
"23:mMax = 99999"
,
# Z maximum mass "TimeShower:mMaxGamma = 0" ] # Z/gamma* combination scale
"TimeShower:mMaxGamma = 0"
]
# Z/gamma* combination scale
#--------------------------------------------------------------
# Dipole option Pythia8
#--------------------------------------------------------------
genSeq
.
Pythia8
.
Commands
+=
[
"SpaceShower:dipoleRecoil = on"
]
# #---------------------------------------------------------------------------------------------------
# # Generator Filters
# #---------------------------------------------------------------------------------------------------
from
GeneratorFilters.GeneratorFiltersConf
import
ParentChildFilter
filtSeq
+=
ParentChildFilter
(
"hbbFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
5
])
filtSeq
+=
ParentChildFilter
(
"hZZFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
23
])
filtSeq
.
Expression
=
"hbbFilter and hZZFilter"
\ No newline at end of file
508xxx/508676/mc.MGPy8EG_NNPDF30NLO_hhml_yy0l_l1cvv1cv1.py
0 → 100644
View file @
bebd61ff
import
MadGraphControl.MadGraph_NNPDF30NLO_Base_Fragment
from
MadGraphControl.MadGraphUtils
import
*
beamEnergy
=-
999
if
hasattr
(
runArgs
,
'ecmEnergy'
):
beamEnergy
=
runArgs
.
ecmEnergy
/
2.
else
:
raise
RuntimeError
(
"No center of mass energy found."
)
#---------------------------------------------------------------------------------------------------
# Due to the low filter efficiency, the number of generated events are set to safefactor times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor
=
30
nevents
=
5000
*
safefactor
if
runArgs
.
maxEvents
>
0
:
nevents
=
runArgs
.
maxEvents
*
safefactor
mode
=
0
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
=
{}
parameters
[
'NEW'
]
=
{
'CV'
:
'1.0000000'
,
# CV
'C2V'
:
'1.0000000'
,
# C2V
'C3'
:
'1.0000000'
}
# C3
#---------------------------------------------------------------------------------------------------
# Setting higgs mass to 125 GeV for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
[
'MASS'
]
=
{
'25'
:
'1.250000e+02'
}
#MH
#---------------------------------------------------------------------------------------------------
# Setting some parameters for run_card.dat
#---------------------------------------------------------------------------------------------------
extras
=
{
'lhe_version'
:
'2.0'
,
'cut_decays'
:
'F'
,
'scale'
:
'125'
,
'dsqrt_q2fact1'
:
'125'
,
'dsqrt_q2fact2'
:
'125'
,
'nevents'
:
int
(
nevents
)}
#---------------------------------------------------------------------------------------------------
# Generating non-resonant VBF-Only HH process with MadGraph
# Parameters are set above
#---------------------------------------------------------------------------------------------------
process
=
"""
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
import model /cvmfs/atlas.cern.ch/repo/sw/Generators/madgraph/models/latest/HHVBF_UFO
generate p p > h h j j $$ z w+ w- / a j QED=4
output -f"""
process_dir
=
new_process
(
process
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the run_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new run_card.dat from an existing one
# Using the values given in "extras" above for the selected parameters when setting up the run_card
# If not set in "extras", default values are used
#---------------------------------------------------------------------------------------------------
modify_run_card
(
process_dir
=
process_dir
,
runArgs
=
runArgs
,
settings
=
extras
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the param_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new param_card.dat from an existing one
# Used values given in "parameters" for CV, C2V, C3
# Higgs mass is set to 125 GeV by "higgsMass"
#---------------------------------------------------------------------------------------------------
modify_param_card
(
process_dir
=
process_dir
,
params
=
parameters
)
#---------------------------------------------------------------------------------------------------
# Printing cards
#---------------------------------------------------------------------------------------------------
print_cards
()
#---------------------------------------------------------------------------------------------------
# Generate events
#---------------------------------------------------------------------------------------------------
generate
(
process_dir
=
process_dir
,
runArgs
=
runArgs
)
#---------------------------------------------------------------------------------------------------
# Move output files into the appropriate place, with the appropriate name
#---------------------------------------------------------------------------------------------------
arrange_output
(
process_dir
=
process_dir
,
runArgs
=
runArgs
,
lhe_version
=
2
,
saveProcDir
=
True
)
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig
.
generators
=
[
"MadGraph"
,
"Pythia8"
]
evgenConfig
.
description
=
"SM diHiggs production, decay to multi-lepton, with Powheg-Box-V2, at NLO + full top mass."
evgenConfig
.
keywords
=
[
"hh"
,
"SM"
,
"SMHiggs"
,
"nonResonant"
,
"VBF"
,
"multilepton"
]
evgenConfig
.
contact
=
[
'Varsha Senthilkumar <varsha.senthilkumar@cern.ch>'
]
evgenConfig
.
nEventsPerJob
=
10000
#evgenConfig.tune = "MMHT2014"
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include
(
"Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py"
)
include
(
"Pythia8_i/Pythia8_MadGraph.py"
)
# configure Pythia8
genSeq
.
Pythia8
.
Commands
+=
[
"25:oneChannel = on 0.353 100 24 -24 "
,
# WW decay
"25:addChannel = on 0.043 100 23 23 "
,
# ZZ decay
"25:addChannel = on 0.104 100 15 -15 "
,
# tautau decay
"25:addChannel = on 0.5 100 22 22 "
,
# yy decay
"24:mMin = 0"
,
# W minimum mass
"24:mMax = 99999"
,
# W maximum mass
"23:mMin = 0"
,
# Z minimum mass
"23:mMax = 99999"
,
# Z maximum mass "TimeShower:mMaxGamma = 0" ] # Z/gamma* combination scale
"TimeShower:mMaxGamma = 0"
]
# Z/gamma* combination scale
#--------------------------------------------------------------
# Dipole option Pythia8
#--------------------------------------------------------------
genSeq
.
Pythia8
.
Commands
+=
[
"SpaceShower:dipoleRecoil = on"
]
# #---------------------------------------------------------------------------------------------------
# # Generator Filters
# #---------------------------------------------------------------------------------------------------
from
GeneratorFilters.GeneratorFiltersConf
import
ParentChildFilter
filtSeq
+=
ParentChildFilter
(
"hyyFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
22
])
filtSeq
+=
ParentChildFilter
(
"hXXFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
15
,
23
,
24
])
from
GeneratorFilters.GeneratorFiltersConf
import
MultiElecMuTauFilter
filtSeq
+=
MultiElecMuTauFilter
(
"LepOneFilter"
)
filtSeq
.
LepOneFilter
.
IncludeHadTaus
=
True
filtSeq
.
LepOneFilter
.
NLeptons
=
1
filtSeq
.
LepOneFilter
.
MinPt
=
7000
filtSeq
.
LepOneFilter
.
MinVisPtHadTau
=
15000
filtSeq
.
LepOneFilter
.
MaxEta
=
3
filtSeq
.
Expression
=
"hyyFilter and hXXFilter and not LepOneFilter"
\ No newline at end of file
508xxx/508677/mc.MGPy8EG_NNPDF30NLO_hhml_yy1l_l1cvv1cv1.py
0 → 100644
View file @
bebd61ff
import
MadGraphControl.MadGraph_NNPDF30NLO_Base_Fragment
from
MadGraphControl.MadGraphUtils
import
*
beamEnergy
=-
999
if
hasattr
(
runArgs
,
'ecmEnergy'
):
beamEnergy
=
runArgs
.
ecmEnergy
/
2.
else
:
raise
RuntimeError
(
"No center of mass energy found."
)
#---------------------------------------------------------------------------------------------------
# Due to the low filter efficiency, the number of generated events are set to safefactor times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor
=
30
nevents
=
5000
*
safefactor
if
runArgs
.
maxEvents
>
0
:
nevents
=
runArgs
.
maxEvents
*
safefactor
mode
=
0
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
=
{}
parameters
[
'NEW'
]
=
{
'CV'
:
'1.0000000'
,
# CV
'C2V'
:
'1.0000000'
,
# C2V
'C3'
:
'1.0000000'
}
# C3
#---------------------------------------------------------------------------------------------------
# Setting higgs mass to 125 GeV for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
[
'MASS'
]
=
{
'25'
:
'1.250000e+02'
}
#MH
#---------------------------------------------------------------------------------------------------
# Setting some parameters for run_card.dat
#---------------------------------------------------------------------------------------------------
extras
=
{
'lhe_version'
:
'2.0'
,
'cut_decays'
:
'F'
,
'scale'
:
'125'
,
'dsqrt_q2fact1'
:
'125'
,
'dsqrt_q2fact2'
:
'125'
,
'nevents'
:
int
(
nevents
)}
#---------------------------------------------------------------------------------------------------
# Generating non-resonant VBF-Only HH process with MadGraph
# Parameters are set above
#---------------------------------------------------------------------------------------------------
process
=
"""
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
import model /cvmfs/atlas.cern.ch/repo/sw/Generators/madgraph/models/latest/HHVBF_UFO
generate p p > h h j j $$ z w+ w- / a j QED=4
output -f"""
process_dir
=
new_process
(
process
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the run_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new run_card.dat from an existing one
# Using the values given in "extras" above for the selected parameters when setting up the run_card
# If not set in "extras", default values are used
#---------------------------------------------------------------------------------------------------
modify_run_card
(
process_dir
=
process_dir
,
runArgs
=
runArgs
,
settings
=
extras
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the param_card
# https://gitlab.cern.ch/atlas/athena/-/tree/master/Generators/MadGraphControl
# Build a new param_card.dat from an existing one
# Used values given in "parameters" for CV, C2V, C3
# Higgs mass is set to 125 GeV by "higgsMass"
#---------------------------------------------------------------------------------------------------
modify_param_card
(
process_dir
=
process_dir
,
params
=
parameters
)
#---------------------------------------------------------------------------------------------------
# Printing cards
#---------------------------------------------------------------------------------------------------
print_cards
()
#---------------------------------------------------------------------------------------------------
# Generate events
#---------------------------------------------------------------------------------------------------
generate
(
process_dir
=
process_dir
,
runArgs
=
runArgs
)
#---------------------------------------------------------------------------------------------------
# Move output files into the appropriate place, with the appropriate name
#---------------------------------------------------------------------------------------------------
arrange_output
(
process_dir
=
process_dir
,
runArgs
=
runArgs
,
lhe_version
=
2
,
saveProcDir
=
True
)
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig
.
generators
=
[
"MadGraph"
,
"Pythia8"
]
evgenConfig
.
description
=
"SM diHiggs production, decay to multi-lepton, with Powheg-Box-V2, at NLO + full top mass."
evgenConfig
.
keywords
=
[
"hh"
,
"SM"
,
"SMHiggs"
,
"nonResonant"
,
"VBF"
,
"multilepton"
]
evgenConfig
.
contact
=
[
'Varsha Senthilkumar <varsha.senthilkumar@cern.ch>'
]
evgenConfig
.
nEventsPerJob
=
10000
#evgenConfig.tune = "MMHT2014"
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include
(
"Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py"
)
include
(
"Pythia8_i/Pythia8_MadGraph.py"
)
# configure Pythia8
genSeq
.
Pythia8
.
Commands
+=
[
"25:oneChannel = on 0.353 100 24 -24 "
,
# WW decay
"25:addChannel = on 0.043 100 23 23 "
,
# ZZ decay
"25:addChannel = on 0.104 100 15 -15 "
,
# tautau decay
"25:addChannel = on 0.5 100 22 22 "
,
# yy decay
"24:mMin = 0"
,
# W minimum mass
"24:mMax = 99999"
,
# W maximum mass
"23:mMin = 0"
,
# Z minimum mass
"23:mMax = 99999"
,
# Z maximum mass "TimeShower:mMaxGamma = 0" ] # Z/gamma* combination scale
"TimeShower:mMaxGamma = 0"
]
# Z/gamma* combination scale
#--------------------------------------------------------------
# Dipole option Pythia8
#--------------------------------------------------------------
genSeq
.
Pythia8
.
Commands
+=
[
"SpaceShower:dipoleRecoil = on"
]
# #---------------------------------------------------------------------------------------------------
# # Generator Filters
# #---------------------------------------------------------------------------------------------------
from
GeneratorFilters.GeneratorFiltersConf
import
ParentChildFilter
filtSeq
+=
ParentChildFilter
(
"hyyFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
22
])
filtSeq
+=
ParentChildFilter
(
"hXXFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
15
,
23
,
24
])
from
GeneratorFilters.GeneratorFiltersConf
import
MultiElecMuTauFilter
filtSeq
+=
MultiElecMuTauFilter
(
"LepOneFilter"
)
filtSeq
.
LepOneFilter
.
IncludeHadTaus
=
True
filtSeq
.
LepOneFilter
.
NLeptons
=
1
filtSeq
.
LepOneFilter
.
MinPt
=
7000
filtSeq
.
LepOneFilter
.
MinVisPtHadTau
=
15000
filtSeq
.
LepOneFilter
.
MaxEta
=
3
filtSeq
+=
MultiElecMuTauFilter
(
"LepTwoFilter"
)
filtSeq
.
LepTwoFilter
.
IncludeHadTaus
=
True
filtSeq
.
LepTwoFilter
.
NLeptons
=
2
filtSeq
.
LepTwoFilter
.
MinPt
=
7000
filtSeq
.
LepTwoFilter
.
MinVisPtHadTau
=
15000
filtSeq
.
LepTwoFilter
.
MaxEta
=
3
filtSeq
.
Expression
=
"hyyFilter and hXXFilter and LepOneFilter and not LepTwoFilter"
\ No newline at end of file
508xxx/508678/mc.MGPy8EG_NNPDF30NLO_hhml_yy2l_l1cvv1cv1.py
0 → 100644
View file @
bebd61ff
import
MadGraphControl.MadGraph_NNPDF30NLO_Base_Fragment
from
MadGraphControl.MadGraphUtils
import
*
beamEnergy
=-
999
if
hasattr
(
runArgs
,
'ecmEnergy'
):
beamEnergy
=
runArgs
.
ecmEnergy
/
2.
else
:
raise
RuntimeError
(
"No center of mass energy found."
)
#---------------------------------------------------------------------------------------------------
# Due to the low filter efficiency, the number of generated events are set to safefactor times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor
=
30
nevents
=
5000
*
safefactor
if
runArgs
.
maxEvents
>
0
:
nevents
=
runArgs
.
maxEvents
*
safefactor
mode
=
0
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat