Skip to content

Commit

Permalink
Merge pull request #3649 from stahlleiton/POWHEG_TTbar_2023
Browse files Browse the repository at this point in the history
Add TTbar, single top, ZZ, WW and WZ POWHEG input cards for PbPb 5.36 TeV
  • Loading branch information
menglu21 authored Feb 28, 2024
2 parents 53af187 + dc4be50 commit b153f04
Show file tree
Hide file tree
Showing 7 changed files with 362 additions and 0 deletions.
Original file line number Diff line number Diff line change
@@ -0,0 +1,78 @@
! TTbar production parameters

numevts NEVENTS
iseed SEED
ih1 1 ! hadron 1
ih2 1 ! hadron 2


! To be set only if using LHA pdfs
lhans1 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208
lhans2 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208
ebeam1 2680 ! energy of beam 1
ebeam2 2680 ! energy of beam 2
qmass 172.5 ! mass of heavy quark in GeV
facscfact 1 ! factorization scale factor: mufact=muref*facscfact
renscfact 1 ! renormalization scale factor: muren=muref*renscfact

hdamp 237.8775

topdecaymode 22222 ! an integer of 5 digits that are either 0, or 2, representing in
! the order the maximum number of the following particles(antiparticles)
! in the final state: e mu tau up charm
! For example
! 22222 All decays (up to 2 units of everything)
! 20000 both top go into e l nu (with the appropriate signs)
! 10011 one top goes into electron (or positron), the other into (any) hadrons,
! or one top goes into charm, the other into up
! 00022 Fully hadronic
! 00002 Fully hadronic with two charms
! 00011 Fully hadronic with a single charm
! 00012 Fully hadronic with at least one charm


