]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C | |
5 | C Modified for HIJING program | |
6 | c | |
7 | c modification July 22, 1997 In pyremnn put an upper limit | |
8 | c on the total pt kick the parton can accumulate via multiple | |
9 | C scattering. Set the upper limit to be the sqrt(s)/2, | |
10 | c this is fix cronin bug for Pb+Pb events at SPS energy. | |
11 | c | |
12 | C | |
13 | C Last modification Oct. 1993 to comply with non-vax | |
14 | C machines' compiler | |
15 | C | |
16 | C | |
17 | SUBROUTINE LU1ENT_HIJING(IP,KF,PE,THE,PHI) | |
18 | ||
19 | C...Purpose: to store one parton/particle in commonblock LUJETS. | |
20 | #include "lujets_hijing.inc" | |
21 | #include "ludat1_hijing.inc" | |
22 | #include "ludat2_hijing.inc" | |
23 | ||
24 | C...Standard checks. | |
25 | MSTU(28)=0 | |
26 | IF(MSTU(12).GE.1) CALL LULIST_HIJING(0) | |
27 | IPA=MAX(1,IABS(IP)) | |
28 | IF(IPA.GT.MSTU(4)) CALL LUERRM_HIJING(21, | |
29 | &'(LU1ENT_HIJING:) writing outside LUJETS_HIJING memory') | |
30 | KC=LUCOMP_HIJING(KF) | |
31 | IF(KC.EQ.0) CALL LUERRM_HIJING(12 | |
32 | $ ,'(LU1ENT_HIJING:) unknown flavour code') | |
33 | ||
34 | C...Find mass. Reset K, P and V vectors. | |
35 | PM=0. | |
36 | IF(MSTU(10).EQ.1) PM=P(IPA,5) | |
37 | IF(MSTU(10).GE.2) PM=ULMASS_HIJING(KF) | |
38 | DO 100 J=1,5 | |
39 | K(IPA,J)=0 | |
40 | P(IPA,J)=0. | |
41 | 100 V(IPA,J)=0. | |
42 | ||
43 | C...Store parton/particle in K and P vectors. | |
44 | K(IPA,1)=1 | |
45 | IF(IP.LT.0) K(IPA,1)=2 | |
46 | K(IPA,2)=KF | |
47 | P(IPA,5)=PM | |
48 | P(IPA,4)=MAX(PE,PM) | |
49 | PA=SQRT(P(IPA,4)**2-P(IPA,5)**2) | |
50 | P(IPA,1)=PA*SIN(THE)*COS(PHI) | |
51 | P(IPA,2)=PA*SIN(THE)*SIN(PHI) | |
52 | P(IPA,3)=PA*COS(THE) | |
53 | ||
54 | C...Set N. Optionally fragment/decay. | |
55 | N=IPA | |
56 | IF(IP.EQ.0) CALL LUEXEC_HIJING | |
57 | ||
58 | RETURN | |
59 | END |