





title: SPheno files



permalink: /SPheno_files/










[Category:SPheno](/Category:SPheno "wikilink") All routines generated by



SARAH are strongly inspired by the intrinsic routines of SPheno to ensure that they interact nicely with all other SPheno functions. We give here some details of the different files which are written by



SARAH as well as of the main functions for the case that the user wants to do some changes.






### <span>BranchingRatios_\[Model\].f90</span>






Contains the routines for calculating the branching ratios for all particles. Checks are done, if three body decays for fermions are necessary or not.






### <span>Couplings_\[Model\].f90</span>






All vertices calculated by



SARAH are written to this files. Also functions for calculating different subsets of these couplings are generated:






 <span>AllCouplings</span>: all couplings for the mass eigenstates (used for the calculation of decays)



 <span>CouplingsForSMfermion</span>: couplings involved in the one loop self energy of SM fermions



 <span>CouplingsForVectorBosons</span>: couplings involved in the one loop self energy of SM vector bosons



 <span>CouplingsForTadpoles</span>: couplings involved in the calculation of the one loop tadpoles



 <span>CouplingsForLoop</span>: couplings needed for the one loop self energies of all particles






For all couplings involved in loop calculation, the mixing matrices of external particles are replaced by the identity matrix.






### <span>CouplingsForDecays_\[Model\].f90</span>






Calculates the running couplings at the mass scale of each decaying particle. The SUSY parameters are obtained by running the RGEs from the SUSY scale to the desired energy scale. The calculation of the SM gauge and Yukawa coupling is different for particles with masses above or below the SUSY scale. For lighter particles, the SM RGEs are used to run from *M*<sub>*Z*</sub> to the mass scale, for heavier particles the SUSY RGEs are taken to run from *M*<sub>*S**U**S**Y*</sub>.



In addition, the running treelevel masses are calculated. These are used to get effective couplings of the Higgs bosons to two gluons or photons. Furthermore, also the couplings of a light Higgs to a real and virtual vector boson (*W*, *Z*) are calculated.






### <span>InputOutput_\[Model\].f90</span>






This file contains the routines for reading the <span>LesHouches.in.\[Model\]</span> input file and writing the output <span>SPheno.spc.\[Model\]</span>.






##### Input






The following blocks are changed in comparison to the standard LesHouches/<span>SPheno</span> input, see also sec. \[sec:SPHENOINPUT\].






 <span>Block MODSEL</span>



 <span>Block MINPAR</span>



 <span>Block EXTPAR</span>



 <span>Block MODSEL</span>



 <span>Block SPHENOINPUT</span>






##### Output






Routines for writing the following information to the LesHouches spectrum file are generated:






 Parameters at the SUSY scale:



As block names for the output, the entries in <span>LesHouches</span> in the parameter file of



SARAH are used. If this entry is missing, the name for the block is generated automatically.



 Masses at the SUSY scale:



The PDG of the particles file of SARAH is used. If this is missing or 0, the mass is not written in order to exclude unphysical states.



 Optionally, decay widths and branching ratios for the particle in <span>ListDecayParticles</span> and <span>ListDecayParticles3B</span>



 Optionally, results for low energy constraints



 Optionally, effective couplings of the Higgs to SM particles which can be used for instance by <span>HiggsBounds</span>



 Optionally, the GUT values of all parameters



 Optionally, the loop contributions of all particles to the effective coupling of the Higgs to two photons or gluons






In addition, that file contains a routine to write the parameters to a WHIZARD specific output file






### <span>LoopCouplings_\[Model\].f90</span>






Routines for calculating the following couplings:






 Corrections to SM gauge couplings due to all particles heavier than *M*<sub>*Z*</sub>, i.e. normally SUSY particles and top quark (<span>AlphaSDR, AlphaEwDR</span>)



 Effective couplings of scalar/pseudoscalar Higgs to SM fermions (e.g. `CouphhtoFe`).



 `DeltaVB`: oneloop correction to *G*<sub>*F*</sub> calculated from <math>\\mu






\\rightarrow e \\sum_{i j} \\nu_i \\bar{\\nu}_j</math>






### <span>LoopMasses_\[Model\].f90</span>






