Page 1 of 1

Vertices output in SPheno

Posted: 14. Nov 2017, 13:02
by samadrita
Hi,

I want to get some desired vertices information in the SPheno output file.
SARAH provides all the vertices in terms of mixing matrices and diagonalising matrices in its MakeTeX[ ] output.
But how can we get the numerical values of these vertices in the SPheno output files.
Is there any specific flag to write the same?

Regards,
Samadrita

Re: Vertices output in SPheno

Posted: 14. Nov 2017, 13:55
by FStaub
Hi,

no, there is no option to write numerical values in the SPheno output.
The easiest option to get values for vertices is to read in the SPheno spc in Mathematica after running Start["MODEL"], ie. something like

Code: Select all

ReadSpectrumFile["SPheno.spc.MSSM", ""];
sub = MakeSubNum /@ AllLesHouchesInput;
(* => brings the content of the spc in form of a substitution list *)

subLH = Table[
  ParameterDefinitions[[i, 1]] -> (LesHouches /. 
      ParameterDefinitions[[i, 2]]) /. {List[a_, b_] :> 
     ToExpression[ToUpperCase[ToString[a]]][b]}, {i, 1, 
   Length[ParameterDefinitions]}]
(* => substitutes couplings by the Les Houches entries *

(* example: all quartic CP-even Higgs vertices *)   
TestVertex = Vertex[{hh, hh, hh, hh}][[2, 1]];
Flatten[Table[{gt1, gt2, gt3, gt4, TestVertex}, {gt1, 1, 2}, {gt2, 1, 2}, {gt3, 1, 2}, {gt4, 1, 2}], 3] /. subLH //. sub[[1]]
Cheers,
Florian

Re: Vertices output in SPheno

Posted: 15. Nov 2017, 07:41
by samadrita
Dear Florian,

Thanks for your speedy reply.

However, I want a further clarification in this regard.
It was written in [https://arxiv.org/abs/1109.5147] that, "it is also possible to use some self defined function. The Fortran code of that
function has to included in the array SelfDefinedFunctions in SPheno.m. Later on it will be written to Model Data.f90".

So if we write the Mathematica expressions of vertices as some functions inside SPheno.m under some new block, then is it possible to write the numerical values of those vertices in SPheno.spc.[$Model] every time we change LesHouches.in.[$Model] ?
If yes then how to write that part inside SPheno.spc ?

Thanking you,
Samadrita

Re: Vertices output in SPheno

Posted: 15. Nov 2017, 14:11
by FStaub
Hi,

the sentence you quote refers explicitly to the definition of boundary conditions. That's not related to the topic here.

If you want to see vertices in the spc output, you need to hack the InputOutput_$MODEL.f90 file. You don't need to write the vertices again, but you can call the functions from Couplings_$MODEL.f90 (after loading this module) which contain the expressions. There is no way to do this automatically, but some knowledge of Fortran might be needed.

Cheers,
Florian

Re: Vertices output in SPheno

Posted: 16. Nov 2017, 06:35
by samadrita
Thanks for the clarification..
We will see what we can do with the Fortran code.

Best Regards,
Samadrita