mc.PhPy8_A14_NNPDF31_ttbb_4FS_hd15mt_dilep_v2_valid.py 2.78 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
evgenConfig.description = "POWHEG-BOX-RES/OpenLoops+Pythia8+EvtGen ttbb (4FS), mur=1/2*[mT(top)*mT(tbar)*mT(b)*mT(bbar)]**(1/4), muf=1/2*[mT(top)+mT(tbar)+mT(b)+mT(bbar)+mT(gluon)], dileptonic channel, hdamp=1.5 mTop, A14 NNPDF23 LO tune, decays with Powheg"
evgenConfig.keywords = [ 'SM', 'top', 'ttbar', 'bbbar', '2lepton']
evgenConfig.contact = ["lars.ferencz@cern.ch"]
evgenConfig.generators += [ 'Powheg' ]
evgenConfig.nEventsPerJob = 1000

#--------------------------------------------------------------
# Powheg
#--------------------------------------------------------------
# Load ATLAS defaults for the Powheg ttbb process
include('PowhegControl/PowhegControl_ttbb_Common.py')

### integration was done with ATHENA_PROC_NUMBER=10
PowhegConfig.ncall1       = "400000"
PowhegConfig.ncall2       = "400000"
PowhegConfig.ncall2rm     = "600000"
PowhegConfig.nubound      = "100000"

### muR/muF recommendation from LHC Higgs WG
PowhegConfig.runningscales = 3

### bornzerodamp
PowhegConfig.bornzerodamp = 5 

### hdamp settings
PowhegConfig.hdamp      = 1.5*172.5  # fixed hdamp off
PowhegConfig.dynhdamp   = 0   # 0=off, 1=on
PowhegConfig.dynhdampPF = 0. # prefactor of dynhdamp, typically 0.25, 0.5 or 1.0

### Fold parameter reducing the negative eventweight fraction
PowhegConfig.foldcsi      = "5"
PowhegConfig.foldphi      = "5"
PowhegConfig.foldy        = "5"
PowhegConfig.for_reweighting = 1

### decay handled by PowHeg
PowhegConfig.decay_mode   = "t t~ > b l+ vl b~ l- vl~"

#Dynamical scale variations
PowhegConfig.mu_F         = [ 1.0, 2.0, 0.5, 1.0, 1.0, 0.5, 2.0, 0.5, 2.0, 0.25, 0.5,  0.25, 0.25, 1.0  ] # List of factorisation scales which pairs with renormalisation scale below
PowhegConfig.mu_R         = [ 1.0, 1.0, 1.0, 2.0, 0.5, 0.5, 2.0, 2.0, 0.5, 0.5,  0.25, 0.25, 1.0,  0.25 ] # List of renormalisation scales


### Set specific PDF set
PowhegConfig.PDF          = [320900, 260400, 13191, 25410, 92000, 265400, 266400, 320500 ]


# --------------------------------------------------------------
# Generate events
# --------------------------------------------------------------
PowhegConfig.generate()

#--------------------------------------------------------------
# Pythia8 showering with the A14 NNPDF2.3 tune
#--------------------------------------------------------------
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("Pythia8_i/Pythia8_Powheg_Main31.py")
### Change pTdef value
genSeq.Pythia8.Commands += [ 'Powheg:pTdef = 1' ] 
genSeq.Pythia8.Commands += [ 'Powheg:pTHard = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 4' ]
genSeq.Pythia8.Commands += [ 'Powheg:veto = 1' ]
genSeq.Pythia8.Commands += [ 'Powheg:vetoCount = 3' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTemt  = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:emitted = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:MPIveto = 0' ]