Madspin crash for ttbar decays

Hi,

Could you try with 2.1.2_beta2?

I think that version should solve your problem.

Cheers,

Olivier
On Jul 2, 2014, at 10:26 AM, Inês Ochoa <email address hidden> wrote:

> New question #251088 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/251088
>
> Hi,
>
> I'm trying to use MadSpin to decay p p > t t~ events, but the program is crashing. (I pasted the content of the MS_debug file below). I've changed a few values on the param_card but it was done according to the calculator in http://madgraph.hep.uiuc.edu/Calculators/sm/sm_ufo_calc.html. It's not clear to me from the error message what could be the problem...
>
> Thanks!
> Inês
>
> ----------------------------------------------------------------------------------------------------------------------------------------------------------
> set seed 486380505
> import tt_8TeV/Events/run_01/events.lhe
> decay t > w+ b, w+ > mu+ vm
> decay t~ > w- b~, w- > j j
> launch
> Traceback (most recent call last):
> File "/home/inesochoa/MG5_aMC_v2_1_0/madgraph/interface/extended_cmd.py", line 872, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/inesochoa/MG5_aMC_v2_1_0/madgraph/interface/extended_cmd.py", line 865, in onecmd_orig
> return func(arg, **opt)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/interface_madspin.py", line 459, in do_launch
> self.options)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1966, in __init__
> self.generate_all_matrix_element()
> File "/home/inesochoa/MG5_aMC_v2_1_0/madgraph/various/misc.py", line 88, in f_with_no_logger
> out = f(self, *args, **opt)
> File "/home/inesochoa/MG5_aMC_v2_1_0/tests/unit_tests/various/test_aloha.py", line 50, in deco_f_set
> out = f(*args, **opt)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 2503, in generate_all_matrix_element
> self.get_branching_ratio()
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 2740, in get_branching_ratio
> width.extract_br(self.list_branches, self.mgcmd)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1302, in extract_br
> self.launch_width_evaluation(to_decay, mgcmd)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1567, in launch_width_evaluation
> self.extract_br_from_card(pjoin(self.path_me, 'param_card.dat'))
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1652, in extract_br_from_card
> self.extract_br_for_antiparticle()
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1527, in extract_br_for_antiparticle
> self.width_value[anti_res]=self.width_value[res]
> KeyError: 'w-'
> Value of current Options:
> zeromass_for_max_weight : 5
> seed : 486380505
> BW_effect : 1
> max_weight_ps_point : 400
> nb_sigma : 5.64021693405
> Nevents_for_max_weigth : 139
> BW_cut : 15.0
> max_weight : -1
> curr_dir : /home/inesochoa/MG5_aMC_v2_1_0/tt_8TeV/Events/run_01/../..
> ms_dir : None
> max_running_process : 100
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

Thanks for the suggestion! I tried repeating the same generation with the same cards using version 2.1.2_beta2, but I get exactly the same crash when running MadSpin...

Cheers
Inês

Hi Ines,

I have try to reproduce your error but I didn’t find a way to.

Could you give me all the information that I need to reproduce your problem.
The critical point is probably your card.

You can send me attachment to my email address: <email address hidden>

Cheers,

Olivier

On Jul 2, 2014, at 1:56 PM, Inês Ochoa <email address hidden> wrote:

> Question #251088 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/251088
>
> Inês Ochoa posted a new comment:
> Hi Olivier,
>
> Thanks for the suggestion! I tried repeating the same generation with
> the same cards using version 2.1.2_beta2, but I get exactly the same
> crash when running MadSpin...
>
> Cheers
> Inês
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Ines,

Thanks for the file. The good news is that I’m able to reproduce your problem now.
So i’m starting to investigate the problem.

Cheers,

Olivier
On Jul 2, 2014, at 2:41 PM, olivier mattelaer <email address hidden> wrote:

> Hi Ines,
>
> I have try to reproduce your error but I didn’t find a way to.
>
> Could you give me all the information that I need to reproduce your problem.
> The critical point is probably your card.
>
> You can send me attachment to my email address: <email address hidden>
>
> Cheers,
>
> Olivier
>
>
> On Jul 2, 2014, at 1:56 PM, Inês Ochoa <email address hidden> wrote:
>
>> Question #251088 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/251088
>>
>> Inês Ochoa posted a new comment:
>> Hi Olivier,
>>
>> Thanks for the suggestion! I tried repeating the same generation with
>> the same cards using version 2.1.2_beta2, but I get exactly the same
>> crash when running MadSpin...
>>
>> Cheers
>> Inês
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadGraph5_aMC@NLO.
>

Hi Ines,

The reason for the bug is that some of the partial width are define (like the one for the W) but not all of them (like the one for the top). This make the code to crash.
So if you add the line of the branching ratio for the top, it works nicely.
I have made a patch to handle this situation in the next version of the code which will be release this week.

Cheers and thanks,

Olivier

Hi Olivier,

I didn't have this line before though - can you send me an example? I tried adding

# BR NDA ID1 ID2 ...
   9.98e-01 2 5 24

but I get an error as well.

Thank you for your patience!

Inês

Hi,

I add the line
 1.0 2 5 24 # 1.5
or
 1.0 2 5 -24 # 1.5
(not sure which one)
and it was working

Cheers,

Olivier

On Jul 3, 2014, at 12:11 PM, Inês Ochoa <email address hidden> wrote:

