ValueError Problem

Asked by teddym

Hi
When I generate some process using add-on models hgg_plugin and taudecay_UFO, I encounter following problems:

INFO: Idle: 0, Running: 0, Completed: 84 [ 10h 10m ]
INFO: Combining runs
Error when reading /home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/SubProcesses/P3_gg_hgg_h_tamtap_tam_pimpi0vl_tap_pippi0vl/G10a0/results.dat
Command "multi_run 10 run_default " interrupted with error:
ValueError : empty string for float()

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Kentarou Mawatari Edit question
Last query:
Last reply:
Revision history for this message
teddym (niepanchongsheng) said :
#1

Following is the message in the log file

multi_run 10 run_default
Traceback (most recent call last):
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/extended_cmd.py", line 1009, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/extended_cmd.py", line 964, in onecmd_orig
    return func(arg, **opt)
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/madevent_interface.py", line 2357, in do_multi_run
    self.exec_cmd('generate_events %s_%s -f' % (main_name, i), postcmd=False)
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/extended_cmd.py", line 1035, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/extended_cmd.py", line 964, in onecmd_orig
    return func(arg, **opt)
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/madevent_interface.py", line 1994, in do_generate_events
    self.exec_cmd('refine %s' % nb_event, postcmd=False)
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/extended_cmd.py", line 1035, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/extended_cmd.py", line 964, in onecmd_orig
    return func(arg, **opt)
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/madevent_interface.py", line 2851, in do_refine
    combine_runs.CombineRuns(self.me_dir)
 File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/combine_runs.py", line 78, in __init__
    self.sum_multichannel(channel)
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/combine_runs.py", line 101, in sum_multichannel
    filepath=pjoin(path, 'results.dat'))
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/sum_html.py", line 409, in add_results
    oneresult.read_results(filepath)
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/sum_html.py", line 300, in read_results
    data = [secure_float(d) for d in line.split()]
  File "/home/satya/Simulations/Higgs_tau_tau/MG5_aMC_v2_4_2/HiggsTauTau/VBF/pp_hjj_matched_4pi2v/bin/internal/sum_html.py", line 297, in secure_float
    return float(m.group(1))*10**(float(m.group(2)))
ValueError: empty string for float()

Revision history for this message
teddym (niepanchongsheng) said :
#2

How should I solve this problem?

Thanks

Ted

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#3

What is in the log file of
/VBF/pp_hjj_matched_4pi2v/SubProcesses/P3_gg_hgg_h_tamtap_tam_pimpi0vl_tap_pippi0vl/G10a0/

Cheers,

Olivier

Revision history for this message
teddym (niepanchongsheng) said :
#4

Hi
There are something like:

"""
Relative number of events:
  0.4997E+00 0.0000E+00
  0.5003E+00 0.0000E+00
 Events:
      639628 0
      640372 0
 Error done early 30556 156883
 Error done early 30557 156883
 Error done early 30558 156883
 Error done early 30559 156883
 Error done early 30560 156883
 Error done early 30561 156883
.....
......

 Error done early 156721 156883
At line 164 of file rw_events.f (unit = 27, file = '/tmp/gfortrantmpUaogqr')
Fortran runtime error: Bad address
 Error done early 156722 156883
......
.....
 Error done early 156881 156883
 Error done early 156882 156883
 Error done early 156883 156883
 Found 156883 events.
 Wrote 2 events.
 Actual xsec 1.89065836577810694E-003
 Correct abs xsec 1.89065836577810694E-003
 Event xsec 1.89065836577810694E-003
 Events wgts > 1: 0
 % Cross section > 1: 0.0000000000000000 0.0000000000000000

ls status:
events.lhe
ftn25
ftn26
input_sg.txt
log.txt
moffset.dat
results.dat
run_default_0_app.log

"""

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#5

Hi,

Looks like you have multiple problem.
Some linked to your file system (not able to handle the file size that we wrote) some of them related with the running of the code.

The most important problem seems to be:

> Found 156883 events.
> Wrote 2 events.

So the unweighting efficiency (the event generation) has extremely bad efficiency.
This is typically related to a divergency of your matrix element.
i.e. that you have forget to put a cut to remove some of the LO divergency of the matrix-element.

All the other problem are likely to be solve automatically when you would have solve that one.

Cheers,

Olivier

> On Aug 1, 2016, at 17:42, teddym <email address hidden> wrote:
>
> Question #310949 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/310949
>
> Status: Needs information => Open
>
> teddym gave more information on the question:
> Hi
> There are something like:
>
> """
> Relative number of events:
> 0.4997E+00 0.0000E+00
> 0.5003E+00 0.0000E+00
> Events:
> 639628 0
> 640372 0
> Error done early 30556 156883
> Error done early 30557 156883
> Error done early 30558 156883
> Error done early 30559 156883
> Error done early 30560 156883
> Error done early 30561 156883
> .....
> ......
>
> Error done early 156721 156883
> At line 164 of file rw_events.f (unit = 27, file = '/tmp/gfortrantmpUaogqr')
> Fortran runtime error: Bad address
> Error done early 156722 156883
> ......
> .....
> Error done early 156881 156883
> Error done early 156882 156883
> Error done early 156883 156883
> Found 156883 events.
> Wrote 2 events.
> Actual xsec 1.89065836577810694E-003
> Correct abs xsec 1.89065836577810694E-003
> Event xsec 1.89065836577810694E-003
> Events wgts > 1: 0
> % Cross section > 1: 0.0000000000000000 0.0000000000000000
>
> ls status:
> events.lhe
> ftn25
> ftn26
> input_sg.txt
> log.txt
> moffset.dat
> results.dat
> run_default_0_app.log
>
> """
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
teddym (niepanchongsheng) said :
#6

Hi:
   Now the situation is that:

When I generate just p p > h j j, h > ta- ta+, MG works well. But when I add model taudecay_UFO to decay taus ( ta > pi pi vt)based on above process, then

    1. It will take very long time to improve the cross section. ( nearly 20h)
    2. after this and just before combine_events, it will throw "ValueError : empty string for float()" and exit.

Also, I use FeynRule and write two additional effective vertex ( tau-rho-v and rho-pi-pi ) into standard model, and generate UFO model files. And using this new generated model files, there is no such problem. MG works also very well.

So I doubt that there maybe some issues in taudecay_UFO or maybe there is something I need to pay attention when using taudecay_UFO?

Revision history for this message
Kentarou Mawatari (kentarou-mawatari) said :
#7

Do you start from e.g. 'sm-lepton_masses'?
The tau width may cause the problem, I guess.

Kentarou

Revision history for this message
teddym (niepanchongsheng) said :
#8

Dear Kentarou

      Yes, I was starting from sm-lepton_masses.

      Then what should I use for tau-width? right now I use 2.27E-12. Thanks!

Best!

Revision history for this message
Kentarou Mawatari (kentarou-mawatari) said :
#9

With the latest MG5aMC (v.2.4.3), I got 10,000 events without any problem on my laptop for the VBF process:

import model sm-lepton_masses
add model taudecay_UFO
generate p p > j j h $$ w+ w- z, (h > ta- ta+, ta- > vt pi- pi0, ta+ > vt~ pi+ pi0)

Do you want to generate gluon-fusion H+2jet events, instead of VBF?

Kentarou

Can you help with this problem?

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

To post a message you must log in.