pyrvi1.f

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

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