About NLO model file loop_qcd_qed_sm

Dear Subhasish,

What it is that you are trying to compute? Your current setup does not make sense to me.

best,
Rikkert

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://answers.launchpad.net/mg5amcnlo/+question/294100
>
> 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://answers.launchpad.net/mg5amcnlo/+question/294100
>
> You received this question notification because you asked the question.
>

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

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://answers.launchpad.net/mg5amcnlo/+question/294100
>
> 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.

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://server06.fynu.ucl.ac.be/projects/madgraph *
#* and *
#* http://amcatnlo.cern.ch *
#* *
#************************************************************
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#************************************************************
launch auto
Traceback (most recent call last):
  File
"/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/madgraph/interface/extended_cmd.py",
line 908, in onecmd
    return self.onecmd_orig(line, **opt)
  File
"/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/madgraph/interface/extended_cmd.py",
line 897, in onecmd_orig
    return func(arg, **opt)
  File
"/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/madgraph/interface/amcatnlo_run_interface.py",
line 1210, in do_launch
    evt_file = self.run(mode, options)
  File
"/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/madgraph/interface/amcatnlo_run_interface.py",
line 1365, in run
    jobs_to_collect,integration_step,mode,mode_dict[mode])
  File
"/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/madgraph/interface/amcatnlo_run_interface.py",
line 1654, in collect_the_results

jobs_to_run_new=self.update_jobs_to_run(req_acc,integration_step,jobs_to_run,fixed_order)
  File
"/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/madgraph/interface/amcatnlo_run_interface.py",
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:
              text_editor : None
      notification_center : True
                    pjfry :
/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/PJFry/lib
       cluster_local_path : /cvmfs/cp3.uclouvain.be/madgraph/
       group_subprocesses : Auto
ignore_six_quark_processes : False
    loop_optimized_output : True
    cluster_status_update : (600, 30)
             pythia8_path : None
                hwpp_path : None
                    golem :
/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/golem95/lib
          pythia-pgs_path :
/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/pythia-pgs
                  td_path : None
             delphes_path :
/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/Delphes
              thepeg_path : None
             cluster_type : condor
      exrootanalysis_path :
/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/ExRootAnalysis
                      OLP : MadLoop
                 applgrid : applgrid-config
               eps_viewer : None
                  fastjet :
/home/shibananda/RESEARCH/PACKAGES_USE/FASTJET-3.1.3/bin/fastjet-config
                 run_mode : 2
              web_browser : None
   automatic_html_opening : False
        cluster_temp_path : None
             cluster_size : 100
            cluster_queue : None
             syscalc_path :
/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3/SysCalc
         madanalysis_path : None
                   lhapdf :
/home/shibananda/RESEARCH/PACKAGES_USE/LHAPDF616/bin/lhapdf-config
             stdout_level : 20
            f2py_compiler : None
               hepmc_path : None
                  amcfast : amcfast-config
       cluster_retry_wait : 300
      output_dependencies : external
         fortran_compiler : None
                  nb_core : 4
         loop_color_flows : False
              auto_update : 7
         cluster_nb_retry : 1
                 mg5_path :
/home/shibananda/RESEARCH/EVENT_GENERATORS/MG5_aMC_v2_3_3
                  timeout : 60
                    gauge : Feynman
      complex_mass_scheme : False
             cpp_compiler : None
   max_npoint_for_channel : 0
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.3.3 2015-10-25 *
#* *
#* 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 loop_optimized_output 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~
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/Antiparticle distinction (0=own anti)
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/Antiparticle distinction (0=own anti)
#***********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# 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
                        ! number of points and iter. below)
# 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=default))
#***********************************************************************
# 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
                                  ! scale by this factor
#***********************************************************************
# 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/set_scales.f) *
#***********************************************************************
 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_scale_choice ! Choose one of the preselected 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://amcatnlo.cern.ch/FxFx_merging.htm for details. *
# 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*Gamma) *
#***********************************************************************
 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/cuts.f) *
#***********************************************************************
 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/cuts.f) *
#***********************************************************************
  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://amcfast.hepforge.org)*
#***********************************************************************
 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://answers.launchpad.net/mg5amcnlo/+question/294100
>
> 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://answers.launchpad.net/mg5amcnlo/+question/294100
> >
> > 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://answers.launchpad.net/mg5amcnlo/+question/294100/+confirm?answer_id=3
>
> 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/294100
>
> You received this question notification because you asked the question.
>

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

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://answers.launchpad.net/mg5amcnlo/+question/294100
>
> 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://answers.launchpad.net/mg5amcnlo/+question/294100/+confirm?answer_id=5
>
> 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/294100
>
> You received this question notification because you asked the question.
>