! Parameters for the generation of spin correlations in t tbar decays
tdec/wmass 80.377 ! https://pdg.lbl.gov/2023/reviews/rpp2023-rev-phys-constants.pdf
tdec/wwidth 2.085 ! https://pdg.lbl.gov/2023/tables/rpp2023-sum-gauge-higgs-bosons.pdf
tdec/bmass 4.78 ! ! https://pdg.lbl.gov/2023/tables/rpp2023-sum-quarks.pdf (pole mass)
tdec/twidth 1.330 ! Eq. (72.1) http://pdg.lbl.gov/2017/reviews/rpp2017-rev-top-quark.pdf
tdec/elbranching 0.1086 ! https://pdg.lbl.gov/2023/tables/rpp2023-sum-gauge-higgs-bosons.pdf
tdec/emass 0.000510998950 ! https://pdg.lbl.gov/2023/reviews/rpp2023-rev-phys-constants.pdf
tdec/mumass 0.1056583755 ! https://pdg.lbl.gov/2023/tables/rpp2023-sum-leptons.pdf
tdec/taumass 1.77686 ! https://pdg.lbl.gov/2023/tables/rpp2023-sum-leptons.pdf
tdec/dmass 0.100 !
tdec/umass 0.100 !
tdec/smass 0.200 !
tdec/cmass 1.67 ! https://pdg.lbl.gov/2023/tables/rpp2023-sum-quarks.pdf (pole mass)
tdec/sin2cabibbo 0.0506 ! 0.22500*0.22500, CKM_Vus in Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
! Parameters to allow-disallow use of stored data
use-old-grid 1 ! if 1 use old grid if file pwggrids.dat is present (# 1: regenerate)
use-old-ubound 1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; # 1: regenerate

ncall1 10000 ! number of calls for initializing the integration grid
itmx1 5 ! number of iterations for initializing the integration grid
ncall2 100000 ! number of calls for computing the integral and finding upper bound
itmx2 5 ! number of iterations for computing the integral and finding upper bound
foldcsi 1 ! number of folds on x integration
foldy 1 ! number of folds on y integration
foldphi 1 ! number of folds on phi integration
nubound 100000 ! number of bbarra calls to setup norm of upper bounding function
iymax 1 ! <= 10, normalization of upper bounding function in iunorm X iunorm square in y, log(m2qq)

xupbound 2 ! increase upper bound for radiation generation

pdfreweight 0
storeinfo_rwgt 0 ! store weight information
dampreweight 0 ! h_damp reweighting (mt/2, mt, mt*2)
withnegweights 1 ! default 0











Original file line number Diff line number Diff line change
@@ -0,0 +1,35 @@
! WW boson production parameters


numevts NEVENTS ! number of events to be generated
ih1 1 ! hadron 1 (1 for protons, -1 for antiprotons)
ih2 1 ! hadron 2 (1 for protons, -1 for antiprotons)
ebeam1 2680 ! energy of beam 1
ebeam2 2680 ! energy of beam 2
! To be set only if using LHA pdfs
lhans1 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208
lhans2 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208

! Parameters to allow or not the use of stored data
use-old-grid 1 ! if 1 use old grid if file pwggrids.dat is present (<> 1 regenerate)
use-old-ubound 1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; <> 1 regenerate
ncall1 400000 ! number of calls for initializing the integration grid
itmx1 5 ! number of iterations for initializing the integration grid
ncall2 1000000 ! number of calls for computing the integral and finding upper bound
itmx2 5 ! number of iterations for computing the integral and finding upper bound
foldcsi 1 ! number of folds on csi integration
foldy 1 ! number of folds on y integration
foldphi 1 ! number of folds on phi integration
nubound 1000000 ! number of bbarra calls to setup norm of upper bounding function

! OPTIONAL PARAMETERS
leptonic 1
renscfact 1 ! (default 1d0) ren scale factor: muren = muref * renscfact
facscfact 1 ! (default 1d0) fac scale factor: mufact = muref * facscfact

iseed SEED ! initialize random number sequence

pdfreweight 0 ! PDF reweighting
storeinfo_rwgt 0
withnegweights 1

Original file line number Diff line number Diff line change
@@ -0,0 +1,35 @@
! WZ boson production parameters
mllmin 4d0 ! default 0.1 GeV this is minimum invar mass for Z leptons

numevts NEVENTS
ih1 1 ! hadron 1 (1 for protons, -1 for antiprotons)
ih2 1 ! hadron 2 (1 for protons, -1 for antiprotons)
ebeam1 2680d0 ! energy of beam 1
ebeam2 2680d0 ! energy of beam 2
! To be set only if using LHA pdfs
lhans1 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208
lhans2 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208

! Parameters to allow or not the use of stored data
use-old-grid 1 ! if 1 use old grid if file pwggrids.dat is present (<> 1 regenerate)
use-old-ubound 1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; <> 1 regenerate
ncall1 400000 ! number of calls for initializing the integration grid
itmx1 4 ! number of iterations for initializing the integration grid
ncall2 1000000 ! number of calls for computing the integral and finding upper bound
itmx2 5 ! number of iterations for computing the integral and finding upper bound
foldcsi 1 ! number of folds on csi integration
foldy 1 ! number of folds on y integration
foldphi 1 ! number of folds on phi integration
nubound 1000000 ! number of bbarra calls to setup norm of upper bounding function

! OPTIONAL PARAMETERS
lll 1

renscfact 1d0 ! (default 1d0) ren scale factor: muren = muref * renscfact
facscfact 1d0 ! (default 1d0) fac scale factor: mufact = muref * facscfact

iseed SEED

pdfreweight 0
storeinfo_rwgt 0
withnegweights 1
Original file line number Diff line number Diff line change
@@ -0,0 +1,34 @@
! ZZ boson production parameters
mllmin 4d0 ! default 0.1 GeV this is minimum invar mass for Z leptons

numevts NEVENTS ! number of events to be generated
ih1 1 ! hadron 1 (1 for protons, -1 for antiprotons)
ih2 1 ! hadron 2 (1 for protons, -1 for antiprotons)
ebeam1 2680d0 ! energy of beam 1
ebeam2 2680d0 ! energy of beam 2
! To be set only if using LHA pdfs
lhans1 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208
lhans2 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208

! Parameters to allow or not the use of stored data
use-old-grid 1 ! if 1 use old grid if file pwggrids.dat is present (<> 1 regenerate)
use-old-ubound 1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; <> 1 regenerate
ncall1 1000000 ! number of calls for initializing the integration grid
itmx1 4 ! number of iterations for initializing the integration grid
ncall2 1000000 ! number of calls for computing the integral and finding upper bounds
itmx2 5 ! number of iterations for computing the integral and finding upper bound
foldcsi 2 ! number of folds on csi integration
foldy 5 ! number of folds on y integration
foldphi 2 ! number of folds on phi integration
nubound 1000000 ! number of bbarra calls to setup norm of upper bounding function

! OPTIONAL PARAMETERS
leptons-nu 1
renscfact 1d0 ! (default 1d0) ren scale factor: muren = muref * renscfact
facscfact 1d0 ! (default 1d0) fac scale factor: mufact = muref * facscfact

iseed SEED ! initialize random number sequence

pdfreweight 0 ! PDF reweighting
storeinfo_rwgt 0
withnegweights 1
Original file line number Diff line number Diff line change
@@ -0,0 +1,34 @@
! ZZ boson production parameters
mllmin 4d0 ! default 0.1 GeV this is minimum invar mass for Z leptons

numevts NEVENTS ! number of events to be generated
ih1 1 ! hadron 1 (1 for protons, -1 for antiprotons)
ih2 1 ! hadron 2 (1 for protons, -1 for antiprotons)
ebeam1 2680d0 ! energy of beam 1
ebeam2 2680d0 ! energy of beam 2
! To be set only if using LHA pdfs
lhans1 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208
lhans2 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208

! Parameters to allow or not the use of stored data
use-old-grid 1 ! if 1 use old grid if file pwggrids.dat is present (<> 1 regenerate)
use-old-ubound 1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; <> 1 regenerate
ncall1 1000000 ! number of calls for initializing the integration grid
itmx1 4 ! number of iterations for initializing the integration grid
ncall2 1000000 ! number of calls for computing the integral and finding upper bounds
itmx2 5 ! number of iterations for computing the integral and finding upper bound
foldcsi 2 ! number of folds on csi integration
foldy 5 ! number of folds on y integration
foldphi 2 ! number of folds on phi integration
nubound 1000000 ! number of bbarra calls to setup norm of upper bounding function

! OPTIONAL PARAMETERS
leptonic 1
renscfact 1d0 ! (default 1d0) ren scale factor: muren = muref * renscfact
facscfact 1d0 ! (default 1d0) fac scale factor: mufact = muref * facscfact

iseed SEED ! initialize random number sequence

pdfreweight 0 ! PDF reweighting
storeinfo_rwgt 0
withnegweights 1
Original file line number Diff line number Diff line change
@@ -0,0 +1,73 @@
! ST-wtchannel production parameters

! GENERAL POWHEG PARAMETERS

numevts NEVENTS ! number of events to be generated
ih1 1 ! hadron 1 (1 for protons, -1 for antiprotons)
ih2 1 ! hadron 2 (1 for protons, -1 for antiprotons)
ebeam1 2680d0 ! energy of beam 1
ebeam2 2680d0 ! energy of beam 2

iseed SEED ! initialize random number sequence

! To be set only if using LHA pdfs
lhans1 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208
lhans2 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208

! Parameters to allow or not the use of stored data
use-old-grid 1 ! if 1 use old grid if file pwggrids.dat is present (<> 1 regenerate)
use-old-ubound 1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; <> 1 regenerate

ncall1 50000 ! number of calls for initializing the integration grid
itmx1 5 ! number of iterations for initializing the integration grid
ncall2 50000 ! number of calls for computing the integral and finding upper bound
itmx2 5 ! number of iterations for computing the integral and finding upper bound
foldcsi 1 ! number of folds on csi integration
foldy 1 ! number of folds on y integration
foldphi 1 ! number of folds on phi integration
nubound 50000 ! number of bbarra calls to setup norm of upper bounding function
xupbound 2d0 ! increase upper bound for radiation generation


withdamp 1
hdamp 237.8775


! mandatory production parameters
ttype -1 ! 1 for t, -1 for tbar


! mandatory parameters used in decay generation
topdecaymode 11111 ! decay mode: the 5 digits correspond to the following
! top-decay channels (l,mu,tau,u,c)
! 0 means close, 1 open
wdecaymode 11111 ! decay mode: the 5 digits correspond to the following
! primary-w-decay channels (l,mu,tau,u,c)
! 0 means close, 1 open
tdec/elbranching 0.1086 ! https://pdg.lbl.gov/2023/tables/rpp2023-sum-gauge-higgs-bosons.pdf

! optional production parameters
! (defaults defined in init_couplings.f)
topmass 172.5 ! https://pdg.lbl.gov/2023/tables/rpp2023-sum-quarks.pdf (pole mass)
wmass 80.377 ! https://pdg.lbl.gov/2023/reviews/rpp2023-rev-phys-constants.pdf
sthw2 0.23121 ! https://pdg.lbl.gov/2023/reviews/rpp2023-rev-phys-constants.pdf
alphaem_inv 137.035999084 ! https://pdg.lbl.gov/2023/reviews/rpp2023-rev-phys-constants.pdf



CKM_Vud 0.97435 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vus 0.22500 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vub 0.00369 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vcd 0.22486 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vcs 0.97349 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vcb 0.04182 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vtd 0.00857 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vts 0.04110 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vtb 0.999118 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf



pdfreweight 0 ! PDF reweighting
dampreweight 0 ! h_damp reweighting (mt/2, mt, mt*2)
storeinfo_rwgt 0 ! store weight information
withnegweights 1 ! default
Original file line number Diff line number Diff line change
@@ -0,0 +1,73 @@
! ST-wtchannel production parameters

! GENERAL POWHEG PARAMETERS

numevts NEVENTS ! number of events to be generated
ih1 1 ! hadron 1 (1 for protons, -1 for antiprotons)
ih2 1 ! hadron 2 (1 for protons, -1 for antiprotons)
ebeam1 2680d0 ! energy of beam 1
ebeam2 2680d0 ! energy of beam 2

iseed SEED ! initialize random number sequence

! To be set only if using LHA pdfs
lhans1 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208
lhans2 904400 ! nPDF: EPPS21nlo_CT18Anlo_Pb208

! Parameters to allow or not the use of stored data
use-old-grid 1 ! if 1 use old grid if file pwggrids.dat is present (<> 1 regenerate)
use-old-ubound 1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; <> 1 regenerate

ncall1 50000 ! number of calls for initializing the integration grid
itmx1 5 ! number of iterations for initializing the integration grid
ncall2 50000 ! number of calls for computing the integral and finding upper bound
itmx2 5 ! number of iterations for computing the integral and finding upper bound
foldcsi 1 ! number of folds on csi integration
foldy 1 ! number of folds on y integration
foldphi 1 ! number of folds on phi integration
nubound 50000 ! number of bbarra calls to setup norm of upper bounding function
xupbound 2d0 ! increase upper bound for radiation generation


withdamp 1
hdamp 237.8775


! mandatory production parameters
ttype 1 ! 1 for t, -1 for tbar


! mandatory parameters used in decay generation
topdecaymode 11111 ! decay mode: the 5 digits correspond to the following
! top-decay channels (l,mu,tau,u,c)
! 0 means close, 1 open
wdecaymode 11111 ! decay mode: the 5 digits correspond to the following
! primary-w-decay channels (l,mu,tau,u,c)
! 0 means close, 1 open
tdec/elbranching 0.1086 ! https://pdg.lbl.gov/2023/tables/rpp2023-sum-gauge-higgs-bosons.pdf

! optional production parameters
! (defaults defined in init_couplings.f)
topmass 172.5 ! https://pdg.lbl.gov/2023/tables/rpp2023-sum-quarks.pdf (pole mass)
wmass 80.377 ! https://pdg.lbl.gov/2023/reviews/rpp2023-rev-phys-constants.pdf
sthw2 0.23121 ! https://pdg.lbl.gov/2023/reviews/rpp2023-rev-phys-constants.pdf
alphaem_inv 137.035999084 ! https://pdg.lbl.gov/2023/reviews/rpp2023-rev-phys-constants.pdf



CKM_Vud 0.97435 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vus 0.22500 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vub 0.00369 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vcd 0.22486 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vcs 0.97349 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vcb 0.04182 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vtd 0.00857 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vts 0.04110 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf
CKM_Vtb 0.999118 ! Table (12.27) https://pdg.lbl.gov/2023/reviews/rpp2023-rev-ckm-matrix.pdf



pdfreweight 0 ! PDF reweighting
dampreweight 0 ! h_damp reweighting (mt/2, mt, mt*2)
storeinfo_rwgt 0 ! store weight information
withnegweights 1 ! default

0 comments on commit b153f04

Please sign in to comment.