MadGraph5_aMC@NLO

Bug #1955926
Comment #4

Comment 4 for bug 1955926

Revision history for this message

Olivier Mattelaer (olivier-mattelaer) wrote on 2021-12-29: Re: [Bug 1955926] MG Output Error for H > a a Reaction using @NLO Extension of Model SM

Hi,

So the problem is within this expression:

I execute UVGC_153_116 = ( ( ( 0 if mdl_MW==0 else (mdl_ee__exp__4*mdl_complexi*mdl_vev__exp__2)/(16.*cmath.pi**2*mdl_MW__exp__2*mdl_sw__exp__4) - (mdl_ee__exp__4*mdl_complexi*mdl_vev__exp__2)/(24.*cmath.pi*mdl_MW__exp__2*mdl_sw__exp__4*mdl_sqrt__3) ) if mdl_MH==mdl_MW else ( 0 if mdl_MW==0 else ( 0 if mdl_MH__exp__2<=4*mdl_MW__exp__2 else -(mdl_ee__exp__4*mdl_complexi*mdl_MW__exp__2*mdl_vev__exp__2*reglog(-(mdl_MU_R__exp__2*(-(mdl_MH__exp__2/mdl_MU_R__exp__2) + (2*mdl_MW__exp__2)/mdl_MU_R__exp__2 + (mdl_MH*cmath.sqrt(mdl_MH__exp__2/mdl_MU_R__exp__2 - (4*mdl_MW__exp__2)/mdl_MU_R__exp__2))/MU_R))/(2.*mdl_MW__exp__2)))/(8.*MU_R*cmath.pi**2*mdl_MH__exp__3*mdl_sw__exp__4*cmath.sqrt(mdl_MH__exp__2/mdl_MU_R__exp__2 - (4*mdl_MW__exp__2)/mdl_MU_R__exp__2)) ) + ( -(mdl_ee__exp__4*mdl_complexi*mdl_MW__exp__2*mdl_vev__exp__2*reglog((mdl_MU_R__exp__2*(-(mdl_MH__exp__2/mdl_MU_R__exp__2) + (2*mdl_MW__exp__2)/mdl_MU_R__exp__2 + (mdl_MH*cmath.sqrt(mdl_MH__exp__2/mdl_MU_R__exp__2 - (4*mdl_MW__exp__2)/mdl_MU_R__exp__2))/MU_R))/(2.*mdl_MW__exp__2)))/(8.*MU_R*cmath.pi**2*mdl_MH__exp__3*mdl_sw__exp__4*cmath.sqrt(mdl_MH__exp__2/mdl_MU_R__exp__2 - (4*mdl_MW__exp__2)/mdl_MU_R__exp__2)) if mdl_MH__exp__2<4*mdl_MW__exp__2 else 0 ) + (mdl_ee__exp__4*mdl_complexi*mdl_vev__exp__2)/(16.*cmath.pi**2*mdl_MH__exp__2*mdl_sw__exp__4) ) + ( ( 0 if mdl_MW==0 else (mdl_ee__exp__4*mdl_complexi*mdl_vev__exp__2)/(16.*cmath.pi**2*mdl_MH__exp__2*mdl_sw__exp__4) ) if mdl_MH__exp__2==4*mdl_MW__exp__2 else 0 ) ) + ( (mdl_ee__exp__4*mdl_complexi*mdl_vev__exp__2)/(16.*cmath.pi**2*mdl_MH__exp__2*mdl_sw__exp__4) if mdl_MW==0 else 0 ) if mdl_MH else ( 0 if mdl_MW==0 else -(mdl_ee__exp__4*mdl_complexi*mdl_vev__exp__2)/(96.*cmath.pi**2*mdl_MW__exp__2*mdl_sw__exp__4) ) )

given the error, the issue is within the sub-expression
mdl_MH__exp__2<=4*mdl_MW__exp__2

