MC15.345060.PowhegPythia8EvtGen_Hto4l_CT10_AZNLOCTEQ6L1_ggH125.py

Go to the documentation of this file.

#--------------------------------------------------------------
# Powheg ggH_quark_mass_effects setup
#--------------------------------------------------------------
 
# Must save input gen file names before it is overwritten
input_powheg_files = runArgs.inputGeneratorFile   
 
 
include('PowhegControl/PowhegControl_ggF_H_Common.py')
 
 
# Temporarily, we setup Powheg to use some variables to be the same in the Hto4l config. These
# should eventually be taken from the incoming LHE file
 
# Set Powheg variables, overriding defaults (RDS: defaults are ok. 2019/02)
# Note: width_H will be overwritten in case of CPS.
# PowhegConfig.mass.H  = 125.
# PowhegConfig.width.H = 0.00407
 
# PowhegConfig.generateRunCard()
# PowhegConfig.generateEvents()
# PowhegConfig.generate(run_card_only=True)
PowhegConfig.generate(create_run_card_only=True)
 
#--------------------------------------------------------------
# configure Hto4l
#--------------------------------------------------------------
include('Hto4lControl/Hto4lPowhegCommon.py')
 
# set name of LHE file - used below in merger
 
print "runArgs", runArgs
 
# powhegInputFile                          = "TXT.10248297._000705.events"   # for interactive tests
# powhegInputFile                          = runArgs.inputGeneratorFile    # for prod jobs
# powhegInputFile                          = "XXX.TTT" # runArgs.inputGeneratorFile    # for prod jobs
 
# runargs.maxEvents number of generated events
 
 
#--------------------------------------------------------------
# configure Hto4lPowhegMerger
#--------------------------------------------------------------
 
include('Hto4lControl/Hto4lPowhegMergeCommon.py')
 
# input powheg list of files comes from the input generator file names
hto4lPowhegMerger.input_powheg_file_name = input_powheg_files
 
# setup inputs for the individual Hto4l jobs - argument must correspond to the number of files per
# job on the grid submission:
hto4lPowhegMerger.setUpInputs(8)
 
# Generate a file for 4e
hto4lConfig4e.channel      = "eeee"
hto4lConfig4e.lheInput     = hto4lPowhegMerger.input_powheg_to_hto4l4e_file_name
hto4lConfig4e.lheOutput    = hto4lPowhegMerger.output_hto4l4e_file_name
hto4lConfig4e.nEvents      = int( runArgs.maxEvents / 4 + 0.5 ) # 1/4 for this channel
hto4lConfig4e.maxDiffWidth = 3.2052962671115008E-002 # from first run with value set to 0
hto4lConfig4e.random_seed  = 123*runArgs.randomSeed # separate random seed for 4e
 
hto4lConfig4e.generateRunCard()
hto4lConfig4e.generateEvents()
 
# Generate a file for 4mu
hto4lConfig4mu.channel      = "mmmm"
hto4lConfig4mu.lheInput     = hto4lPowhegMerger.input_powheg_to_hto4l4mu_file_name
hto4lConfig4mu.lheOutput    = hto4lPowhegMerger.output_hto4l4mu_file_name
hto4lConfig4mu.nEvents      = int( runArgs.maxEvents / 4 + 0.5 ) # 1/4 for this channel
hto4lConfig4mu.maxDiffWidth = 1.3396654719615330E-002
hto4lConfig4mu.random_seed  = 368*runArgs.randomSeed # separate random seed for 4mu
 
hto4lConfig4mu.generateRunCard()
hto4lConfig4mu.generateEvents()
 
# Generate a file for 2e2mu
hto4lConfig2e2mu.channel      = "eemm"
hto4lConfig2e2mu.lheInput     = hto4lPowhegMerger.input_powheg_to_hto4l2e2mu_file_name
hto4lConfig2e2mu.lheOutput    = hto4lPowhegMerger.output_hto4l2e2mu_file_name
hto4lConfig2e2mu.nEvents      = int( 1.1*(runArgs.maxEvents / 2) ) # 1/2 for this channel + extra 10% for events failing pythia checks
hto4lConfig2e2mu.maxDiffWidth = 0.12975814265233268
hto4lConfig2e2mu.random_seed  = 278*runArgs.randomSeed # separate random seed for 2e2mu
 
hto4lConfig2e2mu.generateRunCard()
hto4lConfig2e2mu.generateEvents()
 
 
#--------------------------------------------------------------
# run final LHE file merging with Hto4lPowhegMerger
#--------------------------------------------------------------
 
hto4lPowhegMerger.merge()
 
 
print "runArgs", runArgs
 
#--------------------------------------------------------------
# Pythia8 showering
# note: Main31 is set in Pythia8_AZNLO_CTEQ6L1_EvtGen_Common.py
#--------------------------------------------------------------
include('MC15JobOptions/Pythia8_AZNLO_CTEQ6L1_EvtGen_Common.py')
 
#--------------------------------------------------------------
# Pythia8 main31 update
#--------------------------------------------------------------
# genSeq.Pythia8.UserModes += [ 'Main31:NFinal = 2' ]
 
if "UserHooks" in genSeq.Pythia8.__slots__.keys():
 
  genSeq.Pythia8.Commands  += ['Powheg:NFinal = 2']
 
else:
 
  genSeq.Pythia8.UserModes += [ 'Main31:NFinal = 2' ]
 
#--------------------------------------------------------------
# Higgs at Pythia8
#--------------------------------------------------------------
genSeq.Pythia8.Commands += ["TimeShower:QEDshowerByL = off",
                            "TimeShower:QEDshowerByGamma = off",
                            '25:onMode = off',     # decay of Higgs
                            '23:onMode = off' ]    # decay of Z
 
 
# genSeq.Pythia8.Commands += ["TimeShower:QEDshowerByOther = off"]
 
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = "POWHEG+HTO4L+PYTHIA8+EVTGEN, ggH H->llll+X mh=125 GeV CPS"
evgenConfig.keywords    = [ "Higgs", "SMHiggs", "ZZ" ]
evgenConfig.contact     = [ 'R.D.Schaffer@cern.ch' ]