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
infrastructure
MC15JobOptions
Commits
fe11ddbe
Commit
fe11ddbe
authored
Feb 11, 2021
by
Cyril Becot
Browse files
New tag
parent
1ac8e90a
Changes
11
Hide whitespace changes
Inline
Side-by-side
ChangeLog
View file @
fe11ddbe
2021-02-11 Cyril Becot <cyril.becot@cern.ch>
* add 300768
* add 451892-451898
* tag MC15JobOptions-01-04-43/
2021-02-04 Joany Manjarres <joany@cern.ch>
* tag MC15JobOptions-01-04-42
* add 399670-883
...
...
share/DSID300xxx/MC15.300768.Pythia8B_Photospp_A14_CTEQ6L1_AntiBs_Jpsimu3p5mu3p5_phi.py
0 → 100644
View file @
fe11ddbe
##############################################################
# Job options for Pythia8B_i generation of antiBs->J/psi(mu4m4)phi(KK)
##############################################################
evgenConfig
.
description
=
"Signal antiBs->J/psi(mu3p5mu3p5)phi(KK)"
evgenConfig
.
keywords
=
[
"exclusive"
,
"Bs"
,
"Jpsi"
,
"2muon"
]
evgenConfig
.
minevents
=
500
include
(
"MC15JobOptions/nonStandard/Pythia8B_A14_CTEQ6L1_Common.py"
)
include
(
'MC15JobOptions/nonStandard/Pythia8B_Photospp.py'
)
include
(
"MC15JobOptions/Pythia8B_exclusiveAntiB_Common.py"
)
genSeq
.
Pythia8B
.
Commands
+=
[
'PhaseSpace:pTHatMin = 11.'
]
genSeq
.
Pythia8B
.
QuarkPtCut
=
9.0
genSeq
.
Pythia8B
.
AntiQuarkPtCut
=
0.0
genSeq
.
Pythia8B
.
QuarkEtaCut
=
2.5
genSeq
.
Pythia8B
.
AntiQuarkEtaCut
=
102.5
genSeq
.
Pythia8B
.
RequireBothQuarksPassCuts
=
True
# define pseudo-J/psi for exclusive decay
genSeq
.
Pythia8B
.
Commands
+=
[
'999443:all = myJ/psi void 3 0 0 3.09692 0.00009 3.09602 3.09782 0'
]
# name antiName spinType chargeType colType m0 mWidth mMin mMax tau0
# force its decay to mu+ mu-
genSeq
.
Pythia8B
.
Commands
+=
[
'999443:onMode = off'
]
genSeq
.
Pythia8B
.
Commands
+=
[
'999443:oneChannel = on 1. 0 -13 13'
]
# onMode bRatio meMode products
# force B+ -> myJ/psi K+
genSeq
.
Pythia8B
.
Commands
+=
[
'531:addChannel = 3 1. 0 999443 333'
]
# onMode bRatio meMode products
genSeq
.
Pythia8B
.
Commands
+=
[
'531:tau0 = 0.4527'
]
genSeq
.
Pythia8B
.
SignalPDGCodes
=
[
-
531
,
999443
,
-
13
,
13
,
333
,
321
,
-
321
]
genSeq
.
Pythia8B
.
SignalPtCuts
=
[
0.0
,
0.0
,
0.0
,
0.0
,
0.0
,
0.8
,
0.8
]
# no cuts on muons here -- apply trigger-like selection below
genSeq
.
Pythia8B
.
SignalEtaCuts
=
[
102.5
,
102.5
,
102.5
,
102.5
,
102.5
,
2.6
,
2.6
]
# no cuts on muons here -- apply trigger-like selection below
genSeq
.
Pythia8B
.
NHadronizationLoops
=
4
genSeq
.
Pythia8B
.
TriggerPDGCode
=
13
genSeq
.
Pythia8B
.
TriggerStatePtCut
=
[
3.5
]
genSeq
.
Pythia8B
.
TriggerStateEtaCut
=
2.6
genSeq
.
Pythia8B
.
MinimumCountPerCut
=
[
2
]
share/DSID451xxx/MC15.451892.MadGraphHerwig7EvtGen_PDF23LO_X270tohh_bbyy.py
0 → 100644
View file @
fe11ddbe
from
MadGraphControl.MadGraphUtils
import
*
mode
=
0
mhh
=
2.700000e+02
#---------------------------------------------------------------------------------------------------
# Setting mHH and WHH for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
=
{
'1560'
:
'2.700000e+02'
,
#MHH
'1561'
:
'1.000000e-02'
}
#WHH
#---------------------------------------------------------------------------------------------------
# Setting higgs mass to 125 GeV for param_card.dat
#---------------------------------------------------------------------------------------------------
higgsMass
=
{
'25'
:
'1.250000e+02'
}
#MH
#---------------------------------------------------------------------------------------------------
# Setting some parameters for run_card.dat: Parton shower is HERWIG7
#---------------------------------------------------------------------------------------------------
extras
=
{
'lhe_version'
:
'2.0'
,
'cut_decays'
:
'F'
,
'pdlabel'
:
"'nn23lo1'"
,
'fixed_ren_scale'
:
'True'
,
#
'fixed_fac_scale'
:
'True'
,
#
'scale'
:
'135'
,
#
'dsqrt_q2fact1'
:
'135'
,
#
'dsqrt_q2fact2'
:
'135'
,
#
# 'parton_shower':'Herwig7',
'ptj'
:
'0'
,
'ptb'
:
'0'
,
'pta'
:
'0'
,
'ptjmax'
:
'-1'
,
'ptbmax'
:
'-1'
,
'ptamax'
:
'-1'
,
'etaj'
:
'-1'
,
'etab'
:
'-1'
,
'etaa'
:
'-1'
,
'etajmin'
:
'0'
,
'etabmin'
:
'0'
,
'etaamin'
:
'0'
,
'mmaa'
:
'0'
,
'mmaamax'
:
'-1'
,
'mmbb'
:
'0'
,
'mmbbmax'
:
'-1'
,
'drjj'
:
'0'
,
'drbb'
:
'0'
,
'draa'
:
'0'
,
'drbj'
:
'0'
,
'draj'
:
'0'
,
'drab'
:
'0'
,
'drjjmax'
:
'-1'
,
'drbbmax'
:
'-1'
,
'draamax'
:
'-1'
,
'drbjmax'
:
'-1'
,
'drajmax'
:
'-1'
,
'drabmax'
:
'-1'
}
#---------------------------------------------------------------------------------------------------
# Generating di-higgs through Heavy Higgs resonance with MadGraph
#---------------------------------------------------------------------------------------------------
fcard
=
open
(
'proc_card_mg5.dat'
,
'w'
)
if
(
runArgs
.
runNumber
==
451892
):
fcard
.
write
(
"""
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 HeavyHiggsTHDM
generate p p > hh > h h
output -f"""
)
fcard
.
write
(
'set Mhh %s
\n
'
%
(
mhh
))
fcard
.
close
()
else
:
raise
RuntimeError
(
"runNumber %i not recognised in these jobOptions."
%
runArgs
.
runNumber
)
beamEnergy
=-
999
if
hasattr
(
runArgs
,
'ecmEnergy'
):
beamEnergy
=
runArgs
.
ecmEnergy
/
2.
else
:
raise
RuntimeError
(
"No center of mass energy found."
)
#---------------------------------------------------------------------------------------------------
# Filter efficiency is low
# Thus, setting the number of generated events to 100 times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor
=
100
nevents
=
10000
*
safefactor
if
runArgs
.
maxEvents
>
0
:
nevents
=
runArgs
.
maxEvents
*
safefactor
runName
=
'run_01'
process_dir
=
new_process
(
card_loc
=
'proc_card_mg5.dat'
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the param_card
# Build a new param_card.dat from an existing one
# Used values given in "parameters" for MHH and WHH, if not set there, default values are used
# Higgs mass is set to 125 GeV by "higgsMass"
#---------------------------------------------------------------------------------------------------
build_param_card
(
param_card_old
=
'param_card.HeavyScalar.dat'
,
param_card_new
=
'param_card_new.dat'
,
masses
=
higgsMass
,
extras
=
parameters
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the run_card
# 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
#---------------------------------------------------------------------------------------------------
build_run_card
(
run_card_old
=
get_default_runcard
(
proc_dir
=
process_dir
),
run_card_new
=
'run_card.dat'
,
xqcut
=
0
,
nevts
=
nevents
,
rand_seed
=
runArgs
.
randomSeed
,
beamEnergy
=
beamEnergy
,
extras
=
extras
)
print_cards
()
generate
(
run_card_loc
=
'run_card.dat'
,
param_card_loc
=
'param_card_new.dat'
,
mode
=
mode
,
proc_dir
=
process_dir
,
run_name
=
runName
)
outputDS
=
arrange_output
(
run_name
=
runName
,
proc_dir
=
process_dir
,
outputDS
=
runName
+
'._00001.events.tar.gz'
)
runArgs
.
inputGeneratorFile
=
outputDS
#--------------------------------------------------------------
# Herwig 7 (H7UE) showering setup
#--------------------------------------------------------------
# initialize Herwig7 generator configuration for showering of LHE files
include
(
"MC15JobOptions/Herwig7_LHEF.py"
)
# configure Herwig7
Herwig7Config
.
add_commands
(
"set /Herwig/Partons/RemnantDecayer:AllowTop Yes"
)
Herwig7Config
.
me_pdf_commands
(
order
=
"LO"
,
name
=
"NNPDF23_lo_as_0130_qed"
)
Herwig7Config
.
tune_commands
()
Herwig7Config
.
lhef_mg5amc_commands
(
lhe_filename
=
runArgs
.
inputGeneratorFile
,
me_pdf_order
=
"LO"
)
evgenConfig
.
tune
=
"H7.1-Default"
# add EvtGen
include
(
"MC15JobOptions/Herwig71_EvtGen.py"
)
#HW7 settings and Higgs BR
Herwig7Config
.
add_commands
(
"""
#set /Herwig/Shower/KinematicsReconstructor:ReconstructionOption General
#set /Herwig/Shower/KinematicsReconstructor:InitialInitialBoostOption LongTransBoost
set /Herwig/Shower/LtoLGammaSudakov:pTmin 0.000001
set /Herwig/Shower/QtoGammaQSudakov:Alpha /Herwig/Shower/AlphaQED
cd /Herwig/Particles
create /ThePEG/ParticleData XH
setup XH 1560 XH MASS WIDTH MAX_WIDTH 0 0 0 1 0
#set /Herwig/EventHandlers/LHEReader:Decayer /Herwig/Decays/Mambo
do /Herwig/Particles/h0:SelectDecayModes h0->gamma,gamma; h0->b,bbar;
do /Herwig/Particles/h0:PrintDecayModes
set /Herwig/Particles/h0/h0->gamma,gamma;:BranchingRatio 0.5
set /Herwig/Particles/h0/h0->b,bbar;:BranchingRatio 0.5
"""
)
# run Herwig7
Herwig7Config
.
run
()
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig
.
generators
=
[
"MadGraph"
,
"Herwig7"
,
"EvtGen"
]
if
(
runArgs
.
runNumber
==
451892
):
evgenConfig
.
description
=
"Di-Higgs production through 270 GeV Heavy Higgs resonance which decays to bbyy."
evgenConfig
.
keywords
=
[
"hh"
,
"BSM"
,
"BSMHiggs"
,
"resonance"
,
"diphoton"
,
"bottom"
]
evgenConfig
.
contact
=
[
'Florian Beisiegel <florian.beisiegel@cern.ch>'
]
evgenConfig
.
inputfilecheck
=
runName
#runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
#---------------------------------------------------------------------------------------------------
# Filter for bbyy
#---------------------------------------------------------------------------------------------------
from
GeneratorFilters.GeneratorFiltersConf
import
ParentChildFilter
filtSeq
+=
ParentChildFilter
(
"HbbFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
5
])
filtSeq
+=
ParentChildFilter
(
"HyyFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
22
])
filtSeq
.
Expression
=
"HbbFilter and HyyFilter"
share/DSID451xxx/MC15.451893.MadGraphHerwig7EvtGen_PDF23LO_X290tohh_bbyy.py
0 → 100644
View file @
fe11ddbe
from
MadGraphControl.MadGraphUtils
import
*
mode
=
0
mhh
=
2.900000e+02
#---------------------------------------------------------------------------------------------------
# Setting mHH and WHH for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
=
{
'1560'
:
'2.900000e+02'
,
#MHH
'1561'
:
'1.000000e-02'
}
#WHH
#---------------------------------------------------------------------------------------------------
# Setting higgs mass to 125 GeV for param_card.dat
#---------------------------------------------------------------------------------------------------
higgsMass
=
{
'25'
:
'1.250000e+02'
}
#MH
#---------------------------------------------------------------------------------------------------
# Setting some parameters for run_card.dat: Parton shower is HERWIG7
#---------------------------------------------------------------------------------------------------
extras
=
{
'lhe_version'
:
'2.0'
,
'cut_decays'
:
'F'
,
'pdlabel'
:
"'nn23lo1'"
,
'fixed_ren_scale'
:
'True'
,
#
'fixed_fac_scale'
:
'True'
,
#
'scale'
:
'145'
,
#
'dsqrt_q2fact1'
:
'145'
,
#
'dsqrt_q2fact2'
:
'145'
,
#
# 'parton_shower':'Herwig7',
'ptj'
:
'0'
,
'ptb'
:
'0'
,
'pta'
:
'0'
,
'ptjmax'
:
'-1'
,
'ptbmax'
:
'-1'
,
'ptamax'
:
'-1'
,
'etaj'
:
'-1'
,
'etab'
:
'-1'
,
'etaa'
:
'-1'
,
'etajmin'
:
'0'
,
'etabmin'
:
'0'
,
'etaamin'
:
'0'
,
'mmaa'
:
'0'
,
'mmaamax'
:
'-1'
,
'mmbb'
:
'0'
,
'mmbbmax'
:
'-1'
,
'drjj'
:
'0'
,
'drbb'
:
'0'
,
'draa'
:
'0'
,
'drbj'
:
'0'
,
'draj'
:
'0'
,
'drab'
:
'0'
,
'drjjmax'
:
'-1'
,
'drbbmax'
:
'-1'
,
'draamax'
:
'-1'
,
'drbjmax'
:
'-1'
,
'drajmax'
:
'-1'
,
'drabmax'
:
'-1'
}
#---------------------------------------------------------------------------------------------------
# Generating di-higgs through Heavy Higgs resonance with MadGraph
#---------------------------------------------------------------------------------------------------
fcard
=
open
(
'proc_card_mg5.dat'
,
'w'
)
if
(
runArgs
.
runNumber
==
451893
):
fcard
.
write
(
"""
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 HeavyHiggsTHDM
generate p p > hh > h h
output -f"""
)
fcard
.
write
(
'set Mhh %s
\n
'
%
(
mhh
))
fcard
.
close
()
else
:
raise
RuntimeError
(
"runNumber %i not recognised in these jobOptions."
%
runArgs
.
runNumber
)
beamEnergy
=-
999
if
hasattr
(
runArgs
,
'ecmEnergy'
):
beamEnergy
=
runArgs
.
ecmEnergy
/
2.
else
:
raise
RuntimeError
(
"No center of mass energy found."
)
#---------------------------------------------------------------------------------------------------
# Filter efficiency is low
# Thus, setting the number of generated events to 100 times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor
=
100
nevents
=
10000
*
safefactor
if
runArgs
.
maxEvents
>
0
:
nevents
=
runArgs
.
maxEvents
*
safefactor
runName
=
'run_01'
process_dir
=
new_process
(
card_loc
=
'proc_card_mg5.dat'
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the param_card
# Build a new param_card.dat from an existing one
# Used values given in "parameters" for MHH and WHH, if not set there, default values are used
# Higgs mass is set to 125 GeV by "higgsMass"
#---------------------------------------------------------------------------------------------------
build_param_card
(
param_card_old
=
'param_card.HeavyScalar.dat'
,
param_card_new
=
'param_card_new.dat'
,
masses
=
higgsMass
,
extras
=
parameters
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the run_card
# 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
#---------------------------------------------------------------------------------------------------
build_run_card
(
run_card_old
=
get_default_runcard
(
proc_dir
=
process_dir
),
run_card_new
=
'run_card.dat'
,
xqcut
=
0
,
nevts
=
nevents
,
rand_seed
=
runArgs
.
randomSeed
,
beamEnergy
=
beamEnergy
,
extras
=
extras
)
print_cards
()
generate
(
run_card_loc
=
'run_card.dat'
,
param_card_loc
=
'param_card_new.dat'
,
mode
=
mode
,
proc_dir
=
process_dir
,
run_name
=
runName
)
outputDS
=
arrange_output
(
run_name
=
runName
,
proc_dir
=
process_dir
,
outputDS
=
runName
+
'._00001.events.tar.gz'
)
runArgs
.
inputGeneratorFile
=
outputDS
#--------------------------------------------------------------
# Herwig 7 (H7UE) showering setup
#--------------------------------------------------------------
# initialize Herwig7 generator configuration for showering of LHE files
include
(
"MC15JobOptions/Herwig7_LHEF.py"
)
# configure Herwig7
Herwig7Config
.
add_commands
(
"set /Herwig/Partons/RemnantDecayer:AllowTop Yes"
)
Herwig7Config
.
me_pdf_commands
(
order
=
"LO"
,
name
=
"NNPDF23_lo_as_0130_qed"
)
Herwig7Config
.
tune_commands
()
Herwig7Config
.
lhef_mg5amc_commands
(
lhe_filename
=
runArgs
.
inputGeneratorFile
,
me_pdf_order
=
"LO"
)
evgenConfig
.
tune
=
"H7.1-Default"
# add EvtGen
include
(
"MC15JobOptions/Herwig71_EvtGen.py"
)
#HW7 settings and Higgs BR
Herwig7Config
.
add_commands
(
"""
#set /Herwig/Shower/KinematicsReconstructor:ReconstructionOption General
#set /Herwig/Shower/KinematicsReconstructor:InitialInitialBoostOption LongTransBoost
set /Herwig/Shower/LtoLGammaSudakov:pTmin 0.000001
set /Herwig/Shower/QtoGammaQSudakov:Alpha /Herwig/Shower/AlphaQED
cd /Herwig/Particles
create /ThePEG/ParticleData XH
setup XH 1560 XH MASS WIDTH MAX_WIDTH 0 0 0 1 0
#set /Herwig/EventHandlers/LHEReader:Decayer /Herwig/Decays/Mambo
do /Herwig/Particles/h0:SelectDecayModes h0->gamma,gamma; h0->b,bbar;
do /Herwig/Particles/h0:PrintDecayModes
set /Herwig/Particles/h0/h0->gamma,gamma;:BranchingRatio 0.5
set /Herwig/Particles/h0/h0->b,bbar;:BranchingRatio 0.5
"""
)
# run Herwig7
Herwig7Config
.
run
()
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig
.
generators
=
[
"MadGraph"
,
"Herwig7"
,
"EvtGen"
]
if
(
runArgs
.
runNumber
==
451893
):
evgenConfig
.
description
=
"Di-Higgs production through 290 GeV Heavy Higgs resonance which decays to bbyy."
evgenConfig
.
keywords
=
[
"hh"
,
"BSM"
,
"BSMHiggs"
,
"resonance"
,
"diphoton"
,
"bottom"
]
evgenConfig
.
contact
=
[
'Florian Beisiegel <florian.beisiegel@cern.ch>'
]
evgenConfig
.
inputfilecheck
=
runName
#runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
#---------------------------------------------------------------------------------------------------
# Filter for bbyy
#---------------------------------------------------------------------------------------------------
from
GeneratorFilters.GeneratorFiltersConf
import
ParentChildFilter
filtSeq
+=
ParentChildFilter
(
"HbbFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
5
])
filtSeq
+=
ParentChildFilter
(
"HyyFilter"
,
PDGParent
=
[
25
],
PDGChild
=
[
22
])
filtSeq
.
Expression
=
"HbbFilter and HyyFilter"
share/DSID451xxx/MC15.451894.MadGraphHerwig7EvtGen_PDF23LO_X3125tohh_bbyy.py
0 → 100644
View file @
fe11ddbe
from
MadGraphControl.MadGraphUtils
import
*
mode
=
0
mhh
=
3.125000e+02
#---------------------------------------------------------------------------------------------------
# Setting mHH and WHH for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters
=
{
'1560'
:
'3.125000e+02'
,
#MHH
'1561'
:
'1.000000e-02'
}
#WHH
#---------------------------------------------------------------------------------------------------
# Setting higgs mass to 125 GeV for param_card.dat
#---------------------------------------------------------------------------------------------------
higgsMass
=
{
'25'
:
'1.250000e+02'
}
#MH
#---------------------------------------------------------------------------------------------------
# Setting some parameters for run_card.dat: Parton shower is HERWIG7
#---------------------------------------------------------------------------------------------------
extras
=
{
'lhe_version'
:
'2.0'
,
'cut_decays'
:
'F'
,
'pdlabel'
:
"'nn23lo1'"
,
'fixed_ren_scale'
:
'True'
,
#
'fixed_fac_scale'
:
'True'
,
#
'scale'
:
'156.25'
,
#
'dsqrt_q2fact1'
:
'156.25'
,
#
'dsqrt_q2fact2'
:
'156.25'
,
#
# 'parton_shower':'Herwig7',
'ptj'
:
'0'
,
'ptb'
:
'0'
,
'pta'
:
'0'
,
'ptjmax'
:
'-1'
,
'ptbmax'
:
'-1'
,
'ptamax'
:
'-1'
,
'etaj'
:
'-1'
,
'etab'
:
'-1'
,
'etaa'
:
'-1'
,
'etajmin'
:
'0'
,
'etabmin'
:
'0'
,
'etaamin'
:
'0'
,
'mmaa'
:
'0'
,
'mmaamax'
:
'-1'
,
'mmbb'
:
'0'
,
'mmbbmax'
:
'-1'
,
'drjj'
:
'0'
,
'drbb'
:
'0'
,
'draa'
:
'0'
,
'drbj'
:
'0'
,
'draj'
:
'0'
,
'drab'
:
'0'
,
'drjjmax'
:
'-1'
,
'drbbmax'
:
'-1'
,
'draamax'
:
'-1'
,
'drbjmax'
:
'-1'
,
'drajmax'
:
'-1'
,
'drabmax'
:
'-1'
}
#---------------------------------------------------------------------------------------------------
# Generating di-higgs through Heavy Higgs resonance with MadGraph
#---------------------------------------------------------------------------------------------------
fcard
=
open
(
'proc_card_mg5.dat'
,
'w'
)
if
(
runArgs
.
runNumber
==
451894
):
fcard
.
write
(
"""
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 HeavyHiggsTHDM
generate p p > hh > h h
output -f"""
)
fcard
.
write
(
'set Mhh %s
\n
'
%
(
mhh
))
fcard
.
close
()
else
:
raise
RuntimeError
(
"runNumber %i not recognised in these jobOptions."
%
runArgs
.
runNumber
)
beamEnergy
=-
999
if
hasattr
(
runArgs
,
'ecmEnergy'
):
beamEnergy
=
runArgs
.
ecmEnergy
/
2.
else
:
raise
RuntimeError
(
"No center of mass energy found."
)
#---------------------------------------------------------------------------------------------------
# Filter efficiency is low
# Thus, setting the number of generated events to 100 times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor
=
100
nevents
=
10000
*
safefactor
if
runArgs
.
maxEvents
>
0
:
nevents
=
runArgs
.
maxEvents
*
safefactor
runName
=
'run_01'
process_dir
=
new_process
(
card_loc
=
'proc_card_mg5.dat'
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the param_card
# Build a new param_card.dat from an existing one
# Used values given in "parameters" for MHH and WHH, if not set there, default values are used
# Higgs mass is set to 125 GeV by "higgsMass"
#---------------------------------------------------------------------------------------------------
build_param_card
(
param_card_old
=
'param_card.HeavyScalar.dat'
,
param_card_new
=
'param_card_new.dat'
,
masses
=
higgsMass
,
extras
=
parameters
)
#---------------------------------------------------------------------------------------------------
# Using the helper function from MadGraphControl for setting up the run_card
# 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
#---------------------------------------------------------------------------------------------------
build_run_card
(
run_card_old
=
get_default_runcard
(
proc_dir
=
process_dir
),
run_card_new
=
'run_card.dat'
,
xqcut
=
0
,
nevts
=
nevents
,
rand_seed
=
runArgs
.
randomSeed
,
beamEnergy
=
beamEnergy
,
extras
=
extras
)
print_cards
()
generate
(
run_card_loc
=
'run_card.dat'
,
param_card_loc
=
'param_card_new.dat'
,
mode
=
mode
,
proc_dir
=
process_dir
,
run_name
=
runName
)
outputDS
=
arrange_output
(
run_name
=
runName
,
proc_dir
=
process_dir
,
outputDS
=
runName
+
'._00001.events.tar.gz'
)
runArgs
.
inputGeneratorFile
=
outputDS
#--------------------------------------------------------------
# Herwig 7 (H7UE) showering setup
#--------------------------------------------------------------
# initialize Herwig7 generator configuration for showering of LHE files
include
(
"MC15JobOptions/Herwig7_LHEF.py"
)
# configure Herwig7
Herwig7Config
.
add_commands
(
"set /Herwig/Partons/RemnantDecayer:AllowTop Yes"
)
Herwig7Config
.
me_pdf_commands
(
order
=
"LO"
,
name
=
"NNPDF23_lo_as_0130_qed"
)
Herwig7Config
.
tune_commands
()
Herwig7Config
.
lhef_mg5amc_commands
(
lhe_filename
=
runArgs
.
inputGeneratorFile
,
me_pdf_order
=
"LO"
)
evgenConfig
.
tune
=
"H7.1-Default"
# add EvtGen
include
(
"MC15JobOptions/Herwig71_EvtGen.py"
)
#HW7 settings and Higgs BR
Herwig7Config
.
add_commands
(
"""
#set /Herwig/Shower/KinematicsReconstructor:ReconstructionOption General
#set /Herwig/Shower/KinematicsReconstructor:InitialInitialBoostOption LongTransBoost
set /Herwig/Shower/LtoLGammaSudakov:pTmin 0.000001
set /Herwig/Shower/QtoGammaQSudakov:Alpha /Herwig/Shower/AlphaQED
cd /Herwig/Particles
create /ThePEG/ParticleData XH
setup XH 1560 XH MASS WIDTH MAX_WIDTH 0 0 0 1 0
#set /Herwig/EventHandlers/LHEReader:Decayer /Herwig/Decays/Mambo
do /Herwig/Particles/h0:SelectDecayModes h0->gamma,gamma; h0->b,bbar;
do /Herwig/Particles/h0:PrintDecayModes
set /Herwig/Particles/h0/h0->gamma,gamma;:BranchingRatio 0.5
set /Herwig/Particles/h0/h0->b,bbar;:BranchingRatio 0.5
"""
)
# run Herwig7
Herwig7Config
.
run
()
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig
.
generators
=
[
"MadGraph"
,
"Herwig7"
,
"EvtGen"
]
if
(
runArgs
.
runNumber
==
451894
):
evgenConfig
.
description
=
"Di-Higgs production through 312.5 GeV Heavy Higgs resonance which decays to bbyy."
evgenConfig
.
keywords
=
[
"hh"
,
"BSM"
,
"BSMHiggs"
,
"resonance"
,
"diphoton"
,
"bottom"
]
evgenConfig
.
contact
=
[
'Florian Beisiegel <florian.beisiegel@cern.ch>'
]
evgenConfig
.
inputfilecheck
=
runName
#runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'