From: morsch Date: Thu, 10 Jun 2010 13:57:41 +0000 (+0000) Subject: Correct ran function. X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=6082051fc9b7a31238e8d0ec5f2ad5919d7d026f;p=u%2Fmrichter%2FAliRoot.git Correct ran function. --- diff --git a/HIJING/hijing1_36/hijing.F b/HIJING/hijing1_36/hijing.F index 790d8e5e5a7..9df2e9215ac 100644 --- a/HIJING/hijing1_36/hijing.F +++ b/HIJING/hijing1_36/hijing.F @@ -189,9 +189,9 @@ C 5 R=HIRND(1) c if(IHNT2(1).EQ.2) then - rnd1=max(RAN0(NSEED),1.0e-20) - rnd2=max(RAN0(NSEED),1.0e-20) - rnd3=max(RAN0(NSEED),1.0e-20) + rnd1=max(RLU_HIJING(NSEED),1.0e-20) + rnd2=max(RLU_HIJING(NSEED),1.0e-20) + rnd3=max(RLU_HIJING(NSEED),1.0e-20) R=-0.5*(log(rnd1)*4.38/2.0+log(rnd2)*0.85/2.0 & +4.38*0.85*log(rnd3)/(4.38+0.85)) endif @@ -246,9 +246,9 @@ C****************************** 15 R=HIRND(2) c if(IHNT2(3).EQ.2) then - rnd1=max(RAN0(NSEED),1.0e-20) - rnd2=max(RAN0(NSEED),1.0e-20) - rnd3=max(RAN0(NSEED),1.0e-20) + rnd1=max(RLU_HIJING(NSEED),1.0e-20) + rnd2=max(RLU_HIJING(NSEED),1.0e-20) + rnd3=max(RLU_HIJING(NSEED),1.0e-20) R=-0.5*(log(rnd1)*4.38/2.0+log(rnd2)*0.85/2.0 & +4.38*0.85*log(rnd3)/(4.38+0.85)) endif @@ -645,7 +645,7 @@ C ****** identify the mother particle V3MIN1=RARB/GAMMA0 V3MIN2=1.0/MAX(1.0,5.08*ABS(PATT(I,3))) VATT_MIN=MAX(V3MIN1,V3MIN2) - VATT03=VATT03+(0.5-RLU(0))*VATT_MIN + VATT03=VATT03+(0.5-RLU_HIJING(0))*VATT_MIN amt2=P(I,1)**2+P(I,2)**2+P(I,5)**2 IF(amt2.GT.0.0) THEN tauf=0.2*2.0*P(I,3)/amt2 @@ -733,7 +733,7 @@ C of the produced hadrons V3MIN1=RARB/GAMMA0 V3MIN2=1.0/MAX(1.0,5.08*ABS(PATT(I,3))) VATT_MIN=MAX(V3MIN1,V3MIN2) - VATT03=VATT03+(0.5-RLU(0))*VATT_MIN + VATT03=VATT03+(0.5-RLU_HIJING(0))*VATT_MIN amt2=P(I,1)**2+P(I,2)**2+P(I,5)**2 IF(amt2.GT.0.0) THEN tauf=0.2*2.0*P(I,3)/amt2