> Question #251088 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/251088
>
> Inês Ochoa posted a new comment:
> Hi Olivier,
>
> I didn't have this line before though - can you send me an example? I
> tried adding
>
> # BR NDA ID1 ID2 ...
> 9.98e-01 2 5 24
>
> but I get an error as well.
>
> Thank you for your patience!
>
> Inês
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

I'm afraid now I only get the crash ZeroDivisionError : float division by zero.
I've pasted the MS_debug file below.

Thanks!

----------------------------------------------------------------------------------------------------------------------------------------------------------

set seed 25905577
import tt_8TeV/Events/run_01/events.lhe.gz
decay t > w+ b, w+ > mu+ vm
decay t~ > w- b~, w- > mu- vm~
launch
launch
Traceback (most recent call last):
  File "/home/inesochoa/MG5_aMC_v2_1_0/madgraph/interface/extended_cmd.py", line 872, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/inesochoa/MG5_aMC_v2_1_0/madgraph/interface/extended_cmd.py", line 865, in onecmd_orig
    return func(arg, **opt)
  File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/interface_madspin.py", line 459, in do_launch
    self.options)
  File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1966, in __init__
    self.generate_all_matrix_element()
  File "/home/inesochoa/MG5_aMC_v2_1_0/madgraph/various/misc.py", line 88, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/home/inesochoa/MG5_aMC_v2_1_0/tests/unit_tests/various/test_aloha.py", line 50, in deco_f_set
    out = f(*args, **opt)
  File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 2503, in generate_all_matrix_element
    self.get_branching_ratio()
  File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 2740, in get_branching_ratio
    width.extract_br(self.list_branches, self.mgcmd)
  File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1287, in extract_br
    self.extract_br_from_banner(self.banner)
  File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1612, in extract_br_from_banner
    return self.extract_br_from_card(param_card)
  File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1634, in extract_br_from_card
    relative_diff=abs(recalculated_width-width_in_the_banner)/recalculated_width
ZeroDivisionError: float division by zero
Value of current Options:
  zeromass_for_max_weight : 5
                     seed : 25905577
                BW_effect : 1
      max_weight_ps_point : 400
                 nb_sigma : 4.5
   Nevents_for_max_weigth : 75
                   BW_cut : 15.0
               max_weight : -1
                 curr_dir : /home/inesochoa/MG5_aMC_v2_1_0/tt_8TeV/Events/run_01/../..
                   ms_dir : None
      max_running_process : 100

Ok,

Just wait the next version, it should work out of the box in that one.
If you have bazaar installed on your computer, you can already install it via the command
bzr branch lp:~maddevelopers/mg5amc/2.1.2
otherwise, this version will be available via classic download in less than 24h.

Cheers,

Olivier

On Jul 3, 2014, at 2:26 PM, Inês Ochoa <email address hidden> wrote:

> Question #251088 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/251088
>
> Inês Ochoa posted a new comment:
> Hi Olivier,
>
> I'm afraid now I only get the crash ZeroDivisionError : float division by zero.
> I've pasted the MS_debug file below.
>
> Thanks!
>
> ----------------------------------------------------------------------------------------------------------------------------------------------------------
>
> set seed 25905577
> import tt_8TeV/Events/run_01/events.lhe.gz
> decay t > w+ b, w+ > mu+ vm
> decay t~ > w- b~, w- > mu- vm~
> launch
> launch
> Traceback (most recent call last):
> File "/home/inesochoa/MG5_aMC_v2_1_0/madgraph/interface/extended_cmd.py", line 872, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/inesochoa/MG5_aMC_v2_1_0/madgraph/interface/extended_cmd.py", line 865, in onecmd_orig
> return func(arg, **opt)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/interface_madspin.py", line 459, in do_launch
> self.options)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1966, in __init__
> self.generate_all_matrix_element()
> File "/home/inesochoa/MG5_aMC_v2_1_0/madgraph/various/misc.py", line 88, in f_with_no_logger
> out = f(self, *args, **opt)
> File "/home/inesochoa/MG5_aMC_v2_1_0/tests/unit_tests/various/test_aloha.py", line 50, in deco_f_set
> out = f(*args, **opt)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 2503, in generate_all_matrix_element
> self.get_branching_ratio()
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 2740, in get_branching_ratio
> width.extract_br(self.list_branches, self.mgcmd)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1287, in extract_br
> self.extract_br_from_banner(self.banner)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1612, in extract_br_from_banner
> return self.extract_br_from_card(param_card)
> File "/home/inesochoa/MG5_aMC_v2_1_0/MadSpin/decay.py", line 1634, in extract_br_from_card
> relative_diff=abs(recalculated_width-width_in_the_banner)/recalculated_width
> ZeroDivisionError: float division by zero
> Value of current Options:
> zeromass_for_max_weight : 5
> seed : 25905577
> BW_effect : 1
> max_weight_ps_point : 400
> nb_sigma : 4.5
> Nevents_for_max_weigth : 75
> BW_cut : 15.0
> max_weight : -1
> curr_dir : /home/inesochoa/MG5_aMC_v2_1_0/tt_8TeV/Events/run_01/../..
> ms_dir : None
> max_running_process : 100
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

That's fine, less than 24h works for me! Thanks a lot, I'll try again tomorrow with the new version and will let you know.

Cheers
Inês

Thanks Olivier Mattelaer, that solved my question.