p p > t t~

Hi,

In your process, you do not have any final state jet in the final state.
(you "only" have one top and an anti-top).

So your cut
 3000.0 = ptj1min ! minimum pt for the leading jet in pt
 3000.0 = ptj2min ! minimum pt for the second jet in pt

Does not make any sense and likely cut 100% of your process.

You either have to apply such cut after parton shower (and/or decay) but not at this level of the computation.

Cheers,

Olivier

PS: In principle such cut is not even displayed in the run_card for such process for that precise reason that it does not make sense to specify cut on "jet" for your process.

> On 1 Aug 2021, at 07:45, Cherry Parba <email address hidden> wrote:
>
> New question #698211 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/698211
>
> Good day Madgraph expert,
>
> I am trying to simulate, p p > t t~ at center of mass energy of 100 TeV and i want to put some cuts. I want that atleast 2 jets in an event has pt > 3000 GeV, |eta| < 3.0 and deltaEta(eta1,eta2) < 2.4 .
>
> But in does not work,
>
> Survey return zero cross section.
> Typical reasons are the following:
> 1) A massive s-channel particle has a width set to zero.
> 2) The pdf are zero for at least one of the initial state particles
> or you are using maxjetflavor=4 for initial state b:s.
> 3) The cuts are too strong.
> Please check/correct your param_card and/or your run_card.
>
>
> This is the banner text file that i used,
>
> <LesHouchesEvents version="3.0">
> <header>
> <!--
> #*********************************************************************
> # *
> # MadGraph5_aMC@NLO *
> # *
> # Going Beyond *
> # *
> # http://madgraph.hep.uiuc.edu *
> # http://madgraph.phys.ucl.ac.be *
> # http://amcatnlo.cern.ch *
> # *
> # The MadGraph5_aMC@NLO team *
> # *
> #....................................................................*
> # *
> # This file contains all the information necessary to reproduce *
> # the events generated: *
> # *
> # 1. software version *
> # 2. proc_card : code generation info including model *
> # 3. param_card : model primary parameters in the LH format *
> # 4. run_card : running parameters (collider and cuts) *
> # 5. pythia_card : present only if pythia has been run *
> # 6. pgs_card : present only if pgs has been run *
> # 7. delphes_cards : present only if delphes has been run *
> # *
> # *
> #*********************************************************************
> -->
> <MGVersion>
> 3.1.1
> </MGVersion>
> <MG5ProcCard>
> <![CDATA[
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.1.1 2021-05-28 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #************************************************************
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set max_t_for_channel 99
> set loop_optimized_output True
> set low_mem_multicore_nlo_generation False
> set default_unset_couplings 99
> set include_lepton_initiated_processes False
> set zerowidth_tchannel True
> set nlo_mixed_expansion True
> set loop_color_flows False
> set gauge unitary
> set complex_mass_scheme False
> set max_npoint_for_channel 0
> 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~
> define l+ = e+ mu+
> define l- = e- mu-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> define p = p b b~
> define j = j b b~
> generate p p > t t~
> output ttbar1_background
> ]]>
> </MG5ProcCard>
> <MGRunCard>
> <![CDATA[
> #*********************************************************************
> # 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 *
> # *
> # To display more options, you can type the command: *
> # update full_run_card *
> #*********************************************************************
> #
> #*********************************************************************
> # Tag name for the run (one word) *
> #*********************************************************************
> tag_1 = run_tag ! name of the run
> #*********************************************************************
> # Number of events and rnd seed *
> # Warning: Do not generate more than 1M events in a single run *
> #*********************************************************************
> 1071 = 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, 4=photon from muon *
> #*********************************************************************
> 1 = lpp1 ! beam 1 type
> 1 = lpp2 ! beam 2 type
> 50000.0 = ebeam1 ! beam 1 total energy in GeV
> 50000.0 = ebeam2 ! beam 2 total energy in GeV
> # To see polarised beam options: type "update beam_pol"
>
> #*********************************************************************
> # PDF CHOICE: this automatically fixes also alpha_s and its evol. *
> #*********************************************************************
> nn23lo1 = pdlabel ! PDF set
> 230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
> # To see heavy ion options: type "update ion_pdf"
> #*********************************************************************
> # 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
> #*********************************************************************
> # Type and output format
> #*********************************************************************
> False = gridpack !True = setting up the grid pack
> -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
> average = event_norm ! average/sum. Normalization of the weight in the LHEF
> # To see MLM/CKKW merging options: type "update MLM" or "update CKKW"
>
> #*********************************************************************
> #
> #*********************************************************************
> # Phase-Space Optimization strategy (basic options)
> #*********************************************************************
> 0 = nhel ! using helicities importance sampling or not.
> ! 0: sum over helicity, 1: importance sampling
> 2 = sde_strategy ! default integration strategy (hep-ph/2021.xxxxx)
> ! 1 is old strategy (using amp square)
> ! 2 is new strategy (using only the denominator)
> # To see advanced option for Phase-Space optimization: type "update psoptim"
> #*********************************************************************
> # Generation bias, check the wiki page below for more information: *
> # 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
> #*********************************************************************
> None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
> {} = bias_parameters ! Specifies the parameters of the module.
> #
> #*******************************
> # Parton level cuts definition *
> #*******************************
> #
> #
> #*********************************************************************
> # BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
> #*********************************************************************
> 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
> #*********************************************************************
> # Standard Cuts *
> #*********************************************************************
> # Minimum and maximum pt's (for max, -1 means no cut) *
> #*********************************************************************
> 20.0 = ptj ! minimum pt for the jets
> 10.0 = ptl ! minimum pt for the charged leptons
> 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
> -1.0 = ptjmax ! maximum pt for the jets
> -1.0 = ptlmax ! maximum pt for the charged leptons
> -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
> {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
>
> #*********************************************************************
> # 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
> {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})
>
> #*********************************************************************
> # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
> #*********************************************************************
> 3.0 = etaj ! max rap for the jets
> 2.5 = etal ! max rap for the charged leptons
> 0.0 = etalmin ! main rap for the charged leptons
> {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
> #*********************************************************************
> # Minimum and maximum DeltaR distance *
> #*********************************************************************
> 0.4 = drjj ! min distance between jets
> 0.4 = drjl ! min distance between jet and lepton
> -1.0 = drjjmax ! max distance between jets
> -1.0 = drjlmax ! max distance between jet and lepton
> #*********************************************************************
> # Minimum and maximum invariant mass for pairs *
> #*********************************************************************
> 0.0 = mmjj ! min invariant mass of a jet pair
> -1.0 = mmjjmax ! max invariant mass of a jet pair
> {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
> {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
> ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
> #*********************************************************************
> # 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 / neutrino *
> # for pair of lepton includes only same flavor, opposite charge
> #*********************************************************************
> #*********************************************************************
> # Inclusive cuts *
> #*********************************************************************
> 0.0 = xptj ! minimum pt for at least one jet
> #*********************************************************************
> # Control the pt's of the jets sorted by pt *
> #*********************************************************************
> 3000.0 = ptj1min ! minimum pt for the leading jet in pt
> 3000.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 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
> #*********************************************************************
> # WBF cuts *
> #*********************************************************************
> 0.0 = xetamin ! minimum rapidity for two jets in the WBF case
> 2.4 = deltaeta ! minimum rapidity for two jets in the WBF case
> #*********************************************************************
> # maximal pdg code for quark to be considered as a light jet *
> # (otherwise b cuts are applied) *
> #*********************************************************************
> 5 = maxjetflavor ! Maximum jet pdg code
> #*********************************************************************
> #
> #*********************************************************************
> # Store info for systematics studies *
> # WARNING: Do not use for interference type of computation *
> #*********************************************************************
> True = use_syst ! Enable systematics studies
> #
> systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
> ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
> # Syscalc is deprecated but to see the associate options type'update syscalc'
> ]]>
> </MGRunCard>
> <slha>
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
> ######################################################################
> ## ##
> ## Width set on Auto will be computed following the information ##
> ## present in the decay.py files of the model. ##
> ## See arXiv:1402.1178 for more details. ##
> ## ##
> ######################################################################
>
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> Block mass
> 5 4.700000e+00 # MB
> 6 1.730000e+02 # MT
> 15 1.777000e+00 # MTA
> 23 9.118800e+01 # MZ
> 25 1.250000e+02 # MH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> 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
> 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 8.041900e+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.325070e+02 # aEWM1
> 2 1.166390e-05 # Gf
> 3 1.180000e-01 # aS (Note that Parameter not used if you use a PDF set)
>
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> Block yukawa
> 5 4.700000e+00 # ymb
> 6 1.730000e+02 # ymt
> 15 1.777000e+00 # ymtau
>
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.491500e+00 # WT
> DECAY 23 2.441404e+00 # WZ
> DECAY 24 2.047600e+00 # WW
> DECAY 25 6.382339e-03 # WH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> 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
> </slha>
> </header>
> </LesHouchesEvents>
>
> would you mind, helping me to fixed the cuts .
>
>
>
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Ah so that is why, when i do , p p > t t~, (t > b w+, w+ > mu vm), (t ~ > b~ w-, w- > j j) it works. So at generator level, to make it sense, i will cut particles that is originally involve on the process.