zb jet and z c jet

Dear Miranda,

I’m confuse but I do not see why the cross-section should be the same since
1) your quark are massless (and therefore the mass of the b is the same as the mass of the c)
2) you only have QCD diagram and therefore b and c are the same.
3) you have the same cut

Cheers,

Olivier

On 22 Jul 2015, at 08:12, miranda <email address hidden> wrote:

> New question #269530 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/269530
>
> Dear authors,
> I want to calculate cross section of Z bozon plus heavy quark jets QQ~ (p p~ > z Q Q~).
> In fact I want to compare the result of Zcc~ and Zbb~ with the following run_card.dat. I expect to have R_cc~/bb~ ∝ 2 but I have the same value for cross section of Zcc~ and Zbb~
> For this aim
> Import model sm-no_b_mass
> Define p = p b b~
> Define j = c c~ (In case of Zbb~ I define j = b b~)
> Generate p p > z j j
> I have kinematic region 20<pT<80 and \eta< 2.5 as you can see in run card.
> Also in case of bb~ I set maxjetflavor=5
> I'm confused. Where is my problem?
> Could you please help me with this issue ?
>
> Thank you in advance for your consideration
>
>
> #*********************************************************************
> # 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 *
> #*********************************************************************
> .false. = 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. *
> #*********************************************************************
> 10000 = 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
> 980 = ebeam1 ! beam 1 total energy in GeV
> 980 = ebeam2 ! beam 2 total energy in GeV
> #*********************************************************************
> # Beam polarization from -100 (left-handed) to 100 (right-handed) *
> #*********************************************************************
> 0 = polbeam1 ! beam polarization for beam 1
> 0 = polbeam2 ! beam polarization for beam 2
> #*********************************************************************
> # PDF CHOICE: this automatically fixes also alpha_s and its evol. *
> #*********************************************************************
> 'lhapdf' = pdlabel ! PDF set
> 10550 = lhaid
> #*********************************************************************
> # Renormalization and factorization scales *
> #*********************************************************************
> F = fixed_ren_scale ! if .true. use fixed ren scale
> F = fixed_fac_scale ! if .true. use fixed fac scale
> 91.1880 = scale ! fixed ren scale
> 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
> 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
> 1 = scalefact ! scale factor for event-by-event scales
> #*********************************************************************
> # 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
> 2 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
> 1 = alpsfact ! scale factor for QCD emission vx
> F = chcluster ! cluster only according to channel diag
> T = pdfwgt ! for ickkw=1, perform pdf reweighting
> 5 = asrwgtflavor ! highest quark flavor for a_s reweight
> T = clusinfo ! include clustering tag in output
> #*********************************************************************
> #**********************************************************
> #
> #**********************************************************
> # Automatic ptj and mjj cuts if xqcut > 0
> # (turn off for VBF and single top processes)
> #**********************************************************
> T = auto_ptj_mjj ! Automatic setting of ptj and mjj
> #**********************************************************
> #
> #**********************************
> # BW cutoff (M+/-bwcutoff*Gamma)
> #**********************************
> 15 = bwcutoff ! (M+/-bwcutoff*Gamma)
> #**********************************************************
> # Apply pt/E/eta/dr/mij cuts on decay products or not
> # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
> #**********************************************************
> T = 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) *
> #*********************************************************************
> 20 = ptj ! minimum pt for the jets
> 0 = ptb ! minimum pt for the b
> 0 = pta ! minimum pt for the photons
> 0 = ptl ! minimum pt for the charged leptons
> 0 = misset ! minimum missing Et (sum of neutrino's momenta)
> 0 = ptheavy ! minimum pt for one heavy final state
> 0 = ptonium ! minimum pt for the quarkonium states
> 90 = ptjmax ! maximum pt for the jets
> -1 = ptbmax ! maximum pt for the b
> -1 = ptamax ! maximum pt for the photons
> -1 = ptlmax ! maximum pt for the charged leptons
> -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
> #*********************************************************************
> # Minimum and maximum E's (in the lab frame) *
> #*********************************************************************
> 0 = ej ! minimum E for the jets
> 0 = eb ! minimum E for the b
> 0 = ea ! minimum E for the photons
> 0 = el ! minimum E for the charged leptons
> -1 = ejmax ! maximum E for the jets
> -1 = ebmax ! maximum E for the b
> -1 = eamax ! maximum E for the photons
> -1 = elmax ! maximum E for the charged leptons
> #*********************************************************************
> # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
> #*********************************************************************
> 2.5 = etaj ! max rap for the jets
> 0 = etab ! max rap for the b
> 0 = etaa ! max rap for the photons
> 0 = etal ! max rap for the charged leptons
> 0 = etaonium ! max rap for the quarkonium states
> 0 = etajmin ! min rap for the jets
> 0 = etabmin ! min rap for the b
> 0 = etaamin ! min rap for the photons
> 0 = etalmin ! main rap for the charged leptons
> #*********************************************************************
> # Minimum and maximum DeltaR distance *
> #*********************************************************************
> 0.5 = drjj ! min distance between jets
> 0.0 = drbb ! min distance between b's
> 0.0 = drll ! min distance between leptons
> 0.0 = draa ! min distance between gammas
> 0 = drbj ! min distance between b and jet
> 0.0 = draj ! min distance between gamma and jet
> 0.0 = drjl ! min distance between jet and lepton
> 0 = drab ! min distance between gamma and b
> 0 = drbl ! min distance between b and lepton
> 0.0 = dral ! min distance between gamma and lepton
> -1 = drjjmax ! max distance between jets
> -1 = drbbmax ! max distance between b's
> -1 = drllmax ! max distance between leptons
> -1 = draamax ! max distance between gammas
> -1 = drbjmax ! max distance between b and jet
> -1 = drajmax ! max distance between gamma and jet
> -1 = drjlmax ! max distance between jet and lepton
> -1 = drabmax ! max distance between gamma and b
> -1 = drblmax ! max distance between b and lepton
> -1 = dralmax ! maxdistance between gamma and lepton
> #*********************************************************************
> # Minimum and maximum invariant mass for pairs *
> #*********************************************************************
> 0 = mmjj ! min invariant mass of a jet pair
> 0 = mmbb ! min invariant mass of a b pair
> 0 = mmaa ! min invariant mass of gamma gamma pair
> 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
> -1 = mmjjmax ! max invariant mass of a jet pair
> -1 = mmbbmax ! max invariant mass of a b pair
> -1 = mmaamax ! max invariant mass of gamma gamma pair
> -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
> #*********************************************************************
> # Minimum and maximum invariant mass for all letpons *
> #*********************************************************************
> 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
> -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
> #*********************************************************************
> # Minimum and maximum pt for 4-momenta sum of leptons *
> #*********************************************************************
> 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
> -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
> #*********************************************************************
> # Inclusive cuts *
> #*********************************************************************
> 0 = xptj ! minimum pt for at least one jet
> 0 = xptb ! minimum pt for at least one b
> 0 = xpta ! minimum pt for at least one photon
> 0 = xptl ! minimum pt for at least one charged lepton
> #*********************************************************************
> # Control the pt's of the jets sorted by pt *
> #*********************************************************************
> 0 = ptj1min ! minimum pt for the leading jet in pt
> 0 = ptj2min ! minimum pt for the second jet in pt
> 0 = ptj3min ! minimum pt for the third jet in pt
> 0 = ptj4min ! minimum pt for the fourth jet in pt
> -1 = ptj1max ! maximum pt for the leading jet in pt
> -1 = ptj2max ! maximum pt for the second jet in pt
> -1 = ptj3max ! maximum pt for the third jet in pt
> -1 = 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 = ptl1min ! minimum pt for the leading lepton in pt
> 0 = ptl2min ! minimum pt for the second lepton in pt
> 0 = ptl3min ! minimum pt for the third lepton in pt
> 0 = ptl4min ! minimum pt for the fourth lepton in pt
> -1 = ptl1max ! maximum pt for the leading lepton in pt
> -1 = ptl2max ! maximum pt for the second lepton in pt
> -1 = ptl3max ! maximum pt for the third lepton in pt
> -1 = ptl4max ! maximum pt for the fourth lepton in pt
> #*********************************************************************
> # Control the Ht(k)=Sum of k leading jets *
> #*********************************************************************
> 0 = htjmin ! minimum jet HT=Sum(jet pt)
> -1 = htjmax ! maximum jet HT=Sum(jet pt)
> 0 = ihtmin !inclusive Ht for all partons (including b)
> -1 = ihtmax !inclusive Ht for all partons (including b)
> 0 = ht2min ! minimum Ht for the two leading jets
> 0 = ht3min ! minimum Ht for the three leading jets
> 0 = ht4min ! minimum Ht for the four leading jets
> -1 = ht2max ! maximum Ht for the two leading jets
> -1 = ht3max ! maximum Ht for the three leading jets
> -1 = 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 = 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 = xetamin ! minimum rapidity for two jets in the WBF case
> 0 = deltaeta ! minimum rapidity for two jets in the WBF case
> #*********************************************************************
> # KT DURHAM CUT *
> #*********************************************************************
> -1 = 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 = xqcut ! minimum kt jet measure between partons
> #*********************************************************************
> #
> #*********************************************************************
> # Store info for systematics studies *
> # WARNING: If use_syst is T, matched Pythia output is *
> # meaningful ONLY if plotted taking matchscale *
> # reweighting into account! *
> #*********************************************************************
> F = use_syst ! Enable systematics studies
> #
> #**************************************
> # Parameter of the systematics study
> # will be use by SysCalc (if install)
> #**************************************
> #
> 0.5 1 2 = sys_scalefact # Central scale factors
> 0.5 1 2 = sys_alpsfact # \alpha_s emission scale factors
> 30 50 = sys_matchscale # variation of merging scale
> # PDF sets and number of members (0 or none for all members).
> CT10nlo.LHgrid = sys_pdf # matching scales
> # MSTW2008nlo68cl.LHgrid 1 = sys_pdf
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

I truly appreciate your help in solving my problem.

Best Regards
Miranda

On Wed, Jul 22, 2015 at 4:16 PM, Olivier Mattelaer <
<email address hidden>> wrote:

> Your question #269530 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/269530
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Dear Miranda,
>
> I’m confuse but I do not see why the cross-section should be the same since
> 1) your quark are massless (and therefore the mass of the b is the same as
> the mass of the c)
> 2) you only have QCD diagram and therefore b and c are the same.
> 3) you have the same cut
>
> Cheers,
>
> Olivier
>
> On 22 Jul 2015, at 08:12, miranda <email address hidden>
> wrote:
>
> > New question #269530 on MadGraph5_aMC@NLO:
> > https://answers.launchpad.net/mg5amcnlo/+question/269530
> >
> > Dear authors,
> > I want to calculate cross section of Z bozon plus heavy quark jets QQ~
> (p p~ > z Q Q~).
> > In fact I want to compare the result of Zcc~ and Zbb~ with the following
> run_card.dat. I expect to have R_cc~/bb~ ∝ 2 but I have the same value for
> cross section of Zcc~ and Zbb~
> > For this aim
> > Import model sm-no_b_mass
> > Define p = p b b~
> > Define j = c c~ (In case of Zbb~ I define j = b b~)
> > Generate p p > z j j
> > I have kinematic region 20<pT<80 and \eta< 2.5 as you can see in run
> card.
> > Also in case of bb~ I set maxjetflavor=5
> > I'm confused. Where is my problem?
> > Could you please help me with this issue ?
> >
> > Thank you in advance for your consideration
> >
> >
> > #*********************************************************************
> > # 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 *
> > #*********************************************************************
> > .false. = 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. *
> > #*********************************************************************
> > 10000 = 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
> > 980 = ebeam1 ! beam 1 total energy in GeV
> > 980 = ebeam2 ! beam 2 total energy in GeV
> > #*********************************************************************
> > # Beam polarization from -100 (left-handed) to 100 (right-handed) *
> > #*********************************************************************
> > 0 = polbeam1 ! beam polarization for beam 1
> > 0 = polbeam2 ! beam polarization for beam 2
> > #*********************************************************************
> > # PDF CHOICE: this automatically fixes also alpha_s and its evol. *
> > #*********************************************************************
> > 'lhapdf' = pdlabel ! PDF set
> > 10550 = lhaid
> > #*********************************************************************
> > # Renormalization and factorization scales *
> > #*********************************************************************
> > F = fixed_ren_scale ! if .true. use fixed ren scale
> > F = fixed_fac_scale ! if .true. use fixed fac scale
> > 91.1880 = scale ! fixed ren scale
> > 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
> > 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
> > 1 = scalefact ! scale factor for event-by-event scales
> > #*********************************************************************
> > # 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
> > 2 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
> > 1 = alpsfact ! scale factor for QCD emission vx
> > F = chcluster ! cluster only according to channel diag
> > T = pdfwgt ! for ickkw=1, perform pdf reweighting
> > 5 = asrwgtflavor ! highest quark flavor for a_s reweight
> > T = clusinfo ! include clustering tag in output
> > #*********************************************************************
> > #**********************************************************
> > #
> > #**********************************************************
> > # Automatic ptj and mjj cuts if xqcut > 0
> > # (turn off for VBF and single top processes)
> > #**********************************************************
> > T = auto_ptj_mjj ! Automatic setting of ptj and mjj
> > #**********************************************************
> > #
> > #**********************************
> > # BW cutoff (M+/-bwcutoff*Gamma)
> > #**********************************
> > 15 = bwcutoff ! (M+/-bwcutoff*Gamma)
> > #**********************************************************
> > # Apply pt/E/eta/dr/mij cuts on decay products or not
> > # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
> > #**********************************************************
> > T = 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) *
> > #*********************************************************************
> > 20 = ptj ! minimum pt for the jets
> > 0 = ptb ! minimum pt for the b
> > 0 = pta ! minimum pt for the photons
> > 0 = ptl ! minimum pt for the charged leptons
> > 0 = misset ! minimum missing Et (sum of neutrino's momenta)
> > 0 = ptheavy ! minimum pt for one heavy final state
> > 0 = ptonium ! minimum pt for the quarkonium states
> > 90 = ptjmax ! maximum pt for the jets
> > -1 = ptbmax ! maximum pt for the b
> > -1 = ptamax ! maximum pt for the photons
> > -1 = ptlmax ! maximum pt for the charged leptons
> > -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
> > #*********************************************************************
> > # Minimum and maximum E's (in the lab frame) *
> > #*********************************************************************
> > 0 = ej ! minimum E for the jets
> > 0 = eb ! minimum E for the b
> > 0 = ea ! minimum E for the photons
> > 0 = el ! minimum E for the charged leptons
> > -1 = ejmax ! maximum E for the jets
> > -1 = ebmax ! maximum E for the b
> > -1 = eamax ! maximum E for the photons
> > -1 = elmax ! maximum E for the charged leptons
> > #*********************************************************************
> > # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
> > #*********************************************************************
> > 2.5 = etaj ! max rap for the jets
> > 0 = etab ! max rap for the b
> > 0 = etaa ! max rap for the photons
> > 0 = etal ! max rap for the charged leptons
> > 0 = etaonium ! max rap for the quarkonium states
> > 0 = etajmin ! min rap for the jets
> > 0 = etabmin ! min rap for the b
> > 0 = etaamin ! min rap for the photons
> > 0 = etalmin ! main rap for the charged leptons
> > #*********************************************************************
> > # Minimum and maximum DeltaR distance *
> > #*********************************************************************
> > 0.5 = drjj ! min distance between jets
> > 0.0 = drbb ! min distance between b's
> > 0.0 = drll ! min distance between leptons
> > 0.0 = draa ! min distance between gammas
> > 0 = drbj ! min distance between b and jet
> > 0.0 = draj ! min distance between gamma and jet
> > 0.0 = drjl ! min distance between jet and lepton
> > 0 = drab ! min distance between gamma and b
> > 0 = drbl ! min distance between b and lepton
> > 0.0 = dral ! min distance between gamma and lepton
> > -1 = drjjmax ! max distance between jets
> > -1 = drbbmax ! max distance between b's
> > -1 = drllmax ! max distance between leptons
> > -1 = draamax ! max distance between gammas
> > -1 = drbjmax ! max distance between b and jet
> > -1 = drajmax ! max distance between gamma and jet
> > -1 = drjlmax ! max distance between jet and lepton
> > -1 = drabmax ! max distance between gamma and b
> > -1 = drblmax ! max distance between b and lepton
> > -1 = dralmax ! maxdistance between gamma and lepton
> > #*********************************************************************
> > # Minimum and maximum invariant mass for pairs *
> > #*********************************************************************
> > 0 = mmjj ! min invariant mass of a jet pair
> > 0 = mmbb ! min invariant mass of a b pair
> > 0 = mmaa ! min invariant mass of gamma gamma pair
> > 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
> > -1 = mmjjmax ! max invariant mass of a jet pair
> > -1 = mmbbmax ! max invariant mass of a b pair
> > -1 = mmaamax ! max invariant mass of gamma gamma pair
> > -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
> > #*********************************************************************
> > # Minimum and maximum invariant mass for all letpons *
> > #*********************************************************************
> > 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
> > -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
> > #*********************************************************************
> > # Minimum and maximum pt for 4-momenta sum of leptons *
> > #*********************************************************************
> > 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
> > -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
> > #*********************************************************************
> > # Inclusive cuts *
> > #*********************************************************************
> > 0 = xptj ! minimum pt for at least one jet
> > 0 = xptb ! minimum pt for at least one b
> > 0 = xpta ! minimum pt for at least one photon
> > 0 = xptl ! minimum pt for at least one charged lepton
> > #*********************************************************************
> > # Control the pt's of the jets sorted by pt *
> > #*********************************************************************
> > 0 = ptj1min ! minimum pt for the leading jet in pt
> > 0 = ptj2min ! minimum pt for the second jet in pt
> > 0 = ptj3min ! minimum pt for the third jet in pt
> > 0 = ptj4min ! minimum pt for the fourth jet in pt
> > -1 = ptj1max ! maximum pt for the leading jet in pt
> > -1 = ptj2max ! maximum pt for the second jet in pt
> > -1 = ptj3max ! maximum pt for the third jet in pt
> > -1 = 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 = ptl1min ! minimum pt for the leading lepton in pt
> > 0 = ptl2min ! minimum pt for the second lepton in pt
> > 0 = ptl3min ! minimum pt for the third lepton in pt
> > 0 = ptl4min ! minimum pt for the fourth lepton in pt
> > -1 = ptl1max ! maximum pt for the leading lepton in pt
> > -1 = ptl2max ! maximum pt for the second lepton in pt
> > -1 = ptl3max ! maximum pt for the third lepton in pt
> > -1 = ptl4max ! maximum pt for the fourth lepton in pt
> > #*********************************************************************
> > # Control the Ht(k)=Sum of k leading jets *
> > #*********************************************************************
> > 0 = htjmin ! minimum jet HT=Sum(jet pt)
> > -1 = htjmax ! maximum jet HT=Sum(jet pt)
> > 0 = ihtmin !inclusive Ht for all partons (including b)
> > -1 = ihtmax !inclusive Ht for all partons (including b)
> > 0 = ht2min ! minimum Ht for the two leading jets
> > 0 = ht3min ! minimum Ht for the three leading jets
> > 0 = ht4min ! minimum Ht for the four leading jets
> > -1 = ht2max ! maximum Ht for the two leading jets
> > -1 = ht3max ! maximum Ht for the three leading jets
> > -1 = 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 = 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 = xetamin ! minimum rapidity for two jets in the WBF case
> > 0 = deltaeta ! minimum rapidity for two jets in the WBF case
> > #*********************************************************************
> > # KT DURHAM CUT *
> > #*********************************************************************
> > -1 = 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 = xqcut ! minimum kt jet measure between partons
> > #*********************************************************************
> > #
> > #*********************************************************************
> > # Store info for systematics studies *
> > # WARNING: If use_syst is T, matched Pythia output is *
> > # meaningful ONLY if plotted taking matchscale *
> > # reweighting into account! *
> > #*********************************************************************
> > F = use_syst ! Enable systematics studies
> > #
> > #**************************************
> > # Parameter of the systematics study
> > # will be use by SysCalc (if install)
> > #**************************************
> > #
> > 0.5 1 2 = sys_scalefact # Central scale factors
> > 0.5 1 2 = sys_alpsfact # \alpha_s emission scale factors
> > 30 50 = sys_matchscale # variation of merging scale
> > # PDF sets and number of members (0 or none for all members).
> > CT10nlo.LHgrid = sys_pdf # matching scales
> > # MSTW2008nlo68cl.LHgrid 1 = sys_pdf
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadGraph5_aMC@NLO.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/mg5amcnlo/+question/269530/+confirm?answer_id=0
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/mg5amcnlo/+question/269530
>
> You received this question notification because you asked the question.
>