Commit 74447202 authored by Cyril Pascal Becot's avatar Cyril Pascal Becot
Browse files

See ChangeLog


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@814157 4525493e-7705-40b1-a816-d608a930855b
parent b21f39c5
2018-05-15 Cyril Becot <cyril.becot@cern.ch>
* add 406228-406235
* add 300460-300461
* up 395244-395258
* add susycontrol/param_card.SM.TT.charginoRPV.dat
* up 411034-411039
* add 396090-396091
* add 345912
* tag MC15JobOptions-00-09-93
2018-05-15 E. M. Lobodzinska <ewelina@mail.desy.de>
* fix merge_lhe.py obo Timothee Theveneaux-Pelzer - AGENE-1576
* Blacklist pcache affected by the skipevents problem -AGENE-1517
......
#######################################################################
# Job options for Pythia8B_i generation of Lambda_b->J/psi(mumu)Lambda
#######################################################################
evgenConfig.description = "Signal Bd->J/psi(mumu)Kstar0(Kpi)"
evgenConfig.keywords = ["exclusive","Jpsi","2muon","Bd"]
evgenConfig.minevents = 200
include("MC15JobOptions/nonStandard/Pythia8B_A14_CTEQ6L1_Common.py")
include('MC15JobOptions/nonStandard/Pythia8B_Photospp.py')
include("MC15JobOptions/Pythia8B_exclusiveB_Common.py")
genSeq.Pythia8B.Commands += ['PhaseSpace:pTHatMin = 11.']
genSeq.Pythia8B.QuarkPtCut = 0.0
genSeq.Pythia8B.AntiQuarkPtCut = 9.0
genSeq.Pythia8B.QuarkEtaCut = 102.5
genSeq.Pythia8B.AntiQuarkEtaCut = 2.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
genSeq.Pythia8B.Commands += ['511:addChannel = 2 1. 0 999443 313'] # onMode bRatio meMode products
genSeq.Pythia8B.SignalPDGCodes = [ 511, 999443, -13, 13, 313, 321, -211 ]
genSeq.Pythia8B.SignalPtCuts = [ 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0 ] # no muon eta cut here - see trigger cuts below
genSeq.Pythia8B.SignalEtaCuts = [ 102.5, 102.5, 102.5, 102.5, 102.5, 2.5, 2.5 ] # no muon eta cut here - see trigger cuts below
genSeq.Pythia8B.OutputLevel = INFO
genSeq.Pythia8B.NHadronizationLoops = 150
genSeq.Pythia8B.TriggerPDGCode = 13
genSeq.Pythia8B.TriggerStatePtCut = [4.0]
genSeq.Pythia8B.TriggerStateEtaCut = 2.5
genSeq.Pythia8B.TriggerOppositeCharges = True
genSeq.Pythia8B.MinimumCountPerCut = [2]
include('MC15JobOptions/BSignalFilter.py')
filtSeq.BSignalFilter.LVL1MuonCutOn = False
filtSeq.BSignalFilter.LVL2MuonCutOn = False
filtSeq.BSignalFilter.Cuts_Final_mu_switch = True
filtSeq.BSignalFilter.B_PDGCode = 511
filtSeq.BSignalFilter.InvMass_switch = True
filtSeq.BSignalFilter.InvMass_PartId1 = 321
filtSeq.BSignalFilter.InvMass_PartId2 = -211
filtSeq.BSignalFilter.InvMassMax = 1019.460 + 11.
filtSeq.BSignalFilter.InvMassMin = 1019.460 - 11.
filtSeq.BSignalFilter.TotalInvMass_switch = True
filtSeq.BSignalFilter.TotalInvMassMax = 5850.
filtSeq.BSignalFilter.TotalInvMassMin = 4950.
filtSeq.BSignalFilter.InvMass_PartFakeMass1 = 493.677
filtSeq.BSignalFilter.InvMass_PartFakeMass2 = 493.677
#######################################################################
# Job options for Pythia8B_i generation of Lambda_b->J/psi(mumu)Lambda
#######################################################################
evgenConfig.description = "Signal Bd->J/psi(mumu)Kpi"
evgenConfig.keywords = ["exclusive","Jpsi","2muon","Bd"]
evgenConfig.minevents = 200
include("MC15JobOptions/nonStandard/Pythia8B_A14_CTEQ6L1_Common.py")
include('MC15JobOptions/nonStandard/Pythia8B_Photospp.py')
include("MC15JobOptions/Pythia8B_exclusiveB_Common.py")
genSeq.Pythia8B.Commands += ['PhaseSpace:pTHatMin = 11.']
genSeq.Pythia8B.QuarkPtCut = 0.0
genSeq.Pythia8B.AntiQuarkPtCut = 9.0
genSeq.Pythia8B.QuarkEtaCut = 102.5
genSeq.Pythia8B.AntiQuarkEtaCut = 2.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
genSeq.Pythia8B.Commands += ['511:addChannel = 2 1. 0 999443 321 -211'] # onMode bRatio meMode products
genSeq.Pythia8B.SignalPDGCodes = [ 511, 999443, -13, 13, 321, -211 ]
genSeq.Pythia8B.SignalPtCuts = [ 0.0, 0.0, 0.0, 0.0, 1.0, 1.0 ] # no muon eta cut here - see trigger cuts below
genSeq.Pythia8B.SignalEtaCuts = [ 102.5, 102.5, 102.5, 102.5, 2.5, 2.5 ] # no muon eta cut here - see trigger cuts below
genSeq.Pythia8B.OutputLevel = INFO
genSeq.Pythia8B.NHadronizationLoops = 80
genSeq.Pythia8B.TriggerPDGCode = 13
genSeq.Pythia8B.TriggerStatePtCut = [4.0]
genSeq.Pythia8B.TriggerStateEtaCut = 2.5
genSeq.Pythia8B.TriggerOppositeCharges = True
genSeq.Pythia8B.MinimumCountPerCut = [2]
include('MC15JobOptions/BSignalFilter.py')
filtSeq.BSignalFilter.LVL1MuonCutOn = False
filtSeq.BSignalFilter.LVL2MuonCutOn = False
filtSeq.BSignalFilter.Cuts_Final_mu_switch = True
filtSeq.BSignalFilter.B_PDGCode = 511
filtSeq.BSignalFilter.InvMass_switch = True
filtSeq.BSignalFilter.InvMass_PartId1 = 321
filtSeq.BSignalFilter.InvMass_PartId2 = -211
filtSeq.BSignalFilter.InvMassMax = 1019.460 + 11.
filtSeq.BSignalFilter.InvMassMin = 1019.460 - 11.
filtSeq.BSignalFilter.TotalInvMass_switch = True
filtSeq.BSignalFilter.TotalInvMassMax = 5850.
filtSeq.BSignalFilter.TotalInvMassMin = 4950.
filtSeq.BSignalFilter.InvMass_PartFakeMass1 = 493.677
filtSeq.BSignalFilter.InvMass_PartFakeMass2 = 493.677
from MadGraphControl.MadGraphUtils import *
mode=0
mcid = 345912
#---------------------------------------------------------------------------------------------------
# Setting parameters for param_card.dat
#---------------------------------------------------------------------------------------------------
parameters = {}
parameters['NEW'] = {'CV': '1.0000000', # CV
'C2V': '1.0000000', # C2V
'C3': '2.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) which decays to bbtautau (hadhad)."
evgenConfig.keywords = ["BSM", "BSMHiggs", "nonResonant", "VBF", "bottom"]
evgenConfig.contact = ['Danilo Ferreira de Lima <dferreir@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 = [211,213,215,311,321,323,10232,10323,20213,20232,20323,30213,100213,100323,1000213]
filtSeq.XtoVVDecayFilterExtended.PDGChild2 = [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"
include('MC15JobOptions/MadGraphControl_SimplifiedModel_C1N2N1_GGMHino_Filter.py')
evgenConfig.minevents = 2000
evgenConfig.minevents = 1000
include('MC15JobOptions/MadGraphControl_SimplifiedModel_C1N2N1_GGMHino_Filter.py')
evgenConfig.minevents = 1000
include('MC15JobOptions/MadGraphControl_SimplifiedModel_C1N2N1_GGMHino_Filter.py')
evgenConfig.minevents = 1000
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
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