Commit 013e9008 authored by MC JO registration's avatar MC JO registration
Browse files

add 346330-1, 346635-8, 450869-73, 450792-826, 312374-403

parent 9e543120
2019-07-09 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* tagging MC15JobOptions-01-00-77
* add 346330-1, 346635-8, 450869-73, 450792-826, 312374-403
* update MadGraphControl_VLL.py
2019-07-09 E.M. Lobodzinska<ewelina2mail.desy.de>
* tag MC15JobOptions-01-00-76
* add 369108
* tag MC15JobOptions-01-00-76
* add 369108
2019-07-05 Joany Manjarres <joany@cern.ch>
* tag MC15JobOptions-01-00-75
......@@ -53,7 +58,7 @@
2019-07-01 E. M. Lobodzinska <ewelina@mail.desy.de>
* tag MC15JobOptions-01-00-65
* 450694-6
* 450694-6
2019-07-01 Thorsten Kuhl <Thorsten.Kuhl@cern.ch>
* tagging MC15JobOptions-01-00-64
......
......@@ -18,24 +18,47 @@ def taup_decay_width(m_taup):
def vtaup_decay_width(m_taup):
epsilon=0.04
mW=79.956
mW=80.385
rw=mW**2 / m_taup**2
decaywidth = m_taup * (1-rw)**2 * (2+1/rw) * epsilon**2 * rw / (32*math.pi)
return decaywidth
# Define constants
minevents=5000
mode=0
dsid_to_mtaup = { 999993 : 100,
999998 : 130,
999999 : 200,
310515 : 300,
999996 : 400,
999995 : 500,
999994 : 600,
}
dsid_to_mtaup_nlep = {
312403 : [ 130, 2 ],
312402 : [ 130, 3 ],
312401 : [ 130, 4 ],
312400 : [ 200, 2 ],
312399 : [ 200, 3 ],
312398 : [ 200, 4 ],
312397 : [ 300, 2 ],
312396 : [ 300, 3 ],
312395 : [ 300, 4 ],
312394 : [ 400, 2 ],
312393 : [ 400, 3 ],
312392 : [ 400, 4 ],
312391 : [ 500, 2 ],
312390 : [ 500, 3 ],
312389 : [ 500, 4 ],
312388 : [ 600, 2 ],
312387 : [ 600, 3 ],
312386 : [ 600, 4 ],
312385 : [ 700, 2 ],
312384 : [ 700, 3 ],
312383 : [ 700, 4 ],
312382 : [ 800, 2 ],
312381 : [ 800, 3 ],
312380 : [ 800, 4 ],
312379 : [ 900, 2 ],
312378 : [ 900, 3 ],
312377 : [ 900, 4 ],
312376 : [ 1000, 2 ],
312375 : [ 1000, 3 ],
312374 : [ 1000, 4 ],
}
# Get arguments
DSID=-1
......@@ -62,18 +85,15 @@ if hasattr(runArgs,'maxEvents'):
else:
nevents=5000
mtaup = dsid_to_mtaup.get(DSID)
if mtaup == None:
try:
mtaup = dsid_to_mtaup_nlep.get(DSID)[0]
nlep = dsid_to_mtaup_nlep.get(DSID)[1]
except (KeyError, IndexError) as e:
raise RuntimeError("VLL joboptions run with invalid runNumber: {}".format(DSID))
evt_multiplier = 3.2
if mtaup <= 100:
evt_multiplier = 15.0
elif mtaup <= 130:
evt_multiplier = 10.0
elif mtaup <= 300:
evt_multiplier = 5.0
evt_multiplier *= 1.05
evt_multiplier = 25.05
if mtaup <= 300:
evt_multiplier *= 2
# Write MG5 process card
proc_card = open('proc_card_mg5.dat','w')
......@@ -99,10 +119,10 @@ process_dir = new_process()
# Set up the run card with mostly default values
run_card_extras = {
'pdlabel' : 'lhapdf',
'lhaid' : 247000,
# 'pdlabel' : 'lhapdf',
# 'lhaid' : 247000,
'pdlabel' : 'nn23lo1',
'sys_pdf' : 'cteq6l1',
'event_norm' : 'sum',
'scale' : str(mtaup),
'dsqrt_q2fact1' : str(mtaup),
'dsqrt_q2fact2' : str(mtaup)
......@@ -147,20 +167,74 @@ if 'ATHENA_PROC_NUMBER' in os.environ:
opts.nprocs = 0
# Shower
evgenConfig.description = 'MadGraph5+Pythia8 for Vectorlike leptons of mass {}'.format(mtaup)
evgenConfig.description = 'MadGraph5+Pythia8 for Vectorlike leptons of mass {} GeV, {} lepton final state'.format(mtaup,nlep)
evgenConfig.contact = ['Daniel Wilbern <daniel.wilbern@cern.ch>']
evgenConfig.keywords += ['exotic', 'multilepton', 'tau']
evgenConfig.generators += ['MadGraph','Pythia8']
evgenConfig.inputfilecheck = output_DS_prefix
evgenConfig.minevents = minevents
runArgs.inputGeneratorFile=output_DS
include('MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py')
include('MC15JobOptions/Pythia8_MadGraph.py')
# Filter
include('MC15JobOptions/MultiLeptonFilter.py')
filtSeq.MultiLeptonFilter.Ptcut = 12000.
filtSeq.MultiLeptonFilter.Etacut = 2.7
filtSeq.MultiLeptonFilter.NLeptons = 2
# Filters
from GeneratorFilters.GeneratorFiltersConf import MultiLeptonWithParentFilter
if nlep == 2:
filtSeq += MultiLeptonWithParentFilter('LightFilter')
LightFilter = filtSeq.LightFilter
LightFilter.PDGOrigin = [9000005, 9000006]
LightFilter.PDGIntermediate = [23, 24, 25]
LightFilter.MinPt = 18000.
LightFilter.MaxEta = 2.8
LightFilter.NLeptonsMin = 2
LightFilter.NLeptonsMax = 2
LightFilter.IncludeLepTaus = True
LightFilter.IncludeHadTaus = False
LightFilter.VetoHadTaus = True
elif nlep == 3:
filtSeq += MultiLeptonWithParentFilter('LeptonFilter')
LeptonFilter = filtSeq.LeptonFilter
LeptonFilter.PDGOrigin = [9000005, 9000006]
LeptonFilter.PDGIntermediate = [23, 24, 25]
LeptonFilter.MinPt = 15000.
LeptonFilter.MaxEta = 2.8
LeptonFilter.NLeptonsMin = 3
LeptonFilter.NLeptonsMax = 3
LeptonFilter.IncludeLepTaus = True
LeptonFilter.IncludeHadTaus = True
LeptonFilter.VetoHadTaus = False
filtSeq += MultiLeptonWithParentFilter('LightFilter')
LightFilter = filtSeq.LightFilter
LightFilter.PDGOrigin = [9000005, 9000006]
LightFilter.PDGIntermediate = [23, 24, 25]
LightFilter.MinPt = 18000.
LightFilter.MaxEta = 2.8
LightFilter.NLeptonsMin = 2
LightFilter.NLeptonsMax = 3
LightFilter.IncludeLepTaus = True
LightFilter.IncludeHadTaus = False
LightFilter.VetoHadTaus = False
elif nlep == 4:
filtSeq += MultiLeptonWithParentFilter('LeptonFilter')
LeptonFilter = filtSeq.LeptonFilter
LeptonFilter.PDGOrigin = [9000005, 9000006]
LeptonFilter.PDGIntermediate = [23, 24, 25]
LeptonFilter.MinPt = 15000.
LeptonFilter.MaxEta = 2.8
LeptonFilter.NLeptonsMin = 4
LeptonFilter.NLeptonsMax = 10
LeptonFilter.IncludeLepTaus = True
LeptonFilter.IncludeHadTaus = True
LeptonFilter.VetoHadTaus = False
filtSeq += MultiLeptonWithParentFilter('LightFilter')
LightFilter = filtSeq.LightFilter
LightFilter.PDGOrigin = [9000005, 9000006]
LightFilter.PDGIntermediate = [23, 24, 25]
LightFilter.MinPt = 18000.
LightFilter.MaxEta = 2.8
LightFilter.NLeptonsMin = 2
LightFilter.NLeptonsMax = 10
LightFilter.IncludeLepTaus = True
LightFilter.IncludeHadTaus = False
LightFilter.VetoHadTaus = False
Supports Markdown
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