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