Tau Decay (1 >3 non-resonant decays)

Asked by Elham E Khoda

Hello,

I have two prblems:

1. In reference to answer #249080: It seems that 1>3 decay in MadSpin was in the to-do list. I was trying generate tau decaying to pi-, pi0 ,vt using MG5_aMC_v2_2_2 keeping MadSpin on. But MadSpin was unable to generate these 1 > 3 body decays. Is it already updated to some other version? Or is there any way such that MadSpin can actually deal with this kind of phase space?

2. Till now I was working with MG5_aMC_v2_2_2. And taudecay_UFO was working properly. Today I was trying to work with MG5_aMC_v2_2_3. But I am unable to add model taudecay_UFO with sm/sm-lepton_masses. It is showing the following message:

WARNING: QED has different hierarchy use the minimal value (2, 1) => 1
Command "add model taudecay_UFO " interrupted in sub-command:
"add model taudecay_UFO" with error:
UnboundLocalError : local variable 'identify_pid' referenced before assignment

I will be highly grateful if you can answer these questions for me.

Thanking you,
Elham E Khoda

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

1. Do not know the exact status on this but never seen a version close to a beta so far.
On the other hand, the next version of the code is likely to have a partial spin correlation mode which will be able to handle 1>3 body decay.
Not sure if it will be include or not at this stage (I have to finish to implement it first).

2.
Thanks for reporting this.
Here is the patch to make it work: It will be include in our next release.

=== modified file 'models/usermod.py'
--- models/usermod.py 2015-02-03 11:27:35 +0000
+++ models/usermod.py 2015-03-09 13:09:34 +0000
@@ -832,8 +832,8 @@

         # Adding automatically identification for anti-particle if needed
         # + define identify_pid which keep tracks of the pdg_code identified
+ identify_pid = {}
         if identify_particles:
- identify_pid = {}
             for new, old in identify_particles.items():
                 new_part = next((p for p in model.all_particles if p.name==new), None)
                 old_part = next((p for p in self.particles if p.name==old), None)

Cheers,

Olivier

On 09 Mar 2015, at 15:26, Elham E Khoda <email address hidden> wrote:

> New question #263408 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/263408
>
> Hello,
>
> I have two prblems:
>
> 1. In reference to answer #249080: It seems that 1>3 decay in MadSpin was in the to-do list. I was trying generate tau decaying to pi-, pi0 ,vt using MG5_aMC_v2_2_2 keeping MadSpin on. But MadSpin was unable to generate these 1 > 3 body decays. Is it already updated to some other version? Or is there any way such that MadSpin can actually deal with this kind of phase space?
>
>
> 2. Till now I was working with MG5_aMC_v2_2_2. And taudecay_UFO was working properly. Today I was trying to work with MG5_aMC_v2_2_3. But I am unable to add model taudecay_UFO with sm/sm-lepton_masses. It is showing the following message:
>
> WARNING: QED has different hierarchy use the minimal value (2, 1) => 1
> Command "add model taudecay_UFO " interrupted in sub-command:
> "add model taudecay_UFO" with error:
> UnboundLocalError : local variable 'identify_pid' referenced before assignment
>
>
> I will be highly grateful if you can answer these questions for me.
>
>
> Thanking you,
> Elham E Khoda
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Elham E Khoda (elhamekhoda) said :
#2

Dear Oliver,

Thank you very much for your help.
I can now add models on top of sm/sm-lepton_masses/mssm after doing the modification you suggested.

Thanking you,
Elham

Revision history for this message
Elham E Khoda (elhamekhoda) said :
#3

Thanks Olivier Mattelaer, that solved my question.