Commit 23506a03 authored by Franziska Iltzsche's avatar Franziska Iltzsche
Browse files

update 345069 - 345977, add 346030


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@813807 4525493e-7705-40b1-a816-d608a930855b
parent 169e28e3
2018-04-19 Franziska Iltzsche <franziska.speiser@cern.ch>
* update 345069 - 345977
* add 346030
2018-04-18 Josh McFayden <mcfayden@cern.ch>
* tagging MC15JobOptions-00-09-63
* Fix for csv file for 345940-2
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 999999
mcid = 345968
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 999999
mcid = 345969
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 999999
mcid = 345970
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 999999
mcid = 345971
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 999999
mcid = 345972
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 999999
mcid = 345973
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 999999
mcid = 345974
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 999999
mcid = 345975
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 999999
mcid = 345976
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 999999
mcid = 345977
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
......
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 346030
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters = {}
parameters['NEW'] = {'CV': '1.0000000', # CV
'C2V': '0.0000000', # C2V
'C3': '0.0000000'} # C3
#---------------------------------------------------------------------------------------------------
# Setting higgs mass to 125 GeV for param_card.dat
#---------------------------------------------------------------------------------------------------
higgsMass={'25':'1.250000e+02'} #MH
#---------------------------------------------------------------------------------------------------
# Setting some parameters for run_card.dat
#---------------------------------------------------------------------------------------------------
extras = { 'lhe_version':'2.0',
'cut_decays':'F',
'pdlabel':"'nn23lo1'",
'scale':'125',
'dsqrt_q2fact1':'125',
'dsqrt_q2fact2':'125'}
#---------------------------------------------------------------------------------------------------
# Generating di-higgs through Heavy Higgs resonance with MadGraph
# Due to the parameters settings above, the resonant contribution is eliminated
# Thus, the non-resonant hh production with different trilinear Higgs couplings is enabled
#---------------------------------------------------------------------------------------------------
fcard = open('proc_card_mg5.dat','w')
if (runArgs.runNumber == mcid):
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 /cvmfs/atlas.cern.ch/repo/sw/Generators/madgraph/models/latest/HHVBF_UFO
generate p p > h h j j QCD=0
output -f""")
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 ~1%
# 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=30
nevents=5000*safefactor
if runArgs.maxEvents > 0:
nevents=runArgs.maxEvents*safefactor
process_dir = new_process()
#---------------------------------------------------------------------------------------------------
# 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',
nevts=nevents,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,xqcut=0,extras=extras)
#---------------------------------------------------------------------------------------------------
# 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='%s/Cards/param_card.dat' % process_dir,param_card_new='param_card_new.dat',masses=higgsMass,params=parameters)
print_cards()
runName='run_01'
#process_dir = new_process()
generate(run_card_loc='run_card.dat',param_card_loc='param_card_new.dat',mode=mode,proc_dir=process_dir,run_name=runName)
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz')
#---------------------------------------------------------------------------------------------------
# Pythia8 Showering with A14_NNPDF23LO
#---------------------------------------------------------------------------------------------------
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_LHEF.py")
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig.generators = ["MadGraph", "Pythia8", "EvtGen"]
if (runArgs.runNumber == mcid):
evgenConfig.description = "Non-resonant di-Higgs production through vector-boson-fusion (VBF) with trilinear coupling of lambda01 (SMx1) which decays to bbtautau (lephad)."
evgenConfig.keywords = ["BSM", "BSMHiggs", "nonResonant", "VBF", "bottom"]
evgenConfig.contact = ['Danilo Ferreira de Lima <dferreir@cern.ch>', ' Puja Saha <puja.saha@cern.ch>']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
#---------------------------------------------------------------------------------------------------
# Decaying hh to bbtautau with Pythia8
# Preferred way to select particular decay channels without modifying the Higgs
# width or cross section is to turn off all decays but those of interest, not to
# modify the branching ratios
#---------------------------------------------------------------------------------------------------
genSeq.Pythia8.Commands += [ "25:onMode=off", # turn off decay modes
"25:onIfAny=5 15" ] # allow modes with b or tau
#---------------------------------------------------------------------------------------------------
# Generator Filters
#---------------------------------------------------------------------------------------------------
from GeneratorFilters.GeneratorFiltersConf import ParentChildFilter
filtSeq += ParentChildFilter("hbbFilter", PDGParent = [25], PDGChild = [5])
filtSeq += ParentChildFilter("hTauTauFilter", PDGParent = [25], PDGChild = [15])
include("MC15JobOptions/XtoVVDecayFilterExtended.py")
filtSeq.XtoVVDecayFilterExtended.PDGGrandParent = 25
filtSeq.XtoVVDecayFilterExtended.PDGParent = 15
filtSeq.XtoVVDecayFilterExtended.StatusParent = 2
filtSeq.XtoVVDecayFilterExtended.PDGChild1 = [11, 13]
filtSeq.XtoVVDecayFilterExtended.PDGChild2 = [24,211,213,215,311,321,323,10232,10323,20213,20232,20323,30213,100213,100323,1000213]
#filtSeq.XtoVVDecayFilterExtended.PDGChild2 = [111,130,211,221,223,310,311,321,323]
#---------------------------------------------------------------------------------------------------
# Filter for 2 leptons (inc tau(had)) with pt cuts on e/mu and tau(had)
#---------------------------------------------------------------------------------------------------
from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
filtSeq += MultiElecMuTauFilter("TauPtFilter")
filtSeq.TauPtFilter.IncludeHadTaus = True
filtSeq.TauPtFilter.NLeptons = 2
filtSeq.TauPtFilter.MinPt = 13000.
filtSeq.TauPtFilter.MinVisPtHadTau = 15000.
filtSeq.TauPtFilter.MaxEta = 5.
filtSeq.Expression = "hbbFilter and hTauTauFilter and XtoVVDecayFilterExtended and TauPtFilter"
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment