COSMOS v7.655  COSMOSv7655
(AirShowerMC)
codeConv.f
Go to the documentation of this file.
1 ! convert cosmos code once into kf code
2 ! and kf code is converted into target code
3  integer code, subcode, charge, kf, ccode
4  character*3 inc, outc
5 
6  integer idtrafo
7  write(0,*) ' cosmos code -> kf code -> target code'
8  write(0,*) ' '
9  write(0,*) 'Enter target code by selecting from'
10  write(0,*) " pdg qgs nxs ghe sib cor"
11  write(0,*) " here epos =nxs cor=Corsika=Geant3"
12 !.....'nxs' = epos
13 !.....'pdg' = PDG 1996
14 !.....'qgs' = QGSJet
15 !.....'ghe' = Gheisha
16 !.....'sib' = Sibyll
17 !.....'cor' = Corsika (GEANT)
18  inc = 'pdg'
19  read(*, *) outc
20  write(0,*)'cos--> ', inc, "-->",outc
21  do while (.true.)
22  write(0,*) 'Enter code, subcode, charge'
23  read(*,*,end=100) code, subcode, charge
24  write(0,*) code, subcode, charge
25  call ccos2kf(code, subcode, charge, kf)
26  write(0,*) ' kf =',kf
27  ccode=idtrafo(inc, outc, kf)
28  write(0,*) ' target code =', ccode
29  write(0,*)
30  enddo
31  100 continue
32  end
33 
34 !------------------------------------------------------------------------------
35  integer function idtrafo(code1,code2,idi)
36 !------------------------------------------------------------------------------
37 !.....tranforms id of code1 (=idi) into id of code2 (=idtrafocx)
38 !.....supported codes:
39 !.....'nxs' = epos
40 !.....'pdg' = PDG 1996
41 !.....'qgs' = QGSJet
42 !.....'ghe' = Gheisha
43 !.....'sib' = Sibyll
44 !.....'cor' = Corsika (GEANT)
45 
46 ! --- ighenex(I)=EPOS CODE CORRESPONDING TO GHEISHA CODE I ---
47  common /ighnx/ ighenex(35)
48  data ighenex/
49  $ 10, 11, -12, 12, -14, 14, 120, 110,
50  $ -120, 130, 20, -20, -130, 1120, -1120, 1220,
51  $ -1220, 2130, -2130, 1130, 1230, 2230, -1130, -1230,
52  $ -2230, 1330, 2330, -1330, -2330, 17, 18, 19,
53  $ 3331, -3331, 30/
54 
55 ! --- DATA STMTS. FOR GEANT/GHEISHA PARTICLE CODE CONVERSIONS ---
56 ! --- KIPART(I)=GHEISHA CODE CORRESPONDING TO GEANT CODE I ---
57 ! --- IKPART(I)=GEANT CODE CORRESPONDING TO GHEISHA CODE I ---
58  dimension kipart(48),ikpart(35)
59  DATA kipart/
60  $ 1, 3, 4, 2, 5, 6, 8, 7,
61  $ 9, 12, 10, 13, 16, 14, 15, 11,
62  $ 35, 18, 20, 21, 22, 26, 27, 33,
63  $ 17, 19, 23, 24, 25, 28, 29, 34,
64  $ 35, 35, 35, 35, 35, 35, 35, 35,
65  $ 35, 35, 35, 35, 30, 31, 32, 35/
66 
67  DATA ikpart/
68  $ 1, 4, 2, 3, 5, 6, 8, 7,
69  $ 9, 11, 16, 10, 12, 14, 15, 13,
70  $ 25, 18, 26, 19, 20, 21, 27, 28,
71  $ 29, 22, 23, 30, 31, 45, 46, 47,
72  $ 24, 32, 48/
73 !-------------------------------------------------------------------------------
74 
75  character*3 code1,code2
76  parameter(ncode=5,nidt=334)
77  integer idt(ncode,nidt)
78  double precision drangen,dummy
79 
80  data ((idt(i,j),i=1,ncode),j= 1,18)/
81  * 1,2,99,99,99 !u quark
82  * , 2,1,99,99,99 !d
83  * , 3,3,99,99,99 !s
84  * , 4,4,99,99,99 !c
85  * , 5,5,99,99,99 !b
86  * , 6,6,99,99,99 !t
87  * , 10,22,99,1,1 !gamma
88  * , 9 ,21,99,99,99 !gluon
89  * , 12,11,11,4,3 !e-
90  * , -12,-11,-11,3,2 !e+
91  * , 11,12,99,2,15 !nu_e-
92  * , -11,-12,99,-2,16 !nu_e+
93  * , 14,13,99,6,5 !mu-
94  * , -14,-13,99,5,4 !mu+
95  * , 13,14,99,2,17 !nu_mu-
96  * , -13,-14,99,-2,18 !nu_mu+
97  * , 16,15,99,99,19 !tau-
98  * , 15,16,99,99,20 / !nu_tau-
99  data ((idt(i,j),i=1,ncode),j= 19,40)/
100  * 110,111,0,8,6 !pi0
101  * , 120,211,1,7,7 !pi+
102  * , -120,-211,-1,9,8 !pi-
103  * , 220,221,10,99,23 !eta
104  * , 130,321,4,10,9 !k+
105  * , -130,-321,-4,13,10 !k-
106  * , 230,311,5,11,21 !k0
107  * , -230,-311,-5,12,22 !k0b
108  * , 20,310,5,11,12 !kshort
109  * , -20,-310,-5,12,11 !klong
110  * , 330,331,99,99,24 !etaprime
111  * , 111,113,99,99,27 !rho0
112  * , 121,213,99,99,25 !rho+
113  * , -121,-213,99,99,26 !rho-
114  * , 221,223,99,99,32 !omega
115  * , 131,323,99,99,28 !k*+
116  * , -131,-323,99,99,29 !k*-
117  * , 231,313,99,99,30 !k*0
118  * , -231,-313,99,99,31 !k*0b
119  * , 331,333,99,99,33 !phi
120  $ , -140,421,8,99,99 !D0(1.864)
121  $ , 240,-411,7,99,99 / !D(1.869)-
122  data ((idt(i,j),i=1,ncode),j= 41,59)/
123  * 1120,2212,2,14,13 !proton
124  * , 1220,2112,3,16,14 !neutron
125  * , 2130,3122,6,18,39 !lambda
126  * , 1130,3222,99,20,34 !sigma+
127  * , 1230,3212,99,21,35 !sigma0
128  * , 2230,3112,99,22,36 !sigma-
129  * , 1330,3322,99,26,37 !xi0
130  * , 2330,3312,99,27,38 !xi-
131  * , 1111,2224,99,99,40 !delta++
132  * , 1121,2214,99,99,41 !delta+
133  * , 1221,2114,99,99,42 !delta0
134  * , 2221,1114,99,99,43 !delta-
135  * , 1131,3224,99,99,44 !sigma*+
136  * , 1231,3214,99,99,45 !sigma*0
137  * , 2231,3114,99,99,46 !sigma*-
138  * , 1331, 3324,99,99,47 !xi*0
139  * , 2331, 3314,99,99,48 !xi*-
140  * , 3331, 3334,99,24,49 !omega-
141  $ , 2140, 4122,9,99,99 / !LambdaC(2.285)+
142  data ((idt(i,j),i=1,ncode),j= 60,64)/
143  $ 17,99,99,30,1002 ! Deuteron
144  $ ,18,99,99,31,1003 ! Triton
145  $ ,19,99,99,32,1004 ! Alpha
146  $ ,0,99,99,0,0 ! Air
147  * ,99,99,99,99,99 / ! unknown
148  data ((idt(i,j),i=1,ncode),j= 65,79)/
149  $ 1112,32224,99,99,99 ! Delta(1600)++
150  $ , 1112, 2222,99,99,99 ! Delta(1620)++
151  $ , 1113,12224,99,99,99 ! Delta(1700)++
152  $ , 1114,12222,99,99,99 ! Delta(1900)++
153  $ , 1114, 2226,99,99,99 ! Delta(1905)++
154  $ , 1114,22222,99,99,99 ! Delta(1910)++
155  $ , 1114,22224,99,99,99 ! Delta(1920)++
156  $ , 1114,12226,99,99,99 ! Delta(1930)++
157  $ , 1114, 2228,99,99,99 ! Delta(1950)++
158  $ , 2222,31114,99,99,99 ! Delta(1600)-
159  $ , 2222, 1112,99,99,99 ! Delta(1620)-
160  $ , 2223,11114,99,99,99 ! Delta(1700)-
161  $ , 2224,11112,99,99,99 ! Delta(1900)-
162  $ , 2224, 1116,99,99,99 ! Delta(1905)-
163  $ , 2224,21112,99,99,99 / ! Delta(1910)-
164  data ((idt(i,j),i=1,ncode),j= 80,94)/
165  $ 2224,21114,99,99,99 ! Delta(1920)-
166  $ ,2224,11116,99,99,99 ! Delta(1930)-
167  $ ,2224, 1118,99,99,99 ! Delta(1950)-
168  $ ,1122,12212,99,99,99 ! N(1440)+
169  $ ,1123, 2124,99,99,99 ! N(1520)+
170  $ ,1123,22212,99,99,99 ! N(1535)+
171  $ ,1124,32214,99,99,99 ! Delta(1600)+
172  $ ,1124, 2122,99,99,99 ! Delta(1620)+
173  $ ,1125,32212,99,99,99 ! N(1650)+
174  $ ,1125, 2216,99,99,99 ! N(1675)+
175  $ ,1125,12216,99,99,99 ! N(1680)+
176  $ ,1126,12214,99,99,99 ! Delta(1700)+
177  $ ,1127,22124,99,99,99 ! N(1700)+
178  $ ,1127,42212,99,99,99 ! N(1710)+
179  $ ,1127,32124,99,99,99 / ! N(1720)+
180  data ((idt(i,j),i=1,ncode),j= 95,109)/
181  $ 1128,12122,99,99,99 ! Delta(1900)+
182  $ ,1128, 2126,99,99,99 ! Delta(1905)+
183  $ ,1128,22122,99,99,99 ! Delta(1910)+
184  $ ,1128,22214,99,99,99 ! Delta(1920)+
185  $ ,1128,12126,99,99,99 ! Delta(1930)+
186  $ ,1128, 2218,99,99,99 ! Delta(1950)+
187  $ ,1222,12112,99,99,99 ! N(1440)0
188  $ ,1223, 1214,99,99,99 ! N(1520)0
189  $ ,1223,22112,99,99,99 ! N(1535)0
190  $ ,1224,32114,99,99,99 ! Delta(1600)0
191  $ ,1224, 1212,99,99,99 ! Delta(1620)0
192  $ ,1225,32112,99,99,99 ! N(1650)0
193  $ ,1225, 2116,99,99,99 ! N(1675)0
194  $ ,1225,12116,99,99,99 ! N(1680)0
195  $ ,1226,12114,99,99,99 / ! Delta(1700)0
196  data ((idt(i,j),i=1,ncode),j= 110,124)/
197  $ 1227,21214,99,99,99 ! N(1700)0
198  $ ,1227,42112,99,99,99 ! N(1710)0
199  $ ,1227,31214,99,99,99 ! N(1720)0
200  $ ,1228,11212,99,99,99 ! Delta(1900)0
201  $ ,1228, 1216,99,99,99 ! Delta(1905)0
202  $ ,1228,21212,99,99,99 ! Delta(1910)0
203  $ ,1228,22114,99,99,99 ! Delta(1920)0
204  $ ,1228,11216,99,99,99 ! Delta(1930)0
205  $ ,1228, 2118,99,99,99 ! Delta(1950)0
206  $ ,1233,13122,99,99,99 ! Lambda(1405)0
207  $ ,1234, 3124,99,99,99 ! Lambda(1520)0
208  $ ,1235,23122,99,99,99 ! Lambda(1600)0
209  $ ,1235,33122,99,99,99 ! Lambda(1670)0
210  $ ,1235,13124,99,99,99 ! Lambda(1690)0
211  $ ,1236,13212,99,99,99 / ! Sigma(1660)0
212  data ((idt(i,j),i=1,ncode),j= 125,139)/
213  $ 1236,13214,99,99,99 ! Sigma(1670)0
214  $ ,1237,23212,99,99,99 ! Sigma(1750)0
215  $ ,1237, 3216,99,99,99 ! Sigma(1775)0
216  $ ,1238,43122,99,99,99 ! Lambda(1800)0
217  $ ,1238,53122,99,99,99 ! Lambda(1810)0
218  $ ,1238, 3126,99,99,99 ! Lambda(1820)0
219  $ ,1238,13126,99,99,99 ! Lambda(1830)0
220  $ ,1238,23124,99,99,99 ! Lambda(1890)0
221  $ ,1239,13216,99,99,99 ! Sigma(1915)0
222  $ ,1239,23214,99,99,99 ! Sigma(1940)0
223  $ ,1132,13222,99,99,99 ! Sigma(1660)+
224  $ ,1132,13224,99,99,99 ! Sigma(1670)+
225  $ ,1133,23222,99,99,99 ! Sigma(1750)+
226  $ ,1133,3226,99,99,99 ! Sigma(1775)+
227  $ ,1134,13226,99,99,99 / ! Sigma(1915)+
228  data ((idt(i,j),i=1,ncode),j= 140,146)/
229  $ 1134,23224,99,99,99 ! Sigma(1940)+
230  $ ,2232,13112,99,99,99 ! Sigma(1660)-
231  $ ,2232,13114,99,99,99 ! Sigma(1670)-
232  $ ,2233,23112,99,99,99 ! Sigma(1750)-
233  $ ,2233,3116,99,99,99 ! Sigma(1775)-
234  $ ,2234,13116,99,99,99 ! Sigma(1915)-
235  $ ,2234,23114,99,99,99 / ! Sigma(1940)-
236  data ((idt(i,j),i=1,ncode),j= 147,159)/
237  $ 5,7,99,99,99 ! quark b'
238  $ ,6,8,99,99,99 ! quark t'
239  $ ,16,17,99,99,99 ! lepton tau'
240  $ ,15,18,99,99,99 ! lepton nu' tau
241  $ ,90,23,99,99,99 ! Z0
242  $ ,80,24,99,99,99 ! W+
243  $ ,81,25,99,99,99 ! h0
244  $ ,85,32,99,99,99 ! Z'0
245  $ ,86,33,99,99,99 ! Z''0
246  $ ,87,34,99,99,99 ! W'+
247  $ ,82,35,99,99,99 ! H0
248  $ ,83,36,99,99,99 ! A0
249  $ ,84,37,99,99,99 / ! H+
250  data ((idt(i,j),i=1,ncode),j= 160,184)/
251  $ 1200,2101,99,99,99 ! diquark ud_0
252  $ ,2300,3101,99,99,99 ! diquark sd_0
253  $ ,1300,3201,99,99,99 ! diquark su_0
254  $ ,2400,4101,99,99,99 ! diquark cd_0
255  $ ,1400,4201,99,99,99 ! diquark cu_0
256  $ ,3400,4301,99,99,99 ! diquark cs_0
257  $ ,2500,5101,99,99,99 ! diquark bd_0
258  $ ,1500,5201,99,99,99 ! diquark bu_0
259  $ ,3500,5301,99,99,99 ! diquark bs_0
260  $ ,4500,5401,99,99,99 ! diquark bc_0
261  $ ,2200,1103,99,99,99 ! diquark dd_1
262  $ ,1200,2103,99,99,99 ! diquark ud_1
263  $ ,1100,2203,99,99,99 ! diquark uu_1
264  $ ,2300,3103,99,99,99 ! diquark sd_1
265  $ ,1300,3203,99,99,99 ! diquark su_1
266  $ ,3300,3303,99,99,99 ! diquark ss_1
267  $ ,2400,4103,99,99,99 ! diquark cd_1
268  $ ,1400,4203,99,99,99 ! diquark cu_1
269  $ ,3400,4303,99,99,99 ! diquark cs_1
270  $ ,4400,4403,99,99,99 ! diquark cc_1
271  $ ,2500,5103,99,99,99 ! diquark bd_1
272  $ ,1500,5203,99,99,99 ! diquark bu_1
273  $ ,3500,5303,99,99,99 ! diquark bs_1
274  $ ,4500,5403,99,99,99 ! diquark bc_1
275  $ ,5500,5503,99,99,99 / ! diquark bb_1
276  data ((idt(i,j),i=1,ncode),j= 185,188)/
277  $ 800000091,91,99,99,99 ! parton system in cluster fragmentation (pythia)
278  $ ,800000092,92,99,99,99 ! parton system in string fragmentation (pythia)
279  $ ,800000093,93,99,99,99 ! parton system in independent system (pythia)
280  $ ,800000094,94,99,99,99 / ! CMshower (pythia)
281  data ((idt(i,j),i=1,ncode),j= 189,208)/
282  $ -340,431,99,99,99 ! Ds+
283  $ ,340,-431,99,99,99 ! Ds-
284  $ ,-241,413,99,99,99 ! D*+
285  $ ,241,-413,99,99,99 ! D*-
286  $ ,-141,423,99,99,99 ! D*0
287  $ ,141,-423,99,99,99 ! D*0b
288  $ ,-341,433,99,99,99 ! Ds*+
289  $ ,341,-433,99,99,99 ! Ds*-
290  $ ,250,511,99,99,99 ! B0
291  $ ,150,521,99,99,99 ! B+
292  $ ,350,531,99,99,99 ! B0s+
293  $ ,450,541,99,99,99 ! Bc+
294  $ ,251,513,99,99,99 ! B*0
295  $ ,151,523,99,99,99 ! B*+
296  $ ,351,533,99,99,99 ! B*0s+
297  $ ,451,543,99,99,99 ! B*c+
298  $ ,440,441,99,99,99 ! etac
299  $ ,441,443,99,99,99 ! J/psi
300  $ ,550,551,99,99,99 ! etab
301  $ ,551,553,99,99,99 / ! Upsilon
302  data ((idt(i,j),i=1,ncode),j= 209,264)/
303  $ 2240,4112,99,99,99 ! sigmac0
304  $ ,1240,4212,99,99,99 ! sigmac+
305  $ ,1140,4222,99,99,99 ! sigmac++
306  $ ,2241,4114,99,99,99 ! sigma*c0
307  $ ,1241,4214,99,99,99 ! sigma*c+
308  $ ,1141,4224,99,99,99 ! sigma*c++
309  $ ,3240,4132,99,99,99 ! Xic0
310  $ ,2340,4312,99,99,99 ! Xi'c0
311  $ ,3140,4232,99,99,99 ! Xic+
312  $ ,1340,4322,99,99,99 ! Xi'c+
313  $ ,3340,4332,99,99,99 ! omegac0
314  $ ,2341,4314,99,99,99 ! Xi*c0
315  $ ,1341,4324,99,99,99 ! Xi*c+
316  $ ,3341,4334,99,99,99 ! omega*c0
317  $ ,2440,4412,99,99,99 ! dcc
318  $ ,2441,4414,99,99,99 ! dcc*
319  $ ,1440,4422,99,99,99 ! ucc
320  $ ,1441,4424,99,99,99 ! ucc*
321  $ ,3440,4432,99,99,99 ! scc
322  $ ,3441,4434,99,99,99 ! scc*
323  $ ,4441,4444,99,99,99 ! ccc*
324  $ ,2250,5112,99,99,99 ! sigmab-
325  $ ,2150,5122,99,99,99 ! lambdab0
326  $ ,3250,5132,99,99,99 ! sdb
327  $ ,4250,5142,99,99,99 ! cdb
328  $ ,1250,5212,99,99,99 ! sigmab0
329  $ ,1150,5222,99,99,99 ! sigmab+
330  $ ,3150,5232,99,99,99 ! sub
331  $ ,4150,5242,99,99,99 ! cub
332  $ ,2350,5312,99,99,99 ! dsb
333  $ ,1350,5322,99,99,99 ! usb
334  $ ,3350,5332,99,99,99 ! ssb
335  $ ,4350,5342,99,99,99 ! csb
336  $ ,2450,5412,99,99,99 ! dcb
337  $ ,1450,5422,99,99,99 ! ucb
338  $ ,3450,5432,99,99,99 ! scb
339  $ ,4450,5442,99,99,99 ! ccb
340  $ ,2550,5512,99,99,99 ! dbb
341  $ ,1550,5522,99,99,99 ! ubb
342  $ ,3550,5532,99,99,99 ! sbb
343  $ ,3550,5542,99,99,99 ! scb
344  $ ,2251,5114,99,99,99 ! sigma*b-
345  $ ,1251,5214,99,99,99 ! sigma*b0
346  $ ,1151,5224,99,99,99 ! sigma*b+
347  $ ,2351,5314,99,99,99 ! dsb*
348  $ ,1351,5324,99,99,99 ! usb*
349  $ ,3351,5334,99,99,99 ! ssb*
350  $ ,2451,5414,99,99,99 ! dcb*
351  $ ,1451,5424,99,99,99 ! ucb*
352  $ ,3451,5434,99,99,99 ! scb*
353  $ ,4451,5444,99,99,99 ! ccb*
354  $ ,2551,5514,99,99,99 ! dbb*
355  $ ,1551,5524,99,99,99 ! ubb*
356  $ ,3551,5534,99,99,99 ! sbb*
357  $ ,4551,5544,99,99,99 ! cbb*
358  $ ,5551,5554,99,99,99 / ! bbb*
359  data ((idt(i,j),i=1,ncode),j= 265,295)/
360  $ 123,10213,99,99,99 ! b1
361  $ ,122,10211,99,99,99 ! a0+
362  $ ,233,10313,99,99,99 ! K0_1
363  $ ,232,10311,99,99,99 ! K*0_1
364  $ ,133,10323,99,99,99 ! K+_1
365  $ ,132,10321,99,99,99 ! K*+_1
366  $ ,143,10423,99,99,99 ! D0_1
367  $ ,132,10421,99,99,99 ! D*0_1
368  $ ,243,10413,99,99,99 ! D+_1
369  $ ,242,10411,99,99,99 ! D*+_1
370  $ ,343,10433,99,99,99 ! D+s_1
371  $ ,342,10431,99,99,99 ! D*0s+_1
372  $ ,223,10113,99,99,99 ! b_10
373  $ ,222,10111,99,99,99 ! a_00
374  $ ,113,10223,99,99,99 ! h_10
375  $ ,112,10221,99,99,99 ! f_00
376  $ ,333,10333,99,99,99 ! h'_10
377  $ ,332,10331,99,99,99 ! f'_00
378  $ ,443,10443,99,99,99 ! h_1c0
379  $ ,442,10441,99,99,99 ! Xi_0c0
380  $ ,444,10443,99,99,99 ! psi'
381  $ ,253,10513,99,99,99 ! db_10
382  $ ,252,10511,99,99,99 ! db*_00
383  $ ,153,10523,99,99,99 ! ub_10
384  $ ,152,10521,99,99,99 ! ub*_00
385  $ ,353,10533,99,99,99 ! sb_10
386  $ ,352,10531,99,99,99 ! sb*_00
387  $ ,453,10543,99,99,99 ! cb_10
388  $ ,452,10541,99,99,99 ! cb*_00
389  $ ,553,10553,99,99,99 ! Upsilon'
390  $ ,552,10551,99,99,99 / ! Upsilon'*
391  data ((idt(i,j),i=1,ncode),j= 296,325)/
392  $ 124,20213,99,99,99 ! a_1+
393  $ ,125,215,99,99,99 ! a_2+
394  $ ,234,20313,99,99,99 ! K*0_1
395  $ ,235,315,99,99,99 ! K*0_2
396  $ ,134,20323,99,99,99 ! K*+_1
397  $ ,135,325,99,99,99 ! K*+_2
398  $ ,144,20423,99,99,99 ! D*_10
399  $ ,135,425,99,99,99 ! D*_20
400  $ ,244,20413,99,99,99 ! D*_1+
401  $ ,245,415,99,99,99 ! D*_2+
402  $ ,344,20433,99,99,99 ! D*_1s+
403  $ ,345,435,99,99,99 ! D*_2s+
404  $ ,224,20113,99,99,99 ! a_10
405  $ ,225,115,99,99,99 ! a_20
406  $ ,114,20223,99,99,99 ! f_10
407  $ ,115,225,99,99,99 ! f_20
408  $ ,334,20333,99,99,99 ! f'_10
409  $ ,335,335,99,99,99 ! f'_20
410  $ ,444,20443,99,99,99 ! Xi_1c0
411  $ ,445,445,99,99,99 ! Xi_2c0
412  $ ,254,20513,99,99,99 ! db*_10
413  $ ,255,515,99,99,99 ! db*_20
414  $ ,154,20523,99,99,99 ! ub*_10
415  $ ,155,525,99,99,99 ! ub*_20
416  $ ,354,20533,99,99,99 ! sb*_10
417  $ ,355,535,99,99,99 ! sb*_20
418  $ ,454,20543,99,99,99 ! cb*_10
419  $ ,455,545,99,99,99 ! cb*_20
420  $ ,554,20553,99,99,99 ! bb*_10
421  $ ,555,555,99,99,99 / ! bb*_20
422  data ((idt(i,j),i=1,ncode),j= 326,nidt)/
423  $ 11099,9900110,99,99,99 ! diff pi0 state
424  $ ,12099,9900210,99,99,99 ! diff pi+ state
425  $ ,22099,9900220,99,99,99 ! diff omega state
426  $ ,33099,9900330,99,99,99 ! diff phi state
427  $ ,44099,9900440,99,99,99 ! diff J/psi state
428  $ ,112099,9902210,99,99,99 ! diff proton state
429  $ ,122099,9902110,99,99,99 ! diff neutron state
430  $ ,800000110,110,99,99,99 ! Reggeon
431  $ ,800000990,990,99,99,99 / ! Pomeron
432 
433 
434 
435 
436  nidtmx=64
437  id1=idi
438  if(code1.eq.'nxs')then
439  i=1
440  elseif(code1.eq.'pdg')then
441  i=2
442  elseif(code1.eq.'qgs')then
443  i=3
444  elseif(code1.eq.'ghe')then
445  id1=ighenex(id1)
446  i=1
447  elseif(code1.eq.'sib')then
448  i=5
449  elseif(code1.eq.'cor')then
450  id1=kipart(id1)
451  id1=ighenex(id1)
452  i=1
453  else
454  stop "unknown code in idtrafo"
455  endif
456  if(code2.eq.'nxs')then
457  j=1
458  ji=j
459  if(i.eq.5.and.id1.gt.1004)then !nucleus from Sibyll
460  idtrafo=(id1-1000)*100
461  return
462  elseif(id1.eq.130.and.i.eq.2)then
463  idtrafo=-20
464  return
465  endif
466  if(i.eq.2) nidtmx=nidt
467  elseif(code2.eq.'pdg')then
468  j=2
469  ji=j
470  if(id1.eq.-20.and.i.eq.1)then
471  idtrafo=130
472  return
473  endif
474  if(i.eq.1) nidtmx=nidt
475  elseif(code2.eq.'qgs')then
476  j=3
477  ji=j
478  elseif(code2.eq.'ghe')then
479  j=4
480  ji=j
481  elseif(code2.eq.'sib')then
482  j=5
483  ji=j
484  elseif(code2.eq.'cor')then
485  j=4
486  ji=6
487  else
488  stop "unknown code in idtrafo"
489  endif
490  iad1=abs(id1)
491  isi=sign(1,id1)
492 
493  do n=1,nidtmx
494  if(iad1.eq.abs(idt(i,n)))then
495  m=1
496  if(n+m.lt.nidt)then
497  do while(abs(idt(i,n+m)).eq.iad1)
498  m=m+1
499  enddo
500  endif
501  mm=0
502  if(m.gt.1)then
503  if(m.eq.2.and.idt(i,n)*idt(i,n+1).lt.0)then
504  if(id1.eq.idt(i,n+1))mm=1
505  isi=1
506  else
507  mm=int(drangen(dummy)*dble(m))
508  endif
509  endif
510  idtrafo=idt(j,n+mm)*isi
511  if(abs(idtrafo).eq.99)call utstop('New particle not allowed ')
512  if(idtrafo.lt.0.and.j.eq.4)then !gheisha id always >0
513  iadtr=abs(idtrafo)
514  if(iadtr.ge.20.and.iadtr.le.22)then
515  idtrafo=iadtr+3
516  elseif(iadtr.eq.26.or.iadtr.eq.27)then
517  idtrafo=iadtr+2
518  elseif(iadtr.ge.14)then
519  idtrafo=iadtr+1
520  else
521  idtrafo=iadtr
522  endif
523  endif
524  if(ji.eq.6)idtrafo=ikpart(idtrafo)
525  return
526  end if
527  enddo
528 
529  !print *, 'idtrafo:',code1,'->', code2,id1,' not found. '
530  !stop'(2403200818) '
531  idtrafo=0
532 
533  end
534  real*8 function drangen(u)
535  real*8 u
536  drangen= rndc(u)
537  end
538  subroutine utstop(msg)
539  character*(*) msg
540  write(0,*) msg
541  end
542 
543 
544 
parameter
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
i
nodes i
Definition: cstdatmos0.fNew2.h:93
true
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon true
Definition: cblkElemag.h:7
idtrafo
integer function idtrafo(code1, code2, idi)
Definition: codeConv.f:36
rndc
subroutine rndc(u)
Definition: rnd.f:91
ccos2kf
subroutine ccos2kf(code, subcode, chg, kf)
Definition: ckf2cos.f:216
j
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
Definition: cblkHeavy.h:36
code
*Zfirst p fm *Zfirst p Zfirst p code
Definition: ZavoidUnionMap.h:1
charge
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p charge
Definition: ZavoidUnionMap.h:1
m
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
Definition: ZlatfitBD.h:35
utstop
subroutine utstop(msg)
Definition: codeConv.f:539
drangen
real *8 function drangen(u)
Definition: codeConv.f:535
n
integer n
Definition: Zcinippxc.h:1
subcode
*Zfirst p fm *Zfirst p Zfirst p Zfirst p subcode
Definition: ZavoidUnionMap.h:1