← Back to branch summary

~alifson/chiralityflow/trunk

« back to all changes in this revision

Viewing changes to vendor/IREGI/src/qcdloop/ql/qlbox1.f

  • Committer: andrew.lifson at lu
  • Date: 2021-09-01 15:34:39 UTC
  • Revision ID: andrew.lifson@thep.lu.se-20210901153439-7fasjhav4cp4m88r
testing a new repository of a madgraph folder

Show diffs side-by-side

added added

removed removed

Lines of Context:
vendor/IREGI/src/qcdloop/ql/qlbox1.f
 
1
      subroutine qlbox1(Y,musq,Ires)
 
2
c    $I_4^{D=4-2 \epsilon}(0,0,0,0;s_{12},s_{23};0,0,0,0)$}
 
3
c----%\cite{Bern:1993kr}
 
4
c----\bibitem{Bern:1993kr}
 
5
c----  Z.~Bern, L.~J.~Dixon and D.~A.~Kosower,
 
6
c----  %``Dimensionally regulated pentagon integrals,''
 
7
c----  Nucl.\ Phys.\ B {\bf 412}, 751 (1994)
 
8
c----  [arXiv:hep-ph/9306240].
 
9
c----  %%CITATION = HEP-PH 9306240;%%
 
10
c----  Eqn (I.11)
 
11
c----Cayley matrix
 
12
c                [                 s12        ]
 
13
c                [   0      0    - ---    0   ]
 
14
c                [                  2         ]
 
15
c                [                            ]
 
16
c                [                        s23 ]
 
17
c                [   0      0      0    - --- ]
 
18
c                [                         2  ]
 
19
c           Y1 = [                            ]
 
20
c                [   s12                      ]
 
21
c                [ - ---    0      0      0   ]
 
22
c                [    2                       ]
 
23
c                [                            ]
 
24
c                [          s23               ]
 
25
c                [   0    - ---    0      0   ]
 
26
c                [           2                ]
 
27
      implicit none
 
28
      include 'qlconstants.f'
 
29
      double precision musq,Y(4,4),si,ta
 
30
      double complex Ires(-2:0),qllnrat,fac
 
31
 
 
32
      si=two*Y(1,3)
 
33
      ta=two*Y(2,4)
 
34
      fac=dcmplx(1d0/(si*ta))     
 
35
      Ires(-2)=fac*ctwo*ctwo
 
36
      Ires(-1)=fac*ctwo*(
 
37
     . -qllnrat(ta,musq)-qllnrat(si,musq))
 
38
      Ires( 0)=fac*(qllnrat(ta,musq)**2+qllnrat(si,musq)**2
 
39
     . -qllnrat(ta,si)**2-dcmplx(pi**2))
 
40
      return
 
41
      end
 
42