]> git.uio.no Git - u/mrichter/AliRoot.git/blame - EPOS/epos167/epos-dky-165.f
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Fix for a case of TPC-TOF combination in Bayesian PID
[u/mrichter/AliRoot.git] / EPOS / epos167 / epos-dky-165.f
CommitLineData
9ef1c2d9 1c-----------------------------------------------------------------------
2 subroutine decayall(n)
3c-----------------------------------------------------------------------
4c decay of objects n to nptl, including their children
5c-----------------------------------------------------------------------
6 include 'epos.inc'
7 common/cttaus/tpro,zpro,ttar,ztar,ttaus,detap,detat
8 double precision tpro,zpro,ttar,ztar,ttaus,detap,detat
9 ttaus=1
10 np1=n
11 1 np2=nptl
12 do ip=np1,np2
13 if(istptl(ip).eq.0)then !consider last generation particles
14 call hdecas(ip,iret)
15 if(iret.eq.1)stop'error in hdecas detected in decay'
16 endif
17 enddo
18 np1=np2+1
19 if(np1.le.nptl)goto1
20 end
21
22
23c-----------------------------------------------------------------------
24 subroutine hdecas(i,iret)
25c-----------------------------------------------------------------------
26c decay of object i (main decay routine)
27c-----------------------------------------------------------------------
28
29 include 'epos.inc'
30 double precision tpro,zpro,ttar,ztar,ttaus,detap,detat,zor,tor
31 common/cttaus/tpro,zpro,ttar,ztar,ttaus,detap,detat
32 double precision ttaux,ttauz
33 integer jcdu(nflav,2)
34 common/cdfptl/idfptl(mxptl)
35
36 iret=0
37 nptlb=nptl
38
39c no last generation -> no decay
40
41 if(istptl(i).ne.0)return
42
43 if(nptl.gt.mxptl-10)then
44 call alist('end&',1,nptl)
45 call utstop('hdecas: mxptl too small&')
46 endif
47c entry
48
49 call utpri('hdecas',ish,ishini,5)
50 ttauz=ttaus
51
52
53c small droplet decay
54
55 if(iabs(idptl(i)).gt.1e8)then
56 stop'hdecas: no longer supported (2). '
57 endif
58
59c ordinary decay
60
61 call idmass(111,amrho0)
62 call idmass(221,amomeg)
63 ioi=iorptl(i)
64 if(ioi.gt.0)then
65 if(.not.(iabs(idptl(ioi)).lt.10000
66 *.and.jorptl(i).eq.0))then
67 if(idptl(i).eq.111)idptl(i)=221
68 if(idptl(i).eq.221.and.rangen().gt.0.5)idptl(i)=111
69 endif
70 endif
71
72 ida=iabs(idptl(i))
73 if(.not.(iappl.eq.7.and.i.eq.1))then
74 if(mod(ndecay ,10).eq.1
75 *.and.ida.ne.0.and.ida.lt.10000)goto1000
76 if(mod(ndecay/10 ,10).eq.1.and.ida.eq. 20)goto1000
77 if(mod(ndecay/100 ,10).eq.1.and.ida.eq.2130)goto1000
78 if(mod(ndecay/1000 ,10).eq.1.and.ida.eq.1130)goto1000
79 if(mod(ndecay/1000 ,10).eq.1.and.ida.eq.2230)goto1000
80 if(mod(ndecay/10000 ,10).eq.1.and.ida.eq.2330)goto1000
81 if(mod(ndecay/10000 ,10).eq.1.and.ida.eq.1330)goto1000
82 if(mod(ndecay/100000 ,10).eq.1.and.ida.eq.3331)goto1000
83 if(mod(ndecay/1000000,10).eq.1.and.ida.eq. 110)goto1000
84
85 if(nrnody.gt.0)then
86 do nod=1,nrnody
87 if(idptl(i).eq.nody(nod))goto1000
88 enddo
89 endif
90
91 endif
92
93 call hdecay(i,iret)
94 if(iret.eq.1)goto1000
95 if(nptl.le.nptlb)then
96 iret=-1
97 goto 1000
98 endif
99
1002000 continue ! ---successful decay---
101
102 istptl(i)=1
103 ifrptl(1,i)=nptlb+1
104 ifrptl(2,i)=nptl
105
106 t=tivptl(2,i)
107 x=xorptl(1,i)+(t-xorptl(4,i))*pptl(1,i)/pptl(4,i)
108 y=xorptl(2,i)+(t-xorptl(4,i))*pptl(2,i)/pptl(4,i)
109 z=xorptl(3,i)+(t-xorptl(4,i))*pptl(3,i)/pptl(4,i)
110 call jtaux(t,z,ttaux)
111 ttaus=ttaux
112 if( ttaus.gt.0d0 ) then
113 call jtauin
114 call jtaus(z,ttest,sz)
115 if (abs(t-ttest).gt.1e-5.and.ish.ge.1) then
116 call utmsg('hdecas')
117 write(ifch,*)'***** t /= ttest'
118 write(ifch,*)t,ttest,i,z,t,xorptl(3,i),xorptl(4,i)
119 $ ,pptl(3,i),pptl(4,i)
120 call utmsgf
121 endif
122 endif
123
124c loop over decay products
125
126 do 20 n=nptlb+1,nptl
127 iorptl(n)=i
128 jorptl(n)=0
129 istptl(n)=0
130 ifrptl(1,n)=0
131 ifrptl(2,n)=0
132 if(idfptl(i).eq.0)then
133 idfptl(n)=0
134 else
135 idfptl(n)=1
136 endif
137 rad=0
138 phi=0
139 ti=t
140 zi=z
141 xorptl(1,n)=x + rad*cos(phi)
142 xorptl(2,n)=y + rad*sin(phi)
143 xorptl(3,n)=zi
144 xorptl(4,n)=ti
145 io=n
1461 io=iorptl(io)
147 if(ish.ge.4)write(ifch,*)'io = ',io,' origin: ',iorptl(io)
148 if(iorptl(io).gt.0)goto1
149 if(ish.ge.4)write(ifch,*)'origin: ',io,idptl(io)
150 zor=xorptl(3,io)
151 tor=xorptl(4,io)
152 call idquac(io,nq,ndummy1,ndummy2,jcdu)
153 r=rangen()
154 tauran=-taurea*alog(r)
155 call jtaix(n,tauran,zor,tor,zis,tis)
156 tivptl(1,n)=amax1(ti,tis)
157 call idtau(idptl(n),pptl(4,n),pptl(5,n),taugm)
158 r=rangen()
159 tivptl(2,n)=t+taugm*(-alog(r))
160 ityptl(n)=ityptl(i)
161 radptl(n)=0.
162 dezptl(n)=0.
163 itsptl(n)=itsptl(i)
164 rinptl(n)=rinptl(i)
16520 continue
166
167 if(iabs(idptl(nptlb+1)).le.6) then
168 call gakli2(0,0)
169 write (*,*) 'nptlb+1,nptl:',nptlb+1,nptl
170 istptl(nptlb+1)=1
171 do n=nptlb+2,nptl
172 istptl(n)=20
173 enddo
174 call gakfra(iret)
175 if(iret.eq.1)goto1000
176 call gakli2(0,0)
177 endif
178
1791000 continue
180 ttaus=ttauz
181 call jtauin
182 call utprix('hdecas',ish,ishini,5)
183 return
184 end
185
186c-----------------------------------------------------------------------
187 subroutine hdecay(ip,iret)
188c-----------------------------------------------------------------------
189c decays particle ip from /cptl/
190c for ip being no resonance: call hnbaaa
191c for ip being resonance: standard resonance decay procedure
192c-----------------------------------------------------------------------
193 include 'epos.inc'
194 double precision tpro,zpro,ttar,ztar,ttaus,detap,detat
195 common/cttaus/tpro,zpro,ttar,ztar,ttaus,detap,detat
196 common/wco/wmass2,wgam2
197 parameter (mxlook=10000,mxdky=2000)
198 common/dkytab/look(mxlook),cbr(mxdky),mode(5,mxdky)
199 dimension pgen(5,10),rnd(10),u(3),beta(3)
200 1 ,reduce(10)
201 dimension prest(4,10),kno(10)
202 data reduce/1.,1.,2.,5.,15.,60.,250.,1500.,1.2E4,1.2E5/
203 data twome/1.022006e-3/
204
205c fctn definitions
206 dot(i1,i2)=prest(4,i1)*prest(4,i2)-prest(1,i1)*prest(1,i2)
207 *-prest(2,i1)*prest(2,i2)-prest(3,i1)*prest(3,i2)
208c charged w propagator.
209 wprop(z)=(z-wmass2**2)**2+(wmass2*wgam2)**2
210
211 call utpri('hdecay',ish,ishini,5)
212
213 ipp=ip
214 iret=0
215 nptlb=nptl
216
217 if(ish.ge.5)write(ifch,*)'ip,id,mass: ',ip,idptl(ip),pptl(5,ip)
218c no k_long decay
219c ---------------
220c if(idptl(ip).eq.-20)goto1000
221
222c select decay mode
223c -----------------
224 ntry=0
2252 ntry=ntry+1
226 if(ntry.gt.100)then
227 if(ish.ge.1)then
228 call utmsg('hdecay')
229 write(ifch,*)'***** decay not possible. iret = 1.'
230 call utmsgf
231 endif
232 iret=1
233 goto1000
234 endif
235 idlv1=idptl(ip)
236 amss=pptl(5,ip)
237
238c Decay of deuteron
239
240 if(idlv1.eq.17)then
241 amss=1.01*amss
242 naddptl=2
243 call idmass(1120,amnew)
244 pptl(5,nptl+1)=amnew
245 idptl(nptl+1)=1120
246 sum=amnew
247 call idmass(1220,amnew)
248 pptl(5,nptl+2)=amnew
249 idptl(nptl+2)=1220
250 sum=sum+amnew
251 goto 111
252 endif
253
254c Decay of triton
255
256 if(idlv1.eq.18)then
257 amss=1.01*amss
258 naddptl=3
259 call idmass(1120,amnew)
260 pptl(5,nptl+1)=amnew
261 idptl(nptl+1)=1120
262 sum=amnew
263 call idmass(1220,amnew)
264 pptl(5,nptl+2)=amnew
265 idptl(nptl+2)=1220
266 sum=sum+amnew
267 call idmass(1220,amnew)
268 pptl(5,nptl+3)=amnew
269 idptl(nptl+3)=1220
270 sum=sum+amnew
271 goto 111
272 endif
273
274c Decay of alpha
275
276 if(idlv1.eq.19)then
277 amss=1.01*amss
278 naddptl=4
279 call idmass(1120,amnew)
280 pptl(5,nptl+1)=amnew
281 idptl(nptl+1)=1120
282 sum=amnew
283 call idmass(1220,amnew)
284 pptl(5,nptl+2)=amnew
285 idptl(nptl+2)=1220
286 sum=sum+amnew
287 call idmass(1120,amnew)
288 pptl(5,nptl+3)=amnew
289 idptl(nptl+3)=1120
290 sum=sum+amnew
291 call idmass(1220,amnew)
292 pptl(5,nptl+4)=amnew
293 idptl(nptl+4)=1220
294 sum=sum+amnew
295 goto 111
296 endif
297
2981 continue !select one of the decay channel
299 ipoint=look(iabs(idlv1))-1
300 if(idlv1.eq.-20)ipoint=look(320)-1
301 if(ipoint.lt.0) goto1000
302 try=rangen()
303100 ipoint=ipoint+1
304 if(ish.ge.4)write(ifch,*)'ipoint,cbr,try',ipoint,cbr(ipoint),try
305 if(try.gt.cbr(ipoint)) goto100
306 naddptl=0
307 sum=0.
308c nstart=nptl+1 !?????????????????unused
309 do 110 i=1,5 !store id and mass of products
310 if(mode(i,ipoint).eq.0) goto 110
311 if(nptl+naddptl+1.gt.mxptl) goto 9999
312 if(iabs( mode(1,ipoint)) .le. 6.and.i.eq.2)then !decay into quark ???
313 call vedi(mode(1,ipoint),mode(2,ipoint),k3,idlv1)
314 idptl(new)=idlv1
315 call idmass(idlv1,amnew)
316 pptl(5,new)=amnew
317 sum=pptl(5,new)
318 else !decay into particles
319 naddptl=naddptl+1
320 new=nptl+naddptl
321 idptl(new)=mode(i,ipoint)
322 idlv1=idptl(new)
323 call idmass(idlv1,pptl(5,new))
324 sum=sum+pptl(5,new)
325 endif
326 110 continue
327 111 continue
328 if(naddptl.ne.1.and.sum.ge.amss)goto 2
329 112 naddptl1=naddptl-1
330 do 120 j=1,5
331 pgen(j,1)=pptl(j,ip)
332120 continue
333 pgen(5,1)=amss !needed because of deuteron, triton and alpha decay and OK
334
335 pgen(5,naddptl)=pptl(5,nptl+naddptl)
336 if(naddptl.eq.1) goto 700 !one body decay
337 if(naddptl.eq.2) goto 400 !two body decay
338
339 if(ish.ge.4)write(ifch,*)'>= 3 body decay'
340
341c use kroll-wada distribution for pi0 and eta dalitz decays.
342c ----------------------------------------------
343 if(.not.((idptl(ip).eq.110.or.idptl(ip).eq.220).and.
344 1iabs(idptl(nptl+2)).eq.12)) goto 130
345 ntry=0 !decay of pi0 or eta into electron
346125 ntry=ntry+1
347 if(ntry.gt.10)then
348 if(ish.ge. 0)then
349 call utmsg('hdecay')
350 write(ifch,*)'***** ntry > 10. iret = 1.'
351 write(ifch,*)'***** amee,ree,wtee',amee,ree,wtee
352 call utmsgf
353 endif
354 iret=1
355 goto1000
356 endif
357 amee=twome*(pptl(5,ip)/twome)**rangen()
358 ree=(twome/amee)**2
359 wtee=(1.-(amee/pptl(5,ip))**2)**3*sqrt(1.-ree)*(1.+.5*ree)
360 if(wtee.lt.rangen()) goto125
361 pgen(5,2)=amee
362 goto400
363130 continue
364
365c calculate maximum phase-space weight
366c ------------------------------------
367 wtmax=1./reduce(naddptl)
368 sum1=pgen(5,1)
369 sum2=sum-pptl(5,nptl+1)
370 do 200 i=1,naddptl1
371 wtmax=wtmax*utpcm(sum1,sum2,pptl(5,nptl+i))
372 sum1=sum1-pptl(5,nptl+i)
373 sum2=sum2-pptl(5,nptl+i+1)
374200 continue
375
376c generate uniform naddptl-body phase space
377c --------------------------------------
378 ntry=0
379300 ntry=ntry+1
380 if(ntry.gt.10000)then
381 if(ish.ge. 0)then
382 call utmsg('hdecay')
383 write(ifch,*)'***** infinite loop (2). iret = 1.'
384 write(ifch,*)'***** ip,idptl(ip),pptl(5,ip):'
385 *,ip,idptl(ip),pptl(5,ip)
386 write(ifch,*)'***** wt,wtmax:',wt,wtmax
387 write(ifch,*)'***** i,pgen(5,i),pptl(5,nptl+i),idptl(nptl+i):'
388 do i=1,naddptl
389 write(ifch,*)i,pgen(5,i),pptl(5,nptl+i),idptl(nptl+i)
390 enddo
391 call utmsgf
392 endif
393 iret=1
394 goto1000
395 endif
396 rnd(1)=1.
397 do 310 i=2,naddptl1
398 rnew=rangen()
399 i1=i-1
400 do 320 jj1=1,i1
401 j=i-jj1
402 jsave=j+1
403 if(rnew.le.rnd(j)) goto310
404 rnd(jsave)=rnd(j)
405320 continue
406310 rnd(jsave)=rnew
407 rnd(naddptl)=0.
408 wt=1.
409 sum1=sum
410 do 330 i=2,naddptl
411 sum1=sum1-pptl(5,nptl+i-1)
412 pgen(5,i)=sum1+rnd(i)*(pgen(5,1)-sum)
413 a=pgen(5,i-1)
414 b=pgen(5,i)
415 c=pptl(5,nptl+i-1)
416 wt=wt*utpcm(a,b,c)
417330 continue
418 if(wt.lt.rangen()*wtmax) goto300
419
420c carry out two-body decays in pgen frames
421c ----------------------------------------
422400 continue
423 if(ish.ge.4)write(ifch,*)'2 body decay'
424 do 410 i=1,naddptl1
425 qcm=utpcm(pgen(5,i),pgen(5,i+1),pptl(5,nptl+i))
426 u(3)=2.*rangen()-1.
427 phi=2.*pi*rangen()
428 u(1)=sqrt(1.-u(3)**2)*cos(phi)
429 u(2)=sqrt(1.-u(3)**2)*sin(phi)
430 do 420 j=1,3
431 pptl(j,nptl+i)=qcm*u(j)
432 pgen(j,i+1)=-pptl(j,nptl+i)
433420 continue
434 pptl(4,nptl+i)=sqrt(qcm**2+pptl(5,nptl+i)**2)
435 pgen(4,i+1)=sqrt(qcm**2+pgen(5,i+1)**2)
436410 continue
437 do 430 j=1,4
438 pptl(j,nptl+naddptl)=pgen(j,naddptl)
439430 continue
440
441c boost pgen frames to lab frame
442c also save momenta in rest frame (last frame)
443c -------------------------------------------------
444 do 500 ii=1,naddptl1
445 i=naddptl-ii
446 do 510 j=1,3
447 beta(j)=pgen(j,i)/pgen(4,i)
448510 continue
449 gamma=pgen(4,i)/pgen(5,i)
450 do 520 k=i,naddptl
451 k1=nptl+k
452 bp=beta(1)*pptl(1,k1)+beta(2)*pptl(2,k1)+beta(3)*pptl(3,k1)
453 do 530 j=1,3
454 prest(j,k)=pptl(j,k1)
455 pptl(j,k1)=pptl(j,k1)+gamma*beta(j)*(pptl(4,k1)
456 1+bp*gamma/(gamma+1.))
457530 continue
458 prest(4,k)=pptl(4,k1)
459 pptl(4,k1)=gamma*(pptl(4,k1)+bp)
460 if(pptl(4,k1).lt.1.d-5)then
461 pptl(4,k1)=sqrt(pptl(1,k1)*pptl(1,k1)+pptl(2,k1)*pptl(2,k1)
462 & +pptl(3,k1)*pptl(3,k1))
463 endif
464520 continue
465500 continue
466
467c matrix elements
468c ---------------
469 if(iabs(idptl(ip)).eq.14)then !muon decay
470 goto 650
471 elseif(naddptl.eq.3)then
472 if(idptl(ip).eq.221.or.idptl(ip).eq.331)then !omeg and phi decay
473 goto 610
474 elseif(iabs(idptl(ip)).eq.130.or. !Kl and K decay
475 1 idptl(ip).eq.-20)then
476 if(iabs(idptl(nptl+2)).lt.20)then !semi-leptonic
477 goto 630
478 else !hadronic
479 goto 640
480 endif
481 elseif(iabs(idptl(nptl+1)).lt.20.and. !other semi-leptonic decay
482 1 idptl(nptl+1).ne.10)then
483 goto 620
484 elseif(iabs(idptl(nptl+2)).le.6)then
485 goto 605 !decay into quark
486 else
487 goto 800
488 endif
489 else
490 goto 800
491 endif
492
493 605 wt=pptl(5,ip)*pptl(5,nptl+1)*dot(2,3)
494 IF(wt.LT.rangen()*pptl(5,ip)**4/16.) goto 300
495 ams=sqrt(dot(2,2)+dot(3,3)+2.*dot(2,3))
496 kno(1)=idptl(nptl+2)
497 kno(2)=idptl(nptl+3)
498 if(ammin(kno(1),kno(2)).gt.ams)then
499 call vedi(kno(1),kno(2),iddum,idlv2)
500 idptl(nptl+2)=idlv2
501 call idmass(idlv2,amnew2)
502 pptl(5,nptl+2)=amnew2
503 naddptl=2
504 goto 112
505 endif
506c......multiplicity
507 PS =sqrt(dot(2,2))
508 psq=sqrt(dot(3,3))
509c PSP=PS !!???????????????unused
510 np=0 !!!!?????
511 nq=2
512 CNDE=4.5*LOG(MAX((ams-PS-PSQ)/0.7,1.1))
513c IF(MMAT.EQ.12) CNDE=CNDE+PARJ(63)
514 769 NTRY=NTRY+1
515 IF(NTRY.GT.1000) THEN
516 write(*,*)'hdecay caught in infinite loop'
517 goto 1000
518 ENDIF
519 GAUSS=SQRT(-2.*CNDE*LOG(MAX(1E-10,rangen())))*
520 & SIN(2.*pi*rangen())
521 ND=0.5+0.5*NP+0.25*NQ+CNDE+GAUSS
522 IF(ND.LT.NP+NQ/2.OR.ND.LT.2.OR.ND.GT.10) GOTO 769
523
524
525c......choose hadrons
526
527
528 132 kno(3)=kno(1)
529 kno(4)=kno(2)
530
531 CONTINUE
532 IF(ND.EQ.NP+NQ/2) GOTO 773
533 DO I=nptl+2,nptl+2+nd-nq/2-1
534 JT=2+1+INT((NQ-1) * rangen() )
535 CALL vedi(kno(JT),0,KFL2,idlv3)
536 idptl(i)=idlv3
537c IF(K(I,2).EQ.0) GOTO 769
538 kno(JT)=-KFL2
539 enddo
540 773 CONTINUE
541 CALL vedi(kno(3),kno(4),KFLDMP,idlv4)
542 idptl(nptl+2+nd-nq/2)=idlv4
543 sum=0.
544 do i=nptl+2,nptl+2+nd-nq/2
545 call idmass(idptl(i),am)
546 pptl(5,i)=am
547 sum=sum+am
548 enddo
549 if(sum.gt.ams) goto 769
550c......goto phase space dis....
551 ip=nptl+2+nd-nq/2+1
552 do j=1,4
553 pptl(j,ip)=pptl(j,ipp)-pptl(j,nptl+1)
554 enddo
555 pptl(5,ip)=ams
556 idptl(ip)=sign(80,idptl(ipp))
557 nptl=nptl+1
558 naddptl=nd
559 goto 112
560
561
562c omeg and phi decay
563c use vectors in rest frame
564c ------------------------------
565610 wt=(pptl(5,nptl+1)*pptl(5,nptl+2)*pptl(5,nptl+3))**2
566 1-(pptl(5,nptl+1)*dot(2,3))**2
567 2-(pptl(5,nptl+2)*dot(1,3))**2
568 3-(pptl(5,nptl+3)*dot(1,2))**2
569 4+2.*dot(1,2)*dot(2,3)*dot(1,3)
570 if(wt.lt.rangen()*pptl(5,ip)**6/108.) goto300
571 goto800
572
573c semileptonic and quark decays
574c use vectors in rest frame, where ip has (m,0,0,0)
575c include w propagator
576c ------------------------------------------------------
577620 wt=(pptl(5,ip)*prest(4,2))*dot(1,3)
578 s12=pptl(5,nptl+1)**2+pptl(5,nptl+2)**2+2.*dot(1,2)
579 s12max=pptl(5,ip)**2
580 wt=wt*wprop(s12max)/wprop(s12)
581 if(wt.lt.rangen()*pptl(5,ip)**4/16.) goto 300
582 goto 800
583
584c semileptonic kaon decays
585c use vectors in rest frame, where ip has (m,0,0,0)
586c include form factor FML
587c ------------------------------------------------------
588630 if(iabs(idptl(ip)).eq.130)then
589 if(iabs(idptl(nptl+2)).eq.12)then
590 ncha=1 !K -> Pi0 + e + Nu
591 else
592 ncha=2 !K -> Pi0 + Mu + Nu
593 endif
594 else
595 if(iabs(idptl(nptl+2)).eq.12)then
596 ncha=3 !K0 -> Pi + e + Nu
597 else
598 ncha=4 !K0 -> Pi + Mu + Nu
599 endif
600 endif
601
602 wt=FML(ncha,pptl(5,ip),pptl(5,nptl+1),pptl(5,nptl+2)
603 & ,prest(4,1),prest(4,2),prest(4,3))
604 if(wt.lt.rangen()) goto 300
605 goto 800
606
607c hadronic kaon decays
608c use vectors in rest frame, where ip has (m,0,0,0)
609c include form factor FM
610c ------------------------------------------------------
611640 if(iabs(idptl(ip)).eq.130)then
612 if(iabs(idptl(nptl+3)).eq.120)then
613 ncha=1 !K -> 3 Pi
614 else
615 ncha=2 !K -> Pi + 2 Pi0
616 endif
617 else
618 if(iabs(idptl(nptl+1)).eq.110)then
619 ncha=3 !K0 -> 3 Pi0
620 else
621 ncha=4 !K0 -> 2 Pi + Pi0
622 endif
623 endif
624 S0=(pptl(5,ip)**2+pptl(5,nptl+1)**2+pptl(5,nptl+2)**2
625 & +pptl(5,nptl+3)**2)/3.d0
626 S1=pptl(5,ip)**2+pptl(5,nptl+1)**2-2.*prest(4,1)*pptl(5,ip)
627 S2=pptl(5,ip)**2+pptl(5,nptl+2)**2-2.*prest(4,2)*pptl(5,ip)
628 S3=pptl(5,ip)**2+pptl(5,nptl+3)**2-2.*prest(4,3)*pptl(5,ip)
629 wt=FM(ncha,S0,S1,S2,S3)
630 if(wt.lt.rangen()) goto 300
631 goto 800
632
633c muon decays
634c use vectors in rest frame, where ip has (m,0,0,0)
635c include form factor FMU
636c ------------------------------------------------------
637650 xxx=2.*prest(4,1)/pptl(5,ip) !reduced energy of electron
638 if(xxx.gt.1.) goto 300
639 wt=FMU(xxx)
640 rrr=rangen()
641 if(wt.lt.rrr) goto 300
642 goto 800
643
644c one-particle decays
645c -------------------
646700 continue
647 do 710 j=1,5
648 pptl(j,nptl+1)=pptl(j,ip)
649710 continue
650
651c swap particles and antiparticles if idptl(ip)<0
652c -----------------------------------------------
653 800 continue
654 if(iabs(idptl(ip)).eq.80)then
655 nptl=nptl-1
656 naddptl=naddptl+1
657 endif
658 if(idptl(ipp).ge.0.or.iabs(idptl(ipp)).eq.20) goto 900
659 do 810 i=1,naddptl
660 idabs=iabs(idptl(nptl+i))
661 ifl1=idabs/1000
662 ifl2=mod(idabs/100,10)
663 ifl3=mod(idabs/10,10)
664 if(ifl1.eq.0.and.ifl2.ne.0.and.ifl2.eq.ifl3) goto 810
665 if(idabs.eq.9.or.idabs.eq.10.or.idabs.eq.20) goto 810
666 if(idabs.eq.29.or.idabs.eq.30.or.idabs.eq.40) goto 810
667 idptl(nptl+i)=-idptl(nptl+i)
668 810 continue
669
670 900 continue
671 nptl=nptl+naddptl
672 if(nptl.gt.mxptl)call utstop('hdecay: nptl>mxptl&')
673c nqk=0 !???????????????????unused
674 if(iabs(idptl(nptl)).lt.10.or.mod(idptl(nptl),100).eq.0)then
675c call utstop('hdecay: decay ptcl is parton&')
676 endif
677
678c print
679c -----
680
681 if(ish.ge.3)then
682 write(ifch,140)sngl(ttaus)
683 140 format(/' ----------------------------'/
684 *' decay at tau =',f6.2/
685 *' ----------------------------')
686 write(ifch,*)'decaying object:'
687 call alist('&',ip,ip)
688 write(ifch,*)'decay products:'
689 call alist('&',nptlb+1,nptl)
690 endif
691 if(ish.ge.5)then
692 write(ifch,*)'momentum sum:'
693 do kk=1,5
694 pptl(kk,nptl+1)=0
695 do ii=nptlb+1,nptl
696 pptl(kk,nptl+1)=pptl(kk,nptl+1)+pptl(kk,ii)
697 enddo
698 enddo
699 call alist('&',nptl+1,nptl+1)
700 endif
701
702c exit
703c ----
704
705 1000 continue
706 ip=ipp
707 if(iret.ne.0.and.ish.ge.1)write(ifmt,'(a)')'hdecay: redo event'
708 call utprix('hdecay',ish,ishini,5)
709 return
710
711 9999 call utstop('hdecay: mxptl too small&')
712 end
713
714c---------------------------------------------------------------------
715 subroutine vedi(k1,k2,k3,id)
716c---------------------------------------------------------------------
717 include 'epos.inc'
718 if(k2.eq.0)then
719 if(rangen().lt.pdiqua.and.iabs(k1).lt.6)then
720 ifl1=int(rangen()/pud)+1
721 ifl2=int(rangen()/pud)+1
722 k3=-min(ifl1,ifl2)*1000-max(ifl1,ifl2)*100
723 else
724 k3=int(rangen()/pud)+1
725 endif
726 if(k1.gt.0.and.k1.le.6)k3=-k3
727 if(k1.lt.-1000)k3=-k3
728 else
729 k3=k2
730 endif
731 id=idsp(k1,k3)
732 if(iabs(id).le.999) then
733 ids=max(mod(iabs(id)/100,10),mod(iabs(id)/10,10))
734 if(ids.le.2)then
735 idr=sign(iabs(id)+int(rangen()+0.5),id)
736 elseif(ids.eq.3)then
737 idr=sign(iabs(id)+int(rangen()+0.6),id)
738 else
739 idr=sign(iabs(id)+int(rangen()+0.75),id)
740 endif
741 else
742 idr=sign(iabs(id)+int(0.5+rangen()),id)
743 endif
744 id=idr
745 if(ish.ge.5)write(ifch,*) 'Flavor:',k1,k2,k3,id
746 end
747
748c-----------------------------------------------------------------------
749 subroutine hdecin(lprint)
750c-----------------------------------------------------------------------
751c sets up /dkytab/
752c-----------------------------------------------------------------------
753 include 'epos.inc'
754 common/wco/wmass2,wgam2
755 dimension imode(6)
756 character*8 idlabl,lmode(6),lres
757 character*8 iblank
758 logical lprint
759 parameter (mxlook=10000,mxdky=2000)
760 common/dkytab/look(mxlook),cbr(mxdky),mode(5,mxdky)
761 common/nodcay/nodcay,noeta,nopi0,nonunu,noevol,nohadr
762 logical nodcay,noeta,nopi0,nonunu,noevol,nohadr
763 parameter (ndectb=1192)
764 real dectab(7,ndectb)
765
766 data ((dectab(i,j),i=1,7),j= 1, 18)/
767 * 110., .98850, 10., 10., 0., 0., 0.
768 *, 110.,1.00000, 10., 12., -12., 0., 0.
769 *, 220., .38000, 10., 10., 0., 0., 0.
770 *, 220., .71000, 110., 110., 110., 0., 0.
771 *, 220., .94600, 120.,-120., 110., 0., 0.
772 *, 220., .99500, 120.,-120., 10., 0., 0.
773 *, 220.,1.00000, 10., 12., -12., 0., 0.
774 *, 330., .44100, 220., 120.,-120., 0., 0.
775 *, 330., .66100, 220., 110., 110., 0., 0.
776 *, 330., .95900, 111., 10., 0., 0., 0.
777 *, 330., .98000, 221., 10., 0., 0., 0.
778 *, 330.,1.00000, 10., 10., 0., 0., 0.
779 *, 121.,1.00000, 120., 110., 0., 0., 0.
780 *, 111., .99989, 120.,-120., 0., 0., 0.
781 *, 111., .99993, 12., -12., 0., 0., 0.
782 *, 111.,1.00000, 14., -14., 0., 0., 0.
783 *, 221., .89900, 120.,-120., 110., 0., 0.
784 *, 221., .91200, 120.,-120., 0., 0., 0./
785 data ((dectab(i,j),i=1,7),j= 19, 36)/
786 * 221., .99992, 110., 10., 0., 0., 0.
787 *, 221.,1.00000, 12., -12., 0., 0., 0.
788 *, 331., .48600, 130.,-130., 0., 0., 0.
789 *, 331., .83700, 230.,-230., 0., 0., 0.
790 *, 331., .98400, 120.,-120., 110., 0., 0.
791 *, 331., .99944, 220., 10., 0., 0., 0.
792 *, 331., .99975, 12., -12., 0., 0., 0.
793 *, 331.,1.00000, 14., -14., 0., 0., 0.
794 *, 230., .50000, 20., 0., 0., 0., 0.
795 *, 230.,1.00000, -20., 0., 0., 0., 0.
796 *, 131., .66670, 230., 120., 0., 0., 0.
797 *, 131.,1.00000, 130., 110., 0., 0., 0.
798 *, 231., .66670, 130.,-120., 0., 0., 0.
799 *, 231.,1.00000, 230., 110., 0., 0., 0.
800 *, 240., .11000, 12., -11., 230., 0., 0.
801 *, 240., .17000, 12., -11., 231., 0., 0.
802 *, 240., .28000, 14., -13., 230., 0., 0.
803 *, 240., .34000, 14., -13., 231., 0., 0./
804 data ((dectab(i,j),i=1,7),j= 37, 54)/
805 * 240., .37800, 230.,-120., 0., 0., 0.
806 *, 240., .56300, 230.,-121., 0., 0., 0.
807 *, 240., .60800, 231.,-120., 0., 0., 0.
808 *, 240., .62100, 230.,-120., 110., 0., 0.
809 *, 240., .71000, 130.,-120.,-120., 0., 0.
810 *, 240., .80100, 230.,-120.,-120., 120., 0.
811 *, 240., .87900, 130.,-120.,-120., 110., 0.
812 *, 240., .95400, 230.,-120., 110., 110., 0.
813 *, 240., .96600, 230.,-130., 0., 0., 0.
814 *, 240., .97600, 331.,-120., 0., 0., 0.
815 *, 240., .98800,-130., 231., 0., 0., 0.
816 *, 240.,1.00000,-131., 230., 0., 0., 0.
817 *, 140., .04500, -12., 11., 130., 0., 0.
818 *, 140., .07500, -12., 11., 131., 0., 0.
819 *, 140., .12000, -14., 13., 130., 0., 0.
820 *, 140., .15000, -14., 13., 131., 0., 0.
821 *, 140., .20300, 130.,-120., 0., 0., 0.
822 *, 140., .22700, 230., 110., 0., 0., 0./
823 data ((dectab(i,j),i=1,7),j= 55, 72)/
824 * 140., .24700, 230., 220., 0., 0., 0.
825 *, 140., .28900, 230., 221., 0., 0., 0.
826 *, 140., .45100, 130.,-121., 0., 0., 0.
827 *, 140., .53600, 131.,-120., 0., 0., 0.
828 *, 140., .56200, 231., 110., 0., 0., 0.
829 *, 140., .57600, 230., 111., 0., 0., 0.
830 *, 140., .58700, 130.,-120., 110., 0., 0.
831 *, 140., .60300, 230.,-120., 120., 0., 0.
832 *, 140., .72700, 130.,-120.,-120., 120., 0.
833 *, 140., .87600, 230.,-120., 120., 110., 0.
834 *, 140., .96900, 130.,-120., 110., 110., 0.
835 *, 140.,1.00000, 230., 110., 110., 110., 0.
836 *, 340., .03250, 12., -11., 220., 0., 0.
837 *, 340., .06500, 12., -11., 331., 0., 0.
838 *, 340., .09750, 14., -13., 220., 0., 0.
839 *, 340., .13000, 14., -13., 331., 0., 0.
840 *, 340., .17900,-130., 230., 0., 0., 0.
841 *, 340., .22800,-120., 220., 0., 0., 0./
842 data ((dectab(i,j),i=1,7),j= 73, 90)/
843 * 340., .33800,-131., 230., 0., 0., 0.
844 *, 340., .44800,-130., 231., 0., 0., 0.
845 *, 340., .55800,-120., 331., 0., 0., 0.
846 *, 340., .57500,-130., 230., 110., 0., 0.
847 *, 340., .59200,-230., 230.,-120., 0., 0.
848 *, 340., .69400,-130., 230.,-120., 120., 0.
849 *, 340., .79600,-130., 230., 110., 110., 0.
850 *, 340., .89800,-130., 130.,-120., 110., 0.
851 *, 340.,1.00000,-230., 230.,-120., 110., 0.
852 *, 241., .64000, 140.,-120., 0., 0., 0.
853 *, 241., .92000, 240., 110., 0., 0., 0.
854 *, 241.,1.00000, 240., 10., 0., 0., 0.
855 *, 141., .55000, 140., 110., 0., 0., 0.
856 *, 141.,1.00000, 140., 10., 0., 0., 0.
857 *, 341.,1.00000, 340., 10., 0., 0., 0.
858 *, 441., .07400, 12., -12., 0., 0., 0.
859 *, 441., .14800, 14., -14., 0., 0., 0.
860 *, 441., .15210,-121., 120., 0., 0., 0./
861 data ((dectab(i,j),i=1,7),j= 91,108)/
862 * 441., .15620, 111., 110., 0., 0., 0.
863 *, 441., .16020, 121.,-120., 0., 0., 0.
864 *, 441., .16300,-121., 111., 120., 0., 0.
865 *, 441., .16580, 121.,-121., 110., 0., 0.
866 *, 441., .16860, 121., 111.,-120., 0., 0.
867 *, 441., .28740, 120.,-120., 130.,-130., 0.
868 *, 441., .40620, 110., 110., 130.,-130., 0.
869 *, 441., .52500, 120.,-120., 120.,-120., 0.
870 *, 441., .64380, 120.,-120., 110., 110., 0.
871 *, 441., .76260, 110., 110., 110., 110., 0.
872 *, 441., .88130, 120.,-120., 230.,-230., 0.
873 *, 441.,1.00000, 110., 110., 230., 230., 0.
874 *, 150., .06000, -12., 11., 140., 0., 0.
875 *, 150., .12000, -12., 11., 141., 0., 0.
876 *, 150., .18000, -14., 13., 140., 0., 0.
877 *, 150., .24000, -14., 13., 141., 0., 0.
878 *, 150., .25500, -16., 15., 140., 0., 0.
879 *, 150., .27000, -16., 15., 141., 0., 0./
880 data ((dectab(i,j),i=1,7),j=109,122)/
881 * 150., .28050, 140., 120., 0., 0., 0.
882 *, 150., .29100, 140., 121., 0., 0., 0.
883 *, 150., .30150, 141., 120., 0., 0., 0.
884 *, 150., .31200, 141., 121., 0., 0., 0.
885 *, 150., .32650, 140.,-340., 0., 0., 0.
886 *, 150., .34100, 140.,-341., 0., 0., 0.
887 *, 150., .35550, 141.,-340., 0., 0., 0.
888 *, 150., .37000, 141.,-341., 0., 0., 0.
889 *, 150., 0.820 , 1., -4., 1., -2., 0.
890 *, 150., 0.920 , 1., -2., 1., -4., 0.
891 *, 150., 0.975 , 1., -4., 4., -3., 0.
892 *, 150., 0.985 , 1., -3., 4., -4., 0.
893 *, 150., 0.995 , 1., -1., 1., -2., 0.
894 *, 150., 1. , 1., -1., 4., -3., 0./
895 data ((dectab(i,j),i=1,7),j=123,142)
896 */ 250., .06000, -12., 11., 240., 0., 0.
897 *, 250., .12000, -12., 11., 241., 0., 0.
898 *, 250., .18000, -14., 13., 240., 0., 0.
899 *, 250., .24000, -14., 13., 241., 0., 0.
900 *, 250., .25500, -16., 15., 240., 0., 0.
901 *, 250., .27000, -16., 15., 241., 0., 0.
902 *, 250., .28050, 240., 120., 0., 0., 0.
903 *, 250., .29100, 240., 121., 0., 0., 0.
904 *, 250., .30150, 241., 120., 0., 0., 0.
905 *, 250., .31200, 241., 121., 0., 0., 0.
906 *, 250., .32650, 240.,-340., 0., 0., 0.
907 *, 250., .34100, 240.,-341., 0., 0., 0.
908 *, 250., .35550, 241.,-340., 0., 0., 0.
909 *, 250., .37000, 241.,-341., 0., 0., 0.
910 *, 250., 0.820 , 2., -4., 1., -2., 0.
911 *, 250., 0.920 , 2., -2., 1., -4., 0.
912 *, 250., 0.975 , 2., -4., 4., -3., 0.
913 *, 250., 0.985 , 2., -3., 4., -4., 0.
914 *, 250., 0.995 , 2., -1., 1., -2., 0.
915 *, 250., 1. , 2., -1., 4., -3., 0./
916 data ((dectab(i,j),i=1,7),j=143,176)/
917 * 238*1. /
918 data ((dectab(i,j),i=1,7),j=177,190)
919 * /350., .06000, 12., -11., 340., 0., 0.
920 *, 350., .12000, 12., -11., 341., 0., 0.
921 *, 350., .18000, 14., -13., 340., 0., 0.
922 *, 350., .24000, 14., -13., 341., 0., 0.
923 *, 350., .25500, 16., -15., 340., 0., 0.
924 *, 350., .27000, 16., -15., 341., 0., 0.
925 *, 350., .28050, 340., 120., 0., 0., 0.
926 *, 350., .29100, 340., 121., 0., 0., 0.
927 *, 350., .30150, 341., 120., 0., 0., 0.
928 *, 350., .31200, 341., 121., 0., 0., 0.
929 *, 350., .32650, 340.,-340., 0., 0., 0.
930 *, 350., .34100, 340.,-341., 0., 0., 0.
931 *, 350., .35550, 341.,-340., 0., 0., 0.
932 *, 350., .37000, 341.,-341., 0., 0., 0./
933 data ((dectab(i,j),i=1,7),j=191,196)/
934 * 350., 0.820 , 3., -4., 1., -2., 0.
935 *, 350., 0.920 , 3., -2., 1., -4., 0.
936 *, 350., 0.975 , 3., -4., 4., -3., 0.
937 *, 350., 0.985 , 3., -3., 4., -4., 0.
938 *, 350., 0.995 , 3., -1., 1., -2., 0.
939 *, 350., 1. , 3., -1., 4., -3., 0./
940 data ((dectab(i,j),i=1,7),j=197,244)/
941 * 336*1. /
942 data ((dectab(i,j),i=1,7),j=245,262)/
943 * 160., .33330, -1., 2., -5., 0., 0.
944 *, 160., .66660, -4., 3., -5., 0., 0.
945 *, 160., .77770, 11., -12., -5., 0., 0.
946 *, 160., .88880, 13., -14., -5., 0., 0.
947 *, 160.,1.00000, -15., 16., -5., 0., 0.
948 *, 260., .33330, -1., 2., -5., 0., 0.
949 *, 260., .66660, -4., 3., -5., 0., 0.
950 *, 260., .77770, -11., 12., -5., 0., 0.
951 *, 260., .88880, -13., 14., -5., 0., 0.
952 *, 260.,1.00000, -15., 16., -5., 0., 0.
953 *, 360., .33330, -1., 2., -5., 0., 0.
954 *, 360., .66660, -4., 3., -5., 0., 0.
955 *, 360., .77770, -11., 12., -5., 0., 0.
956 *, 360., .88880, -13., 14., -5., 0., 0.
957 *, 360.,1.00000, -15., 16., -5., 0., 0.
958 *, 151.,1.00000, 150., 10., 0., 0., 0.
959 *, 251.,1.00000, 250., 10., 0., 0., 0.
960 *, 351.,1.00000, 350., 10., 0., 0., 0./
961 data ((dectab(i,j),i=1,7),j=263,280)/
962 * 161.,1.00000, 160., 10., 0., 0., 0.
963 *, 261.,1.00000, 260., 10., 0., 0., 0.
964 *, 361.,1.00000, 360., 10., 0., 0., 0.
965 *,1230.,1.00000,2130., 10., 0., 0., 0.
966 *,1111.,1.00000,1120., 120., 0., 0., 0.
967 *,1121., .66670,1120., 110., 0., 0., 0.
968 *,1121.,1.00000,1220., 120., 0., 0., 0.
969 *,1221., .66670,1220., 110., 0., 0., 0.
970 *,1221.,1.00000,1120.,-120., 0., 0., 0.
971 *,2221.,1.00000,1220.,-120., 0., 0., 0.
972 *,1131., .88000,2130., 120., 0., 0., 0.
973 *,1131., .94000,1130., 110., 0., 0., 0.
974 *,1131.,1.00000,1230., 120., 0., 0., 0.
975 *,1231., .88000,2130., 110., 0., 0., 0.
976 *,1231., .94000,1130.,-120., 0., 0., 0.
977 *,1231.,1.00000,2230., 120., 0., 0., 0.
978 *,2231., .88000,2130.,-120., 0., 0., 0.
979 *,2231., .94000,1230.,-120., 0., 0., 0./
980 data ((dectab(i,j),i=1,7),j=281,298)/
981 * 2231.,1.00000,2230., 110., 0., 0., 0.
982 *,1331., .66670,2330., 120., 0., 0., 0.
983 *,1331.,1.00000,1330., 110., 0., 0., 0.
984 *,2331., .66670,1330.,-120., 0., 0., 0.
985 *,2331.,1.00000,2330., 110., 0., 0., 0.
986 *, 16., .18000, 12., -11., 15., 0., 0.
987 *, 16., .36000, 14., -13., 15., 0., 0.
988 *, 16., .45100,-120., 15., 0., 0., 0.
989 *, 16., .66000,-121., 15., 0., 0., 0.
990 *, 16., .78000, 110., 110.,-120., 15., 0.
991 *, 16., .83600, 120.,-120.,-120., 15., 0.
992 *, 16.,1.00000, 120., 110.,-120.,-120., 15.
993 *,2140., .03750, -12., 11.,2130., 0., 0.
994 *,2140., .07500, -12., 11.,1231., 0., 0.
995 *,2140., .11250, -14., 13.,2130., 0., 0.
996 *,2140., .15000, -14., 13.,1231., 0., 0.
997 *,2140., .18200,2130., 120., 0., 0., 0.
998 *,2140., .21300,1230., 120., 0., 0., 0./
999 data ((dectab(i,j),i=1,7),j=299,316)/
1000 * 2140., .24400,1120.,-230., 0., 0., 0.
1001 *,2140., .29500,1131., 110., 0., 0., 0.
1002 *,2140., .34600,1231., 120., 0., 0., 0.
1003 *,2140., .39700,1121.,-230., 0., 0., 0.
1004 *,2140., .44800,1111.,-130., 0., 0., 0.
1005 *,2140., .49900,1130., 111., 0., 0., 0.
1006 *,2140., .55000,1230., 121., 0., 0., 0.
1007 *,2140., .60100,1120.,-231., 0., 0., 0.
1008 *,2140., .65800,1120.,-230., 120.,-120., 0.
1009 *,2140., .71500,1120.,-230., 110., 110., 0.
1010 *,2140., .77200,1120.,-130., 120., 110., 0.
1011 *,2140., .82900,1220.,-230., 120., 110., 0.
1012 *,2140., .88600,1220.,-130., 120., 120., 0.
1013 *,2140., .94300,2130., 120., 120.,-120., 0.
1014 *,2140.,1.00000,2130., 120., 110., 110., 0.
1015 *,1140.,1.00000,2140., 120., 0., 0., 0.
1016 *,1240.,1.00000,2140., 110., 0., 0., 0.
1017 *,2240.,1.00000,2140.,-120., 0., 0., 0./
1018 data ((dectab(i,j),i=1,7),j=317,334)/
1019 * 1340., .03750, -12., 11.,1330., 0., 0.
1020 *,1340., .07500, -12., 11.,1331., 0., 0.
1021 *,1340., .11250, -14., 13.,1330., 0., 0.
1022 *,1340., .15000, -14., 13.,1331., 0., 0.
1023 *,1340., .19900,1330., 120., 0., 0., 0.
1024 *,1340., .24800,1231., 130., 0., 0., 0.
1025 *,1340., .28800,1330., 120., 0., 0., 0.
1026 *,1340., .32800,1131.,-230., 0., 0., 0.
1027 *,1340., .36800,1330., 121., 0., 0., 0.
1028 *,1340., .40800,1130.,-230., 0., 0., 0.
1029 *,1340., .44800,1330., 120., 110., 0., 0.
1030 *,1340., .48800,2330., 120., 120., 0., 0.
1031 *,1340., .52800,1130.,-130., 120., 0., 0.
1032 *,1340., .56800,1130.,-230., 110., 0., 0.
1033 *,1340., .60800,1230.,-230., 120., 0., 0.
1034 *,1340., .66400,2130.,-230., 120., 110., 0.
1035 *,1340., .72000,2130.,-130., 120., 120., 0.
1036 *,1340., .77600,1130.,-230., 120.,-120., 0./
1037 data ((dectab(i,j),i=1,7),j=335,352)/
1038 * 1340., .83200,1130.,-230., 110., 110., 0.
1039 *,1340., .88800,1330., 120., 120.,-120., 0.
1040 *,1340., .94400,1330., 120., 110., 110., 0.
1041 *,1340.,1.00000,2330., 120., 120., 110., 0.
1042 *,3140., .03750, -12., 11.,1330., 0., 0.
1043 *,3140., .07500, -12., 11.,1331., 0., 0.
1044 *,3140., .11250, -14., 13.,1330., 0., 0.
1045 *,3140., .15000, -14., 13.,1331., 0., 0.
1046 *,3140., .19900,1330., 120., 0., 0., 0.
1047 *,3140., .24800,1231., 130., 0., 0., 0.
1048 *,3140., .28800,1330., 120., 0., 0., 0.
1049 *,3140., .32800,1131.,-230., 0., 0., 0.
1050 *,3140., .36800,1330., 121., 0., 0., 0.
1051 *,3140., .40800,1130.,-230., 0., 0., 0.
1052 *,3140., .44800,1330., 120., 110., 0., 0.
1053 *,3140., .48800,2330., 120., 120., 0., 0.
1054 *,3140., .52800,1130.,-130., 120., 0., 0.
1055 *,3140., .56800,1130.,-230., 110., 0., 0./
1056 data ((dectab(i,j),i=1,7),j=353,370)/
1057 * 3140., .60800,1230.,-230., 120., 0., 0.
1058 *,3140., .66400,2130.,-230., 120., 110., 0.
1059 *,3140., .72000,2130.,-130., 120., 120., 0.
1060 *,3140., .77600,1130.,-230., 120.,-120., 0.
1061 *,3140., .83200,1130.,-230., 110., 110., 0.
1062 *,3140., .88800,1330., 120., 120.,-120., 0.
1063 *,3140., .94400,1330., 120., 110., 110., 0.
1064 *,3140.,1.00000,2330., 120., 120., 110., 0.
1065 *,2340., .03750, -12., 11.,2330., 0., 0.
1066 *,2340., .07500, -12., 11.,2331., 0., 0.
1067 *,2340., .11250, -14., 13.,2330., 0., 0.
1068 *,2340., .15000, -14., 13.,2331., 0., 0.
1069 *,2340., .17500,2330., 120., 0., 0., 0.
1070 *,2340., .20000,1330., 110., 0., 0., 0.
1071 *,2340., .22500,1130.,-130., 0., 0., 0.
1072 *,2340., .25000,1230.,-230., 0., 0., 0.
1073 *,2340., .29500,2331., 120., 0., 0., 0.
1074 *,2340., .34000,1331., 110., 0., 0., 0./
1075 data ((dectab(i,j),i=1,7),j=371,388)/
1076 * 2340., .38500,1131.,-130., 0., 0., 0.
1077 *,2340., .43000,1231.,-230., 0., 0., 0.
1078 *,2340., .47500,2330., 121., 0., 0., 0.
1079 *,2340., .52000,1330., 111., 0., 0., 0.
1080 *,2340., .56500,1130.,-131., 0., 0., 0.
1081 *,2340., .61000,1230.,-231., 0., 0., 0.
1082 *,2340., .64900,2130.,-230., 120.,-120., 0.
1083 *,2340., .68800,2130.,-230., 110., 110., 0.
1084 *,2340., .72700,2130.,-130., 120., 110., 0.
1085 *,2340., .76600,1130.,-230.,-120., 110., 0.
1086 *,2340., .80500,1130.,-130., 120.,-120., 0.
1087 *,2340., .84400,1130.,-130., 110., 110., 0.
1088 *,2340., .88300,1330., 120.,-120., 110., 0.
1089 *,2340., .92200,1330., 110., 110., 110., 0.
1090 *,2340., .96100,2330., 120., 120.,-120., 0.
1091 *,2340.,1.00000,2330., 120., 110., 110., 0.
1092 *,3240., .03750, -12., 11.,2330., 0., 0.
1093 *,3240., .07500, -12., 11.,2331., 0., 0./
1094 data ((dectab(i,j),i=1,7),j=389,406)/
1095 * 3240., .11250, -14., 13.,2330., 0., 0.
1096 *,3240., .15000, -14., 13.,2331., 0., 0.
1097 *,3240., .17500,2330., 120., 0., 0., 0.
1098 *,3240., .20000,1330., 110., 0., 0., 0.
1099 *,3240., .22500,1130.,-130., 0., 0., 0.
1100 *,3240., .25000,1230.,-230., 0., 0., 0.
1101 *,3240., .29500,2331., 120., 0., 0., 0.
1102 *,3240., .34000,1331., 110., 0., 0., 0.
1103 *,3240., .38500,1131.,-130., 0., 0., 0.
1104 *,3240., .43000,1231.,-230., 0., 0., 0.
1105 *,3240., .47500,2330., 121., 0., 0., 0.
1106 *,3240., .52000,1330., 111., 0., 0., 0.
1107 *,3240., .56500,1130.,-131., 0., 0., 0.
1108 *,3240., .61000,1230.,-231., 0., 0., 0.
1109 *,3240., .64900,2130.,-230., 120.,-120., 0.
1110 *,3240., .68800,2130.,-230., 110., 110., 0.
1111 *,3240., .72700,2130.,-130., 120., 110., 0.
1112 *,3240., .76600,1130.,-230.,-120., 110., 0./
1113 data ((dectab(i,j),i=1,7),j=407,424)/
1114 * 3240., .80500,1130.,-130., 120.,-120., 0.
1115 *,3240., .84400,1130.,-130., 110., 110., 0.
1116 *,3240., .88300,1330., 120.,-120., 110., 0.
1117 *,3240., .92200,1330., 110., 110., 110., 0.
1118 *,3240., .96100,2330., 120., 120.,-120., 0.
1119 *,3240.,1.00000,2330., 120., 110., 110., 0.
1120 *,3340., .07500, -12., 11.,3331., 0., 0.
1121 *,3340., .15000, -12., 11.,3331., 0., 0.
1122 *,3340., .25000,1330.,-230., 0., 0., 0.
1123 *,3340., .31000,3331., 120., 0., 0., 0.
1124 *,3340., .37000,1331.,-230., 0., 0., 0.
1125 *,3340., .43000,1330.,-231., 0., 0., 0.
1126 *,3340., .49000,2330.,-230., 120., 0., 0.
1127 *,3340., .55000,1330.,-230., 110., 0., 0.
1128 *,3340., .61000,1330.,-130., 120., 0., 0.
1129 *,3340., .67500,3331., 120., 120.,-120., 0.
1130 *,3340., .74000,3331., 120., 110., 110., 0.
1131 *,3340., .80500,1330.,-230., 120.,-120., 0./
1132 data ((dectab(i,j),i=1,7),j=425,442)/
1133 * 3340., .87000,1330.,-230., 110., 110., 0.
1134 *,3340., .93500,2330.,-230., 120., 110., 0.
1135 *,3340.,1.00000,2330.,-130., 120., 120., 0.
1136 *,1141.,1.00000,2140., 120., 0., 0., 0.
1137 *,1241.,1.00000,2140., 110., 0., 0., 0.
1138 *,2241.,1.00000,2140.,-120., 0., 0., 0.
1139 *,1341., .66670,2340., 120., 0., 0., 0.
1140 *,1341.,1.00000,1340., 110., 0., 0., 0.
1141 *,2341., .66670,1340.,-120., 0., 0., 0.
1142 *,2341.,1.00000,2340., 110., 0., 0., 0.
1143 *,3341.,1.00000,3340., 110., 0., 0., 0.
1144 *,1150., .06000, 12., -11.,1140., 0., 0.
1145 *,1150., .12000, 12., -11.,1141., 0., 0.
1146 *,1150., .18000, 14., -13.,1140., 0., 0.
1147 *,1150., .24000, 14., -13.,1141., 0., 0.
1148 *,1150., .25500, 16., -15.,1140., 0., 0.
1149 *,1150., .27000, 16., -15.,1141., 0., 0.
1150 *,1150., .28925,1140.,-120., 0., 0., 0./
1151 data ((dectab(i,j),i=1,7),j=443,460)/
1152 * 1150., .30850,1140.,-121., 0., 0., 0.
1153 *,1150., .32775,1141.,-120., 0., 0., 0.
1154 *,1150., .34700,1141.,-121., 0., 0., 0.
1155 *,1150., .35775,1140., 340., 0., 0., 0.
1156 *,1150., .36850,1140., 341., 0., 0., 0.
1157 *,1150., .37925,1141., 340., 0., 0., 0.
1158 *,1150., .39000,1141., 341., 0., 0., 0.
1159 *,1150., .42050,1140.,-120., 110., 0., 0.
1160 *,1150., .45100,1140.,-120., 220., 0., 0.
1161 *,1150., .48150,1140.,-120., 111., 0., 0.
1162 *,1150., .51200,1140.,-120., 221., 0., 0.
1163 *,1150., .54250,1140.,-121., 110., 0., 0.
1164 *,1150., .57300,1140.,-121., 220., 0., 0.
1165 *,1150., .60350,1140.,-121., 111., 0., 0.
1166 *,1150., .63400,1140.,-121., 221., 0., 0.
1167 *,1150., .66450,1141.,-120., 110., 0., 0.
1168 *,1150., .69500,1141.,-120., 220., 0., 0.
1169 *,1150., .72550,1141.,-120., 111., 0., 0./
1170 data ((dectab(i,j),i=1,7),j=461,478)/
1171 * 1150., .75600,1141.,-120., 221., 0., 0.
1172 *,1150., .78650,1141.,-121., 110., 0., 0.
1173 *,1150., .81700,1141.,-121., 220., 0., 0.
1174 *,1150., .84750,1141.,-121., 111., 0., 0.
1175 *,1150., .87800,1141.,-121., 221., 0., 0.
1176 *,1150., .89325,1140.,-130., 230., 0., 0.
1177 *,1150., .90850,1140.,-130., 231., 0., 0.
1178 *,1150., .92375,1140.,-131., 230., 0., 0.
1179 *,1150., .93900,1140.,-131., 231., 0., 0.
1180 *,1150., .95425,1141.,-130., 230., 0., 0.
1181 *,1150., .96950,1141.,-130., 231., 0., 0.
1182 *,1150., .98475,1141.,-131., 230., 0., 0.
1183 *,1150.,1.00000,1141.,-131., 231., 0., 0.
1184 *,1250., .06000, 12., -11.,1240., 0., 0.
1185 *,1250., .12000, 12., -11.,1241., 0., 0.
1186 *,1250., .18000, 14., -13.,1240., 0., 0.
1187 *,1250., .24000, 14., -13.,1241., 0., 0.
1188 *,1250., .25500, 16., -15.,1240., 0., 0./
1189 data ((dectab(i,j),i=1,7),j=479,496)/
1190 * 1250., .27000, 16., -15.,1241., 0., 0.
1191 *,1250., .28925,1240.,-120., 0., 0., 0.
1192 *,1250., .30850,1240.,-121., 0., 0., 0.
1193 *,1250., .32775,1241.,-120., 0., 0., 0.
1194 *,1250., .34700,1241.,-121., 0., 0., 0.
1195 *,1250., .35775,1240., 340., 0., 0., 0.
1196 *,1250., .36850,1240., 341., 0., 0., 0.
1197 *,1250., .37925,1241., 340., 0., 0., 0.
1198 *,1250., .39000,1241., 341., 0., 0., 0.
1199 *,1250., .42050,1240.,-120., 110., 0., 0.
1200 *,1250., .45100,1240.,-120., 220., 0., 0.
1201 *,1250., .48150,1240.,-120., 111., 0., 0.
1202 *,1250., .51200,1240.,-120., 221., 0., 0.
1203 *,1250., .54250,1240.,-121., 110., 0., 0.
1204 *,1250., .57300,1240.,-121., 220., 0., 0.
1205 *,1250., .60350,1240.,-121., 111., 0., 0.
1206 *,1250., .63400,1240.,-121., 221., 0., 0.
1207 *,1250., .66450,1241.,-120., 110., 0., 0./
1208 data ((dectab(i,j),i=1,7),j=497,514)/
1209 * 1250., .69500,1241.,-120., 220., 0., 0.
1210 *,1250., .72550,1241.,-120., 111., 0., 0.
1211 *,1250., .75600,1241.,-120., 221., 0., 0.
1212 *,1250., .78650,1241.,-121., 110., 0., 0.
1213 *,1250., .81700,1241.,-121., 220., 0., 0.
1214 *,1250., .84750,1241.,-121., 111., 0., 0.
1215 *,1250., .87800,1241.,-121., 221., 0., 0.
1216 *,1250., .89325,1240.,-130., 230., 0., 0.
1217 *,1250., .90850,1240.,-130., 231., 0., 0.
1218 *,1250., .92375,1240.,-131., 230., 0., 0.
1219 *,1250., .93900,1240.,-131., 231., 0., 0.
1220 *,1250., .95425,1241.,-130., 230., 0., 0.
1221 *,1250., .96950,1241.,-130., 231., 0., 0.
1222 *,1250., .98475,1241.,-131., 230., 0., 0.
1223 *,1250.,1.00000,1241.,-131., 231., 0., 0.
1224 *,1350., .06000, 12., -11.,1340., 0., 0.
1225 *,1350., .12000, 12., -11.,1341., 0., 0.
1226 *,1350., .18000, 14., -13.,1340., 0., 0./
1227 data ((dectab(i,j),i=1,7),j=515,532)/
1228 * 1350., .24000, 14., -13.,1341., 0., 0.
1229 *,1350., .25500, 16., -15.,1340., 0., 0.
1230 *,1350., .27000, 16., -15.,1341., 0., 0.
1231 *,1350., .28925,1340.,-120., 0., 0., 0.
1232 *,1350., .30850,1340.,-121., 0., 0., 0.
1233 *,1350., .32775,1341.,-120., 0., 0., 0.
1234 *,1350., .34700,1341.,-121., 0., 0., 0.
1235 *,1350., .35775,1340., 340., 0., 0., 0.
1236 *,1350., .36850,1340., 341., 0., 0., 0.
1237 *,1350., .37925,1341., 340., 0., 0., 0.
1238 *,1350., .39000,1341., 341., 0., 0., 0.
1239 *,1350., .42050,1340.,-120., 110., 0., 0.
1240 *,1350., .45100,1340.,-120., 220., 0., 0.
1241 *,1350., .48150,1340.,-120., 111., 0., 0.
1242 *,1350., .51200,1340.,-120., 221., 0., 0.
1243 *,1350., .54250,1340.,-121., 110., 0., 0.
1244 *,1350., .57300,1340.,-121., 220., 0., 0.
1245 *,1350., .60350,1340.,-121., 111., 0., 0./
1246 data ((dectab(i,j),i=1,7),j=533,550)/
1247 * 1350., .63400,1340.,-121., 221., 0., 0.
1248 *,1350., .66450,1341.,-120., 110., 0., 0.
1249 *,1350., .69500,1341.,-120., 220., 0., 0.
1250 *,1350., .72550,1341.,-120., 111., 0., 0.
1251 *,1350., .75600,1341.,-120., 221., 0., 0.
1252 *,1350., .78650,1341.,-121., 110., 0., 0.
1253 *,1350., .81700,1341.,-121., 220., 0., 0.
1254 *,1350., .84750,1341.,-121., 111., 0., 0.
1255 *,1350., .87800,1341.,-121., 221., 0., 0.
1256 *,1350., .89325,1340.,-130., 230., 0., 0.
1257 *,1350., .90850,1340.,-130., 231., 0., 0.
1258 *,1350., .92375,1340.,-131., 230., 0., 0.
1259 *,1350., .93900,1340.,-131., 231., 0., 0.
1260 *,1350., .95425,1341.,-130., 230., 0., 0.
1261 *,1350., .96950,1341.,-130., 231., 0., 0.
1262 *,1350., .98475,1341.,-131., 230., 0., 0.
1263 *,1350.,1.00000,1341.,-131., 231., 0., 0.
1264 *,2150., .06000, 12., -11.,2140., 0., 0./
1265 data ((dectab(i,j),i=1,7),j=551,568)/
1266 * 2150., .12000, 12., -11.,1241., 0., 0.
1267 *,2150., .18000, 14., -13.,2140., 0., 0.
1268 *,2150., .24000, 14., -13.,1241., 0., 0.
1269 *,2150., .25500, 16., -15.,2140., 0., 0.
1270 *,2150., .27000, 16., -15.,1241., 0., 0.
1271 *,2150., .28925,2140.,-120., 0., 0., 0.
1272 *,2150., .30850,2140.,-121., 0., 0., 0.
1273 *,2150., .32775,1241.,-120., 0., 0., 0.
1274 *,2150., .34700,1241.,-121., 0., 0., 0.
1275 *,2150., .35775,2140., 340., 0., 0., 0.
1276 *,2150., .36850,2140., 341., 0., 0., 0.
1277 *,2150., .37925,1241., 340., 0., 0., 0.
1278 *,2150., .39000,1241., 341., 0., 0., 0.
1279 *,2150., .42050,2140.,-120., 110., 0., 0.
1280 *,2150., .45100,2140.,-120., 220., 0., 0.
1281 *,2150., .48150,2140.,-120., 111., 0., 0.
1282 *,2150., .51200,2140.,-120., 221., 0., 0.
1283 *,2150., .54250,2140.,-121., 110., 0., 0./
1284 data ((dectab(i,j),i=1,7),j=569,586)/
1285 * 2150., .57300,2140.,-121., 220., 0., 0.
1286 *,2150., .60350,2140.,-121., 111., 0., 0.
1287 *,2150., .63400,2140.,-121., 221., 0., 0.
1288 *,2150., .66450,1241.,-120., 110., 0., 0.
1289 *,2150., .69500,1241.,-120., 220., 0., 0.
1290 *,2150., .72550,1241.,-120., 111., 0., 0.
1291 *,2150., .75600,1241.,-120., 221., 0., 0.
1292 *,2150., .78650,1241.,-121., 110., 0., 0.
1293 *,2150., .81700,1241.,-121., 220., 0., 0.
1294 *,2150., .84750,1241.,-121., 111., 0., 0.
1295 *,2150., .87800,1241.,-121., 221., 0., 0.
1296 *,2150., .89325,2140.,-130., 230., 0., 0.
1297 *,2150., .90850,2140.,-130., 231., 0., 0.
1298 *,2150., .92375,2140.,-131., 230., 0., 0.
1299 *,2150., .93900,2140.,-131., 231., 0., 0.
1300 *,2150., .95425,1241.,-130., 230., 0., 0.
1301 *,2150., .96950,1241.,-130., 231., 0., 0.
1302 *,2150., .98475,1241.,-131., 230., 0., 0./
1303 data ((dectab(i,j),i=1,7),j=587,604)/
1304 * 2150.,1.00000,1241.,-131., 231., 0., 0.
1305 *,2250., .06000, 12., -11.,2240., 0., 0.
1306 *,2250., .12000, 12., -11.,2241., 0., 0.
1307 *,2250., .18000, 14., -13.,2240., 0., 0.
1308 *,2250., .24000, 14., -13.,2241., 0., 0.
1309 *,2250., .25500, 16., -15.,2240., 0., 0.
1310 *,2250., .27000, 16., -15.,2241., 0., 0.
1311 *,2250., .28925,2240.,-120., 0., 0., 0.
1312 *,2250., .30850,2240.,-121., 0., 0., 0.
1313 *,2250., .32775,2241.,-120., 0., 0., 0.
1314 *,2250., .34700,2241.,-121., 0., 0., 0.
1315 *,2250., .35775,2240., 340., 0., 0., 0.
1316 *,2250., .36850,2240., 341., 0., 0., 0.
1317 *,2250., .37925,2241., 340., 0., 0., 0.
1318 *,2250., .39000,2241., 341., 0., 0., 0.
1319 *,2250., .42050,2240.,-120., 110., 0., 0.
1320 *,2250., .45100,2240.,-120., 220., 0., 0.
1321 *,2250., .48150,2240.,-120., 111., 0., 0./
1322 data ((dectab(i,j),i=1,7),j=605,622)/
1323 * 2250., .51200,2240.,-120., 221., 0., 0.
1324 *,2250., .54250,2240.,-121., 110., 0., 0.
1325 *,2250., .57300,2240.,-121., 220., 0., 0.
1326 *,2250., .60350,2240.,-121., 111., 0., 0.
1327 *,2250., .63400,2240.,-121., 221., 0., 0.
1328 *,2250., .66450,2241.,-120., 110., 0., 0.
1329 *,2250., .69500,2241.,-120., 220., 0., 0.
1330 *,2250., .72550,2241.,-120., 111., 0., 0.
1331 *,2250., .75600,2241.,-120., 221., 0., 0.
1332 *,2250., .78650,2241.,-121., 110., 0., 0.
1333 *,2250., .81700,2241.,-121., 220., 0., 0.
1334 *,2250., .84750,2241.,-121., 111., 0., 0.
1335 *,2250., .87800,2241.,-121., 221., 0., 0.
1336 *,2250., .89325,2240.,-130., 230., 0., 0.
1337 *,2250., .90850,2240.,-130., 231., 0., 0.
1338 *,2250., .92375,2240.,-131., 230., 0., 0.
1339 *,2250., .93900,2240.,-131., 231., 0., 0.
1340 *,2250., .95425,2241.,-130., 230., 0., 0./
1341 data ((dectab(i,j),i=1,7),j=623,640)/
1342 * 2250., .96950,2241.,-130., 231., 0., 0.
1343 *,2250., .98475,2241.,-131., 230., 0., 0.
1344 *,2250.,1.00000,2241.,-131., 231., 0., 0.
1345 *,2350., .06000, 12., -11.,2340., 0., 0.
1346 *,2350., .12000, 12., -11.,2341., 0., 0.
1347 *,2350., .18000, 14., -13.,2340., 0., 0.
1348 *,2350., .24000, 14., -13.,2341., 0., 0.
1349 *,2350., .25500, 16., -15.,2340., 0., 0.
1350 *,2350., .27000, 16., -15.,2341., 0., 0.
1351 *,2350., .28925,2340.,-120., 0., 0., 0.
1352 *,2350., .30850,2340.,-121., 0., 0., 0.
1353 *,2350., .32775,2341.,-120., 0., 0., 0.
1354 *,2350., .34700,2341.,-121., 0., 0., 0.
1355 *,2350., .35775,2340., 340., 0., 0., 0.
1356 *,2350., .36850,2340., 341., 0., 0., 0.
1357 *,2350., .37925,2341., 340., 0., 0., 0.
1358 *,2350., .39000,2341., 341., 0., 0., 0.
1359 *,2350., .42050,2340.,-120., 110., 0., 0./
1360 data ((dectab(i,j),i=1,7),j=641,658)/
1361 * 2350., .45100,2340.,-120., 220., 0., 0.
1362 *,2350., .48150,2340.,-120., 111., 0., 0.
1363 *,2350., .51200,2340.,-120., 221., 0., 0.
1364 *,2350., .54250,2340.,-121., 110., 0., 0.
1365 *,2350., .57300,2340.,-121., 220., 0., 0.
1366 *,2350., .60350,2340.,-121., 111., 0., 0.
1367 *,2350., .63400,2340.,-121., 221., 0., 0.
1368 *,2350., .66450,2341.,-120., 110., 0., 0.
1369 *,2350., .69500,2341.,-120., 220., 0., 0.
1370 *,2350., .72550,2341.,-120., 111., 0., 0.
1371 *,2350., .75600,2341.,-120., 221., 0., 0.
1372 *,2350., .78650,2341.,-121., 110., 0., 0.
1373 *,2350., .81700,2341.,-121., 220., 0., 0.
1374 *,2350., .84750,2341.,-121., 111., 0., 0.
1375 *,2350., .87800,2341.,-121., 221., 0., 0.
1376 *,2350., .89325,2340.,-130., 230., 0., 0.
1377 *,2350., .90850,2340.,-130., 231., 0., 0.
1378 *,2350., .92375,2340.,-131., 230., 0., 0./
1379 data ((dectab(i,j),i=1,7),j=659,720)/
1380 * 434*1./
1381 data ((dectab(i,j),i=1,7),j=721,738)/
1382 * 2350., .93900,2340.,-131., 231., 0., 0.
1383 *,2350., .95425,2341.,-130., 230., 0., 0.
1384 *,2350., .96950,2341.,-130., 231., 0., 0.
1385 *,2350., .98475,2341.,-131., 230., 0., 0.
1386 *,2350.,1.00000,2341.,-131., 231., 0., 0.
1387 *,3150., .06000, 12., -11.,3140., 0., 0.
1388 *,3150., .12000, 12., -11.,1341., 0., 0.
1389 *,3150., .18000, 14., -13.,3140., 0., 0.
1390 *,3150., .24000, 14., -13.,1341., 0., 0.
1391 *,3150., .25500, 16., -15.,3140., 0., 0.
1392 *,3150., .27000, 16., -15.,1341., 0., 0.
1393 *,3150., .28925,3140.,-120., 0., 0., 0.
1394 *,3150., .30850,3140.,-121., 0., 0., 0.
1395 *,3150., .32775,1341.,-120., 0., 0., 0.
1396 *,3150., .34700,1341.,-121., 0., 0., 0.
1397 *,3150., .35775,3140., 340., 0., 0., 0.
1398 *,3150., .36850,3140., 341., 0., 0., 0.
1399 *,3150., .37925,1341., 340., 0., 0., 0./
1400 data ((dectab(i,j),i=1,7),j=739,756)/
1401 * 3150., .39000,1341., 341., 0., 0., 0.
1402 *,3150., .42050,3140.,-120., 110., 0., 0.
1403 *,3150., .45100,3140.,-120., 220., 0., 0.
1404 *,3150., .48150,3140.,-120., 111., 0., 0.
1405 *,3150., .51200,3140.,-120., 221., 0., 0.
1406 *,3150., .54250,3140.,-121., 110., 0., 0.
1407 *,3150., .57300,3140.,-121., 220., 0., 0.
1408 *,3150., .60350,3140.,-121., 111., 0., 0.
1409 *,3150., .63400,3140.,-121., 221., 0., 0.
1410 *,3150., .66450,1341.,-120., 110., 0., 0.
1411 *,3150., .69500,1341.,-120., 220., 0., 0.
1412 *,3150., .72550,1341.,-120., 111., 0., 0.
1413 *,3150., .75600,1341.,-120., 221., 0., 0.
1414 *,3150., .78650,1341.,-121., 110., 0., 0.
1415 *,3150., .81700,1341.,-121., 220., 0., 0.
1416 *,3150., .84750,1341.,-121., 111., 0., 0.
1417 *,3150., .87800,1341.,-121., 221., 0., 0.
1418 *,3150., .89325,3140.,-130., 230., 0., 0./
1419 data ((dectab(i,j),i=1,7),j=757,774)/
1420 * 3150., .90850,3140.,-130., 231., 0., 0.
1421 *,3150., .92375,3140.,-131., 230., 0., 0.
1422 *,3150., .93900,3140.,-131., 231., 0., 0.
1423 *,3150., .95425,1341.,-130., 230., 0., 0.
1424 *,3150., .96950,1341.,-130., 231., 0., 0.
1425 *,3150., .98475,1341.,-131., 230., 0., 0.
1426 *,3150.,1.00000,1341.,-131., 231., 0., 0.
1427 *,3250., .06000, 12., -11.,3240., 0., 0.
1428 *,3250., .12000, 12., -11.,2341., 0., 0.
1429 *,3250., .18000, 14., -13.,3240., 0., 0.
1430 *,3250., .24000, 14., -13.,2341., 0., 0.
1431 *,3250., .25500, 16., -15.,3240., 0., 0.
1432 *,3250., .27000, 16., -15.,2341., 0., 0.
1433 *,3250., .28925,3240.,-120., 0., 0., 0.
1434 *,3250., .30850,3240.,-121., 0., 0., 0.
1435 *,3250., .32775,2341.,-120., 0., 0., 0.
1436 *,3250., .34700,2341.,-121., 0., 0., 0.
1437 *,3250., .35775,3240., 340., 0., 0., 0./
1438 data ((dectab(i,j),i=1,7),j=775,792)/
1439 * 3250., .36850,3240., 341., 0., 0., 0.
1440 *,3250., .37925,2341., 340., 0., 0., 0.
1441 *,3250., .39000,2341., 341., 0., 0., 0.
1442 *,3250., .42050,3240.,-120., 110., 0., 0.
1443 *,3250., .45100,3240.,-120., 220., 0., 0.
1444 *,3250., .48150,3240.,-120., 111., 0., 0.
1445 *,3250., .51200,3240.,-120., 221., 0., 0.
1446 *,3250., .54250,3240.,-121., 110., 0., 0.
1447 *,3250., .57300,3240.,-121., 220., 0., 0.
1448 *,3250., .60350,3240.,-121., 111., 0., 0.
1449 *,3250., .63400,3240.,-121., 221., 0., 0.
1450 *,3250., .66450,2341.,-120., 110., 0., 0.
1451 *,3250., .69500,2341.,-120., 220., 0., 0.
1452 *,3250., .72550,2341.,-120., 111., 0., 0.
1453 *,3250., .75600,2341.,-120., 221., 0., 0.
1454 *,3250., .78650,2341.,-121., 110., 0., 0.
1455 *,3250., .81700,2341.,-121., 220., 0., 0.
1456 *,3250., .84750,2341.,-121., 111., 0., 0./
1457 data ((dectab(i,j),i=1,7),j=793,810)/
1458 * 3250., .87800,2341.,-121., 221., 0., 0.
1459 *,3250., .89325,3240.,-130., 230., 0., 0.
1460 *,3250., .90850,3240.,-130., 231., 0., 0.
1461 *,3250., .92375,3240.,-131., 230., 0., 0.
1462 *,3250., .93900,3240.,-131., 231., 0., 0.
1463 *,3250., .95425,2341.,-130., 230., 0., 0.
1464 *,3250., .96950,2341.,-130., 231., 0., 0.
1465 *,3250., .98475,2341.,-131., 230., 0., 0.
1466 *,3250.,1.00000,2341.,-131., 231., 0., 0.
1467 *,3350., .06000, 12., -11.,3340., 0., 0.
1468 *,3350., .12000, 12., -11.,3341., 0., 0.
1469 *,3350., .18000, 14., -13.,3340., 0., 0.
1470 *,3350., .24000, 14., -13.,3341., 0., 0.
1471 *,3350., .25500, 16., -15.,3340., 0., 0.
1472 *,3350., .27000, 16., -15.,3341., 0., 0.
1473 *,3350., .28925,3340.,-120., 0., 0., 0.
1474 *,3350., .30850,3340.,-121., 0., 0., 0.
1475 *,3350., .32775,3341.,-120., 0., 0., 0./
1476 data ((dectab(i,j),i=1,7),j=811,828)/
1477 * 3350., .34700,3341.,-121., 0., 0., 0.
1478 *,3350., .35775,3340., 340., 0., 0., 0.
1479 *,3350., .36850,3340., 341., 0., 0., 0.
1480 *,3350., .37925,3341., 340., 0., 0., 0.
1481 *,3350., .39000,3341., 341., 0., 0., 0.
1482 *,3350., .42050,3340.,-120., 110., 0., 0.
1483 *,3350., .45100,3340.,-120., 220., 0., 0.
1484 *,3350., .48150,3340.,-120., 111., 0., 0.
1485 *,3350., .51200,3340.,-120., 221., 0., 0.
1486 *,3350., .54250,3340.,-121., 110., 0., 0.
1487 *,3350., .57300,3340.,-121., 220., 0., 0.
1488 *,3350., .60350,3340.,-121., 111., 0., 0.
1489 *,3350., .63400,3340.,-121., 221., 0., 0.
1490 *,3350., .66450,3341.,-120., 110., 0., 0.
1491 *,3350., .69500,3341.,-120., 220., 0., 0.
1492 *,3350., .72550,3341.,-120., 111., 0., 0.
1493 *,3350., .75600,3341.,-120., 221., 0., 0.
1494 *,3350., .78650,3341.,-121., 110., 0., 0./
1495 data ((dectab(i,j),i=1,7),j=829,846)/
1496 * 3350., .81700,3341.,-121., 220., 0., 0.
1497 *,3350., .84750,3341.,-121., 111., 0., 0.
1498 *,3350., .87800,3341.,-121., 221., 0., 0.
1499 *,3350., .89325,3340.,-130., 230., 0., 0.
1500 *,3350., .90850,3340.,-130., 231., 0., 0.
1501 *,3350., .92375,3340.,-131., 230., 0., 0.
1502 *,3350., .93900,3340.,-131., 231., 0., 0.
1503 *,3350., .95425,3341.,-130., 230., 0., 0.
1504 *,3350., .96950,3341.,-130., 231., 0., 0.
1505 *,3350., .98475,3341.,-131., 230., 0., 0.
1506 *,3350.,1.00000,3341.,-131., 231., 0., 0.
1507 *,1160., .33300, 1., -2.,1500., 0., 0.
1508 *,1160., .66700, 4., -3.,1500., 0., 0.
1509 *,1160., .77800, -12., 11.,1500., 0., 0.
1510 *,1160., .88900, -14., 13.,1500., 0., 0.
1511 *,1160.,1.00000, -16., 15.,1500., 0., 0.
1512 *,1260., .33300, 1., -2.,2500., 0., 0.
1513 *,1260., .66700, 4., -3.,2500., 0., 0./
1514 data ((dectab(i,j),i=1,7),j=847,864)/
1515 * 1260., .77800, -12., 11.,2500., 0., 0.
1516 *,1260., .88900, -14., 13.,2500., 0., 0.
1517 *,1260.,1.00000, -16., 15.,2500., 0., 0.
1518 *,2260., .33300, 1., -2.,2500., 0., 0.
1519 *,2260., .66700, 4., -3.,2500., 0., 0.
1520 *,2260., .77800, -12., 11.,2500., 0., 0.
1521 *,2260., .88900, -14., 13.,2500., 0., 0.
1522 *,2260.,1.00000, -16., 15.,2500., 0., 0.
1523 *,2160., .33300, 1., -2.,1500., 0., 0.
1524 *,2160., .66700, 4., -3.,1500., 0., 0.
1525 *,2160., .77800, -12., 11.,1500., 0., 0.
1526 *,2160., .88900, -14., 13.,1500., 0., 0.
1527 *,2160.,1.00000, -16., 15.,1500., 0., 0.
1528 *,1360., .33300, 1., -2.,3500., 0., 0.
1529 *,1360., .66700, 4., -3.,3500., 0., 0.
1530 *,1360., .77800, -12., 11.,3500., 0., 0.
1531 *,1360., .88900, -14., 13.,3500., 0., 0.
1532 *,1360.,1.00000, -16., 15.,3500., 0., 0./
1533 data ((dectab(i,j),i=1,7),j=865,882)/
1534 * 2360., .33300, 1., -2.,3500., 0., 0.
1535 *,2360., .66700, 4., -3.,3500., 0., 0.
1536 *,2360., .77800, -12., 11.,3500., 0., 0.
1537 *,2360., .88900, -14., 13.,3500., 0., 0.
1538 *,2360.,1.00000, -16., 15.,3500., 0., 0.
1539 *,3360., .33300, 1., -2.,3500., 0., 0.
1540 *,3360., .66700, 4., -3.,3500., 0., 0.
1541 *,3360., .77800, -12., 11.,3500., 0., 0.
1542 *,3360., .88900, -14., 13.,3500., 0., 0.
1543 *,3360.,1.00000, -16., 15.,3500., 0., 0.
1544 *,1151.,1.00000,1150., 10., 0., 0., 0.
1545 *,1251.,1.00000,1250., 10., 0., 0., 0.
1546 *,2251.,1.00000,2250., 10., 0., 0., 0.
1547 *,1351.,1.00000,1350., 10., 0., 0., 0.
1548 *,2351.,1.00000,2350., 10., 0., 0., 0.
1549 *,3351.,1.00000,3350., 10., 0., 0., 0.
1550 *,1161.,1.00000,1160., 10., 0., 0., 0.
1551 *,1261.,1.00000,1260., 10., 0., 0., 0./
1552 data ((dectab(i,j),i=1,7),j=883,886)/
1553 * 2261.,1.00000,2260., 10., 0., 0., 0.
1554 *,1361.,1.00000,1360., 10., 0., 0., 0.
1555 *,2361.,1.00000,2360., 10., 0., 0., 0.
1556 *,3361.,1.00000,3360., 10., 0., 0., 0./
1557c *---------------------------------------------
1558c * delta++ resonances
1559c *---------------------------------------------
1560 data ((dectab(i,j),i=1,7),j=887,900)/
1561c *--dl++(1620)---------------------------------
1562 * 1112., .30000,1120., 120., 0., 0., 0.
1563 *,1112., .66000,1111., 110., 0., 0., 0.
1564 *,1112., .90000,1121., 120., 0., 0., 0.
1565 *,1112.,1.00000,1120., 120., 110., 0., 0.
1566c *--dl++(1700)---------------------------------
1567 *,1113., .15000,1120., 120., 0., 0., 0.
1568 *,1113., .51000,1111., 110., 0., 0., 0.
1569 *,1113., .75000,1121., 120., 0., 0., 0.
1570 *,1113.,1.00000,1120., 120., 110., 0., 0.
1571c *--dl++(1925)---------------------------------
1572 *,1114., .28000,1120., 120., 0., 0., 0.
1573 *,1114., .40600,1111., 110., 0., 0., 0.
1574 *,1114., .49000,1121., 120., 0., 0., 0.
1575 *,1114., .69000,1120., 121., 0., 0., 0.
1576 *,1114., .70000,1130., 130., 0., 0., 0.
1577 *,1114.,1.00000,1122., 120., 0., 0., 0./
1578c *---------------------------------------------
1579c * delta- resonances
1580c *---------------------------------------------
1581 data ((dectab(i,j),i=1,7),j=901,914)/
1582c *--dl-(1620)----------------------------------
1583 * 2222., .30000,1220.,-120., 0., 0., 0.
1584 *,2222., .66000,2221., 110., 0., 0., 0.
1585 *,2222., .90000,1221.,-120., 0., 0., 0.
1586 *,2222.,1.00000,1220., 110.,-120., 0., 0.
1587c *--dl-(1700)----------------------------------
1588 *,2223., .15000,1220.,-120., 0., 0., 0.
1589 *,2223., .51000,2221., 110., 0., 0., 0.
1590 *,2223., .75000,1221.,-120., 0., 0., 0.
1591 *,2223.,1.00000,1220., 110.,-120., 0., 0.
1592c *--dl-(1925)----------------------------------
1593 *,2224., .28000,1220.,-120., 0., 0., 0.
1594 *,2224., .40600,2221., 110., 0., 0., 0.
1595 *,2224., .49000,1221.,-120., 0., 0., 0.
1596 *,2224., .69000,1220.,-121., 0., 0., 0.
1597 *,2224., .70000,2230., 230., 0., 0., 0.
1598 *,2224.,1.00000,1222.,-120., 0., 0., 0./
1599c *---------------------------------------------
1600c * n*+ resonances + delta+ resonances
1601c *---------------------------------------------
1602 data ((dectab(i,j),i=1,7),j=915,931)/
1603c *--n*+(1440)----------------------------------
1604 * 1122., .20000,1120., 110., 0., 0., 0.
1605 *,1122., .60000,1220., 120., 0., 0., 0.
1606 *,1122., .68000,1111.,-120., 0., 0., 0.
1607 *,1122., .73000,1121., 110., 0., 0., 0.
1608 *,1122., .76000,1221., 120., 0., 0., 0.
1609 *,1122., .84000,1120., 120.,-120., 0., 0.
1610 *,1122., .87000,1120., 110., 110., 0., 0.
1611 *,1122.,1.00000,1220., 120., 110., 0., 0.
1612c *--n*+(1530)----------------------------------
1613 *,1123., .17000,1120., 110., 0., 0., 0.
1614 *,1123., .51000,1220., 120., 0., 0., 0.
1615 *,1123., .57000,1111.,-120., 0., 0., 0.
1616 *,1123., .61000,1121., 110., 0., 0., 0.
1617 *,1123., .63000,1221., 120., 0., 0., 0.
1618 *,1123., .67000,1120., 120.,-120., 0., 0.
1619 *,1123., .68000,1120., 110., 110., 0., 0.
1620 *,1123., .75000,1220., 120., 110., 0., 0.
1621 *,1123.,1.00000,1120., 220., 0., 0., 0./
1622 data ((dectab(i,j),i=1,7),j=932,948)/
1623c *--dl+(1620)----------------------------------
1624 * 1124., .20000,1120., 110., 0., 0., 0.
1625 *,1124., .30000,1220., 120., 0., 0., 0.
1626 *,1124., .54000,1111.,-120., 0., 0., 0.
1627 *,1124., .58000,1121., 110., 0., 0., 0.
1628 *,1124., .90000,1221., 120., 0., 0., 0.
1629 *,1124., .96000,1120., 120.,-120., 0., 0.
1630 *,1124.,1.00000,1220., 120., 110., 0., 0.
1631c *--n*+(1665)----------------------------------
1632 *,1125., .16700,1120., 110., 0., 0., 0.
1633 *,1125., .49970,1220., 120., 0., 0., 0.
1634 *,1125., .62470,1111.,-120., 0., 0., 0.
1635 *,1125., .70800,1121., 110., 0., 0., 0.
1636 *,1125., .74970,1221., 120., 0., 0., 0.
1637 *,1125., .82080,1120., 120.,-120., 0., 0.
1638 *,1125., .85190,1120., 110., 110., 0., 0.
1639 *,1125., .96300,1220., 120., 110., 0., 0.
1640 *,1125., .97300,1120., 220., 0., 0., 0.
1641 *,1125.,1.00000,2130., 130., 0., 0., 0./
1642 data ((dectab(i,j),i=1,7),j=949,955)/
1643c *--dl+(1700)----------------------------------
1644 * 1126., .10000,1120., 110., 0., 0., 0.
1645 *,1126., .15000,1220., 120., 0., 0., 0.
1646 *,1126., .39000,1111.,-120., 0., 0., 0.
1647 *,1126., .43000,1121., 110., 0., 0., 0.
1648 *,1126., .75000,1221., 120., 0., 0., 0.
1649 *,1126., .91500,1120., 120.,-120., 0., 0.
1650 *,1126.,1.00000,1220., 120., 110., 0., 0./
1651 data ((dectab(i,j),i=1,7),j=956,969)/
1652c *--n*+(1710)----------------------------------
1653 * 1127., .04430,1120., 110., 0., 0., 0.
1654 *,1127., .13290,1220., 120., 0., 0., 0.
1655 *,1127., .23790,1111.,-120., 0., 0., 0.
1656 *,1127., .30790,1121., 110., 0., 0., 0.
1657 *,1127., .34290,1221., 120., 0., 0., 0.
1658 *,1127., .41190,1120., 120.,-120., 0., 0.
1659 *,1127., .48090,1120., 110., 110., 0., 0.
1660 *,1127., .54990,1220., 120., 110., 0., 0.
1661 *,1127., .66070,1120., 220., 0., 0., 0.
1662 *,1127., .72800,2130., 130., 0., 0., 0.
1663 *,1127., .74930,1230., 130., 0., 0., 0.
1664 *,1127., .76000,1130., 230., 0., 0., 0.
1665 *,1127., .84000,1120., 111., 0., 0., 0.
1666 *,1127.,1.00000,1220., 121., 0., 0., 0./
1667 data ((dectab(i,j),i=1,7),j=970,980)/
1668c *--dl+(1925)----------------------------------
1669 * 1128., .18700,1120., 110., 0., 0., 0.
1670 *,1128., .28000,1220., 120., 0., 0., 0.
1671 *,1128., .36400,1111.,-120., 0., 0., 0.
1672 *,1128., .37800,1121., 110., 0., 0., 0.
1673 *,1128., .49000,1221., 120., 0., 0., 0.
1674 *,1128., .62300,1120., 111., 0., 0., 0.
1675 *,1128., .69000,1220., 121., 0., 0., 0.
1676 *,1128., .69350,1130., 230., 0., 0., 0.
1677 *,1128., .69900,1230., 130., 0., 0., 0.
1678 *,1128., .89900,1122., 110., 0., 0., 0.
1679 *,1128.,1.00000,1222., 120., 0., 0., 0./
1680c *---------------------------------------------
1681c * n*0 resonances + delta0 resonances
1682c *---------------------------------------------
1683 data ((dectab(i,j),i=1,7),j=981,997)/
1684c *----------n*0(1440)--------------------------
1685 * 1222., .20000,1220., 110., 0., 0., 0.
1686 *,1222., .60000,1120.,-120., 0., 0., 0.
1687 *,1222., .68000,2221., 120., 0., 0., 0.
1688 *,1222., .73000,1221., 110., 0., 0., 0.
1689 *,1222., .76000,1121.,-120., 0., 0., 0.
1690 *,1222., .84000,1220., 120.,-120., 0., 0.
1691 *,1222., .87000,1220., 110., 110., 0., 0.
1692 *,1222.,1.00000,1120.,-120., 110., 0., 0.
1693c *----------n*0(1530)--------------------------
1694 *,1223., .17000,1220., 110., 0., 0., 0.
1695 *,1223., .51000,1120.,-120., 0., 0., 0.
1696 *,1223., .57000,2221., 120., 0., 0., 0.
1697 *,1223., .61000,1221., 110., 0., 0., 0.
1698 *,1223., .63000,1121.,-120., 0., 0., 0.
1699 *,1223., .67000,1220., 120.,-120., 0., 0.
1700 *,1223., .68000,1220., 110., 110., 0., 0.
1701 *,1223., .75000,1120.,-120., 110., 0., 0.
1702 *,1223.,1.00000,1220., 220., 0., 0., 0./
1703 data ((dectab(i,j),i=1,7),j=998,1014)/
1704c *----------dl0(1620)--------------------------
1705 * 1224., .20000,1220., 110., 0., 0., 0.
1706 *,1224., .30000,1120.,-120., 0., 0., 0.
1707 *,1224., .54000,2221., 120., 0., 0., 0.
1708 *,1224., .58000,1221., 110., 0., 0., 0.
1709 *,1224., .90000,1121.,-120., 0., 0., 0.
1710 *,1224., .96500,1220., 120.,-120., 0., 0.
1711 *,1224.,1.00000,1120.,-120., 110., 0., 0.
1712c *----------n*0(1665)--------------------------
1713 *,1225., .16700,1220., 110., 0., 0., 0.
1714 *,1225., .49970,1120.,-120., 0., 0., 0.
1715 *,1225., .62470,2221., 120., 0., 0., 0.
1716 *,1225., .70800,1221., 110., 0., 0., 0.
1717 *,1225., .74970,1121.,-120., 0., 0., 0.
1718 *,1225., .82080,1220., 120.,-120., 0., 0.
1719 *,1225., .85190,1220., 110., 110., 0., 0.
1720 *,1225., .96300,1120.,-120., 110., 0., 0.
1721 *,1225., .97300,1220., 220., 0., 0., 0.
1722 *,1225.,1.00000,2130., 230., 0., 0., 0./
1723 data ((dectab(i,j),i=1,7),j=1015,1021)/
1724c *----------dl0(1700)--------------------------
1725 * 1226., .10000,1220., 110., 0., 0., 0.
1726 *,1226., .15000,1120.,-120., 0., 0., 0.
1727 *,1226., .39000,2221., 120., 0., 0., 0.
1728 *,1226., .43000,1221., 110., 0., 0., 0.
1729 *,1226., .75000,1121.,-120., 0., 0., 0.
1730 *,1226., .91500,1220., 120.,-120., 0., 0.
1731 *,1226.,1.00000,1120.,-120., 110., 0., 0./
1732 data ((dectab(i,j),i=1,7),j=1022,1035)/
1733c *----------n*0(1710)--------------------------
1734 * 1227., .04430,1220., 110., 0., 0., 0.
1735 *,1227., .13290,1120.,-120., 0., 0., 0.
1736 *,1227., .23790,2221., 120., 0., 0., 0.
1737 *,1227., .30790,1221., 110., 0., 0., 0.
1738 *,1227., .34290,1121.,-120., 0., 0., 0.
1739 *,1227., .41190,1220., 120.,-120., 0., 0.
1740 *,1227., .48090,1220., 110., 110., 0., 0.
1741 *,1227., .54990,1120.,-120., 110., 0., 0.
1742 *,1227., .66070,1220., 220., 0., 0., 0.
1743 *,1227., .72800,2130., 230., 0., 0., 0.
1744 *,1227., .73870,1230., 230., 0., 0., 0.
1745 *,1227., .76000,2230., 130., 0., 0., 0.
1746 *,1227., .92000,1120.,-121., 0., 0., 0.
1747 *,1227.,1.00000,1220., 111., 0., 0., 0./
1748 data ((dectab(i,j),i=1,7),j=1036,1046)/
1749c *----------dl0(1925)--------------------------
1750 * 1228., .18700,1220., 110., 0., 0., 0.
1751 *,1228., .28000,1120.,-120., 0., 0., 0.
1752 *,1228., .36400,2221., 120., 0., 0., 0.
1753 *,1228., .37800,1221., 110., 0., 0., 0.
1754 *,1228., .49000,1121.,-120., 0., 0., 0.
1755 *,1228., .55700,1220., 111., 0., 0., 0.
1756 *,1228., .69000,1120.,-121., 0., 0., 0.
1757 *,1228., .69350,2230., 130., 0., 0., 0.
1758 *,1228., .70000,1230., 230., 0., 0., 0.
1759 *,1228., .80000,1122.,-120., 0., 0., 0.
1760 *,1228.,1.00000,1222., 110., 0., 0., 0./
1761c *---------------------------------------------
1762c * lambda resonances + sigma0 resonances
1763c *---------------------------------------------
1764 data ((dectab(i,j),i=1,7),j=1047,1059)/
1765c *----------lambda(1405)-----------------------
1766 * 1233., .33000,1230., 110., 0., 0., 0.
1767 *,1233., .66000,2230., 120., 0., 0., 0.
1768 *,1233.,1.00000,1130.,-120., 0., 0., 0.
1769c *----------lambda(1520)-----------------------
1770 *,1234., .22500,1120.,-130., 0., 0., 0.
1771 *,1234., .48000,1220.,-230., 0., 0., 0.
1772 *,1234., .62000,1230., 110., 0., 0., 0.
1773 *,1234., .76000,2230., 120., 0., 0., 0.
1774 *,1234., .90000,1130.,-120., 0., 0., 0.
1775 *,1234., .96000,2130., 120.,-120., 0., 0.
1776 *,1234., .99000,2130., 110., 110., 0., 0.
1777 *,1234., .99330,1130.,-120., 110., 0., 0.
1778 *,1234., .99660,2230., 120., 110., 0., 0.
1779 *,1234.,1.00000,1230., 120.,-120., 0., 0./
1780 data ((dectab(i,j),i=1,7),j=1060,1075)/
1781c *----------lambda(1645)-----------------------
1782 * 1235., .10000,1120.,-130., 0., 0., 0.
1783 *,1235., .20000,1220.,-230., 0., 0., 0.
1784 *,1235., .35000,1230., 110., 0., 0., 0.
1785 *,1235., .50000,2230., 120., 0., 0., 0.
1786 *,1235., .65000,1130.,-120., 0., 0., 0.
1787 *,1235., .75000,2130., 120.,-120., 0., 0.
1788 *,1235., .80000,2130., 110., 110., 0., 0.
1789 *,1235., .84500,1130.,-120., 110., 0., 0.
1790 *,1235., .89000,2230., 120., 110., 0., 0.
1791 *,1235., .93500,1230., 120.,-120., 0., 0.
1792 *,1235.,1.00000,2130., 220., 0., 0., 0.
1793c *----------sigma0(1665)-----------------------
1794 *,1236., .10000,1120.,-130., 0., 0., 0.
1795 *,1236., .20000,1220.,-230., 0., 0., 0.
1796 *,1236., .40000,2230., 120., 0., 0., 0.
1797 *,1236., .60000,1130.,-120., 0., 0., 0.
1798 *,1236.,1.00000,2130., 110., 0., 0., 0./
1799 data ((dectab(i,j),i=1,7),j=1076,1084)/
1800c *----------sigma0(1776)-----------------------
1801 * 1237., .17500,1120.,-130., 0., 0., 0.
1802 *,1237., .35000,1220.,-230., 0., 0., 0.
1803 *,1237., .38750,2230., 120., 0., 0., 0.
1804 *,1237., .42500,1130.,-120., 0., 0., 0.
1805 *,1237., .57500,2130., 110., 0., 0., 0.
1806 *,1237., .60000,2231., 120., 0., 0., 0.
1807 *,1237., .62500,1131.,-120., 0., 0., 0.
1808 *,1237., .75000,1234., 110., 0., 0., 0.
1809 *,1237.,1.00000,1230., 220., 0., 0., 0./
1810 data ((dectab(i,j),i=1,7),j=1085,1094)/
1811c *----------lambda(1845)-----------------------
1812 * 1238., .17000,1120.,-130., 0., 0., 0.
1813 *,1238., .34000,1220.,-230., 0., 0., 0.
1814 *,1238., .44000,1230., 110., 0., 0., 0.
1815 *,1238., .54000,2230., 120., 0., 0., 0.
1816 *,1238., .64000,1130.,-120., 0., 0., 0.
1817 *,1238., .70000,1231., 110., 0., 0., 0.
1818 *,1238., .76000,2231., 120., 0., 0., 0.
1819 *,1238., .82000,1131.,-120., 0., 0., 0.
1820 *,1238., .91000,1120.,-131., 0., 0., 0.
1821 *,1238.,1.00000,1220.,-231., 0., 0., 0./
1822 data ((dectab(i,j),i=1,7),j=1095,1106)/
1823c *----------sigma0(1930)-----------------------
1824 * 1239., .07500,1120.,-130., 0., 0., 0.
1825 *,1239., .15000,1220.,-230., 0., 0., 0.
1826 *,1239., .20000,1121.,-130., 0., 0., 0.
1827 *,1239., .25000,1221.,-230., 0., 0., 0.
1828 *,1239., .32500,1120.,-131., 0., 0., 0.
1829 *,1239., .40000,1220.,-231., 0., 0., 0.
1830 *,1239., .47500,2230., 120., 0., 0., 0.
1831 *,1239., .55000,1130.,-120., 0., 0., 0.
1832 *,1239., .70000,2130., 110., 0., 0., 0.
1833 *,1239., .77500,2231., 120., 0., 0., 0.
1834 *,1239., .85000,1131.,-120., 0., 0., 0.
1835 *,1239.,1.00000,1234., 110., 0., 0., 0./
1836c *---------------------------------------------
1837c * sigma+ resonances
1838c *---------------------------------------------
1839 data ((dectab(i,j),i=1,7),j=1107,1118)/
1840c *----------sigma+(1665)-----------------------
1841 * 1132., .20000,1120.,-230., 0., 0., 0.
1842 *,1132., .40000,1130., 110., 0., 0., 0.
1843 *,1132., .60000,1230., 120., 0., 0., 0.
1844 *,1132.,1.00000,2130., 120., 0., 0., 0.
1845c *----------sigma+(1776)-----------------------
1846 *,1133., .35000,1120.,-230., 0., 0., 0.
1847 *,1133., .38750,1130., 110., 0., 0., 0.
1848 *,1133., .42500,1230., 120., 0., 0., 0.
1849 *,1133., .57500,2130., 120., 0., 0., 0.
1850 *,1133., .60000,1131., 110., 0., 0., 0.
1851 *,1133., .62500,1231., 120., 0., 0., 0.
1852 *,1133., .75000,1234., 120., 0., 0., 0.
1853 *,1133.,1.00000,1130., 220., 0., 0., 0./
1854 data ((dectab(i,j),i=1,7),j=1119,1128)/
1855c *----------sigma+(1930)-----------------------
1856 * 1134., .15000,1120.,-230., 0., 0., 0.
1857 *,1134., .22500,1111.,-130., 0., 0., 0.
1858 *,1134., .25000,1121.,-230., 0., 0., 0.
1859 *,1134., .40000,1120.,-231., 0., 0., 0.
1860 *,1134., .47500,1130., 110., 0., 0., 0.
1861 *,1134., .55000,1230., 120., 0., 0., 0.
1862 *,1134., .70000,2130., 120., 0., 0., 0.
1863 *,1134., .77500,1131., 110., 0., 0., 0.
1864 *,1134., .85000,1231., 120., 0., 0., 0.
1865 *,1134.,1.00000,1234., 120., 0., 0., 0./
1866c *---------------------------------------------
1867c * sigma- resonances
1868c *---------------------------------------------
1869 data ((dectab(i,j),i=1,7),j=1129,1140)/
1870c *----------sigma-(1665)-----------------------
1871 * 2232., .20000,1220.,-130., 0., 0., 0.
1872 *,2232., .40000,2230., 110., 0., 0., 0.
1873 *,2232., .60000,1230.,-120., 0., 0., 0.
1874 *,2232.,1.00000,2130.,-120., 0., 0., 0.
1875c *----------sigma-(1776)-----------------------
1876 *,2233., .35000,1220.,-130., 0., 0., 0.
1877 *,2233., .38750,2230., 110., 0., 0., 0.
1878 *,2233., .42500,1230.,-120., 0., 0., 0.
1879 *,2233., .57500,2130.,-120., 0., 0., 0.
1880 *,2233., .60000,2231., 110., 0., 0., 0.
1881 *,2233., .62500,1231.,-120., 0., 0., 0.
1882 *,2233., .75000,1234.,-120., 0., 0., 0.
1883 *,2233.,1.00000,2230., 220., 0., 0., 0./
1884 data ((dectab(i,j),i=1,7),j=1141,1150)/
1885c *----------sigma-(1930)-----------------------
1886 * 2234., .15000,1220.,-130., 0., 0., 0.
1887 *,2234., .17500,1221.,-130., 0., 0., 0.
1888 *,2234., .25000,2221.,-230., 0., 0., 0.
1889 *,2234., .40000,1220.,-131., 0., 0., 0.
1890 *,2234., .47500,2230., 110., 0., 0., 0.
1891 *,2234., .55000,1230.,-120., 0., 0., 0.
1892 *,2234., .70000,2130.,-120., 0., 0., 0.
1893 *,2234., .77500,2231., 110., 0., 0., 0.
1894 *,2234., .85000,1231.,-120., 0., 0., 0.
1895 *,2234.,1.00000,1234.,-120., 0., 0., 0./
1896c *---------------------------------------------
1897c * additional mesonresonances
1898c *---------------------------------------------
1899 data ((dectab(i,j),i=1,7),j=1151,1159)/
1900c *-----------f0(975)---------------------------
1901 * 332., .50000, 120.,-120., 0., 0., 0.
1902 *, 332., .75000, 110., 110., 0., 0., 0.
1903 *, 332., .87500, 130.,-130., 0., 0., 0.
1904 *, 332.,1.00000, 230.,-230., 0., 0., 0.
1905c *-----------a0(980)---------------------------
1906 *, 112., .56000, 110., 220., 0., 0., 0.
1907 *, 112., .78000, 130.,-130., 0., 0., 0.
1908 *, 112.,1.00000, 230.,-230., 0., 0., 0.
1909c *-----------a+(980)---------------------------
1910 *, 122., .60000, 120., 220., 0., 0., 0.
1911 *, 122.,1.00000, 130.,-230., 0., 0., 0./
1912c *---------------------------------------------
1913c * weak baryon decays
1914c *---------------------------------------------
1915 data ((dectab(i,j),i=1,7),j=1160,1169)/
1916c *-----------lambda(1116)----------------------
1917 * 2130.,0.64200,1120.,-120., 0., 0., 0.
1918 *,2130.,1.00000,1220., 110., 0., 0., 0.
1919c *-----------sigma+(1180)----------------------
1920 *,1130.,0.51580,1120., 110., 0., 0., 0.
1921 *,1130.,1.00000,1220., 120., 0., 0., 0.
1922c *-----------sigma-(1180)----------------------
1923 *,2230.,1.00000,1220.,-120., 0., 0., 0.
1924c *---------kaskade-(1360)----------------------
1925 *,2330.,1.00000,2130.,-120., 0., 0., 0.
1926c *---------kaskade0(1360)----------------------
1927 *,1330.,1.00000,2130., 110., 0., 0., 0.
1928c *---------omega-(1680)------------------------
1929 *,3331.,0.68000,2130.,-130., 0., 0., 0.
1930 *,3331.,0.82000,1330.,-120., 0., 0., 0.
1931 *,3331.,1.00000,2330., 110., 0., 0., 0./
1932c *---------------------------------------------
1933c * weak meson decays
1934c *---------------------------------------------
1935 data ((dectab(i,j),i=1,7),j=1170,1171)/
1936c *-----------k0s()--------------------------
1937 * 20., .68610, 120.,-120., 0., 0., 0.
1938 *, 20.,1.00000, 110., 110., 0., 0., 0./
1939 data ((dectab(i,j),i=1,7),j=1172,1192)/
1940c *-----------k0l-------------------------------
1941 * 320., .2113, 110., 110., 110., 0., 0.
1942 *, 320., .2113, 110., 110., 110., 0., 0.
1943 *, 320., .2120, 110., 110., 110., 0., 0.
1944 *, 320., .3380, 120.,-120., 110., 0., 0.
1945 *, 320., .4744, 120., 14., -13., 0., 0.
1946 *, 320., .6108,-120., -14., 13., 0., 0.
1947 *, 320., .8054, 120., 12., -11., 0., 0.
1948 *, 320.,1.0000,-120., -12., 11., 0., 0.
1949c *-----------k+-------------------------------
1950 *, 130., .6352 , -14., 13., 0., 0., 0.
1951 *, 130., .8468 , 120., 110., 0., 0., 0.
1952 *, 130., .9027 , 120., 120.,-120., 0., 0.
1953 *, 130., .92 , 120., 110., 110., 0., 0.
1954 *, 130., .9518 , 110., -14., 13., 0., 0.
1955 *, 130.,1. , 110., -12., 11., 0., 0.
1956c *-----------pi+------------------------------
1957 *, 120., 1. , -14., 13., 0., 0., 0.
1958c *-----------mu-------------------------------
1959 *, 14., 1. , 12., -11., 13., 0., 0.
1960c *-----------etac-------------------------------
1961 *, 440., .32 , 230.,-230., 110., 0., 0.
1962 *, 440., .64 , 220., 110., 110., 0., 0.
1963 *, 440., .76 , 120.,-120., 130.,-130., 0.
1964 *, 440., .88 , 120.,-120., 120.,-120., 0.
1965 *, 440., 1. , 130.,-130., 130.,-130., 0./
1966
1967 call idresi
1968
1969c determine wmass2,wgam2
1970c ----------------------
1971 alfa=1./137.036
1972 gf=1.16570e-5
1973 sin2w=.215
1974 sinw=sqrt(sin2w)
1975c cosw=sqrt(1.-sin2w) !?????????????????unused
1976 amw=sqrt(pi*alfa/(.9304*sqrt(2.)*gf))/sinw
1977 wmass2=amw
1978 call idmass(5,amlep5)
1979 call idmass(6,amlep6)
1980 ngam=12
1981 if(amlep5+amlep6.gt.amw) ngam=9
1982 wgam2=gf*amw**3/(6.*pi*sqrt(2.))*ngam
1983
1984 data iblank/' '/
1985 ird=0
1986 do 1 i=1,mxlook
19871 look(i)=0
1988 do 2 i=1,mxdky
1989 do 3 j=1,5
19903 mode(j,i)=0
19912 cbr(i)=0.
1992 nodcay=.false.
1993 noeta=.false.
1994 nopi0=.false.
1995 nonunu=.false.
1996 noevol=.false.
1997 nohadr=.false.
1998 if(lprint) write(ifch,10)
199910 format('1',30('*')/' *',28x,'*'/
2000 1' *',5x,'isajet decay table',5x,'*'/
2001 2' *',28x,'*'/' ',30('*')//
2002 36x,'part',18x,'decay mode',19x,'cum br',15x,'ident',17x,
2003 4'decay ident')
2004 loop=0
2005 iold=0
2006 if(nodcay) return
2007
2008200 loop=loop+1
2009 if(loop.gt.mxdky) goto9999
2010220 do 210 i=1,5
2011 imode(i)=0
2012 lmode(i)=iblank
2013210 continue
2014 ird=ird+1
2015 if(ird.gt.ndectb)return
2016c if(ird.gt.1171)return ! ??????????????????????????
2017 ires=nint(dectab(1,ird))
2018 br=dectab(2,ird)
2019 do 215 i=1,5
2020215 imode(i)=nint(dectab(2+i,ird))
2021 if(nopi0.and.ires.eq.110) goto220
2022 if(noeta.and.ires.eq.220) goto220
2023 if(ires.eq.iold) goto230
2024 if(ires.lt.0.or.ires.gt.mxlook)
2025 *call utstop('hdecin: ires out of range&')
2026 look(ires)=loop
2027230 iold=ires
2028 cbr(loop)=br
2029 do 240 i=1,5
2030 mode(i,loop)=imode(i)
2031 if(imode(i).ne.0) lmode(i)=idlabl(imode(i))
2032240 continue
2033 lres=idlabl(ires)
2034 if(lprint) write(ifch,20) lres,(lmode(k),k=1,5),
2035 1br,ires,(imode(k),k=1,5)
203620 format(6x,a5,6x,5(a5,2x),3x,f8.5,15x,i5,4x,5(i5,2x))
2037 goto200
2038
20399999 write(ifch,*)'loop=', loop
2040 call utstop('hdecin: loop > mxdky&')
2041
2042 end
2043
2044C -----------------------------------------------
2045 FUNCTION FM(NQ,S0,S1,S2,S3)
2046C -----------------------------------------------
2047C Normalized TRANSITION MATRIX FOR THE DALIZT PLOT DISTRI.
2048C OF K -> 3 PIONS. PARAMETRIZATION OF WEINBERG
2049C AS DESCRIBE IN PARTICLE DATA BOOK.
2050C G IS THE LINEAR COEFFICIENT (SLOPE g)
2051C H IS THE QUADRATIC COEFFICIENT h
2052C D IS THE QUADRATIC COEFFICIENT k
2053C Amax is the maximum of this amplitude (taken from Corsika by D. Heck)
2054C NQ is the decay channel :
2055C 1 - K -> 3 Pi
2056C 2 - K -> Pi + 2 Pi0
2057C 3 - K0 -> 3 Pi0
2058C 4 - K0 -> 2 Pi + Pi0
2059C -----------------------------------------------
2060 DIMENSION G(4),H(4),D(4),Amax(4)
2061 PARAMETER (PIM=139.57E-3)
2062 DATA G/-0.2154,0.594,0.,0.67/
2063 DATA H/0.01,0.035,0.,0.079/
2064 DATA D/-0.01,0.,0.,0.0098/
2065 DATA Amax/1.27,1.84,1.,2.22/
2066
2067 FM=1.+G(NQ)*(S3-S0)/(PIM*PIM)+H(NQ)*((S3-S0)/(PIM*PIM))**2
2068 *+D(NQ)*((S2-S1)/(PIM*PIM))**2
2069 FM=FM/Amax(NQ)
2070
2071 RETURN
2072 END
2073C -----------------------------------------------
2074 FUNCTION FML(N,AM,RM1,RM2,E1S,E2S,E3S)
2075C -----------------------------------------------
2076C Normalized DALITZ PLOT DENSITY (RHO)
2077C OF K -> 1 PION + 2 LEPTONS
2078C AS DESCRIBE IN PARTICLE DATA BOOK.
2079C CLP IS THE LAMBDA + FORM FACTOR COEFFICIENT
2080C CLN IS THE LAMBDA 0 FORM FACTOR COEFFICIENT
2081C EEP IS E'pion
2082C GP IS THE F+(t) FORM FACTOR (t=AM*AM+SM1-2.D0*AM*E1S)
2083C H IS EPS(t)=F-(t)/F+(t) WHERE F- IS CALCULATED FROM F0
2084C Amax is the maximum of this density (taken from Corsika by D. Heck)
2085C N is the decay channel :
2086C 1 - K -> Pi0 + e + Nu
2087C 2 - K -> Pi0 + Mu + Nu
2088C 3 - K0 -> Pi + e + Nu
2089C 4 - K0 -> Pi + Mu + Nu
2090C -----------------------------------------------
2091 DIMENSION CLP(4),CLN(4),Amax(4)
2092 DATA CLP/0.0276,0.031,0.0288,0.034/
2093 DATA CLN/0.0,0.006,0.,0.025/
2094 DATA Amax/1.28e-2,1.194e-2,1.31e-2,1.241e-2/
2095
2096 SM1=RM1*RM1
2097 SM2=RM2*RM2
2098 EEP=0.5D0*(AM*AM+SM1-SM2)/AM-E1S
2099 GP=1.+CLP(N)*(AM*AM+SM1-2.*AM*E1S)/SM1
2100 H=(AM*AM-SM1)/SM1*(CLN(N)-CLP(N))/GP
2101 FML=GP*GP*(AM*(2.*E2S*E3S-AM*EEP)+
2102 *SM2*(0.25*EEP-E3S)+H*SM2*(E3S-0.5*EEP)+
2103 *0.25*H*H*SM2*EEP)
2104 FML=FML/Amax(N)
2105 RETURN
2106 END
2107C -----------------------------------------------
2108 FUNCTION FMU(X)
2109C -----------------------------------------------
2110C PROBABILITY DISTRI. FOR ELECTRON ENERGY FROM MUON DECAY :
2111C MU -> 2NU + E. DESCRIBE IN PARTICLE DATA BOOK.
2112C (SIMPLIFY DIFFERENTIAL DECAY RATE INTEGRATED)
2113C X REDUCED ENERGY OF PARTICLE
2114C -----------------------------------------------
2115
2116 FMU=2.*(3.-2.*X)*X*X
2117
2118 RETURN
2119 END