]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C | |
5 | SUBROUTINE HIJINI | |
bc676b8e | 6 | #define BLANKET_SAVE |
e74335a4 | 7 | #include "hiparnt.inc" |
8 | #include "histrng.inc" | |
9 | #include "hijjet1.inc" | |
10 | #include "hijjet2.inc" | |
11 | #include "hijjet4.inc" | |
12 | SAVE | |
13 | C****************Reset the momentum of initial particles************ | |
14 | C and assign flavors to the proj and targ string * | |
15 | C******************************************************************* | |
16 | NSG=0 | |
17 | NDR=0 | |
18 | IPP=2212 | |
19 | IPT=2212 | |
20 | IF(IHNT2(5).NE.0) IPP=IHNT2(5) | |
21 | IF(IHNT2(6).NE.0) IPT=IHNT2(6) | |
22 | C ********in case the proj or targ is a hadron. | |
23 | C | |
24 | DO 100 I=1,IHNT2(1) | |
25 | PP(I,1)=0.0 | |
26 | PP(I,2)=0.0 | |
27 | PP(I,3)=SQRT(HINT1(1)**2/4.0-HINT1(8)**2) | |
28 | PP(I,4)=HINT1(1)/2 | |
29 | PP(I,5)=HINT1(8) | |
30 | PP(I,6)=0.0 | |
31 | PP(I,7)=0.0 | |
32 | PP(I,8)=0.0 | |
33 | PP(I,9)=0.0 | |
34 | PP(I,10)=0.0 | |
35 | NFP(I,3)=IPP | |
36 | NFP(I,4)=IPP | |
37 | NFP(I,5)=0 | |
38 | NFP(I,6)=0 | |
39 | NFP(I,7)=0 | |
40 | NFP(I,8)=0 | |
41 | NFP(I,9)=0 | |
42 | NFP(I,10)=0 | |
43 | NFP(I,11)=0 | |
44 | NPJ(I)=0 | |
45 | IF(I.GT.ABS(IHNT2(2))) NFP(I,3)=2112 | |
46 | CALL ATTFLV(NFP(I,3),IDQ,IDQQ) | |
47 | NFP(I,1)=IDQ | |
48 | NFP(I,2)=IDQQ | |
49 | NFP(I,15)=-1 | |
50 | IF(ABS(IDQ).GT.1000.OR.(ABS(IDQ*IDQQ).LT.100.AND. | |
51 | & RLU_HIJING(0).LT.0.5)) NFP(I,15)=1 | |
52 | PP(I,14)=ULMASS_HIJING(IDQ) | |
53 | PP(I,15)=ULMASS_HIJING(IDQQ) | |
54 | 100 CONTINUE | |
55 | C | |
56 | DO 200 I=1,IHNT2(3) | |
57 | PT(I,1)=0.0 | |
58 | PT(I,2)=0.0 | |
59 | PT(I,3)=-SQRT(HINT1(1)**2/4.0-HINT1(9)**2) | |
60 | PT(I,4)=HINT1(1)/2.0 | |
61 | PT(I,5)=HINT1(9) | |
62 | PT(I,6)=0.0 | |
63 | PT(I,7)=0.0 | |
64 | PT(I,8)=0.0 | |
65 | PT(I,9)=0.0 | |
66 | PT(I,10)=0.0 | |
67 | NFT(I,3)=IPT | |
68 | NFT(I,4)=IPT | |
69 | NFT(I,5)=0 | |
70 | NFT(I,6)=0 | |
71 | NFT(I,7)=0 | |
72 | NFT(I,8)=0 | |
73 | NFT(I,9)=0 | |
74 | NFT(I,10)=0 | |
75 | NFT(I,11)=0 | |
76 | NTJ(I)=0 | |
77 | IF(I.GT.ABS(IHNT2(4))) NFT(I,3)=2112 | |
78 | CALL ATTFLV(NFT(I,3),IDQ,IDQQ) | |
79 | NFT(I,1)=IDQ | |
80 | NFT(I,2)=IDQQ | |
81 | NFT(I,15)=1 | |
82 | IF(ABS(IDQ).GT.1000.OR.(ABS(IDQ*IDQQ).LT.100.AND. | |
83 | & RLU_HIJING(0).LT.0.5)) NFT(I,15)=-1 | |
84 | PT(I,14)=ULMASS_HIJING(IDQ) | |
85 | PT(I,15)=ULMASS_HIJING(IDQQ) | |
86 | 200 CONTINUE | |
87 | RETURN | |
88 | END |