aMCatNLOError : Some tests failed, run cannot continue.

Asked by jeremy andrea on 2014-04-09

Dear MG/aMC@NLO experts,

I hope this is the right place to ask my question.

I'm willing to calculate NLO cross-section of the tllq SM process using aMC@NLO. I'm putting my proc_card at the end of this message. But at the event generation step (after typing "launch", updating configuration and starting the generation) I'm getting an error that I'm not able to understand. The last few lines dumped on the screen before the crash are :

.
.
.
Collinear tests done for (Born) config 1
 Failures: 1
Collinear test 5 PASSED. Fraction of failures: 0.01

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 6
 FKS partons are: i= 6 j= 2
 with PDGs: i= -5 j= 21

 Enter graph number (iconfig), '0' loops over all graphs
 Fatal error in set_tau_min

Error detected in "launch auto "
write debug file /grid_mnt/opt__sbg__data__safe1/cms/jandrea/aMCatNLO_211/MG5_aMC_v2_0_2/PROCNLO_loop_sm_1/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
aMCatNLOError : Some tests failed, run cannot continue.
 Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.

I copy/pasted the log at the end of this email. I put explicitly the width of the top mass to 0 in the param_card, but it didn't change anything. This is the first time I'm using aMC@NLO, so most probably I did something wrong, and I'm not able to find where my mistake is. Help would be very much appreciated.

Thank you in advance,
Jeremy

-------------------------------------------------------------
The proc_card I used is the following :

import model sm
define p g u c s d b u~ c~ s~ d~ b~
define j = p
define l+ e+ mu+ ta+
define l- e- mu- ta-
define vl ve vm vt
define vl~ ve~ vm~ vt~
generate p p > t l- l+ j / h [QCD] @0
add process p p > t~ l- l+ j / h [QCD] @1

-------------------------------------------------------------
content of the debug/log file :

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.0.2 2014-02-07 *
#* *
#* 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 "/grid_mnt/opt__sbg__data__safe1/cms/jandrea/aMCatNLO_211/MG5_aMC_v2_0_2/madgraph/interface/extended_cmd.py", line 841, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/grid_mnt/opt__sbg__data__safe1/cms/jandrea/aMCatNLO_211/MG5_aMC_v2_0_2/madgraph/interface/extended_cmd.py", line 834, in onecmd_orig
    return func(arg, **opt)
  File "/grid_mnt/opt__sbg__data__safe1/cms/jandrea/aMCatNLO_211/MG5_aMC_v2_0_2/madgraph/interface/amcatnlo_run_interface.py", line 1184, in do_launch
    self.compile(mode, options)
  File "/grid_mnt/opt__sbg__data__safe1/cms/jandrea/aMCatNLO_211/MG5_aMC_v2_0_2/madgraph/interface/amcatnlo_run_interface.py", line 3192, in compile
    self.check_tests(test, this_dir)
  File "/grid_mnt/opt__sbg__data__safe1/cms/jandrea/aMCatNLO_211/MG5_aMC_v2_0_2/madgraph/interface/amcatnlo_run_interface.py", line 3204, in check_tests
    return self.parse_test_mx_log(pjoin(dir, '%s.log' % test))
  File "/grid_mnt/opt__sbg__data__safe1/cms/jandrea/aMCatNLO_211/MG5_aMC_v2_0_2/madgraph/interface/amcatnlo_run_interface.py", line 3216, in parse_test_mx_log
    ' set to 0 in the param_card.dat.')
aMCatNLOError: Some tests failed, run cannot continue.
Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.
Value of current Options:
              text_editor : None
              web_browser : None
        cluster_temp_path : None
       group_subprocesses : False
                  timeout : 60
            cluster_queue : None
         madanalysis_path : None
                   lhapdf : lhapdf-config
ignore_six_quark_processes : False
                 mg5_path : /grid_mnt/opt__sbg__data__safe1/cms/jandrea/aMCatNLO_211/MG5_aMC_v2_0_2
    loop_optimized_output : True
    cluster_status_update : (600, 30)
               hepmc_path : None
             pythia8_path : None
                hwpp_path : None
   automatic_html_opening : False
       cluster_retry_wait : 300
             stdout_level : 20
          pythia-pgs_path : None
                  td_path : None
             delphes_path : None
              thepeg_path : None
             cluster_type : condor
      exrootanalysis_path : None
         fortran_compiler : None
                      OLP : MadLoop
              auto_update : 7
         cluster_nb_retry : 1
               eps_viewer : None
             syscalc_path : None
                  fastjet : None
                    gauge : unitary
      complex_mass_scheme : False
                  nb_core : 16
                 run_mode : 2

Question information

Language:
English Edit question
Status:
Open
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Last query:
2014-04-10
Last reply:
2014-04-10
marco zaro (marco-zaro) said : #1

Dear Jeremy,
can you try doing exactly what you did, but importing the model with massless b quark? That may be the reason of the failing tests.

To import the model you should type

 > import model loop_sm-no_b_mass
 and then the commands you used to generate the process
> define p = ...
> define j= p
> generate ...
> output

