pyrvi3.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pyrvi3.f

 
C*********************************************************************
 
C...PYRVI3
C...Function to integrate true interference contributions
 
      FUNCTION pyrvi3(ID1,ID2,ID3)
 
      IMPLICIT NONE
      DOUBLE PRECISION lo,hi,pyrvi3, pyrvg3, pygaus
      DOUBLE PRECISION res, ab, rm, resm, resw, a, b, rms
      INTEGER id1,id2,id3, idr, idr2, kfr, intres
      LOGICAL mflag,dcmass
      EXTERNAL pyrvg3,pygaus
      common/pyrvnv/ab(2,16,2),rms(0:3),res(6,2),intres(6,3),idr,idr2
     &     ,dcmass,kfr(3)
      common/pyrvpm/rm(0:3),a(2),b(2),resm(2),resw(2),mflag
      SAVE/pyrvnv/,/pyrvpm/
C...Initialize mass and width information
      pyrvi3 = 0d0
      rm(0)  = rms(0)
      rm(1)  = rms(id1)
      rm(2)  = rms(id2)
      rm(3)  = rms(id3)
      resm(1)= res(idr,1)
      resw(1)= res(idr,2)
      resm(2)= res(idr2,1)
      resw(2)= res(idr2,2)
C...A -> B and B -> A for antisparticles
      a(1)   = ab(1+intres(idr,3),intres(idr,1),intres(idr,2))
      b(1)   = ab(2-intres(idr,3),intres(idr,1),intres(idr,2))
      a(2)   = ab(1+intres(idr2,3),intres(idr2,1),intres(idr2,2))
      b(2)   = ab(2-intres(idr2,3),intres(idr2,1),intres(idr2,2))
C...Boundaries and mass flag
      lo     = (rm(1)+rm(2))**2
      hi     = (rm(0)-rm(3))**2
      mflag  = dcmass
      pyrvi3 = pygaus(pyrvg3,lo,hi,1d-3)
      RETURN
      END