diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_customizecards.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_customizecards.dat new file mode 100644 index 000000000000..f0ee62fa058f --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_customizecards.dat @@ -0,0 +1,4 @@ +#put card customizations here +set param_card mass 6 172.5 +set param_card yukawa 6 172.5 +set param_card DECAY 6 auto diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_extramodels.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_extramodels.dat new file mode 100644 index 000000000000..a2a62b1c4be8 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_extramodels.dat @@ -0,0 +1,2 @@ +#Anomalous Top couplings LO model version 1.4 from https://twiki.cern.ch/twiki/bin/view/CMS/TopFCNCmodel +TopFCNC_patch4_UFO.tgz diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_proc_card.dat new file mode 100644 index 000000000000..f015f9ca46db --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_proc_card.dat @@ -0,0 +1,12 @@ +# Model +import model TopFCNC_UFO-kappa_act_ProdConventions_no_c_mass --modelname + +# Multiparticles +define l+ = e+ mu+ ta+ +define l- = e- mu- ta- + +generate p p > t a NP=1, (t > w+ b, w+ > l+ vl ) +add process p p > t~ a NP=1, (t~ > w- b~, w- > l- vl~) + +# Output folder +output TA_Tleptonic_kappa_act_LO -f -nojpeg diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_run_card.dat new file mode 100644 index 000000000000..2c3dbde2bb74 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/TA_Tleptonic_kappa_act_LO/TA_Tleptonic_kappa_act_LO_run_card.dat @@ -0,0 +1,265 @@ +#********************************************************************* +# MadGraph5_aMC@NLO * +# * +# run_card.dat MadEvent * +# * +# This file is used to set the parameters of the run. * +# * +# Some notation/conventions: * +# * +# Lines starting with a '# ' are info or comments * +# * +# mind the format: value = variable ! comment * +#********************************************************************* +# +#******************* +# Running parameters +#******************* +# +#********************************************************************* +# Tag name for the run (one word) * +#********************************************************************* +tag_1 = run_tag ! name of the run +#********************************************************************* +# Run to generate the grid pack * +#********************************************************************* +True = gridpack !True = setting up the grid pack +#********************************************************************* +# Number of events and rnd seed * +# Warning: Do not generate more than 1M events in a single run * +# If you want to run Pythia, avoid more than 50k events in a run. * +#********************************************************************* +100000 = nevents ! Number of unweighted events requested +0 = iseed ! rnd seed (0=assigned automatically=default)) +#********************************************************************* +# Collider type and energy * +# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, * +# 3=photon from electron * +#********************************************************************* +1 = lpp1 ! beam 1 type +1 = lpp2 ! beam 2 type +6500.0 = ebeam1 ! beam 1 total energy in GeV +6500.0 = ebeam2 ! beam 2 total energy in GeV +#********************************************************************* +# Beam polarization from -100 (left-handed) to 100 (right-handed) * +#********************************************************************* +0.0 = polbeam1 ! beam polarization for beam 1 +0.0 = polbeam2 ! beam polarization for beam 2 +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* +lhapdf = pdlabel ! PDF set +$DEFAULT_PDF_SETS = lhaid +$DEFAULT_PDF_MEMBERS = reweight_PDF ! if pdlabel=lhapdf, this is the lhapdf number +#********************************************************************* +# Renormalization and factorization scales * +#********************************************************************* +False = fixed_ren_scale ! if .true. use fixed ren scale +False = fixed_fac_scale ! if .true. use fixed fac scale +91.188 = scale ! fixed ren scale +91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1 +91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2 +-1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices +1.0 = scalefact ! scale factor for event-by-event scales +#********************************************************************* +# Time of flight information. (-1 means not run) +#********************************************************************* +-1.0 = time_of_flight ! threshold below which info is not written +#********************************************************************* +# Matching - Warning! ickkw > 1 is still beta +#********************************************************************* +0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching +1 = highestmult ! for ickkw=2, highest mult group +1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE +1.0 = alpsfact ! scale factor for QCD emission vx +False = chcluster ! cluster only according to channel diag +4 = asrwgtflavor ! highest quark flavor for a_s reweight +True = clusinfo ! include clustering tag in output +3.0 = lhe_version ! Change the way clustering information pass to shower. +#********************************************************************* +#********************************************************** +# +#********************************************************** +# Automatic ptj and mjj cuts if xqcut > 0 +# (turn off for VBF and single top processes) +#********************************************************** +False = auto_ptj_mjj ! Automatic setting of ptj and mjj +#********************************************************** +# +#********************************** +# BW cutoff (M+/-bwcutoff*Gamma) +#********************************** +15.0 = bwcutoff ! (M+/-bwcutoff*Gamma) +#********************************************************** +# Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not +# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) +#************************************************************* +False = cut_decays ! Cut decay products +#************************************************************* +# Number of helicities to sum per event (0 = all helicities) +# 0 gives more stable result, but longer run time (needed for +# long decay chains e.g.). +# Use >=2 if most helicities contribute, e.g. pure QCD. +#************************************************************* +0 = nhel ! Number of helicities used per event +#******************* +# Standard Cuts +#******************* +# +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* +0.0 = ptj ! minimum pt for the jets +0.0 = ptb ! minimum pt for the b +0.0 = pta ! minimum pt for the photons +0.0 = ptl ! minimum pt for the charged leptons +0.0 = misset ! minimum missing Et (sum of neutrino's momenta) +0.0 = ptheavy ! minimum pt for one heavy final state +-1.0 = ptjmax ! maximum pt for the jets +-1.0 = ptbmax ! maximum pt for the b +-1.0 = ptamax ! maximum pt for the photons +-1.0 = ptlmax ! maximum pt for the charged leptons +-1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* +0.0 = ej ! minimum E for the jets +0.0 = eb ! minimum E for the b +0.0 = ea ! minimum E for the photons +0.0 = el ! minimum E for the charged leptons +-1.0 = ejmax ! maximum E for the jets +-1.0 = ebmax ! maximum E for the b +-1.0 = eamax ! maximum E for the photons +-1.0 = elmax ! maximum E for the charged leptons +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* +5.0 = etaj ! max rap for the jets +-1.0 = etab ! max rap for the b +5.0 = etaa ! max rap for the photons +5.0 = etal ! max rap for the charged leptons +0.0 = etajmin ! min rap for the jets +0.0 = etabmin ! min rap for the b +0.0 = etaamin ! min rap for the photons +0.0 = etalmin ! main rap for the charged leptons +#********************************************************************* +# Minimum and maximum DeltaR distance * +#********************************************************************* +0.001 = drjj ! min distance between jets +0.0 = drbb ! min distance between b's +0.001 = drll ! min distance between leptons +0.001 = draa ! min distance between gammas +0.0 = drbj ! min distance between b and jet +0.001 = draj ! min distance between gamma and jet +0.001 = drjl ! min distance between jet and lepton +0.0 = drab ! min distance between gamma and b +0.0 = drbl ! min distance between b and lepton +0.001 = dral ! min distance between gamma and lepton +-1.0 = drjjmax ! max distance between jets +-1.0 = drbbmax ! max distance between b's +-1.0 = drllmax ! max distance between leptons +-1.0 = draamax ! max distance between gammas +-1.0 = drbjmax ! max distance between b and jet +-1.0 = drajmax ! max distance between gamma and jet +-1.0 = drjlmax ! max distance between jet and lepton +-1.0 = drabmax ! max distance between gamma and b +-1.0 = drblmax ! max distance between b and lepton +-1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +# WARNING: for four lepton final state mmll cut require to have * +# different lepton masses for each flavor! * +#********************************************************************* +0.0 = mmjj ! min invariant mass of a jet pair +0.0 = mmbb ! min invariant mass of a b pair +0.0 = mmaa ! min invariant mass of gamma gamma pair +0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair +-1.0 = mmjjmax ! max invariant mass of a jet pair +-1.0 = mmbbmax ! max invariant mass of a b pair +-1.0 = mmaamax ! max invariant mass of gamma gamma pair +-1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair +#********************************************************************* +# Minimum and maximum invariant mass for all letpons * +#********************************************************************* +0.0 = mmnl ! min invariant mass for all letpons (l+- and vl) +-1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl) +#********************************************************************* +# Minimum and maximum pt for 4-momenta sum of leptons * +#********************************************************************* +0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl) +-1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl) +#********************************************************************* +# Inclusive cuts * +#********************************************************************* +0.0 = xptj ! minimum pt for at least one jet +0.0 = xptb ! minimum pt for at least one b +0.0 = xpta ! minimum pt for at least one photon +0.0 = xptl ! minimum pt for at least one charged lepton +#********************************************************************* +# Control the pt's of the jets sorted by pt * +#********************************************************************* +0.0 = ptj1min ! minimum pt for the leading jet in pt +0.0 = ptj2min ! minimum pt for the second jet in pt +0.0 = ptj3min ! minimum pt for the third jet in pt +0.0 = ptj4min ! minimum pt for the fourth jet in pt +-1.0 = ptj1max ! maximum pt for the leading jet in pt +-1.0 = ptj2max ! maximum pt for the second jet in pt +-1.0 = ptj3max ! maximum pt for the third jet in pt +-1.0 = ptj4max ! maximum pt for the fourth jet in pt +0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts +#********************************************************************* +# Control the pt's of leptons sorted by pt * +#********************************************************************* +0.0 = ptl1min ! minimum pt for the leading lepton in pt +0.0 = ptl2min ! minimum pt for the second lepton in pt +0.0 = ptl3min ! minimum pt for the third lepton in pt +0.0 = ptl4min ! minimum pt for the fourth lepton in pt +-1.0 = ptl1max ! maximum pt for the leading lepton in pt +-1.0 = ptl2max ! maximum pt for the second lepton in pt +-1.0 = ptl3max ! maximum pt for the third lepton in pt +-1.0 = ptl4max ! maximum pt for the fourth lepton in pt +#********************************************************************* +# Control the Ht(k)=Sum of k leading jets * +#********************************************************************* +0.0 = htjmin ! minimum jet HT=Sum(jet pt) +-1.0 = htjmax ! maximum jet HT=Sum(jet pt) +0.0 = ihtmin !inclusive Ht for all partons (including b) +-1.0 = ihtmax !inclusive Ht for all partons (including b) +0.0 = ht2min ! minimum Ht for the two leading jets +0.0 = ht3min ! minimum Ht for the three leading jets +0.0 = ht4min ! minimum Ht for the four leading jets +-1.0 = ht2max ! maximum Ht for the two leading jets +-1.0 = ht3max ! maximum Ht for the three leading jets +-1.0 = ht4max ! maximum Ht for the four leading jets +#*********************************************************************** +# Photon-isolation cuts, according to hep-ph/9801442 * +# When ptgmin=0, all the other parameters are ignored * +# When ptgmin>0, pta and draj are not going to be used * +#*********************************************************************** +0.0 = ptgmin ! Min photon transverse momentum +0.4 = R0gamma ! Radius of isolation code +1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442 +1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442 +True = isoEM ! isolate photons from EM energy (photons and leptons) +#********************************************************************* +# WBF cuts * +#********************************************************************* +0.0 = xetamin ! minimum rapidity for two jets in the WBF case +0.0 = deltaeta ! minimum rapidity for two jets in the WBF case +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* +-1.0 = ktdurham +0.4 = dparameter +#********************************************************************* +# maximal pdg code for quark to be considered as a light jet * +# (otherwise b cuts are applied) * +#********************************************************************* +4 = maxjetflavor ! Maximum jet pdg code +#********************************************************************* +# Jet measure cuts * +#********************************************************************* +0.0 = xqcut ! minimum kt jet measure between partons +#********************************************************************* +# +#********************************************************************* diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/create_variations_TLL_Tleptonic.py b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/create_variations_TLL_Tleptonic.py new file mode 100644 index 000000000000..b2b141fc271e --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/create_variations_TLL_Tleptonic.py @@ -0,0 +1,27 @@ +import os + +print "This script will create the cards for ST processes produced by flavour-changing neutral currents." + + +sampledict = {"TLL_Tleptonic_kappa_zct_LO": {"MODELNAME": "TopFCNC_UFO-kappa_zct_ProdConventions_no_c_mass",}, + "TLL_Tleptonic_kappa_zut_LO": {"MODELNAME": "TopFCNC_UFO-kappa_zut_ProdConventions_no_c_mass",} + } + + + + + +for el in sampledict: + os.system("mkdir -p ./" + el) + + for ejC in ["customizecards.dat", "extramodels.dat", "run_card.dat"]: + os.system("cp ./orig_cards/" + ejC + " ./" + el + "/" + el + "_" + ejC) + + # proc card + outproccard = "" + with open("./orig_cards/proc_card.dat", "r") as theF: + for line in theF.readlines(): + outproccard += line.replace("PROCNAME", el).replace("MODELNAME", sampledict[el]["MODELNAME"]) + outprocF = open("./" + el + "/" + el + "_proc_card.dat", "w") + outprocF.write(outproccard) + outprocF.close(); del outprocF diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/customizecards.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/customizecards.dat new file mode 100644 index 000000000000..f0ee62fa058f --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/customizecards.dat @@ -0,0 +1,4 @@ +#put card customizations here +set param_card mass 6 172.5 +set param_card yukawa 6 172.5 +set param_card DECAY 6 auto diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/extramodels.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/extramodels.dat new file mode 100644 index 000000000000..a2a62b1c4be8 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/extramodels.dat @@ -0,0 +1,2 @@ +#Anomalous Top couplings LO model version 1.4 from https://twiki.cern.ch/twiki/bin/view/CMS/TopFCNCmodel +TopFCNC_patch4_UFO.tgz diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/proc_card.dat new file mode 100644 index 000000000000..62105a6b9dd7 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/proc_card.dat @@ -0,0 +1,11 @@ +# Model +import model MODELNAME --modelname +# Multiparticles +define l+ = e+ mu+ ta+ +define l- = e- mu- ta- + +generate p p > t z, z > l+ l-, (t > w+ b, w+ > l+ vl ) NP=1 @10 +add process p p > t~ z, z > l+ l-, (t~ > w- b~, w- > l- vl~) NP=1 @20 + +# Output folder +output PROCNAME -f -nojpeg diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/run_card.dat new file mode 100644 index 000000000000..b1dda0279353 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/ST_FCNC/orig_cards/run_card.dat @@ -0,0 +1,265 @@ +#********************************************************************* +# MadGraph5_aMC@NLO * +# * +# run_card.dat MadEvent * +# * +# This file is used to set the parameters of the run. * +# * +# Some notation/conventions: * +# * +# Lines starting with a '# ' are info or comments * +# * +# mind the format: value = variable ! comment * +#********************************************************************* +# +#******************* +# Running parameters +#******************* +# +#********************************************************************* +# Tag name for the run (one word) * +#********************************************************************* +tag_1 = run_tag ! name of the run +#********************************************************************* +# Run to generate the grid pack * +#********************************************************************* +True = gridpack !True = setting up the grid pack +#********************************************************************* +# Number of events and rnd seed * +# Warning: Do not generate more than 1M events in a single run * +# If you want to run Pythia, avoid more than 50k events in a run. * +#********************************************************************* +100000 = nevents ! Number of unweighted events requested +0 = iseed ! rnd seed (0=assigned automatically=default)) +#********************************************************************* +# Collider type and energy * +# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, * +# 3=photon from electron * +#********************************************************************* +1 = lpp1 ! beam 1 type +1 = lpp2 ! beam 2 type +6500.0 = ebeam1 ! beam 1 total energy in GeV +6500.0 = ebeam2 ! beam 2 total energy in GeV +#********************************************************************* +# Beam polarization from -100 (left-handed) to 100 (right-handed) * +#********************************************************************* +0.0 = polbeam1 ! beam polarization for beam 1 +0.0 = polbeam2 ! beam polarization for beam 2 +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* +lhapdf = pdlabel ! PDF set +$DEFAULT_PDF_SETS = lhaid +$DEFAULT_PDF_MEMBERS = reweight_PDF ! if pdlabel=lhapdf, this is the lhapdf number +#********************************************************************* +# Renormalization and factorization scales * +#********************************************************************* +False = fixed_ren_scale ! if .true. use fixed ren scale +False = fixed_fac_scale ! if .true. use fixed fac scale +91.188 = scale ! fixed ren scale +91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1 +91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2 +-1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices +1.0 = scalefact ! scale factor for event-by-event scales +#********************************************************************* +# Time of flight information. (-1 means not run) +#********************************************************************* +-1.0 = time_of_flight ! threshold below which info is not written +#********************************************************************* +# Matching - Warning! ickkw > 1 is still beta +#********************************************************************* +0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching +1 = highestmult ! for ickkw=2, highest mult group +1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE +1.0 = alpsfact ! scale factor for QCD emission vx +False = chcluster ! cluster only according to channel diag +4 = asrwgtflavor ! highest quark flavor for a_s reweight +True = clusinfo ! include clustering tag in output +3.0 = lhe_version ! Change the way clustering information pass to shower. +#********************************************************************* +#********************************************************** +# +#********************************************************** +# Automatic ptj and mjj cuts if xqcut > 0 +# (turn off for VBF and single top processes) +#********************************************************** +False = auto_ptj_mjj ! Automatic setting of ptj and mjj +#********************************************************** +# +#********************************** +# BW cutoff (M+/-bwcutoff*Gamma) +#********************************** +15.0 = bwcutoff ! (M+/-bwcutoff*Gamma) +#********************************************************** +# Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not +# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) +#************************************************************* +False = cut_decays ! Cut decay products +#************************************************************* +# Number of helicities to sum per event (0 = all helicities) +# 0 gives more stable result, but longer run time (needed for +# long decay chains e.g.). +# Use >=2 if most helicities contribute, e.g. pure QCD. +#************************************************************* +0 = nhel ! Number of helicities used per event +#******************* +# Standard Cuts +#******************* +# +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* +0.0 = ptj ! minimum pt for the jets +0.0 = ptb ! minimum pt for the b +0.0 = pta ! minimum pt for the photons +0.0 = ptl ! minimum pt for the charged leptons +0.0 = misset ! minimum missing Et (sum of neutrino's momenta) +0.0 = ptheavy ! minimum pt for one heavy final state +-1.0 = ptjmax ! maximum pt for the jets +-1.0 = ptbmax ! maximum pt for the b +-1.0 = ptamax ! maximum pt for the photons +-1.0 = ptlmax ! maximum pt for the charged leptons +-1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* +0.0 = ej ! minimum E for the jets +0.0 = eb ! minimum E for the b +0.0 = ea ! minimum E for the photons +0.0 = el ! minimum E for the charged leptons +-1.0 = ejmax ! maximum E for the jets +-1.0 = ebmax ! maximum E for the b +-1.0 = eamax ! maximum E for the photons +-1.0 = elmax ! maximum E for the charged leptons +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* +5.0 = etaj ! max rap for the jets +-1.0 = etab ! max rap for the b +5.0 = etaa ! max rap for the photons +5.0 = etal ! max rap for the charged leptons +0.0 = etajmin ! min rap for the jets +0.0 = etabmin ! min rap for the b +0.0 = etaamin ! min rap for the photons +0.0 = etalmin ! main rap for the charged leptons +#********************************************************************* +# Minimum and maximum DeltaR distance * +#********************************************************************* +0.001 = drjj ! min distance between jets +0.0 = drbb ! min distance between b's +0.001 = drll ! min distance between leptons +0.001 = draa ! min distance between gammas +0.0 = drbj ! min distance between b and jet +0.001 = draj ! min distance between gamma and jet +0.001 = drjl ! min distance between jet and lepton +0.0 = drab ! min distance between gamma and b +0.0 = drbl ! min distance between b and lepton +0.001 = dral ! min distance between gamma and lepton +-1.0 = drjjmax ! max distance between jets +-1.0 = drbbmax ! max distance between b's +-1.0 = drllmax ! max distance between leptons +-1.0 = draamax ! max distance between gammas +-1.0 = drbjmax ! max distance between b and jet +-1.0 = drajmax ! max distance between gamma and jet +-1.0 = drjlmax ! max distance between jet and lepton +-1.0 = drabmax ! max distance between gamma and b +-1.0 = drblmax ! max distance between b and lepton +-1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +# WARNING: for four lepton final state mmll cut require to have * +# different lepton masses for each flavor! * +#********************************************************************* +0.0 = mmjj ! min invariant mass of a jet pair +0.0 = mmbb ! min invariant mass of a b pair +0.0 = mmaa ! min invariant mass of gamma gamma pair +10.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair +-1.0 = mmjjmax ! max invariant mass of a jet pair +-1.0 = mmbbmax ! max invariant mass of a b pair +-1.0 = mmaamax ! max invariant mass of gamma gamma pair +-1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair +#********************************************************************* +# Minimum and maximum invariant mass for all letpons * +#********************************************************************* +0.0 = mmnl ! min invariant mass for all letpons (l+- and vl) +-1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl) +#********************************************************************* +# Minimum and maximum pt for 4-momenta sum of leptons * +#********************************************************************* +0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl) +-1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl) +#********************************************************************* +# Inclusive cuts * +#********************************************************************* +0.0 = xptj ! minimum pt for at least one jet +0.0 = xptb ! minimum pt for at least one b +0.0 = xpta ! minimum pt for at least one photon +0.0 = xptl ! minimum pt for at least one charged lepton +#********************************************************************* +# Control the pt's of the jets sorted by pt * +#********************************************************************* +0.0 = ptj1min ! minimum pt for the leading jet in pt +0.0 = ptj2min ! minimum pt for the second jet in pt +0.0 = ptj3min ! minimum pt for the third jet in pt +0.0 = ptj4min ! minimum pt for the fourth jet in pt +-1.0 = ptj1max ! maximum pt for the leading jet in pt +-1.0 = ptj2max ! maximum pt for the second jet in pt +-1.0 = ptj3max ! maximum pt for the third jet in pt +-1.0 = ptj4max ! maximum pt for the fourth jet in pt +0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts +#********************************************************************* +# Control the pt's of leptons sorted by pt * +#********************************************************************* +0.0 = ptl1min ! minimum pt for the leading lepton in pt +0.0 = ptl2min ! minimum pt for the second lepton in pt +0.0 = ptl3min ! minimum pt for the third lepton in pt +0.0 = ptl4min ! minimum pt for the fourth lepton in pt +-1.0 = ptl1max ! maximum pt for the leading lepton in pt +-1.0 = ptl2max ! maximum pt for the second lepton in pt +-1.0 = ptl3max ! maximum pt for the third lepton in pt +-1.0 = ptl4max ! maximum pt for the fourth lepton in pt +#********************************************************************* +# Control the Ht(k)=Sum of k leading jets * +#********************************************************************* +0.0 = htjmin ! minimum jet HT=Sum(jet pt) +-1.0 = htjmax ! maximum jet HT=Sum(jet pt) +0.0 = ihtmin !inclusive Ht for all partons (including b) +-1.0 = ihtmax !inclusive Ht for all partons (including b) +0.0 = ht2min ! minimum Ht for the two leading jets +0.0 = ht3min ! minimum Ht for the three leading jets +0.0 = ht4min ! minimum Ht for the four leading jets +-1.0 = ht2max ! maximum Ht for the two leading jets +-1.0 = ht3max ! maximum Ht for the three leading jets +-1.0 = ht4max ! maximum Ht for the four leading jets +#*********************************************************************** +# Photon-isolation cuts, according to hep-ph/9801442 * +# When ptgmin=0, all the other parameters are ignored * +# When ptgmin>0, pta and draj are not going to be used * +#*********************************************************************** +0.0 = ptgmin ! Min photon transverse momentum +0.4 = R0gamma ! Radius of isolation code +1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442 +1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442 +True = isoEM ! isolate photons from EM energy (photons and leptons) +#********************************************************************* +# WBF cuts * +#********************************************************************* +0.0 = xetamin ! minimum rapidity for two jets in the WBF case +0.0 = deltaeta ! minimum rapidity for two jets in the WBF case +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* +-1.0 = ktdurham +0.4 = dparameter +#********************************************************************* +# maximal pdg code for quark to be considered as a light jet * +# (otherwise b cuts are applied) * +#********************************************************************* +4 = maxjetflavor ! Maximum jet pdg code +#********************************************************************* +# Jet measure cuts * +#********************************************************************* +0.0 = xqcut ! minimum kt jet measure between partons +#********************************************************************* +# +#********************************************************************* diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/createFCNCcards.py b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/createFCNCcards.py new file mode 100644 index 000000000000..914ac16f581c --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/createFCNCcards.py @@ -0,0 +1,282 @@ +import os + +print "This script will create the cards for all ST & TT processes produced by flavour-changing neutral currents." + +models = {"zeta_zut" : {"1 0.000000e+00 # zzq1" : "1 9.916830e-01 # zzq1"}, + "zeta_zct" : {"2 0.000000e+00 # zzq2" : "2 9.916831e-01 # zzq2"}, + "kappa_zut" : {"1 0.000000e+00 # kzq1" : "1 3.179691e-03 # kzq1"}, + "kappa_zct" : {"2 0.000000e+00 # kzq2" : "2 3.179691e-03 # kzq2"}, + "kappa_aut" : {"1 0.000000e+00 # kaq1" : "1 2.889977e-03 # kaq1"}, + "kappa_act" : {"2 0.000000e+00 # kaq2" : "2 2.889977e-03 # kaq2"}, + "eta_hut" : {"1 0.000000e+00 # ehq1" : "1 9.306701e-01 # ehq1"}, + "eta_hct" : {"2 0.000000e+00 # ehq2" : "2 9.306702e-01 # ehq2"}, +} + + +sampledict = { +# A +"TT_FCNC-aT2AJ_Tleptonic_kappa_act_LO": {"TOPDECAY" : "t > w+ b, w+ > l+ vl NP=0", + "ANTITOPDECAY": "t~ > a c~ NP=1", + "param_card" : models["kappa_act"], + "run_card" : {"0.0 = ptj": "20.0 = ptj", + "0.0 = pta": "20.0 = pta", + "0.0 = ptl": "20.0 = ptl", + "5.0 = etaa": "3.0 = etaa", + "5.0 = etal": "3.0 = etal", + }, + "run_card_extra" : """#********************************************************************* +# Store info for systematics studies * +# WARNING: If use_syst is T, matched Pythia output is * +# meaningful ONLY if plotted taking matchscale * +# reweighting into account! * +#********************************************************************* +True = use_syst ! Enable systematics studies +""", + }, +"TT_FCNC-aT2AJ_Tleptonic_kappa_aut_LO": {"TOPDECAY" : "t > w+ b, w+ > l+ vl NP=0", + "ANTITOPDECAY": "t~ > a u~ NP=1", + "param_card" : models["kappa_aut"], + "run_card" : {"0.0 = ptj": "20.0 = ptj", + "0.0 = pta": "20.0 = pta", + "0.0 = ptl": "20.0 = ptl", + "5.0 = etaa": "3.0 = etaa", + "5.0 = etal": "3.0 = etal", + }, + "run_card_extra" : """#********************************************************************* +# Store info for systematics studies * +# WARNING: If use_syst is T, matched Pythia output is * +# meaningful ONLY if plotted taking matchscale * +# reweighting into account! * +#********************************************************************* +True = use_syst ! Enable systematics studies +""", + }, +"TT_FCNC-T2AJ_aTleptonic_kappa_act_LO": {"TOPDECAY" : "t > a c NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > l- vl~ NP=0", + "param_card" : models["kappa_act"], + "run_card" : {"0.0 = ptj": "20.0 = ptj", + "0.0 = pta": "20.0 = pta", + "0.0 = ptl": "20.0 = ptl", + "5.0 = etaa": "3.0 = etaa", + "5.0 = etal": "3.0 = etal", + }, + "run_card_extra" : """#********************************************************************* +# Store info for systematics studies * +# WARNING: If use_syst is T, matched Pythia output is * +# meaningful ONLY if plotted taking matchscale * +# reweighting into account! * +#********************************************************************* +True = use_syst ! Enable systematics studies +""", + }, +"TT_FCNC-T2AJ_aTleptonic_kappa_aut_LO": {"TOPDECAY" : "t > a u NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > l- vl~ NP=0", + "madspin_card_extra":{"define vl~ = ve~ vm~ vt~":"define vl~ = ve~ vm~ vt~\ndefine inclu = l- l+ j vl vl~"}, + "param_card" : models["kappa_aut"], + "run_card" : {"0.0 = ptj": "20.0 = ptj", + "0.0 = pta": "20.0 = pta", + "0.0 = ptl": "20.0 = ptl", + "5.0 = etaa": "3.0 = etaa", + "5.0 = etal": "3.0 = etal", + }, + "run_card_extra" : """#********************************************************************* +# Store info for systematics studies * +# WARNING: If use_syst is T, matched Pythia output is * +# meaningful ONLY if plotted taking matchscale * +# reweighting into account! * +#********************************************************************* +True = use_syst ! Enable systematics studies +""", + }, + +# Z +"TT_FCNC-aT2ZJ_Tleptonic_kappa_zut_LO": {"TOPDECAY" : "t > w+ b, w+ > l+ vl NP=0", + "ANTITOPDECAY": "t~ > z u~ NP=1", + "do3jets" : False, + "param_card" : models["kappa_zut"], + }, +"TT_FCNC-T2ZJ_aTleptonic_kappa_zut_LO": {"TOPDECAY" : "t > z u NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > l- vl~ NP=0", + "do3jets" : False, + "param_card" : models["kappa_zut"], + }, +"TT_FCNC-aT2ZJ_Tleptonic_kappa_zct_LO": {"TOPDECAY" : "t > w+ b, w+ > l+ vl NP=0", + "ANTITOPDECAY": "t~ > z c~ NP=1", + "do3jets" : False, + "param_card" : models["kappa_zct"], + }, +"TT_FCNC-T2ZJ_aTleptonic_kappa_zct_LO": {"TOPDECAY" : "t > z c NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > l- vl~ NP=0", + "do3jets" : False, + "param_card" : models["kappa_zct"], + }, +"TT_FCNC-aT2ZJ_Thadronic_kappa_zut_LO": {"TOPDECAY" : "t > w+ b, w+ > j j NP=0", + "ANTITOPDECAY": "t~ > z u~ NP=1", + "do3jets" : False, + "param_card" : models["kappa_zut"], + }, +"TT_FCNC-T2ZJ_aThadronic_kappa_zut_LO": {"TOPDECAY" : "t > z u NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > j j NP=0", + "do3jets" : False, + "param_card" : models["kappa_zut"], + }, +"TT_FCNC-aT2ZJ_Thadronic_kappa_zct_LO": {"TOPDECAY" : "t > w+ b, w+ > j j NP=0", + "ANTITOPDECAY": "t~ > z c~ NP=1", + "do3jets" : False, + "param_card" : models["kappa_zct"], + }, +"TT_FCNC-T2ZJ_aThadronic_kappa_zct_LO": {"TOPDECAY" : "t > z c NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > j j NP=0", + "do3jets" : False, + "param_card" : models["kappa_zct"], + }, +"TT_FCNC-aT2ZJ_Tleptonic_zeta_zut_LO" : {"TOPDECAY" : "t > w+ b, w+ > l+ vl NP=0", + "ANTITOPDECAY": "t~ > z u~ NP=1", + "do3jets" : False, + "param_card" : models["zeta_zut"], + }, +"TT_FCNC-T2ZJ_aTleptonic_zeta_zut_LO" : {"TOPDECAY" : "t > z u NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > l- vl~ NP=0", + "do3jets" : False, + "param_card" : models["zeta_zut"], + }, +"TT_FCNC-aT2ZJ_Tleptonic_zeta_zct_LO" : {"TOPDECAY" : "t > w+ b, w+ > l+ vl NP=0", + "ANTITOPDECAY": "t~ > z c~ NP=1", + "do3jets" : False, + "param_card" : models["zeta_zct"], + }, +"TT_FCNC-T2ZJ_aTleptonic_zeta_zct_LO" : {"TOPDECAY" : "t > z c NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > l- vl~ NP=0", + "do3jets" : False, + "param_card" : models["zeta_zct"], + }, +"TT_FCNC-aT2ZJ_Thadronic_zeta_zut_LO" : {"TOPDECAY" : "t > w+ b, w+ > j j NP=0", + "ANTITOPDECAY": "t~ > z u~ NP=1", + "do3jets" : False, + "param_card" : models["zeta_zut"], + }, +"TT_FCNC-T2ZJ_aThadronic_zeta_zut_LO" : {"TOPDECAY" : "t > z u NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > j j NP=0", + "madspin_card_extra":{"define vl~ = ve~ vm~ vt~":"define vl~ = ve~ vm~ vt~\ndefine inclu = l- l+ j vl vl~"}, + "do3jets" : False, + "param_card" : models["zeta_zut"], + }, +"TT_FCNC-aT2ZJ_Thadronic_zeta_zct_LO" : {"TOPDECAY" : "t > w+ b, w+ > j j NP=0", + "ANTITOPDECAY": "t~ > z c~ NP=1", + "do3jets" : False, + "param_card" : models["zeta_zct"], + }, +"TT_FCNC-T2ZJ_aThadronic_zeta_zct_LO" : {"TOPDECAY" : "t > z c NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > j j NP=0", + "do3jets" : False, + "param_card" : models["zeta_zct"], + }, + +## H +"TT_FCNC-aT2HJ_Thadronic_eta_hct_LO" : {"TOPDECAY" : "t > w+ b, w+ > j j NP=0", + "ANTITOPDECAY": "t~ > h c~ NP=1", + "do3jets" : False, + "param_card" : models["eta_hct"], + }, +"TT_FCNC-aT2HJ_Thadronic_eta_hut_LO" : {"TOPDECAY" : "t > w+ b, w+ > j j NP=0", + "ANTITOPDECAY": "t~ > h u~ NP=1", + "do3jets" : False, + "param_card" : models["eta_hut"], + }, +"TT_FCNC-aT2HJ_Tleptonic_eta_hct_LO" : {"TOPDECAY" : "t > w+ b, w+ > l+ vl NP=0", + "ANTITOPDECAY": "t~ > h c~ NP=1", + "do3jets" : False, + "param_card" : models["eta_hct"], + }, +"TT_FCNC-aT2HJ_Tleptonic_eta_hut_LO" : {"TOPDECAY" : "t > w+ b, w+ > l+ vl NP=0", + "ANTITOPDECAY": "t~ > h u~ NP=1", + "do3jets" : False, + "param_card" : models["eta_hut"], + }, +"TT_FCNC-T2HJ_aThadronic_eta_hct_LO" : {"TOPDECAY" : "t > h c NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > j j NP=0", + "do3jets" : False, + "param_card" : models["eta_hct"], + }, +"TT_FCNC-T2HJ_aThadronic_eta_hut_LO" : {"TOPDECAY" : "t > h u NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > j j NP=0", + "do3jets" : False, + "param_card" : models["eta_hut"], + }, +"TT_FCNC-T2HJ_aTleptonic_eta_hct_LO" : {"TOPDECAY" : "t > h c NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > l- vl~ NP=0", + "do3jets" : False, + "param_card" : models["eta_hct"], + }, +"TT_FCNC-T2HJ_aTleptonic_eta_hut_LO" : {"TOPDECAY" : "t > h u NP=1", + "ANTITOPDECAY": "t~ > w- b~, w- > l- vl~ NP=0", + "do3jets" : False, + "param_card" : models["eta_hut"], + }, +} + + +for el in sampledict: + os.system("mkdir -p ./" + el) + + for ejC in ["extramodels.dat"]: + os.system("cp ./orig_cards/" + ejC + " ./" + el + "/" + el + "_" + ejC) + + # proc card + outproccard = "" + with open("./orig_cards/proc_card.dat", "r") as theF: + for line in theF.readlines(): + outproccard += line.replace("PROCNAME", el) + if "NP=0 @12" in line and "do3jets" in sampledict[el]: + if sampledict[el]["do3jets"]: + outproccard += "add process p p > t t~ j j j NP=0 @13\n" + + outprocF = open("./" + el + "/" + el + "_proc_card.dat", "w") + outprocF.write(outproccard) + outprocF.close(); del outprocF + + + # run card + outruncard = "" + with open("./orig_cards/run_card.dat", "r") as theF: + for line in theF.readlines(): + tmpline = line + if "run_card" in sampledict[el]: + for subel,subrep in sampledict[el]["run_card"].iteritems(): + tmpline = tmpline.replace(subel, subrep) + outruncard += tmpline + if "run_card_extra" in sampledict[el]: + outruncard += sampledict[el]["run_card_extra"] + + outrunF = open("./" + el + "/" + el + "_run_card.dat", "w") + outrunF.write(outruncard) + outrunF.close(); del outrunF + + + # param card + outparamcard = "" + with open("./orig_cards/param_card.dat", "r") as theF: + for line in theF.readlines(): + tmpline = line + for subel,subrep in sampledict[el]["param_card"].iteritems(): + tmpline = tmpline.replace(subel, subrep) + outparamcard += tmpline + outparamF = open("./" + el + "/" + el + "_param_card.dat", "w") + outparamF.write(outparamcard) + outparamF.close(); del outparamF + + + # madspin card + outspincard = "" + with open("./orig_cards/madspin_card.dat", "r") as theF: + for line in theF.readlines(): + tmpline = line + if "madspin_card_extra" in sampledict[el]: + for subel,subrep in sampledict[el]["madspin_card_extra"].iteritems(): + tmpline = tmpline.replace(subel, subrep) + + outspincard += tmpline.replace("TOPDECAY", sampledict[el]["TOPDECAY"]).replace("ANTItOPDECAY", sampledict[el]["ANTITOPDECAY"]) + + outspinF = open("./" + el + "/" + el + "_madspin_card.dat", "w") + outspinF.write(outspincard) + outspinF.close(); del outspinF diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/extramodels.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/extramodels.dat new file mode 100644 index 000000000000..a2a62b1c4be8 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/extramodels.dat @@ -0,0 +1,2 @@ +#Anomalous Top couplings LO model version 1.4 from https://twiki.cern.ch/twiki/bin/view/CMS/TopFCNCmodel +TopFCNC_patch4_UFO.tgz diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/madspin_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/madspin_card.dat new file mode 100644 index 000000000000..18f571a751eb --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/madspin_card.dat @@ -0,0 +1,19 @@ +# Settings +set ms_dir ./madspin +set Nevents_for_max_weigth 250 # number of events for the estimate of the max. weight +set max_weight_ps_point 400 # number of PS to estimate the maximum for each event +set max_running_process 1 + +# Multiparticles +define j = g u u~ d d~ s s~ c c~ b b~ +define l+ = e+ mu+ ta+ +define l- = e- mu- ta- +define vl = ve vm vt +define vl~ = ve~ vm~ vt~ + +# SM decay of the quark top +# Anomalous coupling is too small: MadSpin cannot find the anomalous decay and forces SM decay +decay ANTItOPDECAY +decay TOPDECAY + +launch diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/param_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/param_card.dat new file mode 100644 index 000000000000..903d567953bf --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/param_card.dat @@ -0,0 +1,111 @@ +###################################################################### +## PARAM_CARD AUTOMATICALY GENERATED BY MG5 #### +###################################################################### +################################### +## INFORMATION FOR ACEHQ +################################### +BLOCK ACEHQ # + 1 0.000000e+00 # ehq1 + 2 0.000000e+00 # ehq2 + 3 0.000000e+00 # ehq3 +################################### +## INFORMATION FOR ACKAQ +################################### +BLOCK ACKAQ # + 1 0.000000e+00 # kaq1 + 2 0.000000e+00 # kaq2 + 3 0.000000e+00 # kaq3 +################################### +## INFORMATION FOR ACKGQ +################################### +BLOCK ACKGQ # + 1 0.000000e+00 # kgq1 + 2 0.000000e+00 # kgq2 + 3 0.000000e+00 # kgq3 +################################### +## INFORMATION FOR ACKWQ +################################### +BLOCK ACKWQ # + 1 0.000000e+00 # kwq1 + 2 0.000000e+00 # kwq2 + 3 0.000000e+00 # kwq3 +################################### +## INFORMATION FOR ACKZQ +################################### +BLOCK ACKZQ # + 1 0.000000e+00 # kzq1 + 2 0.000000e+00 # kzq2 + 3 0.000000e+00 # kzq3 +################################### +## INFORMATION FOR ACZWQ +################################### +BLOCK ACZWQ # + 1 0.000000e+00 # zwq1 + 2 0.000000e+00 # zwq2 + 3 0.000000e+00 # zwq3 +################################### +## INFORMATION FOR ACZZQ +################################### +BLOCK ACZZQ # + 1 0.000000e+00 # zzq1 + 2 0.000000e+00 # zzq2 + 3 0.000000e+00 # zzq3 +################################### +## INFORMATION FOR CKMBLOCK +################################### +BLOCK CKMBLOCK # + 1 2.277360e-01 # cabi +################################### +## INFORMATION FOR MASS +################################### +BLOCK MASS # + 6 1.725000e+02 # mt + 15 1.777000e+00 # mta + 23 9.118760e+01 # mz + 25 1.250000e+02 # mh + 1 0.000000e+00 # d : 0.0 + 2 0.000000e+00 # u : 0.0 + 3 0.000000e+00 # s : 0.0 + 4 0.000000e+00 # c : 0.0 + 5 0.000000e+00 # b : 0.0 + 11 0.000000e+00 # e- : 0.0 + 12 0.000000e+00 # ve : 0.0 + 13 0.000000e+00 # mu- : 0.0 + 14 0.000000e+00 # vm : 0.0 + 16 0.000000e+00 # vt : 0.0 + 21 0.000000e+00 # g : 0.0 + 22 0.000000e+00 # a : 0.0 + 24 7.982436e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2))) +################################### +## INFORMATION FOR SMINPUTS +################################### +BLOCK SMINPUTS # + 1 1.279000e+02 # aewm1 + 2 1.166370e-05 # gf + 3 1.184000e-01 # as +################################### +## INFORMATION FOR YUKAWA +################################### +BLOCK YUKAWA # + 6 1.725000e+02 # ymt + 15 1.777000e+00 # ymtau +################################### +## INFORMATION FOR DECAY +################################### +DECAY 6 1.31 # wt +DECAY 23 2.495200e+00 # wz +DECAY 24 2.085000e+00 # ww +DECAY 25 4.070000e-03 # wh +DECAY 1 0.000000e+00 # d : 0.0 +DECAY 2 0.000000e+00 # u : 0.0 +DECAY 3 0.000000e+00 # s : 0.0 +DECAY 4 0.000000e+00 # c : 0.0 +DECAY 5 0.000000e+00 # b : 0.0 +DECAY 11 0.000000e+00 # e- : 0.0 +DECAY 12 0.000000e+00 # ve : 0.0 +DECAY 13 0.000000e+00 # mu- : 0.0 +DECAY 14 0.000000e+00 # vm : 0.0 +DECAY 15 0.000000e+00 # ta- : 0.0 +DECAY 16 0.000000e+00 # vt : 0.0 +DECAY 21 0.000000e+00 # g : 0.0 +DECAY 22 0.000000e+00 # a : 0.0 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/proc_card.dat new file mode 100644 index 000000000000..a08fbb4c4986 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/proc_card.dat @@ -0,0 +1,14 @@ +# Model +import model TopFCNC_UFO-all_ProdConventions_no_b_mass --modelname + +# Multiparticles +define j = g u u~ d d~ c c~ s s~ b b~ +define p = g u u~ d d~ c c~ s s~ b b~ + +# Processes +generate p p > t t~ NP=0 @10 +add process p p > t t~ j NP=0 @11 +add process p p > t t~ j j NP=0 @12 + +# Output folder +output PROCNAME -f -nojpeg diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/run_card.dat new file mode 100644 index 000000000000..610d5dde9d8a --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/TT_FCNC/orig_cards/run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# MadGraph5_aMC@NLO * +# * +# run_card.dat MadEvent * +# * +# This file is used to set the parameters of the run. * +# * +# Some notation/conventions: * +# * +# Lines starting with a '# ' are info or comments * +# * +# mind the format: value = variable ! comment * +#********************************************************************* +# +#******************* +# Running parameters +#******************* +# +#********************************************************************* +# Tag name for the run (one word) * +#********************************************************************* +tag_1 = run_tag ! name of the run +#********************************************************************* +# Run to generate the grid pack * +#********************************************************************* +True = gridpack !True = setting up the grid pack +#********************************************************************* +# Number of events and rnd seed * +# Warning: Do not generate more than 1M events in a single run * +# If you want to run Pythia, avoid more than 50k events in a run. * +#********************************************************************* +100000 = nevents ! Number of unweighted events requested +0 = iseed ! rnd seed (0=assigned automatically=default)) +#********************************************************************* +# Collider type and energy * +# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, * +# 3=photon from electron * +#********************************************************************* +1 = lpp1 ! beam 1 type +1 = lpp2 ! beam 2 type +6500.0 = ebeam1 ! beam 1 total energy in GeV +6500.0 = ebeam2 ! beam 2 total energy in GeV +#********************************************************************* +# Beam polarization from -100 (left-handed) to 100 (right-handed) * +#********************************************************************* +0.0 = polbeam1 ! beam polarization for beam 1 +0.0 = polbeam2 ! beam polarization for beam 2 +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* +lhapdf = pdlabel ! PDF set +$DEFAULT_PDF_SETS = lhaid +$DEFAULT_PDF_MEMBERS = reweight_PDF ! if pdlabel=lhapdf, this is the lhapdf number +#********************************************************************* +# Renormalization and factorization scales * +#********************************************************************* +False = fixed_ren_scale ! if .true. use fixed ren scale +False = fixed_fac_scale ! if .true. use fixed fac scale +91.188 = scale ! fixed ren scale +91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1 +91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2 +-1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices +1.0 = scalefact ! scale factor for event-by-event scales +#********************************************************************* +# Time of flight information. (-1 means not run) +#********************************************************************* +-1.0 = time_of_flight ! threshold below which info is not written +#********************************************************************* +# Matching - Warning! ickkw > 1 is still beta +#********************************************************************* +1 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching +1 = highestmult ! for ickkw=2, highest mult group +1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE +1.0 = alpsfact ! scale factor for QCD emission vx +False = chcluster ! cluster only according to channel diag +5 = asrwgtflavor ! highest quark flavor for a_s reweight +True = clusinfo ! include clustering tag in output +3.0 = lhe_version ! Change the way clustering information pass to shower. +#********************************************************************* +#********************************************************** +# +#********************************************************** +# Automatic ptj and mjj cuts if xqcut > 0 +# (turn off for VBF and single top processes) +#********************************************************** +True = auto_ptj_mjj ! Automatic setting of ptj and mjj +#********************************************************** +# +#********************************** +# BW cutoff (M+/-bwcutoff*Gamma) +#********************************** +15.0 = bwcutoff ! (M+/-bwcutoff*Gamma) +#********************************************************** +# Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not +# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) +#************************************************************* +False = cut_decays ! Cut decay products +#************************************************************* +# Number of helicities to sum per event (0 = all helicities) +# 0 gives more stable result, but longer run time (needed for +# long decay chains e.g.). +# Use >=2 if most helicities contribute, e.g. pure QCD. +#************************************************************* +0 = nhel ! Number of helicities used per event +#******************* +# Standard Cuts +#******************* +# +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* +0.0 = ptj ! minimum pt for the jets +0.0 = ptb ! minimum pt for the b +0.0 = pta ! minimum pt for the photons +0.0 = ptl ! minimum pt for the charged leptons +0.0 = misset ! minimum missing Et (sum of neutrino's momenta) +0.0 = ptheavy ! minimum pt for one heavy final state +-1.0 = ptjmax ! maximum pt for the jets +-1.0 = ptbmax ! maximum pt for the b +-1.0 = ptamax ! maximum pt for the photons +-1.0 = ptlmax ! maximum pt for the charged leptons +-1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* +0.0 = ej ! minimum E for the jets +0.0 = eb ! minimum E for the b +0.0 = ea ! minimum E for the photons +0.0 = el ! minimum E for the charged leptons +-1.0 = ejmax ! maximum E for the jets +-1.0 = ebmax ! maximum E for the b +-1.0 = eamax ! maximum E for the photons +-1.0 = elmax ! maximum E for the charged leptons +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* +5.0 = etaj ! max rap for the jets +-1.0 = etab ! max rap for the b +5.0 = etaa ! max rap for the photons +5.0 = etal ! max rap for the charged leptons +0.0 = etajmin ! min rap for the jets +0.0 = etabmin ! min rap for the b +0.0 = etaamin ! min rap for the photons +0.0 = etalmin ! main rap for the charged leptons +#********************************************************************* +# Minimum and maximum DeltaR distance * +#********************************************************************* +0.001 = drjj ! min distance between jets +0.0 = drbb ! min distance between b's +0.001 = drll ! min distance between leptons +0.001 = draa ! min distance between gammas +0.0 = drbj ! min distance between b and jet +0.001 = draj ! min distance between gamma and jet +0.001 = drjl ! min distance between jet and lepton +0.0 = drab ! min distance between gamma and b +0.0 = drbl ! min distance between b and lepton +0.001 = dral ! min distance between gamma and lepton +-1.0 = drjjmax ! max distance between jets +-1.0 = drbbmax ! max distance between b's +-1.0 = drllmax ! max distance between leptons +-1.0 = draamax ! max distance between gammas +-1.0 = drbjmax ! max distance between b and jet +-1.0 = drajmax ! max distance between gamma and jet +-1.0 = drjlmax ! max distance between jet and lepton +-1.0 = drabmax ! max distance between gamma and b +-1.0 = drblmax ! max distance between b and lepton +-1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +# WARNING: for four lepton final state mmll cut require to have * +# different lepton masses for each flavor! * +#********************************************************************* +0.0 = mmjj ! min invariant mass of a jet pair +0.0 = mmbb ! min invariant mass of a b pair +0.0 = mmaa ! min invariant mass of gamma gamma pair +0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair +-1.0 = mmjjmax ! max invariant mass of a jet pair +-1.0 = mmbbmax ! max invariant mass of a b pair +-1.0 = mmaamax ! max invariant mass of gamma gamma pair +-1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair +#********************************************************************* +# Minimum and maximum invariant mass for all letpons * +#********************************************************************* +0.0 = mmnl ! min invariant mass for all letpons (l+- and vl) +-1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl) +#********************************************************************* +# Minimum and maximum pt for 4-momenta sum of leptons * +#********************************************************************* +0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl) +-1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl) +#********************************************************************* +# Inclusive cuts * +#********************************************************************* +0.0 = xptj ! minimum pt for at least one jet +0.0 = xptb ! minimum pt for at least one b +0.0 = xpta ! minimum pt for at least one photon +0.0 = xptl ! minimum pt for at least one charged lepton +#********************************************************************* +# Control the pt's of the jets sorted by pt * +#********************************************************************* +0.0 = ptj1min ! minimum pt for the leading jet in pt +0.0 = ptj2min ! minimum pt for the second jet in pt +0.0 = ptj3min ! minimum pt for the third jet in pt +0.0 = ptj4min ! minimum pt for the fourth jet in pt +-1.0 = ptj1max ! maximum pt for the leading jet in pt +-1.0 = ptj2max ! maximum pt for the second jet in pt +-1.0 = ptj3max ! maximum pt for the third jet in pt +-1.0 = ptj4max ! maximum pt for the fourth jet in pt +0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts +#********************************************************************* +# Control the pt's of leptons sorted by pt * +#********************************************************************* +0.0 = ptl1min ! minimum pt for the leading lepton in pt +0.0 = ptl2min ! minimum pt for the second lepton in pt +0.0 = ptl3min ! minimum pt for the third lepton in pt +0.0 = ptl4min ! minimum pt for the fourth lepton in pt +-1.0 = ptl1max ! maximum pt for the leading lepton in pt +-1.0 = ptl2max ! maximum pt for the second lepton in pt +-1.0 = ptl3max ! maximum pt for the third lepton in pt +-1.0 = ptl4max ! maximum pt for the fourth lepton in pt +#********************************************************************* +# Control the Ht(k)=Sum of k leading jets * +#********************************************************************* +0.0 = htjmin ! minimum jet HT=Sum(jet pt) +-1.0 = htjmax ! maximum jet HT=Sum(jet pt) +0.0 = ihtmin !inclusive Ht for all partons (including b) +-1.0 = ihtmax !inclusive Ht for all partons (including b) +0.0 = ht2min ! minimum Ht for the two leading jets +0.0 = ht3min ! minimum Ht for the three leading jets +0.0 = ht4min ! minimum Ht for the four leading jets +-1.0 = ht2max ! maximum Ht for the two leading jets +-1.0 = ht3max ! maximum Ht for the three leading jets +-1.0 = ht4max ! maximum Ht for the four leading jets +#*********************************************************************** +# Photon-isolation cuts, according to hep-ph/9801442 * +# When ptgmin=0, all the other parameters are ignored * +# When ptgmin>0, pta and draj are not going to be used * +#*********************************************************************** +0.0 = ptgmin ! Min photon transverse momentum +0.4 = R0gamma ! Radius of isolation code +1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442 +1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442 +True = isoEM ! isolate photons from EM energy (photons and leptons) +#********************************************************************* +# WBF cuts * +#********************************************************************* +0.0 = xetamin ! minimum rapidity for two jets in the WBF case +0.0 = deltaeta ! minimum rapidity for two jets in the WBF case +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* +-1.0 = ktdurham +0.4 = dparameter +#********************************************************************* +# maximal pdg code for quark to be considered as a light jet * +# (otherwise b cuts are applied) * +#********************************************************************* +5 = maxjetflavor ! Maximum jet pdg code +#********************************************************************* +# Jet measure cuts * +#********************************************************************* +20.0 = xqcut ! minimum kt jet measure between partons +#*********************************************************************