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