pythiaeRHIC/pythia/pyrvi1.f

0001
0002 C*********************************************************************
0003
0004 C...PYRVI1
0005 C...Function to integrate resonance contributions
0006
0007       FUNCTION PYRVI1(ID1,ID2,ID3)
0008
0009       IMPLICIT NONE
0010       DOUBLE PRECISION LO,HI,PYRVI1,PYRVG1,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 PYRVG1,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       PYRVI1 = 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 C...A->B and B->A for antisparticles
0028       A(1)   = AB(1+INTRES(IDR,3),INTRES(IDR,1),INTRES(IDR,2))
0029       B(1)   = AB(2-INTRES(IDR,3),INTRES(IDR,1),INTRES(IDR,2))
0030 C...Integration boundaries and mass flag
0031       LO     = (RM(1)+RM(2))**2
0032       HI     = (RM(0)-RM(3))**2
0033       MFLAG  = DCMASS
0034       PYRVI1 = PYGAUS(PYRVG1,LO,HI,1D-3)
0035       RETURN
0036       END