pyrvi2.f

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

 
C*********************************************************************
 
C...PYRVI2
C...Function to integrate L-R interference contributions
 
      FUNCTION pyrvi2(ID1,ID2,ID3)
 
      IMPLICIT NONE
      DOUBLE PRECISION lo,hi,pyrvi2, pyrvg2, pygaus
      DOUBLE PRECISION res, ab, rm, resm, resw, a, b, rms
      INTEGER id1,id2,id3, idr, idr2, kfr, intres
      LOGICAL mflag,dcmass
      EXTERNAL pyrvg2,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
      pyrvi2 = 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(idr+1,1)
      resw(2)= res(idr+1,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(idr+1,3),intres(idr+1,1),intres(idr+1,2))
      b(2)   = ab(2-intres(idr+1,3),intres(idr+1,1),intres(idr+1,2))
C...Boundaries and mass flag
      lo     = (rm(1)+rm(2))**2
      hi     = (rm(0)-rm(3))**2
      mflag  = dcmass
      pyrvi2 = pygaus(pyrvg2,lo,hi,1d-3)
      RETURN
      END