New physics, decay chain error

Hi,

The problem is that the model is not valid. The lorentz structure associate to your interaction is not valid:
> Identity(2,1)
is not valid for a FFV, which expect to have a lorentz index in the interaction.

Cheers,

Olivier

On Jul 22, 2014, at 8:31 AM, Jackson Clarke <email address hidden> wrote:

> New question #251994 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/251994
>
> Hi there,
>
> I am trying to implement a simple new physics model in MadGraph5_aMC@NLO v2.1.2.
>
> I first import my own UFO model file, which is the Higgs effective model plus some hidden fermions (hf1, hf2) and a hidden photon (ha).
>
> The process I would like to generate is
> p p > h > hf1 hf1~, hf1 > hf2 ha, hf1~ > hf2~ ha
> where I would like everything on-shell, trying to mimic an ATLAS study 1210.0435.
>
> Everything is fine until I go to generate events. I get the following error:
>
>
> Working on SubProcesses
> P0_gg_hf1hf1x_hf1_hf2ha_hf1x_hf2xha
> INFO: Idle: 0, Running: 1, Completed: 0 [ current time: 16h17 ]
> At line 33 of file FFV2_1.f
> Fortran runtime error: Index '7' of dimension 1 of array 'f1' above upper bound of 6
> INFO: Idle: 0, Running: 0, Completed: 1 [ 5s ]
> INFO: End survey
> Survey return zero cross section...
>
>
> Here is the fortran file FFV2_1.f:
>
>
> C This File is Automatically generated by ALOHA
> C The process calculated in this file is:
> C Identity(2,1)
> C
> SUBROUTINE FFV2_1(F2, V3, COUP, M1, W1,F1)
> IMPLICIT NONE
> COMPLEX*16 CI
> PARAMETER (CI=(0D0,1D0))
> COMPLEX*16 F2(*)
> COMPLEX*16 V3(*)
> REAL*8 P1(0:3)
> REAL*8 M1
> REAL*8 W1
> COMPLEX*16 F1(6)
> COMPLEX*16 DENOM
> COMPLEX*16 COUP
> F1(1) = +F2(1)+V3(1)
> F1(2) = +F2(2)+V3(2)
> P1(0) = -DBLE(F1(1))
> P1(1) = -DBLE(F1(2))
> P1(2) = -DIMAG(F1(2))
> P1(3) = -DIMAG(F1(1))
> DENOM = COUP/(P1(0)**2-P1(1)**2-P1(2)**2-P1(3)**2 - M1 * (M1
> $ -CI* W1))
> F1(3)= DENOM*-CI * V3(3)*(F2(5)*(P1(0)+P1(3))+(F2(6)*(P1(1)
> $ +CI*(P1(2)))-F2(3)*M1))
> F1(4)= DENOM*CI * V3(3)*(F2(5)*(+CI*(P1(2))-P1(1))+(F2(6)*(P1(3)
> $ -P1(0))+F2(4)*M1))
> F1(5)= DENOM*CI * V3(3)*(F2(3)*(P1(3)-P1(0))+(F2(4)*(P1(1)
> $ +CI*(P1(2)))+F2(5)*M1))
> F1(6)= DENOM*-CI * V3(3)*(F2(3)*(+CI*(P1(2))-P1(1))+(F2(4)
> $ *(P1(0)+P1(3))-F2(6)*M1))
> F1(7)= DENOM*-CI * V3(4)*(F2(5)*(P1(0)+P1(3))+(F2(6)*(P1(1)
> $ +CI*(P1(2)))-F2(3)*M1))
> F1(8)= DENOM*CI * V3(4)*(F2(5)*(+CI*(P1(2))-P1(1))+(F2(6)*(P1(3)
> $ -P1(0))+F2(4)*M1))
> F1(9)= DENOM*CI * V3(4)*(F2(3)*(P1(3)-P1(0))+(F2(4)*(P1(1)
> $ +CI*(P1(2)))+F2(5)*M1))
> F1(10)= DENOM*-CI * V3(4)*(F2(3)*(+CI*(P1(2))-P1(1))+(F2(4)
> $ *(P1(0)+P1(3))-F2(6)*M1))
> F1(11)= DENOM*-CI * V3(5)*(F2(5)*(P1(0)+P1(3))+(F2(6)*(P1(1)
> $ +CI*(P1(2)))-F2(3)*M1))
> F1(12)= DENOM*CI * V3(5)*(F2(5)*(+CI*(P1(2))-P1(1))+(F2(6)
> $ *(P1(3)-P1(0))+F2(4)*M1))
> F1(13)= DENOM*CI * V3(5)*(F2(3)*(P1(3)-P1(0))+(F2(4)*(P1(1)
> $ +CI*(P1(2)))+F2(5)*M1))
> F1(14)= DENOM*-CI * V3(5)*(F2(3)*(+CI*(P1(2))-P1(1))+(F2(4)
> $ *(P1(0)+P1(3))-F2(6)*M1))
> F1(15)= DENOM*-CI * V3(6)*(F2(5)*(P1(0)+P1(3))+(F2(6)*(P1(1)
> $ +CI*(P1(2)))-F2(3)*M1))
> F1(16)= DENOM*CI * V3(6)*(F2(5)*(+CI*(P1(2))-P1(1))+(F2(6)
> $ *(P1(3)-P1(0))+F2(4)*M1))
> F1(17)= DENOM*CI * V3(6)*(F2(3)*(P1(3)-P1(0))+(F2(4)*(P1(1)
> $ +CI*(P1(2)))+F2(5)*M1))
> F1(18)= DENOM*-CI * V3(6)*(F2(3)*(+CI*(P1(2))-P1(1))+(F2(4)
> $ *(P1(0)+P1(3))-F2(6)*M1))
> END
>
>
> Generating the production and decays separately works, i.e:
>
> generate p p > h > hf1 hf1~
>
> and
>
> generate hf1 > hf2 ha
> add process hf1~ > hf2~ ha
>
> both work fine on their own.
>
> So what is wrong? Any ideas? I can send through the model file if you'd like.
>
> Thanks,
>
> J
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.