MadSpin and e_ijk color coupling

Dear Simon,

I have discussed with Pierre and we can propose a patch to add the support of the epsilon_ijk.
This is however not tested/validated so please use it care. If it works do not hesitate to inform us and we will then include this patch in the next version.

Cheers,

Olivier

PS The patch:

=== modified file 'MadSpin/decay.py'
--- MadSpin/decay.py 2016-01-14 13:28:22 +0000
+++ MadSpin/decay.py 2016-01-26 09:56:14 +0000
@@ -3687,6 +3687,45 @@
                             d1colup1=maxcol
                             d2colup2=maxcol
                             d2colup1=0
+ elif colord1==-3 and colord2==-3 and colormother == 3:
+ maxcol+=2
+ d1colup1=0
+ d1colup2=maxcol
+ d2colup1=0
+ d2colup2=maxcol-1
+ elif (colord1==-3 and colord2==3 and colormother == 3) or\
+ (colord1==-3 and colord2==3 and colormother == -3):
+ maxcol+=2
+ d1colup1 = 0
+ d1colup2 = maxcol
+ d2colup1 = maxcol-1
+ d2colup2 = 0
+ elif (colord1==3 and colord2==-3 and colormother == 3) or\
+ (colord1==3 and colord2==-3 and colormother == -3):
+ maxcol+=2
+ d1colup1=maxcol
+ d1colup2=0
+ d2colup1=0
+ d2colup2=maxcol-1
+ elif colord1==3 and colord2==3 and colormother == -3:
+ maxcol+=2
+ d1colup1=maxcol
+ d1colup2=0
+ d2colup1=maxcol-1
+ d2colup2=0
+ elif colord2==8 and colord1==8 and colormother ==8:
+ maxcol+=1
+ ran = random.random()
+ if ran> 0.5:
+ d1colup2=colup2
+ d1colup1=maxcol
+ d2colup2=maxcol
+ d2colup1=colup1
+ else:
+ d1colup2=maxcol
+ d1colup1=colup1
+ d2colup2=colup2
+ d2colup1=maxcol
                         else:
                             raise Exception, 'color combination not treated by MadSpin (yet). (%s,%s,%s)' \
                                 % (colord1,colord2,colormother)

 Hi,

 Thx for the patch and the fast answer.
 I will try this patch and will give you a feedback.

 Cheers, Simon.