This file contains the routines for calculating the one loop contributions to the tadpoles and the self energies:






 <span>OneLoopMasses</span>: main routine for calling all other functions



 <span>TadpolesOneLoop</span>: calculates the one loop tadpoles.



 <span>Pi1LoopX</span>: one loop self energy for particle <span>X</span> (scalars or vector bosons (transverse part))



 <span>Sigma1LoopX</span>: calculates *Σ*<sub>*L*</sub>, *Σ*<sub>*R*</sub>, *Σ*<sub>*S*</sub> for fermions



 <span>OneLoopX</span>: calculates the one loop corrected mass for particle <span>X</span> in an iterative way in order to solve



*D**e**t*\[*p*<sub>*i*</sub><sup>2</sup>**1**−*m*<sub>1*L*</sub><sup>2, *X*</sup>(*p*<sup>2</sup>)\] = 0,






for external masses onshell, i.e. *p*<sup>2</sup> = *m*<sub>*T**r**e**e*</sub><sup>2, *X*</sup>.






All these calculation are done in ’t Hooft gauge. Normally, the results are oneloop <span>$\\overline{\\mathrm{DR}}$</span> masses and mixing matrices. For further calculations, the mixing matrices corresponding to an external momentum equal to the heaviest mass eigenstate is taken. However, for SM particles the pole masses are used in the decays as well as in the output.






### <span>LowEnergy_\[Model\].f90</span>






Contains the routines for the calculation of precision observables:






 *b* → *s**γ* (<span>BToQGamma</span>)



 *l*<sub>*i*</sub> → *l*<sub>*j*</sub>*γ* (<span>BrLgammaLp</span>)



 *l*<sub>*i*</sub> → 3*l*<sub>*j*</sub> (with *l* = (*e*, *μ*, *τ*)), (<span>BR1LeptonTo3Leptons</span>)



 the anomalous magnetic moment of leptons (<span>Gminus2</span>)



 electric dipole moments of the charged leptons (<span>LeptonEDM</span>)



 $\\delta\\rho=1\\rho=\\frac{\\Pi_{WW}(0)}{m_W^2}\\frac{\\Pi_{ZZ}(0)}{m_Z^2}$ (*Π*<sub>*Z**Z*</sub>, *Π*<sub>*W**W*</sub> are the selfenergies of the massive vector bosons). (<span>DeltaRho</span>)



 *μ* − *e* conversion in nuclei (Al, Ti, Sr, Sb, Au, Pb) based on the results of (<span>BrLLpHadron</span>)



 *τ* → *l**P*<sup>0</sup> with a pseudoscalar meson *P*<sup>0</sup> (<span>BrLLpHadron</span>) (*π*<sup>0</sup>, *η*, *η*′) based on the results of



 *Z* → *l*<sub>*i*</sub>*l*<sub>*j*</sub> calculated and implemented by Kilian Nickel (based on generic results obtained by FeynArts /<span>FormCalc</span> ) (<span>BrZLLp</span>)



 *B*<sub>*s*, *d*</sub> → *l*<sub>*i*</sub>*l*<sub>*j*</sub> calculated and implemented by Kilian Nickel (based on generic results obtained by FeynArts /<span>FormCalc</span> ) (<span>BrB0LLp</span>)






Note, the calculation is based on the approach and conventions given in the references. However, these expressions are generalized and



SARAH includes new contributions possible in an extended SUSY model. This is done be calculating and implementing the amplitudes for each generic possible diagram for a given process in a general form. The necessity of this is for instance pronounced in Ref. where it has been shown that the *Z*penguins to *l*<sub>*i*</sub> → 3*l*<sub>*j*</sub> often small in the MSSM can dominate in other models like inverse seesaw.






### <span>ModelData_\[Model\].f90</span>






Contains the declaration of global variables and a function for initializing all variables with 0. Only masses of particles, which are integrated out, are initialized with 10<sup>16</sup> GeV. That’s used for the threshold corrections in the first iteration.



In addition, functions for the calculation of GMSB boundary conditions as well as the content of <span>SelfDefinedFunctions</span> is written to that file, see also sec. \[sec:sphenoinputfile\].






### <span>RGEs_\[Model\].f90</span>






That file includes all information about the RGEs. In principle,



SARAH generates three different sets of RGEs for …






 …running of all parameters (gauge couplings, superpotential parameter, soft breaking parameters) from the GUT scale to the SUSY scale



 …running of all parameters and VEVs from SUSY scale to electroweak scale and back



 …running of a minimal set of parameters from electroweak scale to GUT scale: that includes all SM gauge and Yukawa couplings as well as the other parameter the running depends on. Above a threshold scale, the set of parameters consists of all parameter related to the SM gauge and Yukawa couplings.






