hijini.f

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

C            
C
C
        SUBROUTINE hijini
        common/hiparnt/hipr1(100),ihpr2(50),hint1(100),ihnt2(50)
        SAVE  /hiparnt/
        common/histrng/nfp(300,15),pp(300,15),nft(300,15),pt(300,15)
        SAVE  /histrng/
        common/hijjet1/npj(300),kfpj(300,500),pjpx(300,500),
     &                pjpy(300,500),pjpz(300,500),pjpe(300,500),
     &                pjpm(300,500),ntj(300),kftj(300,500),
     &                pjtx(300,500),pjty(300,500),pjtz(300,500),
     &                pjte(300,500),pjtm(300,500)
        SAVE  /hijjet1/
        common/hijjet2/nsg,njsg(900),iasg(900,3),k1sg(900,100),
     &          k2sg(900,100),pxsg(900,100),pysg(900,100),
     &          pzsg(900,100),pesg(900,100),pmsg(900,100)
        SAVE  /hijjet2/
        common/hijjet4/ndr,iadr(900,2),kfdr(900),pdr(900,5)
        SAVE  /hijjet4/
        common/ranseed/nseed
        SAVE  /ranseed/
C****************Reset the momentum of initial particles************
C             and assign flavors to the proj and targ string       *
C*******************************************************************
        nsg=0
        ndr=0
        ipp=2212
        ipt=2212
        IF(ihnt2(5).NE.0) ipp=ihnt2(5)
        IF(ihnt2(6).NE.0) ipt=ihnt2(6)
C               ********in case the proj or targ is a hadron.
C
        DO 100 i=1,ihnt2(1)
        pp(i,1)=0.0
        pp(i,2)=0.0
        pp(i,3)=sqrt(hint1(1)**2/4.0-hint1(8)**2)
        pp(i,4)=hint1(1)/2
        pp(i,5)=hint1(8)
        pp(i,6)=0.0
        pp(i,7)=0.0
        pp(i,8)=0.0
        pp(i,9)=0.0
        pp(i,10)=0.0
        nfp(i,3)=ipp
        nfp(i,4)=ipp
        nfp(i,5)=0
        nfp(i,6)=0
        nfp(i,7)=0
        nfp(i,8)=0
        nfp(i,9)=0
        nfp(i,10)=0
        nfp(i,11)=0
        npj(i)=0
        IF(i.GT.abs(ihnt2(2))) nfp(i,3)=2112
        CALL attflv(nfp(i,3),idq,idqq)
        nfp(i,1)=idq
        nfp(i,2)=idqq
        nfp(i,15)=-1
        IF(abs(idq).GT.1000.OR.(abs(idq*idqq).LT.100.AND.
     &          atl_ran(nseed).LT.0.5)) nfp(i,15)=1
        pp(i,14)=ulmass(idq)
        pp(i,15)=ulmass(idqq)
100     CONTINUE
C
        DO 200 i=1,ihnt2(3)
        pt(i,1)=0.0
        pt(i,2)=0.0
        pt(i,3)=-sqrt(hint1(1)**2/4.0-hint1(9)**2)
        pt(i,4)=hint1(1)/2.0
        pt(i,5)=hint1(9)
        pt(i,6)=0.0
        pt(i,7)=0.0
        pt(i,8)=0.0
        pt(i,9)=0.0
        pt(i,10)=0.0
        nft(i,3)=ipt
        nft(i,4)=ipt
        nft(i,5)=0
        nft(i,6)=0
        nft(i,7)=0
        nft(i,8)=0
        nft(i,9)=0
        nft(i,10)=0
        nft(i,11)=0
        ntj(i)=0
        IF(i.GT.abs(ihnt2(4))) nft(i,3)=2112
        CALL attflv(nft(i,3),idq,idqq)
        nft(i,1)=idq
        nft(i,2)=idqq
        nft(i,15)=1
        IF(abs(idq).GT.1000.OR.(abs(idq*idqq).LT.100.AND.
     &                  atl_ran(nseed).LT.0.5)) nft(i,15)=-1
        pt(i,14)=ulmass(idq)
        pt(i,15)=ulmass(idqq)
200     CONTINUE
        RETURN
        END