FHGetSelf

       FHGetSelf - calculate the renormalized Higgs self-energies

       integer error
       double complex k2
       integer key, dkey, ren
       double complex sig(13), dsig(13)

       integer h0h0, HHHH, A0A0, HmHp
       integer h0HH, h0A0, HHA0
       integer G0G0, h0G0, HHG0, A0G0
       integer GmGp, HmGp
       parameter (h0h0 = 1, HHHH = 2, A0A0 = 3, HmHp = 4)
       parameter (h0HH = 5, h0A0 = 6, HHA0 = 7)
       parameter (G0G0 = 8, h0G0 = 9, HHG0 = 10, A0G0 = 11)
       parameter (GmGp = 12, HmGp = 13)
       #define Key(se) 2**(se-1)

       subroutine FHGetSelf(error, k2, key, sig, dkey, dsig, ren)

       FHGetSelf[k2, key, dkey, ren]

       FHGetSelf  returns  the  renormalized  Higgs  self-energies at momentum-
       squared k2.  The flags and parameters must have been set before with FH-
       SetFlags(1) and FHSetPara(1).

       error (output)
              zero if successful, otherwise the line number in  GetSelf.F  from
              which the error message was emitted

       k2 (input)
              the  momentum-squared  k^2  at which the self-energies are evalu-
              ated.

       key (input)
              a flag determining which of the self-energies are actually evalu-
              ated, e.g. to evaluate the h0-h0 self-energy,  add  Key(h0h0)  to
              key.

       dkey (input)
              a  flag determining which of the derivatives of the self-energies
              are actually evaluated, e.g. to evaluate the  derivative  of  the
              h0-h0 self-energy, add Key(h0h0) to dkey.

       ren (input)
              whether  the unrenormalized (0) or renormalized (1) self-energies
              are output.

       sig(h0h0), sig(HHHH), sig(A0A0), sig(HmHp) (output)
              the h0, HH, A0, and Hp self-energies at k^2 = k2.

       sig(h0HH), sig(h0A0), sig(HHA0) (output)
              the h0-HH, h0-A0, and HH-A0 mixing self-energies at k^2 = k2.

       sig(G0G0), sig(h0G0), sig(HHG0), sig(A0G0) (output)
              the neutral Goldstone self-energies at k^2 = k2.

       sig(GmGp), sig(HmGp) (output)
              the charged Goldstone self-energies at k^2 = k2.

       dsig(i) (output)
              the derivatives of the self-energies with respect to k^2 at k^2 =
              k2, where the index i runs as for the sig(i).

       libFH(1)

                                  24-Apr-2015                      FHGETSELF(1)