So it seems that your W mass is complex and not real.
Do you use some kind of complex mass-scheme or is your benchmark not setup correctly?

Cheers,

Olivier

> On 29 Dec 2021, at 00:37, Brian Keith Ragsdale <email address hidden> wrote:
>
> Nice to hear from you Olivier,
>
> I hope you can access the UFO folder from my drive. Here's the link:
> https://drive.google.com/drive/folders/1JYjoBHBmpxRQnOxq8ZBl7V7HQBtIzHcf?usp=sharing
>
> I'm looking forward to your response!
>
> Regards.
> Brian Keith Ragsdale II
>
> P.S. I've been enjoying yours and other's work in the Tutorials and Help
> directory on the MadGraph website. So amazed by everyone's work and have
> myself a nice collection of reading materials for the next few days
> while I'm out of the office. Cheers!
>
> --
> You received this bug notification because you are subscribed to
> MadGraph5_aMC@NLO.
> https://bugs.launchpad.net/bugs/1955926
>
> Title:
> MG Output Error for H > a a Reaction using @NLO Extension of Model SM
>
> Status in MadGraph5_aMC@NLO:
> New
>
> Bug description:
> Hello,
>
> I used the archetypical SM.fr file to generate an FA folder (GEN, MOD,
> PARS) in FR to obtain CT and R terms with the goal of generating a UFO
> file for the SM@NLO. My model process to determine the success of that
> goal was the Decay Reaction Higgs to DiPhoton. Details of the
> computations are below, followed by the terminal log.
>
> /////////////////////////////////////////////////////////////////////////////
> **FeynArts**
> (3.11)
>
> **FeynRules**
>
> (
> Update notes for FeynRules 2.3.x
>
> 2642 2.3.47 25-05-21 fuks: Mathematica 12.2 and parallelisation are
> fixed
> )
>
> **MadGraph**
> (
> Update notes for MadGraph5_aMC@NLO (in reverse time order)
>
> 2.9.4(28/05/21)
> )
>
> **Mathematica**
> (12.0)
> /////////////////////////////////////////////////////////////////////////////
>
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
> Below is the terminal log with the error highlighted above and below by a series of '-' at the end
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>
> MG5_aMC>import model SMNLO122821/
> INFO: load particles
> INFO: load vertices
> DEBUG: MG5 converter defines UUV3 to P(3,2) + P(3,3)
> DEBUG: MG5 converter defines VVV10 to P(3,1)*Metric(1,2) + P(2,3)*Metric(1,3) + P(1,2)*Metric(2,3)
> DEBUG: MG5 converter defines VVV11 to P(3,2)*Metric(1,2) + P(2,1)*Metric(1,3) + P(1,3)*Metric(2,3)
> DEBUG: MG5 converter defines FFS4 to ProjM(2,1) + ProjP(2,1)
> DEBUG: MG5 converter defines VVVV5 to Metric(1,4)*Metric(2,3) + Metric(1,3)*Metric(2,4)
> DEBUG: MG5 converter defines VVVV6 to Metric(1,3)*Metric(2,4) + Metric(1,2)*Metric(3,4)
> DEBUG: MG5 converter defines VVVV7 to Metric(1,4)*Metric(2,3) + Metric(1,3)*Metric(2,4) + Metric(1,2)*Metric(3,4)
> DEBUG: MG5 converter defines VVVV8 to Metric(1,4)*Metric(2,3) + Metric(1,2)*Metric(3,4)
> DEBUG: MG5 converter defines VVV12 to P(3,1)*Metric(1,2) + P(1,2)*Metric(2,3)
> DEBUG: MG5 converter defines VVV13 to P(3,2)*Metric(1,2) + P(1,3)*Metric(2,3)
> DEBUG: MG5 converter defines VVV14 to P(3,1)*Metric(1,2) + P(2,3)*Metric(1,3)
> DEBUG: MG5 converter defines VVV15 to P(3,2)*Metric(1,2) + P(2,1)*Metric(1,3)
> DEBUG: MG5 converter defines FFV4 to Gamma(3,2,-1)*ProjM(-1,1) + Gamma(3,2,-1)*ProjP(-1,1)
> DEBUG: MG5 converter defines FF6 to ProjM(2,1) + ProjP(2,1)
> DEBUG: MG5 converter defines FF7 to P(-1,1)*Gamma(-1,2,-2)*ProjM(-2,1) + P(-1,1)*Gamma(-1,2,-2)*ProjP(-2,1)
> DEBUG: MG5 converter defines VVV16 to P(2,1)*Metric(1,3) + P(1,2)*Metric(2,3)
> DEBUG: MG5 converter defines VVV17 to P(2,3)*Metric(1,3) + P(1,2)*Metric(2,3)
> DEBUG: model prefixing takes 5.1495606899261475
> This model does not allow Feynman gauge. You will only be able to do tree level QCD loop cmputations with it.
> INFO: Change particles name to pass to MG5 convention
> Kept definitions of multiparticles p / j / l+ / l- / vl / vl~ unchanged
> Defined multiparticle all = g ghg ghg~ u c d s u~ c~ d~ s~ a gha gha~ ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ ghz ghwp ghwm h w- ghz~ ghwp~ ghwm~ ta- ta+
> MG5_aMC>generate h > a a [noborn = QED]
> You will only be able to do tree level and QCD corrections with this model because it does not support Feynman gauge.
> DEBUG: Generating process: h > a a [ noborn = QED ]
> DEBUG: Born diagrams generation skipped by user request.
> DEBUG: Coupling order combinations considered: (QCD,QED,WEIGHTED)
> DEBUG: > No Born contributions for this process.
> DEBUG: > loop : (0,6,W12)
> INFO: Contributing diagrams generated: 0 Born, 9 loops, 4 R2, 1 UV
> 1 processes with 9 diagrams generated in 0.241 s
> Total: 1 processes with 9 diagrams
> MG5_aMC>output haa122821
> INFO: initialize a new directory: haa122821
> INFO: remove old information in haa122821
> INFO: Organizing processes into subprocess groups
> INFO: Generating Helas calls for process: h > a a [ noborn = QED ]
> DEBUG: 1 loop wavefunctions have been reused, for a total of 25 ones
> DEBUG: Computing the loop color basis
> INFO: Processing color information for loop process: h > a a [ noborn = QED ]
> INFO: Creating color matrix loop process: h > a a [ noborn = QED ]
> INFO: Creating files in directory /home/bandit5/maddev/mg5amcnlo/LTS_dev/haa122821/SubProcesses/PV0_0_1_h_aa
> DEBUG: adding in aloha model 17 lorentz struct
> ALOHA: aloha creates VVS1 set of routines with options: L2
> ALOHA: aloha creates FFV1 set of routines with options: L1
> INFO: Computing diagram color coefficients
> INFO: Drawing loop Feynman diagrams for Process: h > a a [ noborn = QED ]
> INFO: Creating files in directory P0_h_aa
> INFO: Generating Feynman diagrams for Process: h > a a [ noborn = QED ]
> INFO: Finding symmetric diagrams for subprocess group h_aa
> Generated helas calls for 1 subprocesses (9 diagrams) in 0.074 s
> ------------------------------------------------------------------------------
> Command "output haa122821" interrupted with error:
> TypeError : '<=' not supported between instances of 'int' and 'complex'
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in 'MG5_debug'.
> Please attach this file to your report.
> ------------------------------------------------------------------------------
>
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/mg5amcnlo/+bug/1955926/+subscriptions
>