For each set of RGEs the following function are generated:






 <span>ParametersToGXX</span>: saves the parameters in a vector of length XX



 <span>GToParametersXX</span>: extracts the parameters from a vector of length XX



 <span>rgeXX</span>: definition of all *β*functions.






SARAH does the following simplification/modifications of the *β*functions before exporting it to Fortran code, in order to increase the speed of the numerical calculation






 All matrix multiplications are replaced by constants, which are calculated at the beginning of each <span>rgeXX</span> routine



 All powers of numbers/parameters are replaced by a constants, which are also calculated first






### <span>Shifts_\[Model\].f90</span>






That file is only created for models including at least one threshold scale with gauge symmetry breaking. It contains the necessary routines to calculate the finite shifts for gauge couplings and gaugino masses at the threshold.






### <span>SPheno_\[Model\].f90</span>






The main program. Calls routines for reading the input, calculating the spectrum, decay widths and low energy observables and writing the output.






### <span>SugraRuns_\[Model\].f90</span>






Main routines for calculating the (s)particle spectrum:






 <span>FirstGuess</span>: calculates a first, approximate spectrum as starting point



 <span>SugraRuns</span>: calculates in an iterative way the spectrum



 <span>BoundaryEW</span>: calculates the starting point of the running for gauge and Yukawa couplings from electroweak data



 <span>BoundaryHS, BoundarySUSY</span>: sets the boundary conditions at the GUT and SUSY scale.



 <span>BoundaryConditionsUp,BoundaryConditionsDown</span>: applies the boundary conditions at the threshold scales when running up/down






### <span>SusyDecays_\[Model\].f90</span>






Calculates all two body decays by using the phase space functions of SPheno . If scalar and pseudo scalar Higgs as well as the quarks are defined in the particles file by the corresponding <span>Description</span> statement, the one loop corrections to the decays *H*, *A* → *q**q̄* from gluons are added, see app. \[app:SPhenoDecays\] for more information.






### <span>SusyMasses_\[Model\].f90</span>






This file contains functions to calculate all running masses at tree level. The mass matrices involving Goldstone bosons are taken in ’t Hooft gauge. Note, that works only, if the Goldstone bosons have been associated with the corresponding vector boson in the <span>particles.m</span> file, see sec. \[sec:gaugefixing\]. The main routines in that file are






 <span>TreeMasses</span>: main function to call the other subroutines



 <span>CalculateX</span>: derives the tree level mass and mixing matrix of particle <span>X</span> by solving the eigensystem of the corresponding mass matrix.






### Three Body Decays, e.g. <span>Glu_\[Model\].f90</span>






Contains all necessary routines to calculate the branching ratios and decay width of a fermion to three other fermions. All possible processes and diagrams are generated by



SARAH and mapped to the phase space function of SPheno .






### Additional Files






 A model specific <span>Makefile</span>



 Template for an input file of the current model: <span>LesHouches.in.\[Model\]</span>, see app. \[sec:LHinput:file\].



 Template for an input file assuming a SUSY scale input: <span>LesHouches.in.\[Model\]_low</span>, see also app. \[sec:LHinput:file\].






LesHouches input file










The LesHouches input file for each model is by default named






LesHouches.in.[Model]






However, it is also possible to give the name of the in and output file as option when running SPheno , see also Ref. :






./[SPheno Directory]/bin/SPheno[Model] [Input file] [Output file]






While the entries of the common block <span>MINPAR</span> and <span>EXTPAR</span> are completely defined by the SPheno input file for



SARAH, some other blocks have distinct values.






### <span>SMINPUTS</span>






The entries of this block accordingly to the SLHA conventions are:






### <span>MODSEL</span>






The block <span>MODSEL</span> is not totally equivalent to the SLHA conventions: since each SPheno module generated by



SARAH handles one specify model, the flags <span>3</span> and <span>4</span> have no effect. In addition, we added the flag <span>2</span> to choose between the different boundary conditions, and flag <span>1</span> is used to choose between a GUT scale and SUSY scale input.






### <span>SPHENOINPUT</span>






The block with SPheno specific commands is similar to the conventions used by the standard SPheno version 3.1 and above. However, some switches have no effect and there are also some new entries which we highlight with in bold.






See also



 


\ No newline at end of file 