Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
pyrghm.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pyrghm.f
1 
2 C*********************************************************************
3 
4 C...PYRGHM
5 C...Auxiliary to PYPOLE.
6 
7  SUBROUTINE pyrghm(MCHI,MA,TANB,MQ,MUR,MD,MTOP,AU,AD,MU,
8  * mhp,hmp,mch,sa,ca,sab,cab,tanba,mglu,deltamt,deltamb)
9  IMPLICIT DOUBLE PRECISION(a-h,l,m,o-z)
10  dimension vh(2,2),m2(2,2),m2p(2,2)
11 C...Parameters.
12  INTEGER mstu,mstj
13  common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
14  SAVE /pydat1/
15 
16  mz = 91.18d0
17  pi = paru(1)
18  v = 174.1d0
19  alpha1 = 0.0101d0
20  alpha2 = 0.0337d0
21  alpha3z = 0.12d0
22  tanba = tanb
23  tanbt = tanb
24 C MBOTTOM(MTOP) = 3. GEV
25  mb = pymrun(5,mtop**2)
26  alpha3 = alpha3z/(1d0 +(11d0 - 10d0/3d0)/4d0/pi*alpha3z*
27  *log(mtop**2/mz**2))
28 C RMTOP= RUNNING TOP QUARK MASS
29  rmtop = mtop/(1d0+4d0*alpha3/3d0/pi)
30  tq = log((mq**2+mtop**2)/mtop**2)
31  tu = log((mur**2 + mtop**2)/mtop**2)
32  td = log((md**2 + mtop**2)/mtop**2)
33 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
34 C
35 C NEW DEFINITION, TGLU.
36 C
37 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
38  tglu = log(mglu**2/mtop**2)
39  sinb = tanb/dsqrt(1d0 + tanb**2)
40  cosb = sinb/tanb
41  IF(ma.GT.mtop)
42  *tanba = tanb*(1d0-3d0/32d0/pi**2*
43  *(rmtop**2/v**2/sinb**2-mb**2/v**2/cosb**2)*
44  *log(ma**2/mtop**2))
45  IF(ma.LT.mtop.OR.ma.EQ.mtop) tanbt = tanba
46  sinb = tanbt/sqrt(1d0 + tanbt**2)
47  cosb = 1d0/dsqrt(1d0 + tanbt**2)
48  g1 = sqrt(alpha1*4d0*pi)
49  g2 = sqrt(alpha2*4d0*pi)
50  g3 = sqrt(alpha3*4d0*pi)
51  hu = rmtop/v/sinb
52  hd = mb/v/cosb
53  CALL pygfxx(ma,tanba,mq,mur,md,mtop,au,ad,mu,mglu,vh,stop1,stop2,
54  *sbot1,sbot2,deltamt,deltamb)
55  IF(mq.GT.mur) tp = tq - tu
56  IF(mq.LT.mur.OR.mq.EQ.mur) tp = tu - tq
57  IF(mq.GT.mur) tdp = tu
58  IF(mq.LT.mur.OR.mq.EQ.mur) tdp = tq
59  IF(mq.GT.md) tpd = tq - td
60  IF(mq.LT.md.OR.mq.EQ.md) tpd = td - tq
61  IF(mq.GT.md) tdpd = td
62  IF(mq.LT.md.OR.mq.EQ.md) tdpd = tq
63 
64  IF(mq.GT.md) dlambda1 = 6d0/96d0/pi**2*g1**2*hd**2*tpd
65  IF(mq.LT.md.OR.mq.EQ.md) dlambda1 = 3d0/32d0/pi**2*
66  * hd**2*(g1**2/3d0+g2**2)*tpd
67 
68  IF(mq.GT.mur) dlambda2 =12d0/96d0/pi**2*g1**2*hu**2*tp
69  IF(mq.LT.mur.OR.mq.EQ.mur) dlambda2 = 3d0/32d0/pi**2*
70  * hu**2*(-g1**2/3d0+g2**2)*tp
71 
72 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
73 C
74 C DLAMBDAP1 AND DLAMBDAP2 ARE THE NEW LOG CORRECTIONS DUE TO
75 C THE PRESENCE OF THE GLUINO MASS. THEY ARE IN GENERAL VERY SMALL,
76 C AND ONLY PRESENT IF THERE IS A HIERARCHY OF MASSES BETWEEN THE
77 C TWO STOPS.
78 C
79 C
80 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
81 
82  dlambdap2 = 0d0
83  IF(mglu.LT.mur.OR.mglu.LT.mq) THEN
84  IF(mq.GT.mur.AND.mglu.GT.mur) THEN
85  dlambdap2 = -4d0/(16d0*pi**2)**2*hu**4*(tq**2-tglu**2)
86  ENDIF
87 
88  IF(mq.GT.mur.AND.mglu.LT.mur) THEN
89  dlambdap2 = -4d0/(16d0*pi**2)**2*hu**4*(tq**2-tu**2)
90  ENDIF
91 
92  IF(mq.GT.mur.AND.mglu.EQ.mur) THEN
93  dlambdap2 = -4d0/(16d0*pi**2)**2*hu**4*(tq**2-tu**2)
94  ENDIF
95 
96  IF(mur.GT.mq.AND.mglu.GT.mq) THEN
97  dlambdap2 = -4d0/(16d0*pi**2)**2*hu**4*(tu**2-tglu**2)
98  ENDIF
99 
100  IF(mur.GT.mq.AND.mglu.LT.mq) THEN
101  dlambdap2 = -4d0/(16d0*pi**2)**2*hu**4*(tu**2-tq**2)
102  ENDIF
103 
104  IF(mur.GT.mq.AND.mglu.EQ.mq) THEN
105  dlambdap2 = -4d0/(16d0*pi**2)**2*hu**4*(tu**2-tq**2)
106  ENDIF
107  ENDIF
108  dlambda3 = 0d0
109  dlambda4 = 0d0
110  IF(mq.GT.md) dlambda3 = -1d0/32d0/pi**2*g1**2*hd**2*tpd
111  IF(mq.LT.md.OR.mq.EQ.md) dlambda3 = 3d0/64d0/pi**2*hd**2*
112  *(g2**2-g1**2/3d0)*tpd
113  IF(mq.GT.mur) dlambda3 = dlambda3 -
114  *1d0/16d0/pi**2*g1**2*hu**2*tp
115  IF(mq.LT.mur.OR.mq.EQ.mur) dlambda3 = dlambda3 +
116  * 3d0/64d0/pi**2*hu**2*(g2**2+g1**2/3d0)*tp
117  IF(mq.LT.mur) dlambda4 = -3d0/32d0/pi**2*g2**2*hu**2*tp
118  IF(mq.LT.md) dlambda4 = dlambda4 - 3d0/32d0/pi**2*g2**2*
119  *hd**2*tpd
120  lambda1 = ((g1**2 + g2**2)/4d0)*
121  * (1d0-3d0*hd**2*(tpd + tdpd)/8d0/pi**2)
122  *+(3d0*hd**4d0/16d0/pi**2) *tpd*(1d0
123  *+ (3d0*hd**2/2d0 + hu**2/2d0
124  *- 8d0*g3**2) * (tpd + 2d0*tdpd)/16d0/pi**2)
125  *+(3d0*hd**4d0/8d0/pi**2) *tdpd*(1d0 + (3d0*hd**2/2d0 + hu**2/2d0
126  *- 8d0*g3**2) * tdpd/16d0/pi**2) + dlambda1
127  lambda2 = ((g1**2 + g2**2)/4d0)*(1d0-3d0*hu**2*
128  *(tp + tdp)/8d0/pi**2)
129  *+(3d0*hu**4d0/16d0/pi**2) *tp*(1d0
130  *+ (3d0*hu**2/2d0 + hd**2/2d0
131  *- 8d0*g3**2) * (tp + 2d0*tdp)/16d0/pi**2)
132  *+(3d0*hu**4d0/8d0/pi**2) *tdp*(1d0 + (3d0*hu**2/2d0 + hd**2/2d0
133  *- 8d0*g3**2) * tdp/16d0/pi**2) + dlambda2 + dlambdap2
134  lambda3 = ((g2**2 - g1**2)/4d0)*(1d0-3d0*
135  *(hu**2)*(tp + tdp)/16d0/pi**2 -3d0*
136  *(hd**2)*(tpd + tdpd)/16d0/pi**2) +dlambda3
137  lambda4 = (- g2**2/2d0)*(1d0
138  *-3d0*(hu**2)*(tp + tdp)/16d0/pi**2
139  *-3d0*(hd**2)*(tpd + tdpd)/16d0/pi**2) +dlambda4
140 
141  lambda5 = 0d0
142  lambda6 = 0d0
143  lambda7 = 0d0
144 
145  m2(1,1) = 2d0*v**2*(lambda1*cosb**2+2d0*lambda6*
146  *cosb*sinb + lambda5*sinb**2) + ma**2*sinb**2
147 
148  m2(2,2) = 2d0*v**2*(lambda5*cosb**2+2d0*lambda7*
149  *cosb*sinb + lambda2*sinb**2) + ma**2*cosb**2
150  m2(1,2) = 2d0*v**2*(lambda6*cosb**2+(lambda3+lambda4)*
151  *cosb*sinb + lambda7*sinb**2) - ma**2*sinb*cosb
152 
153  m2(2,1) = m2(1,2)
154 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
155 CCC THIS IS THE CONTRIBUTION FROM LIGHT CHARGINOS/NEUTRALINOS
156 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
157 
158  mssusy=dsqrt(.5d0*(mq**2+mur**2)+mtop**2)
159 
160  IF(mchi.GT.mssusy) goto 100
161  IF(mchi.LT.mtop) mchi=mtop
162 
163  tchar=log(mssusy**2/mchi**2)
164 
165  deltal12=(9d0/64d0/pi**2*g2**4+5d0/192d0/pi**2*g1**4)*tchar
166  deltal3p4=(3d0/64d0/pi**2*g2**4+7d0/192d0/pi**2*g1**4
167  *+4d0/32d0/pi**2*g1**2*g2**2)*tchar
168 
169  deltam112=2d0*deltal12*v**2*cosb**2
170  deltam222=2d0*deltal12*v**2*sinb**2
171  deltam122=2d0*deltal3p4*v**2*sinb*cosb
172 
173  m2(1,1)=m2(1,1)+deltam112
174  m2(2,2)=m2(2,2)+deltam222
175  m2(1,2)=m2(1,2)+deltam122
176  m2(2,1)=m2(2,1)+deltam122
177 
178  100 CONTINUE
179 
180 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
181 CCC END OF CHARGINOS/NEUTRALINOS
182 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
183 
184  DO 120 i = 1,2
185  DO 110 j = 1,2
186  m2p(i,j) = m2(i,j) + vh(i,j)
187  110 CONTINUE
188  120 CONTINUE
189  trm2p = m2p(1,1) + m2p(2,2)
190  detm2p = m2p(1,1)*m2p(2,2) - m2p(1,2)*m2p(2,1)
191  mh2p = (trm2p - dsqrt(trm2p**2 - 4d0* detm2p))/2d0
192  hm2p = (trm2p + dsqrt(trm2p**2 - 4d0* detm2p))/2d0
193  hmp = dsqrt(hm2p)
194  mch2=ma**2+(lambda5-lambda4)*v**2
195  mch=dsqrt(mch2)
196  IF(mh2p.LT.0.) goto 130
197  mhp = sqrt(mh2p)
198  sin2alpha = 2d0*m2p(1,2)/sqrt(trm2p**2-4d0*detm2p)
199  cos2alpha = (m2p(1,1)-m2p(2,2))/sqrt(trm2p**2-4d0*detm2p)
200  IF(cos2alpha.GE.0.) THEN
201  alpha = asin(sin2alpha)/2d0
202  ELSE
203  alpha = -pi/2d0-asin(sin2alpha)/2d0
204  ENDIF
205  sa = sin(alpha)
206  ca = cos(alpha)
207 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
208 C
209 C HERE THE VALUES OF SAB AND CAB ARE DEFINED, IN ORDER
210 C TO DEFINE THE NEW COUPLINGS OF THE LIGHTEST AND
211 C HEAVY CP-EVEN HIGGS TO THE BOTTOM QUARK.
212 C
213 C
214 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
215  sab = sa*(1d0-deltamb/(1d0+deltamb)*(1d0+ca/sa/tanb))
216  cab = ca*(1d0-deltamb/(1d0+deltamb)*(1d0-sa/ca/tanb))
217  130 CONTINUE
218  RETURN
219  END