Hi,

So the problem is within this expression:

given the error, the issue is within the sub-expression
 mdl_MH__exp__2<=4*mdl_MW__exp__2

So it seems that your W mass is complex and not real.
Do you use some kind of complex mass-scheme or is your benchmark not setup correctly?

Cheers,

Olivier

> On 29 Dec 2021, at 00:37, Brian Keith Ragsdale <1955926@bugs.launchpad.net> wrote:
> 
> Nice to hear from you Olivier,
> 
> I hope you can access the UFO folder from my drive. Here's the link:
> https://drive.google.com/drive/folders/1JYjoBHBmpxRQnOxq8ZBl7V7HQBtIzHcf?usp=sharing
> 
> I'm looking forward to your response!
> 
> Regards.
> Brian Keith Ragsdale II
> 
> P.S. I've been enjoying yours and other's work in the Tutorials and Help
> directory on the MadGraph website. So amazed by everyone's work and have
> myself a nice collection of reading materials for the next few days
> while I'm out of the office. Cheers!
> 
> -- 
> You received this bug notification because you are subscribed to
> MadGraph5_aMC@NLO.
> https://bugs.launchpad.net/bugs/1955926
> 
> Title:
>  MG Output Error for H > a a Reaction using @NLO Extension of Model SM
> 
> Status in MadGraph5_aMC@NLO:
>  New
> 
> Bug description:
>  Hello,
> 
>  I used the archetypical SM.fr file to generate an FA folder (GEN, MOD,
>  PARS) in FR to obtain CT and R terms with the goal of generating a UFO
>  file for the SM@NLO. My model process to determine the success of that
>  goal was the Decay Reaction Higgs to DiPhoton. Details of the
>  computations are below, followed by the terminal log.
> 
>  /////////////////////////////////////////////////////////////////////////////
>  **FeynArts**
>  (3.11)
> 
>  **FeynRules**
> 
>  (
>  Update notes for FeynRules 2.3.x
> 
>  2642 2.3.47 25-05-21 fuks:       Mathematica 12.2 and parallelisation are 
>                                   fixed 
>  )
> 
>  **MadGraph**
>  (
>  Update notes for MadGraph5_aMC@NLO (in reverse time order)
> 
>  2.9.4(28/05/21)
>  )
> 
>  **Mathematica**
>  (12.0)
>  /////////////////////////////////////////////////////////////////////////////
> 
>  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>  Below is the terminal log with the error highlighted above and below by a series of '-' at the end
>  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
> 
>  MG5_aMC>import model SMNLO122821/
>  INFO: load particles
>  INFO: load vertices
>  DEBUG: MG5 converter defines UUV3 to P(3,2) + P(3,3)
>  DEBUG: MG5 converter defines VVV10 to P(3,1)*Metric(1,2) + P(2,3)*Metric(1,3) + P(1,2)*Metric(2,3)
>  DEBUG: MG5 converter defines VVV11 to P(3,2)*Metric(1,2) + P(2,1)*Metric(1,3) + P(1,3)*Metric(2,3)
>  DEBUG: MG5 converter defines FFS4 to ProjM(2,1) + ProjP(2,1)
>  DEBUG: MG5 converter defines VVVV5 to Metric(1,4)*Metric(2,3) + Metric(1,3)*Metric(2,4)
>  DEBUG: MG5 converter defines VVVV6 to Metric(1,3)*Metric(2,4) + Metric(1,2)*Metric(3,4)
>  DEBUG: MG5 converter defines VVVV7 to Metric(1,4)*Metric(2,3) + Metric(1,3)*Metric(2,4) + Metric(1,2)*Metric(3,4)
>  DEBUG: MG5 converter defines VVVV8 to Metric(1,4)*Metric(2,3) + Metric(1,2)*Metric(3,4)
>  DEBUG: MG5 converter defines VVV12 to P(3,1)*Metric(1,2) + P(1,2)*Metric(2,3)
>  DEBUG: MG5 converter defines VVV13 to P(3,2)*Metric(1,2) + P(1,3)*Metric(2,3)
>  DEBUG: MG5 converter defines VVV14 to P(3,1)*Metric(1,2) + P(2,3)*Metric(1,3)
>  DEBUG: MG5 converter defines VVV15 to P(3,2)*Metric(1,2) + P(2,1)*Metric(1,3)
>  DEBUG: MG5 converter defines FFV4 to Gamma(3,2,-1)*ProjM(-1,1) + Gamma(3,2,-1)*ProjP(-1,1)
>  DEBUG: MG5 converter defines FF6 to ProjM(2,1) + ProjP(2,1)
>  DEBUG: MG5 converter defines FF7 to P(-1,1)*Gamma(-1,2,-2)*ProjM(-2,1) + P(-1,1)*Gamma(-1,2,-2)*ProjP(-2,1)
>  DEBUG: MG5 converter defines VVV16 to P(2,1)*Metric(1,3) + P(1,2)*Metric(2,3)
>  DEBUG: MG5 converter defines VVV17 to P(2,3)*Metric(1,3) + P(1,2)*Metric(2,3)
>  DEBUG: model prefixing  takes 5.1495606899261475
>  This model does not allow Feynman gauge. You will only be able to do tree level QCD loop cmputations with it.
>  INFO: Change particles name to pass to MG5 convention
>  Kept definitions of multiparticles p / j / l+ / l- / vl / vl~ unchanged
>  Defined multiparticle all = g ghg ghg~ u c d s u~ c~ d~ s~ a gha gha~ ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ ghz ghwp ghwm h w- ghz~ ghwp~ ghwm~ ta- ta+
>  MG5_aMC>generate h > a a [noborn = QED]
>  You will only be able to do tree level and QCD corrections with this model because it does not support Feynman gauge.
>  DEBUG: Generating process: h > a a [ noborn = QED ]
>  DEBUG: Born diagrams generation skipped by user request.
>  DEBUG: Coupling order combinations considered: (QCD,QED,WEIGHTED)
>  DEBUG:  > No Born contributions for this process.
>  DEBUG:  > loop : (0,6,W12)
>  INFO: Contributing diagrams generated: 0 Born, 9 loops, 4 R2, 1 UV
>  1 processes with 9 diagrams generated in 0.241 s
>  Total: 1 processes with 9 diagrams
>  MG5_aMC>output haa122821
>  INFO: initialize a new directory: haa122821
>  INFO: remove old information in haa122821
>  INFO: Organizing processes into subprocess groups
>  INFO: Generating Helas calls for process: h > a a [ noborn = QED ]
>  DEBUG: 1 loop wavefunctions have been reused, for a total of 25 ones
>  DEBUG: Computing the loop color basis
>  INFO: Processing color information for loop process: h > a a [ noborn = QED ]
>  INFO: Creating color matrix  loop process: h > a a [ noborn = QED ]
>  INFO: Creating files in directory /home/bandit5/maddev/mg5amcnlo/LTS_dev/haa122821/SubProcesses/PV0_0_1_h_aa
>  DEBUG: adding in aloha model 17 lorentz struct
>  ALOHA: aloha creates VVS1 set of routines with options: L2
>  ALOHA: aloha creates FFV1 set of routines with options: L1
>  INFO: Computing diagram color coefficients
>  INFO: Drawing loop Feynman diagrams for Process: h > a a [ noborn = QED ]
>  INFO: Creating files in directory P0_h_aa
>  INFO: Generating Feynman diagrams for Process: h > a a [ noborn = QED ]
>  INFO: Finding symmetric diagrams for subprocess group h_aa
>  Generated helas calls for 1 subprocesses (9 diagrams) in 0.074 s
>  ------------------------------------------------------------------------------
>  Command "output haa122821" interrupted with error:
>  TypeError : '<=' not supported between instances of 'int' and 'complex'
>  Please report this bug on https://bugs.launchpad.net/mg5amcnlo
>  More information is found in 'MG5_debug'.
>  Please attach this file to your report.
>  ------------------------------------------------------------------------------
> 
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/mg5amcnlo/+bug/1955926/+subscriptions
>