pythiaeRHIC/pythia/pyrvi3.f

0001
0002 C*********************************************************************
0003
0004 C...PYRVI3
0005 C...Function to integrate true interference contributions
0006
0007       FUNCTION PYRVI3(ID1,ID2,ID3)
0008
0009       IMPLICIT NONE
0010       DOUBLE PRECISION LO,HI,PYRVI3, PYRVG3, PYGAUS
0011       DOUBLE PRECISION RES, AB, RM, RESM, RESW, A, B, RMS
0012       INTEGER ID1,ID2,ID3, IDR, IDR2, KFR, INTRES
0013       LOGICAL MFLAG,DCMASS
0014       EXTERNAL PYRVG3,PYGAUS
0015       COMMON/PYRVNV/AB(2,16,2),RMS(0:3),RES(6,2),INTRES(6,3),IDR,IDR2
0016      &     ,DCMASS,KFR(3)
0017       COMMON/PYRVPM/RM(0:3),A(2),B(2),RESM(2),RESW(2),MFLAG
0018       SAVE/PYRVNV/,/PYRVPM/
0019 C...Initialize mass and width information
0020       PYRVI3 = 0D0
0021       RM(0)  = RMS(0)
0022       RM(1)  = RMS(ID1)
0023       RM(2)  = RMS(ID2)
0024       RM(3)  = RMS(ID3)
0025       RESM(1)= RES(IDR,1)
0026       RESW(1)= RES(IDR,2)
0027       RESM(2)= RES(IDR2,1)
0028       RESW(2)= RES(IDR2,2)
0029 C...A -> B and B -> A for antisparticles
0030       A(1)   = AB(1+INTRES(IDR,3),INTRES(IDR,1),INTRES(IDR,2))
0031       B(1)   = AB(2-INTRES(IDR,3),INTRES(IDR,1),INTRES(IDR,2))
0032       A(2)   = AB(1+INTRES(IDR2,3),INTRES(IDR2,1),INTRES(IDR2,2))
0033       B(2)   = AB(2-INTRES(IDR2,3),INTRES(IDR2,1),INTRES(IDR2,2))
0034 C...Boundaries and mass flag
0035       LO     = (RM(1)+RM(2))**2
0036       HI     = (RM(0)-RM(3))**2
0037       MFLAG  = DCMASS
0038       PYRVI3 = PYGAUS(PYRVG3,LO,HI,1D-3)
0039       RETURN
0040       END