About NLO model file loop_qcd_qed_sm
Hi,
I am using MG5_v2_3_3 .
I import an inbuilt model : loop_qcd_qed_sm
and generate z > b b~ / no bsm particles [QED].
Thus it generates the process.
output <fn>
While launch, it shows the following error :
>0
INFO: Starting run
INFO: Compiling the code
INFO: Lepton-Lepton collision: Ignoring 'pdlabel' and 'lhaid' in the run_card.
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 4 cores
INFO: Compiling P0_z_bbx_
INFO: P0_z_bbx_
INFO: Checking test output:
INFO: P0_z_bbx_
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: Starting run
INFO: Using 4 cores
INFO: Cleaning previous results
INFO: Doing fixed order NLO
WARNING: No integration channels found for contribution P0_z_bbx_
INFO: Setting up grids
INFO: Idle: 0, Running: 0, Completed: 0 [ current time: 00h34 ]
Command "launch auto " interrupted with error:
ZeroDivisionError : float division by zero
Please report this bug on https:/
More information is found in '/home/
Please attach this file to your report.
INFO:
quit
INFO:
MG5_aMC>
Thanks,
Subhasish
Question information
- Language:
- English Edit question
- Status:
- Open
- Assignee:
- Rikkert Frederix Edit question
- Last query:
- Last reply:
Revision history for this message
![]() |
#1 |
Dear Subhasish,
What it is that you are trying to compute? Your current setup does not make sense to me.
best,
Rikkert
Revision history for this message
![]() |
#2 |
Dear Rikkert,
I was trying to compute the decay width of z > b b~.
Where as I have written in the first mail. "bsm particles" means all other
particles in the model which are not SM particles. May be goldstone bosons
or else. I am using MG5_v2_3_3.
It is the NLO processes which contain only QED-vertices only.
Any other information you want from me ?
Thanks,
Subhasish
On Tue, May 24, 2016 at 1:33 PM, Rikkert Frederix <
<email address hidden>> wrote:
> Your question #294100 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Open => Needs information
>
> Rikkert Frederix requested more information:
> Dear Subhasish,
>
> What it is that you are trying to compute? Your current setup does not
> make sense to me.
>
> best,
> Rikkert
>
> --
> To answer this request for more information, you can either reply to
> this email or enter your reply at the following page:
> https:/
>
> You received this question notification because you asked the question.
>
Revision history for this message
![]() |
#3 |
Dear Subhasish,
Ghosts are not BSM particles, but rather part of the Standard Model in the Feynman gauge. Removing them destroys gauge invariance. Moreover, the model you loaded is a model using a 5 flavour scheme. Hence, the bottom quark is indistinguishable from any of the other light quarks, which makes your process definition ill-defined. Also, you have to be very careful in setting the widths of the particles to keep gauge invariance, in particular in relation with the incoming Z-boson.
However, note that NLO electroweak corrections are not supported by any of the official public versions of MadGraph5_aMC@NLO. They are only available in the development versions of the code, which are available on launchpad website, but they require quite some expertise to use.
best,
Rikkert
Revision history for this message
![]() |
#4 |
Just for completeness,
I have fully forbidden the use of [QED] in the development version (with appropriate error message)
such that the status will be clear for the user.
Cheers,
Olivier
> On May 24, 2016, at 10:32, Rikkert Frederix <email address hidden> wrote:
>
> Question #294100 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Open => Answered
>
> Rikkert Frederix proposed the following answer:
> Dear Subhasish,
>
> Ghosts are not BSM particles, but rather part of the Standard Model in
> the Feynman gauge. Removing them destroys gauge invariance. Moreover,
> the model you loaded is a model using a 5 flavour scheme. Hence, the
> bottom quark is indistinguishable from any of the other light quarks,
> which makes your process definition ill-defined. Also, you have to be
> very careful in setting the widths of the particles to keep gauge
> invariance, in particular in relation with the incoming Z-boson.
>
> However, note that NLO electroweak corrections are not supported by any
> of the official public versions of MadGraph5_aMC@NLO. They are only
> available in the development versions of the code, which are available
> on launchpad website, but they require quite some expertise to use.
>
> best,
> Rikkert
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Revision history for this message
![]() |
#5 |
Yes that's right...
Why am I sending further doubt is because I am getting some decay width,
when I was generating the same process with all the particles taken into
consideration in the loop or outside; but problem with when I discards some
particles.
The debug file I am getting now with the same model and same process but
this time taking j and p without b(b~)- quark.
The following may be useful :
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.3.3 2015-10-25 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* and *
#* http://
#* *
#******
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#******
launch auto
Traceback (most recent call last):
File
"/home/
line 908, in onecmd
return self.onecmd_
File
"/home/
line 897, in onecmd_orig
return func(arg, **opt)
File
"/home/
line 1210, in do_launch
evt_file = self.run(mode, options)
File
"/home/
line 1365, in run
jobs_
File
"/home/
line 1654, in collect_the_results
jobs_to_
File
"/home/
line 1771, in update_jobs_to_run
elif ( req_acc > 0 and err/tot > req_acc*1.2 ) or step == 0:
ZeroDivisionError: float division by zero
Value of current Options:
notificat
/home/shibanand
ignore_
loop_
cluster_
/home/shibanand
/home/shibanand
/home/shibanand
exrootana
/home/shibanand
/home/shibanand
automatic_
/home/shibanand
/home/shibanand
output_
/home/shibanand
complex_
max_
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.3.3 2015-10-25 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set group_subprocesses Auto
set ignore_
set loop_optimized_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_
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~
import model loop_qcd_qed_sm
define j = g u c d s u~ c~ d~ s~
define p = g u c d s u~ c~ d~ s~
generate z > b b~ / ghg ghg~ gha gha~ ghz ghwp ghwm g0 g+ ghz~ ghwp~ \
ghwm~ g- [QED]
output z2bb_no5flavour
#######
## 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 LOOP
#######
Block loop
666 1.000000e+03 # MU_R
#######
## INFORMATION FOR MASS
#######
Block mass
6 1.743000e+02 # MT
23 9.118800e+01 # MZ
24 8.041900e+01 # MW
25 1.200000e+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.000000 # d : 0.0
2 0.000000 # u : 0.0
3 0.000000 # s : 0.0
4 0.000000 # c : 0.0
5 0.000000 # b : 0.0
11 0.000000 # e- : 0.0
12 0.000000 # ve : 0.0
13 0.000000 # mu- : 0.0
14 0.000000 # vm : 0.0
15 0.000000 # ta- : 0.0
16 0.000000 # vt : 0.0
21 0.000000 # g : 0.0
22 0.000000 # a : 0.0
9000002 91.188000 # ghz : MZ
9000003 80.419000 # ghwp : MW
9000004 80.419000 # ghwm : MW
250 91.188000 # g0 : MZ
251 80.419000 # g+ : MW
#######
## INFORMATION FOR SMINPUTS
#######
Block sminputs
1 1.325070e+02 # aEWM1
3 1.180000e-01 # aS
#######
## INFORMATION FOR YUKAWA
#######
Block yukawa
6 1.743000e+02 # ymt
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.354080e+00 # WT
DECAY 23 2.428230e+00 # WZ
DECAY 24 2.028440e+00 # WW
DECAY 25 3.514960e-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.000000 # d : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
DECAY 250 2.428230 # g0 : WZ
DECAY 251 2.028440 # g+ : WW
#======
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#======
Block QNUMBERS 250 # g0
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 251 # g+
1 3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# Lines starting with a hash (#) 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
#******
# Number of LHE events (and their normalization) and the required *
# (relative) accuracy on the Xsec. *
# These values are ignored for fixed order runs *
#******
10000 = nevents ! Number of unweighted events requested
-1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
-1 = nevt_job! Max number of events per job in event generation.
! (-1= no split).
#******
# Normalize the weights of LHE events such that they sum or average to *
# the total cross section *
#******
average = event_norm ! average or sum
#******
# Number of points per itegration channel (ignored for aMC@NLO runs) *
#******
0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
# These numbers are ignored except if req_acc_FO is equal to -1
5000 = npoints_FO_grid ! number of points to setup grids
4 = niters_FO_grid ! number of iter. to setup grids
10000 = npoints_FO ! number of points to compute Xsec
6 = niters_FO ! number of iter. to compute Xsec
#******
# Random number seed *
#******
0 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
#******
0 = lpp1 ! beam 1 type (0 = no PDF)
0 = lpp2 ! beam 2 type (0 = no PDF)
6500.0 = ebeam1 ! beam 1 energy in GeV
6500.0 = ebeam2 ! beam 2 energy in GeV
#******
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#******
nn23nlo = pdlabel ! PDF set
244600 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#******
# Include the NLO Monte Carlo subtr. terms for the following parton *
# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
# WARNING: PYTHIA6PT works only for processes without FSR!!!! *
#******
HERWIG6 = parton_shower
1.0 = shower_scale_factor ! multiply default shower starting
#******
# Renormalization and factorization scales *
# (Default functional form for the non-fixed scales is the sum of *
# the transverse masses of all final state particles and partons. This *
# can be changed in SubProcesses/
#******
False = fixed_ren_scale ! if .true. use fixed ren scale
False = fixed_fac_scale ! if .true. use fixed fac scale
91.118 = muR_ref_fixed ! fixed ren reference scale
91.118 = muF1_ref_fixed ! fixed fact reference scale for pdf1
91.118 = muF2_ref_fixed ! fixed fact reference scale for pdf2
-1 = dynamical_
choices
#******
# Renormalization and factorization scales (advanced and NLO options) *
#******
False = fixed_QES_scale ! if .true. use fixed Ellis-Sexton scale
91.118 = QES_ref_fixed ! fixed Ellis-Sexton reference scale
1.0 = muR_over_ref ! ratio of current muR over reference muR
1.0 = muF1_over_ref ! ratio of current muF1 over reference muF1
1.0 = muF2_over_ref ! ratio of current muF2 over reference muF2
1.0 = QES_over_ref ! ratio of current QES over reference QES
#******
# Reweight flags to get scale dependence and PDF uncertainty *
# For scale dependence: factor rw_scale_up/down around central scale *
# For PDF uncertainty: use LHAPDF with supported set *
#******
True = reweight_scale ! reweight to get scale dependence
0.5 = rw_Rscale_down ! lower bound for ren scale variations
2.0 = rw_Rscale_up ! upper bound for ren scale variations
0.5 = rw_Fscale_down ! lower bound for fact scale variations
2.0 = rw_Fscale_up ! upper bound for fact scale variations
False = reweight_PDF ! reweight to get PDF uncertainty
244601 = PDF_set_min ! First of the error PDF sets
244700 = PDF_set_max ! Last of the error PDF sets
#******
# ickkw parameter: *
# 0: No merging *
# 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
# level. After showering an MLM-type merging should be applied as *
# well. See http://
# 4: UNLOPS merging (with pythia8 only). No interface from within *
# MG5_aMC available, but available in Pythia8. *
# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
#******
0 = ickkw
#******
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15.0 = bwcutoff
#******
# Cuts on the jets *
# Jet clustering is performed by FastJet.
# When matching to a parton shower, these generation cuts should be *
# considerably softer than the analysis cuts. *
# (more specific cuts can be specified in SubProcesses/
#******
1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
0.7 = jetradius ! The radius parameter for the jet algorithm
0.0 = ptj ! Min jet transverse momentum
-1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
#******
# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
# (more specific gen cuts can be specified in SubProcesses/
#******
0.0 = ptl ! Min lepton transverse momentum
-1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
0.0 = drll ! Min distance between opposite sign lepton pairs
0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
0.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
#******
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
#******
0.0 = ptgmin ! Min photon transverse momentum
-1.0 = etagamma ! Max photon abs(pseudo-rap)
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)
#******
# Maximal PDG code for quark to be considered a jet when applying cuts.*
# At least all massless quarks of the model should be included here. *
#******
4 = maxjetflavor
#******
# For aMCfast+APPLGRID use in PDF fitting (http://
#******
0 = iappl ! aMCfast switch (0=OFF, 1=prepare APPLgrids, 2=fill grids)
#******
On Tue, May 24, 2016 at 3:27 PM, Olivier Mattelaer <
<email address hidden>> wrote:
> Your question #294100 on MadGraph5_aMC@NLO changed:
> https:/
>
> Olivier Mattelaer proposed the following answer:
> Just for completeness,
> I have fully forbidden the use of [QED] in the development version (with
> appropriate error message)
> such that the status will be clear for the user.
>
> Cheers,
>
> Olivier
>
> > On May 24, 2016, at 10:32, Rikkert Frederix <
> <email address hidden>> wrote:
> >
> > Question #294100 on MadGraph5_aMC@NLO changed:
> > https:/
> >
> > Status: Open => Answered
> >
> > Rikkert Frederix proposed the following answer:
> > Dear Subhasish,
> >
> > Ghosts are not BSM particles, but rather part of the Standard Model in
> > the Feynman gauge. Removing them destroys gauge invariance. Moreover,
> > the model you loaded is a model using a 5 flavour scheme. Hence, the
> > bottom quark is indistinguishable from any of the other light quarks,
> > which makes your process definition ill-defined. Also, you have to be
> > very careful in setting the widths of the particles to keep gauge
> > invariance, in particular in relation with the incoming Z-boson.
> >
> > However, note that NLO electroweak corrections are not supported by any
> > of the official public versions of MadGraph5_aMC@NLO. They are only
> > available in the development versions of the code, which are available
> > on launchpad website, but they require quite some expertise to use.
> >
> > best,
> > Rikkert
> >
> > --
> > 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:/
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https:/
>
> You received this question notification because you asked the question.
>
Revision history for this message
![]() |
#6 |
Dear Subhasish,
I repeat:
"Note that NLO electroweak corrections are not supported by any of the official public versions of MadGraph5_aMC@NLO. They are only available in the development versions of the code, which are available on launchpad website, but they require quite some expertise to use."
You can keep trying with the official release, but you want get a sensible result. No matter what you try...
Best,
Rikkert
Revision history for this message
![]() |
#7 |
Dear Rikkert,
Thanks you so much for clarification.
Thanks,
Subhasish
On Thu, May 26, 2016 at 12:17 AM, Rikkert Frederix <
<email address hidden>> wrote:
> Your question #294100 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Open => Answered
>
> Rikkert Frederix proposed the following answer:
> Dear Subhasish,
>
> I repeat:
>
> "Note that NLO electroweak corrections are not supported by any of the
> official public versions of MadGraph5_aMC@NLO. They are only available
> in the development versions of the code, which are available on
> launchpad website, but they require quite some expertise to use."
>
> You can keep trying with the official release, but you want get a
> sensible result. No matter what you try...
>
> Best,
> Rikkert
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https:/
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https:/
>
> You received this question notification because you asked the question.
>
Can you help with this problem?
Provide an answer of your own, or ask SUBHASISH BEHERA for more information if necessary.