hiptdi.f

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

C
C
C
  SUBROUTINE hiptdi(PT,PTMAX,IOPT)
  common/hiparnt/hipr1(100),ihpr2(50),hint1(100),ihnt2(50)
  SAVE  /hiparnt/
  common/ranseed/nseed
  SAVE  /ranseed/
  IF(iopt.EQ.2) THEN
    pt=hirnd2(7,0.0,ptmax)
    IF(pt.GT.hipr1(8)) 
     &    pt=hipr1(2)*sqrt(-alog(exp(-hipr1(8)**2/hipr1(2)**2)
     &      -atl_ran(nseed)*(exp(-hipr1(8)**2/hipr1(2)**2)-
     &      exp(-ptmax**2/hipr1(2)**2))))

  ELSE
    pt=hipr1(2)*sqrt(-alog(1.0-atl_ran(nseed)*
     &      (1.0-exp(-ptmax**2/hipr1(2)**2))))
  ENDIF
  ptmax0=max(ptmax,0.01)
  pt=min(ptmax0-0.01,pt)
  RETURN
  END