Commit 752e70ed authored by Giovanni Siragusa's avatar Giovanni Siragusa
Browse files

Update 341004


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@663044 4525493e-7705-40b1-a816-d608a930855b
parent b5283691
2015-04-25 Giovanni Siragusa <siragusa@cern.ch>
* Update 341004
2015-04-24 James Monk <jmonk@cern.ch>
* Add EvtGen, Tauola and Photos fragments for Herwig/Jimmy
* Remove unneeded PowhegJimmy_AUET2_CT10_Common.py and Jimmy_AUET2_CT10_Fragment
......
......@@ -2,10 +2,20 @@ from MadGraphControl.MadGraphUtils import *
mode=0
parameters={'1560':'3.000000e+02',
'1561':'1.000000e-02'}
#---------------------------------------------------------------------------------------------------
# Setting mHH and WHH for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters={'1560':'3.000000e+02', #MHH
'1561':'1.000000e-02'} #WHH
#Defaults for run_card.dat
#---------------------------------------------------------------------------------------------------
# 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'",
......@@ -42,7 +52,9 @@ extras = { 'lhe_version':'2.0',
'drajmax':'-1',
'drabmax':'-1' }
#---------------------------------------------------------------------------------------------------
# Generating di-higgs through Heavy Higgs resonance with MadGraph
#---------------------------------------------------------------------------------------------------
fcard = open('proc_card_mg5.dat','w')
if (runArgs.runNumber == 341004):
fcard.write("""
......@@ -62,10 +74,35 @@ if hasattr(runArgs,'ecmEnergy'):
else:
raise RuntimeError("No center of mass energy found.")
#---------------------------------------------------------------------------------------------------
# Filter efficiency is ~50%
# Thus, setting the number of generated events to 3 times maxEvents,
# to avoid crashing due to not having enough events
# Also putting protection to avoid from crashing when maxEvents=-1
#---------------------------------------------------------------------------------------------------
safefactor=3
nevents=5000*safefactor
if runArgs.maxEvents > 0:
nevents=runArgs.maxEvents*safefactor
#---------------------------------------------------------------------------------------------------
# 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(),run_card_new='run_card.dat',
nevts=runArgs.maxEvents*3,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,extras=extras)
nevts=nevents,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,extras=extras)
build_param_card(param_card_old='param_card.HeavyScalar.dat',param_card_new='param_card_new.dat',extras=parameters)
#---------------------------------------------------------------------------------------------------
# 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)
print_cards()
......@@ -76,23 +113,34 @@ generate(run_card_loc='run_card.dat',param_card_loc='param_card_new.dat',mode=mo
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz')
### Shower
#---------------------------------------------------------------------------------------------------
# Pythia8 Showering with A14_NNPDF23LO
#---------------------------------------------------------------------------------------------------
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
include("MC15JobOptions/Pythia8_LHEF.py")
#---------------------------------------------------------------------------------------------------
# EVGEN Configuration
#---------------------------------------------------------------------------------------------------
evgenConfig.generators = ["MadGraph", "Pythia8", "EvtGen"]
if (runArgs.runNumber == 341004):
evgenConfig.description = "diHiggs production through 260GeV Heavy Higgs resonance, decay to yybb, with Madgraph."
evgenConfig.description = "Di-Higgs production through 300 GeV Heavy Higgs resonance which decays to yybb."
evgenConfig.keywords = ["BSM", "BSMHiggs", "resonance", "diphoton", "bottom"]
evgenConfig.contact = ['Tulin Varol']
evgenConfig.contact = ['Tulin Varol <Tulin.Varol@cern.ch>']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
#---------------------------------------------------------------------------------------------------
# Decaying hh to bbyy with Pythia8
#---------------------------------------------------------------------------------------------------
genSeq.Pythia8.Commands += ["25:oneChannel = on 0.5 100 5 -5 ", #bb decay
"25:addChannel = on 0.5 100 22 22 "] #gammagamma decay
#---------------------------------------------------------------------------------------------------
# Generator Filters
include("XtoVVDecayFilterExtended.py")
#---------------------------------------------------------------------------------------------------
include("MC15JobOptions/XtoVVDecayFilterExtended.py")
filtSeq.XtoVVDecayFilterExtended.PDGGrandParent = 1560
filtSeq.XtoVVDecayFilterExtended.PDGParent = 25
filtSeq.XtoVVDecayFilterExtended.StatusParent = 22
......
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