Problem related with Cot coupling

Asked by Disha Bhatia on 2019-09-30

Hello,

I am facing an extremely strange problem. I am defining an effective model with pseudoscalar coupling to GG and AA. My coupling strength is proportional to Cot[beta]. On running madgraph, I am getting following error which I am not able to follow:
  gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o couplings.o couplings.f
     make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
     make[1]: Entering directory '/home/disha/software_required/MG5_aMC_v2_6_6/bin/PROC_Higgs_Effective_Couplings_UFO_3/Source/CERNLIB'
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o abend.o abend.f
     intparam_definition.inc:55:21:

              MDL_RP1 = COT(DCMPLX(MDL_BETA))*MDL_SIN__BETAP
                          1
     Error: Function ‘cot’ at (1) has no IMPLICIT type
     intparam_definition.inc:55:21:

              MDL_RP1 = COT(DCMPLX(MDL_BETA))*MDL_SIN__BETAP
                          1
     Error: Function ‘cot’ at (1) has no IMPLICIT type
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o dlsqp2.o dlsqp2.f
     <builtin>: recipe for target 'couplings.o' failed
     make[1]: *** [couplings.o] Error 1
     make[1]: Leaving directory '/home/disha/software_required/MG5_aMC_v2_6_6/bin/PROC_Higgs_Effective_Couplings_UFO_3/Source/MODEL'
     makefile:53: recipe for target '../lib/libmodel.a' failed
     make: *** [../lib/libmodel.a] Error 2
     make: *** Waiting for unfinished jobs....

Can you please help me out with this error. I tried converting Cot with Cos/Sin, but it didn't help.

Thank You,
Disha Bhatia

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Disha Bhatia
Solved:
2019-09-30
Last query:
2019-09-30
Last reply:
2019-09-30

Hi,

Cot function is not part of the default f77 routine.
So in order to use it, you need to define it in the file functions_library.py of your UFO model.

like:

cot = Function(name = 'cot',
               arguments = ('z',),
               expression = '1./cmath.tan(z.real)')

or like

cot = Function(name = 'cot',
               arguments = ('z',),
               expression = '1./cmath.tan(z)')

# if you want to allow complex evaluation fo the tan function

This being said, i do not understand why your change is not working.
I would need more details about it in order to be able to comment.

Cheers,

Olivier

> On 30 Sep 2019, at 15:43, Disha Bhatia <email address hidden> wrote:
>
> New question #684493 on MadGraph5_aMC@NLO:
> https://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fanswers.launchpad.net%2Fmg5amcnlo%2F%2Bquestion%2F684493&amp;data=02%7C01%7Colivier.mattelaer%40uclouvain.be%7Ca3acaa38705b467a3ff108d745ac26c9%7C7ab090d4fa2e4ecfbc7c4127b4d582ec%7C0%7C0%7C637054478640626365&amp;sdata=oU9wq8Mogof7tsh2eDuunioa%2B9UQPanYtrRISU35NPs%3D&amp;reserved=0
>
> Hello,
>
> I am facing an extremely strange problem. I am defining an effective model with pseudoscalar coupling to GG and AA. My coupling strength is proportional to Cot[beta]. On running madgraph, I am getting following error which I am not able to follow:
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o couplings.o couplings.f
> make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
> make[1]: Entering directory '/home/disha/software_required/MG5_aMC_v2_6_6/bin/PROC_Higgs_Effective_Couplings_UFO_3/Source/CERNLIB'
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o abend.o abend.f
> intparam_definition.inc:55:21:
>
> MDL_RP1 = COT(DCMPLX(MDL_BETA))*MDL_SIN__BETAP
> 1
> Error: Function ‘cot’ at (1) has no IMPLICIT type
> intparam_definition.inc:55:21:
>
> MDL_RP1 = COT(DCMPLX(MDL_BETA))*MDL_SIN__BETAP
> 1
> Error: Function ‘cot’ at (1) has no IMPLICIT type
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o dlsqp2.o dlsqp2.f
> <builtin>: recipe for target 'couplings.o' failed
> make[1]: *** [couplings.o] Error 1
> make[1]: Leaving directory '/home/disha/software_required/MG5_aMC_v2_6_6/bin/PROC_Higgs_Effective_Couplings_UFO_3/Source/MODEL'
> makefile:53: recipe for target '../lib/libmodel.a' failed
> make: *** [../lib/libmodel.a] Error 2
> make: *** Waiting for unfinished jobs....
>
>
> Can you please help me out with this error. I tried converting Cot with Cos/Sin, but it didn't help.
>
> Thank You,
> Disha Bhatia
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Disha Bhatia (dishabhatia1989) said : #2