Please let me know if this solves your problem.

Best!

Marco

jeremy andrea (jandrea) said : #2

Dear MArco,

thank you very much, it works indeed with the bmassless model.

However, I'm not sure to know if it is consistent with pdfs, as the
b-quarks can be in the initial states and the PDFs consider b-quarks as
massive.

But at least I can get a number. Thank you very much.

Cheers,
Jeremy

Le 10/04/2014 09:16, marco zaro a écrit :
> Your question #246750 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/246750
>
> Status: Open => Answered
>
> marco zaro proposed the following answer:
> Dear Jeremy,
> can you try doing exactly what you did, but importing the model with massless b quark? That may be the reason of the failing tests.
>
> To import the model you should type
>
> > import model loop_sm-no_b_mass
> and then the commands you used to generate the process
>> define p = ...
>> define j= p
>> generate ...
>> output
> Please let me know if this solves your problem.
>
> Best!
>
> Marco
>

marco zaro (marco-zaro) said : #3

Dear Jeremy,
which PDFs are you using?
Most PDF sets are in the five flavour scheme, so they effectively include
the b quark in the proton. With these pdfs (basically all sets on lhapdf
which do not say '4flavour' in their name) it is consistent to use a
massless b also in the hard matrix element.
Let me know if you need more.
Cheers,

Marco

Marco Zaro

On Thu, Apr 10, 2014 at 2:06 PM, jeremy andrea <
<email address hidden>> wrote:

> Question #246750 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/246750
>
> Status: Answered => Open
>
> jeremy andrea is still having a problem:
> Dear MArco,
>
> thank you very much, it works indeed with the bmassless model.
>
> However, I'm not sure to know if it is consistent with pdfs, as the
> b-quarks can be in the initial states and the PDFs consider b-quarks as
> massive.
>
> But at least I can get a number. Thank you very much.
>
> Cheers,
> Jeremy
>
>
> Le 10/04/2014 09:16, marco zaro a écrit :
> > Your question #246750 on MadGraph5_aMC@NLO changed:
> > https://answers.launchpad.net/mg5amcnlo/+question/246750
> >
> > Status: Open => Answered
> >
> > marco zaro proposed the following answer:
> > Dear Jeremy,
> > can you try doing exactly what you did, but importing the model with
> massless b quark? That may be the reason of the failing tests.
> >
> > To import the model you should type
> >
> > > import model loop_sm-no_b_mass
> > and then the commands you used to generate the process
> >> define p = ...
> >> define j= p
> >> generate ...
> >> output
> > Please let me know if this solves your problem.
> >
> > Best!
> >
> > Marco
> >
>
> --
> You received this question notification because you are a direct
> subscriber of the question.
>

jeremy andrea (jandrea) said : #4

Dear Marco,

thank you again for your answers.

I'm using default pdf : cteq6_m.

thank you for confirming that what I'm doing is consistent.

Thank you again for help.
Jeremy

Le 10/04/2014 14:16, marco zaro a écrit :
> Your question #246750 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/246750
>
> Status: Open => Answered
>
> marco zaro proposed the following answer:
> Dear Jeremy,
> which PDFs are you using?
> Most PDF sets are in the five flavour scheme, so they effectively include
> the b quark in the proton. With these pdfs (basically all sets on lhapdf
> which do not say '4flavour' in their name) it is consistent to use a
> massless b also in the hard matrix element.
> Let me know if you need more.
> Cheers,
>
> Marco
>
>
> Marco Zaro
>
>
> On Thu, Apr 10, 2014 at 2:06 PM, jeremy andrea <
> <email address hidden>> wrote:
>
>> Question #246750 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/246750
>>
>> Status: Answered => Open
>>
>> jeremy andrea is still having a problem:
>> Dear MArco,
>>
>> thank you very much, it works indeed with the bmassless model.
>>
>> However, I'm not sure to know if it is consistent with pdfs, as the
>> b-quarks can be in the initial states and the PDFs consider b-quarks as
>> massive.
>>
>> But at least I can get a number. Thank you very much.
>>
>> Cheers,
>> Jeremy
>>
>>
>> Le 10/04/2014 09:16, marco zaro a écrit :
>>> Your question #246750 on MadGraph5_aMC@NLO changed:
>>> https://answers.launchpad.net/mg5amcnlo/+question/246750
>>>
>>> Status: Open => Answered
>>>
>>> marco zaro proposed the following answer:
>>> Dear Jeremy,
>>> can you try doing exactly what you did, but importing the model with
>> massless b quark? That may be the reason of the failing tests.
>>> To import the model you should type
>>>
>>> > import model loop_sm-no_b_mass
>>> and then the commands you used to generate the process
>>>> define p = ...
>>>> define j= p
>>>> generate ...
>>>> output
>>> Please let me know if this solves your problem.
>>>
>>> Best!
>>>
>>> Marco
>>>
>> --
>> You received this question notification because you are a direct
>> subscriber of the question.
>>

Can you help with this problem?

Provide an answer of your own, or ask jeremy andrea for more information if necessary.

To post a message you must log in.