Problem related with Cot coupling

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.

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

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