Hi,

The change worked now. Actually I put the sin[beta] factor directly into
the Lagrangian and kept Cos[beta] as a part of the coupling constant so
that Mathematica doesn't convert it back to Cot[beta].
Now i worked. Sorry should have tried this also first before posting,

Thanks,
Disha

On Mon, Sep 30, 2019 at 4:22 PM Olivier Mattelaer <
<email address hidden>> wrote:

> Your question #684493 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/684493
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi,
>
> Cot function is not part of the default f77 routine.
> So in order to use it, you need to define it in the file
> functions_library.py of your UFO model.
>
> like:
>
> cot = Function(name = 'cot',
> arguments = ('z',),
> expression = '1./cmath.tan(z.real)')
>
> or like
>
> cot = Function(name = 'cot',
> arguments = ('z',),
> expression = '1./cmath.tan(z)')
>
> # if you want to allow complex evaluation fo the tan function
>
> This being said, i do not understand why your change is not working.
> I would need more details about it in order to be able to comment.
>
> Cheers,
>
> Olivier
>
> > On 30 Sep 2019, at 15:43, Disha Bhatia <
> <email address hidden>> wrote:
> >
> > New question #684493 on MadGraph5_aMC@NLO:
> >
> https://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fanswers.launchpad.net%2Fmg5amcnlo%2F%2Bquestion%2F684493&amp;data=02%7C01%7Colivier.mattelaer%40uclouvain.be%7Ca3acaa38705b467a3ff108d745ac26c9%7C7ab090d4fa2e4ecfbc7c4127b4d582ec%7C0%7C0%7C637054478640626365&amp;sdata=oU9wq8Mogof7tsh2eDuunioa%2B9UQPanYtrRISU35NPs%3D&amp;reserved=0
> >
> > Hello,
> >
> > I am facing an extremely strange problem. I am defining an effective
> model with pseudoscalar coupling to GG and AA. My coupling strength is
> proportional to Cot[beta]. On running madgraph, I am getting following
> error which I am not able to follow:
> > gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o
> couplings.o couplings.f
> > make[1]: warning: jobserver unavailable: using -j1. Add '+'
> to parent make rule.
> > make[1]: Entering directory
> '/home/disha/software_required/MG5_aMC_v2_6_6/bin/PROC_Higgs_Effective_Couplings_UFO_3/Source/CERNLIB'
> > gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132
> -c -o abend.o abend.f
> > intparam_definition.inc:55:21:
> >
> > MDL_RP1 = COT(DCMPLX(MDL_BETA))*MDL_SIN__BETAP
> > 1
> > Error: Function ‘cot’ at (1) has no IMPLICIT type
> > intparam_definition.inc:55:21:
> >
> > MDL_RP1 = COT(DCMPLX(MDL_BETA))*MDL_SIN__BETAP
> > 1
> > Error: Function ‘cot’ at (1) has no IMPLICIT type
> > gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132
> -c -o dlsqp2.o dlsqp2.f
> > <builtin>: recipe for target 'couplings.o' failed
> > make[1]: *** [couplings.o] Error 1
> > make[1]: Leaving directory
> '/home/disha/software_required/MG5_aMC_v2_6_6/bin/PROC_Higgs_Effective_Couplings_UFO_3/Source/MODEL'
> > makefile:53: recipe for target '../lib/libmodel.a' failed
> > make: *** [../lib/libmodel.a] Error 2
> > make: *** Waiting for unfinished jobs....
> >
> >
> > Can you please help me out with this error. I tried converting Cot with
> Cos/Sin, but it didn't help.
> >
> > Thank You,
> > Disha Bhatia
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadGraph5_aMC@NLO.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/mg5amcnlo/+question/684493/+confirm?answer_id=0
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/mg5amcnlo/+question/684493
>
> You received this question notification because you asked the question.
>

--
Regards,
Disha Bhatia

Disha Bhatia (dishabhatia1989) said : #3

Hi,

The change worked now. Actually I put the sin[beta] factor directly into
the Lagrangian and kept Cos[beta] as a part of the coupling constant so
that Mathematica doesn't convert it back to Cot[beta].
Now i worked. Sorry should have tried this also first before posting,

Thanks,
Disha