]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | ||
3 | C********************************************************************* | |
4 | ||
5 | SUBROUTINE PYKMAP_HIJING(IVAR,MVAR,VVAR) | |
6 | ||
7 | C...Maps a uniform distribution into a distribution of a kinematical | |
8 | C...variable according to one of the possibilities allowed. It is | |
9 | C...assumed that kinematical limits have been set by a PYKLIM call. | |
10 | #include "ludat2_hijing.inc" | |
11 | #include "pyint1_hijing.inc" | |
12 | #include "pyint2_hijing.inc" | |
13 | ||
14 | C...Convert VVAR to tau variable. | |
15 | ISUB=MINT(1) | |
16 | IF(IVAR.EQ.1) THEN | |
17 | TAUMIN=VINT(11) | |
18 | TAUMAX=VINT(31) | |
19 | IF(MVAR.EQ.3.OR.MVAR.EQ.4) THEN | |
20 | TAURE=VINT(73) | |
21 | GAMRE=VINT(74) | |
22 | ELSEIF(MVAR.EQ.5.OR.MVAR.EQ.6) THEN | |
23 | TAURE=VINT(75) | |
24 | GAMRE=VINT(76) | |
25 | ENDIF | |
26 | IF(MINT(43).EQ.1.AND.(ISET(ISUB).EQ.1.OR.ISET(ISUB).EQ.2)) THEN | |
27 | TAU=1. | |
28 | ELSEIF(MVAR.EQ.1) THEN | |
29 | TAU=TAUMIN*(TAUMAX/TAUMIN)**VVAR | |
30 | ELSEIF(MVAR.EQ.2) THEN | |
31 | TAU=TAUMAX*TAUMIN/(TAUMIN+(TAUMAX-TAUMIN)*VVAR) | |
32 | ELSEIF(MVAR.EQ.3.OR.MVAR.EQ.5) THEN | |
33 | RATGEN=(TAURE+TAUMAX)/(TAURE+TAUMIN)*TAUMIN/TAUMAX | |
34 | TAU=TAURE*TAUMIN/((TAURE+TAUMIN)*RATGEN**VVAR-TAUMIN) | |
35 | ELSE | |
36 | AUPP=ATAN((TAUMAX-TAURE)/GAMRE) | |
37 | ALOW=ATAN((TAUMIN-TAURE)/GAMRE) | |
38 | TAU=TAURE+GAMRE*TAN(ALOW+(AUPP-ALOW)*VVAR) | |
39 | ENDIF | |
40 | VINT(21)=MIN(TAUMAX,MAX(TAUMIN,TAU)) | |
41 | ||
42 | C...Convert VVAR to y* variable. | |
43 | ELSEIF(IVAR.EQ.2) THEN | |
44 | YSTMIN=VINT(12) | |
45 | YSTMAX=VINT(32) | |
46 | IF(MINT(43).EQ.1) THEN | |
47 | YST=0. | |
48 | ELSEIF(MINT(43).EQ.2) THEN | |
49 | IF(ISET(ISUB).LE.2) YST=-0.5*LOG(VINT(21)) | |
50 | IF(ISET(ISUB).GE.3) YST=-0.5*LOG(VINT(26)) | |
51 | ELSEIF(MINT(43).EQ.3) THEN | |
52 | IF(ISET(ISUB).LE.2) YST=0.5*LOG(VINT(21)) | |
53 | IF(ISET(ISUB).GE.3) YST=0.5*LOG(VINT(26)) | |
54 | ELSEIF(MVAR.EQ.1) THEN | |
55 | YST=YSTMIN+(YSTMAX-YSTMIN)*SQRT(VVAR) | |
56 | ELSEIF(MVAR.EQ.2) THEN | |
57 | YST=YSTMAX-(YSTMAX-YSTMIN)*SQRT(1.-VVAR) | |
58 | ELSE | |
59 | AUPP=ATAN(EXP(YSTMAX)) | |
60 | ALOW=ATAN(EXP(YSTMIN)) | |
61 | YST=LOG(TAN(ALOW+(AUPP-ALOW)*VVAR)) | |
62 | ENDIF | |
63 | VINT(22)=MIN(YSTMAX,MAX(YSTMIN,YST)) | |
64 | ||
65 | C...Convert VVAR to cos(theta-hat) variable. | |
66 | ELSEIF(IVAR.EQ.3) THEN | |
67 | RM34=2.*VINT(63)*VINT(64)/(VINT(21)*VINT(2))**2 | |
68 | RSQM=1.+RM34 | |
69 | IF(2.*VINT(71)**2/(VINT(21)*VINT(2)).LT.0.0001) RM34=MAX(RM34, | |
70 | & 2.*VINT(71)**2/(VINT(21)*VINT(2))) | |
71 | CTNMIN=VINT(13) | |
72 | CTNMAX=VINT(33) | |
73 | CTPMIN=VINT(14) | |
74 | CTPMAX=VINT(34) | |
75 | IF(MVAR.EQ.1) THEN | |
76 | ANEG=CTNMAX-CTNMIN | |
77 | APOS=CTPMAX-CTPMIN | |
78 | IF(ANEG.GT.0..AND.VVAR*(ANEG+APOS).LE.ANEG) THEN | |
79 | VCTN=VVAR*(ANEG+APOS)/ANEG | |
80 | CTH=CTNMIN+(CTNMAX-CTNMIN)*VCTN | |
81 | ELSE | |
82 | VCTP=(VVAR*(ANEG+APOS)-ANEG)/APOS | |
83 | CTH=CTPMIN+(CTPMAX-CTPMIN)*VCTP | |
84 | ENDIF | |
85 | ELSEIF(MVAR.EQ.2) THEN | |
86 | RMNMIN=MAX(RM34,RSQM-CTNMIN) | |
87 | RMNMAX=MAX(RM34,RSQM-CTNMAX) | |
88 | RMPMIN=MAX(RM34,RSQM-CTPMIN) | |
89 | RMPMAX=MAX(RM34,RSQM-CTPMAX) | |
90 | ANEG=LOG(RMNMIN/RMNMAX) | |
91 | APOS=LOG(RMPMIN/RMPMAX) | |
92 | IF(ANEG.GT.0..AND.VVAR*(ANEG+APOS).LE.ANEG) THEN | |
93 | VCTN=VVAR*(ANEG+APOS)/ANEG | |
94 | CTH=RSQM-RMNMIN*(RMNMAX/RMNMIN)**VCTN | |
95 | ELSE | |
96 | VCTP=(VVAR*(ANEG+APOS)-ANEG)/APOS | |
97 | CTH=RSQM-RMPMIN*(RMPMAX/RMPMIN)**VCTP | |
98 | ENDIF | |
99 | ELSEIF(MVAR.EQ.3) THEN | |
100 | RMNMIN=MAX(RM34,RSQM+CTNMIN) | |
101 | RMNMAX=MAX(RM34,RSQM+CTNMAX) | |
102 | RMPMIN=MAX(RM34,RSQM+CTPMIN) | |
103 | RMPMAX=MAX(RM34,RSQM+CTPMAX) | |
104 | ANEG=LOG(RMNMAX/RMNMIN) | |
105 | APOS=LOG(RMPMAX/RMPMIN) | |
106 | IF(ANEG.GT.0..AND.VVAR*(ANEG+APOS).LE.ANEG) THEN | |
107 | VCTN=VVAR*(ANEG+APOS)/ANEG | |
108 | CTH=RMNMIN*(RMNMAX/RMNMIN)**VCTN-RSQM | |
109 | ELSE | |
110 | VCTP=(VVAR*(ANEG+APOS)-ANEG)/APOS | |
111 | CTH=RMPMIN*(RMPMAX/RMPMIN)**VCTP-RSQM | |
112 | ENDIF | |
113 | ELSEIF(MVAR.EQ.4) THEN | |
114 | RMNMIN=MAX(RM34,RSQM-CTNMIN) | |
115 | RMNMAX=MAX(RM34,RSQM-CTNMAX) | |
116 | RMPMIN=MAX(RM34,RSQM-CTPMIN) | |
117 | RMPMAX=MAX(RM34,RSQM-CTPMAX) | |
118 | ANEG=1./RMNMAX-1./RMNMIN | |
119 | APOS=1./RMPMAX-1./RMPMIN | |
120 | IF(ANEG.GT.0..AND.VVAR*(ANEG+APOS).LE.ANEG) THEN | |
121 | VCTN=VVAR*(ANEG+APOS)/ANEG | |
122 | CTH=RSQM-1./(1./RMNMIN+ANEG*VCTN) | |
123 | ELSE | |
124 | VCTP=(VVAR*(ANEG+APOS)-ANEG)/APOS | |
125 | CTH=RSQM-1./(1./RMPMIN+APOS*VCTP) | |
126 | ENDIF | |
127 | ELSEIF(MVAR.EQ.5) THEN | |
128 | RMNMIN=MAX(RM34,RSQM+CTNMIN) | |
129 | RMNMAX=MAX(RM34,RSQM+CTNMAX) | |
130 | RMPMIN=MAX(RM34,RSQM+CTPMIN) | |
131 | RMPMAX=MAX(RM34,RSQM+CTPMAX) | |
132 | ANEG=1./RMNMIN-1./RMNMAX | |
133 | APOS=1./RMPMIN-1./RMPMAX | |
134 | IF(ANEG.GT.0..AND.VVAR*(ANEG+APOS).LE.ANEG) THEN | |
135 | VCTN=VVAR*(ANEG+APOS)/ANEG | |
136 | CTH=1./(1./RMNMIN-ANEG*VCTN)-RSQM | |
137 | ELSE | |
138 | VCTP=(VVAR*(ANEG+APOS)-ANEG)/APOS | |
139 | CTH=1./(1./RMPMIN-APOS*VCTP)-RSQM | |
140 | ENDIF | |
141 | ENDIF | |
142 | IF(CTH.LT.0.) CTH=MIN(CTNMAX,MAX(CTNMIN,CTH)) | |
143 | IF(CTH.GT.0.) CTH=MIN(CTPMAX,MAX(CTPMIN,CTH)) | |
144 | VINT(23)=CTH | |
145 | ||
146 | C...Convert VVAR to tau' variable. | |
147 | ELSEIF(IVAR.EQ.4) THEN | |
148 | TAU=VINT(11) | |
149 | TAUPMN=VINT(16) | |
150 | TAUPMX=VINT(36) | |
151 | IF(MINT(43).EQ.1) THEN | |
152 | TAUP=1. | |
153 | ELSEIF(MVAR.EQ.1) THEN | |
154 | TAUP=TAUPMN*(TAUPMX/TAUPMN)**VVAR | |
155 | ELSE | |
156 | AUPP=(1.-TAU/TAUPMX)**4 | |
157 | ALOW=(1.-TAU/TAUPMN)**4 | |
158 | TAUP=TAU/(1.-(ALOW+(AUPP-ALOW)*VVAR)**0.25) | |
159 | ENDIF | |
160 | VINT(26)=MIN(TAUPMX,MAX(TAUPMN,TAUP)) | |
161 | ENDIF | |
162 | ||
163 | RETURN | |
164 | END |