3 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
9 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
115 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
116 +iibar(210),k1(210),k2(210)
121 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
122 COMMON /nuccc/ jt,jtz,jp,jpz,jjproj,jbproj,jjtarg,jbtarg
124 COMMON /cmhico/ cmhis
126 COMMON /resona/ ireso
128 COMMON /trafop/ gamp,bgamp,betp
130 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
132 COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
133 +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
134 +irvs14, irvv11,irvv12,irvv13,irvv14
136 LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
138 COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
139 +ipadis,ishmal,lpauli
141 COMMON /dnun/ nn,np,
nt
143 COMMON /nshmak/ nnshma,npshma,ntshma,nshmac,nshma2
145 COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
146 * bsite(0:1,200),nstatb,nsiteb
148 COMMON /seaqxx/ seaqx,seaqxn
149 COMMON /diffra/ isingd,idiftp,ioudif,iflagd
150 COMMON /evflag/ numev
151 COMMON /neutyy/ neutyp,neudec
152 COMMON /fluctu/ifluct
153 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
154 *idzre(3),izdre(3),idiqrz(7)
155 COMMON /intneu/ndzsu,nzdsu
159 COMMON /bufues/ bnnvv,bnnss,bnnsv,bnnvs,bnncc,
160 * bnndv,bnnvd,bnnds,bnnsd,
162 * bptvv,bptss,bptsv,bptvs,bptcc,bptdv,
165 * beevv,beess,beesv,beevs,beecc,beedv,
169 * ,bnnzd,bnndz,bptzd,bptdz,beezd,beedz
170 COMMON /ncoucs/ bcouvv,bcouss,bcousv,bcouvs,
171 * bcouzz,bcouhh,bcouds,bcousd,
172 * bcoudz,bcouzd,bcoudi,
173 * bcoudv,bcouvd,bcoucc
174 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
175 * anndv,annvd,annds,annsd,
177 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
179 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
182 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
183 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
184 * acouzz,acouhh,acouds,acousd,
185 * acoudz,acouzd,acoudi,
186 * acoudv,acouvd,acoucc
187 common/popcck/pdbck,pdbse,pdbseu,
188 * ijpock,irejck,ick4,ihad4,ick6,ihad6
189 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
190 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
191 *isea43,isea63,irejao
192 COMMON /nstari/nstart
193 COMMON /ncshxx/ncouxh,ncouxt
195 COMMON /nucros/dsigsu,dsigmc,ndsig
196 COMMON /nncms/ gamcm,bgcm,umo,pcm,eproj,pproj
197 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
198 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
199 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
200 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
201 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
202 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
203 COMMON /hdjase/nhse1,nhse2,nhse3,nhase1,nhase2,nhase3
204 COMMON /casadi/casaxx,icasad
205 COMMON /vxsvd/vxsp(50),vxst(50),vxsap(50),vxsat(50),
206 * vxvp(50),vxvt(50),vxdp(50),vxdt(50),
207 * nxsp,nxst,nxsap,nxsat,nxvp,nxvt,nxdp,nxdt
208 dimension vxsss(50,6),vxvvv(50,6),xxxx(50,6)
209 dimension xb(200),bimpp(200)
214 CHARACTER*8 projty,targty
215 COMMON /user1/titled,projty,targty
216 COMMON /user2/cmener,sdfrac,ptlar,istruf,isingx,idubld
217 COMMON /strufu/istrum,istrut
218 COMMON /ptsamp/ isampt
219 common/collis/
s,ijprox,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
220 COMMON /dropjj/dropjt,dropva
221 COMMON /pomtyp/ipim,icon,isig,lmax,mmax,
nmax,difel,difnu
223 COMMON /zentra/ icentr
227 COMMON /taufo/ taufor,ktauge,itauve,incmod
228 common/popcor/pdb,ajsdef
229 COMMON /diquax/amedd,idiqua,idiquu
230 COMMON /colle/nevhad,nvers,ihadrz,nfile
231 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
232 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
233 +prebin,taebin,fermod,etacou
234 COMMON /cronin/cronco,mkcron
235 COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
237 COMMON /secint/isecin
266 OPEN(47,
file=
'GLAUBTAR.DAT',
268 OPEN(37,
file=
'GLAUBCROSSPB.DAT',
306 CALL rluxgo(lux_level,iseed,0,0)
386 CALL
parpt(1,pt1,pt2,ipt,nevt)
388 CALL
dminit(ncases,epn,ppn,ncount,iglaub)
406 IF( iglaub.EQ.1)
THEN
414 WRITE(6,*)
' Printout of important Parameters before DPMJET run.'
415 *,
' Please note for DPMJET input all numbers are floating point!'
416 WRITE(6,*)
'PROJPAR ',ip,ipz
417 WRITE(6,*)
'TARPAR ',it,itz
418 WRITE(6,*)
'MOMENTUM ',ppn
419 WRITE(6,*)
'ENERGY ',epn
420 WRITE(6,*)
'CMENERGY ',umo
421 WRITE(6,*)
'NOFINALE ',ifinal
422 WRITE(6,*)
'EVAPORAT ',ievap
423 WRITE(6,*)
'OUTLEVEL ',ipri,ipev,ippa,ipco,
init,iphkk
424 auauau=rd2out(iseed1,iseed2)
425 WRITE(6,*)
'RANDOMIZ ',iseed1,iseed2,
' Initial RNDM (RM48) seeds'
426 WRITE(6,*)
'STRUCFUN ',istruf+100*istrut
427 WRITE(6,*)
'SAMPT ',isampt
428 WRITE(6,*)
'SELHARD ',0,iophrd, 0,dropjt,ptthr,ptthr2
429 WRITE(6,*)
'SIGMAPOM ',0,isig,ipim+10*icon,imax,mmax,
nmax
430 WRITE(6,*)
'PSHOWER ',ipshow
431 WRITE(6,*)
'CENTRAL ',icentr
432 WRITE(6,*)
'CMHISTO ',cmhis
433 WRITE(6,*)
'SEASU3 ',seasq
434 WRITE(6,*)
'RECOMBIN ',irecom
435 WRITE(6,*)
'SINGDIFF ',isingd
436 WRITE(6,*)
'TAUFOR ',taufor,ktauge,itauve
437 WRITE(6,*)
'POPCORN ',pdb
438 WRITE(6,*)
'POPCORCK ',ijpock,pdbck
439 WRITE(6,*)
'POPCORSE ',pdbse,pdbseu
440 WRITE(6,*)
'CASADIQU ',icasad,casaxx
441 WRITE(6,*)
'DIQUARKS ',idiqua,idiquu,amedd
442 WRITE(6,*)
'HADRONIZ ',ihadrz
443 WRITE(6,*)
'INTPT ',intpt
444 WRITE(6,*)
'PAULI ',lpauli
445 WRITE(6,*)
'FERMI ',fermp,fermod
446 WRITE(6,*)
'CRONINPT ',mkcron,cronco
447 WRITE(6,*)
'SEADISTR ',xseacu+0.95d0,unon,unom,unosea
448 WRITE(6,*)
'SEAQUARK ',seaqx,seaqxn
449 WRITE(6,*)
'SECINTER ',isecin
450 WRITE(6,*)
'XCUTS ',cvq,cdq,csea,ssmima
451 WRITE(6,*)
' Printout of important Parameters before DPMJET run.'
452 *,
' Please note for DPMJET input all numbers are floating point!'
458 ndone=(iiii-1)*ncaset
460 WRITE(6,
'(A,4I5)')
' KKINC call ',iit,iitz,iip,iipz
462 1111
FORMAT(
' NDONE= ',i10)
465 numev = i+(iiii-1)*ncaset
466 IF ((i.EQ.486).OR.(i.EQ.803).OR.(i.EQ.1368).OR.
467 & (i.EQ.1465).OR.(i.EQ.1693).OR.(i.EQ.1808))
THEN
535 CALL
kkinc(epn,iit,iitz,iip,iipz,iiproj,kkmat,
536 * iitarg,nhkkh1,irej)
538 IF(irej.EQ.1)go to 765
615 WRITE(6,
'(A,4I5)')
' KKINC call ',iit,iitz,iip,iipz
617 WRITE(6,
'(2I4,I6,4I4,5F10.2,2I3,I2,I4)')kk,isthkk(kk),idhkk(kk),
618 *jmohkk(1,kk),jmohkk(2,kk),jdahkk(1,kk),jdahkk(2,kk),
619 *(phkk(ll,kk),ll=1,5),idres(kk),idxres(kk),nobam(kk),idbam(kk)
625 WRITE(6,
'(A,4I5)')
' KKINC call ',iit,iitz,iip,iipz
627 WRITE(6,
'(2I4,I6,4I4,5F10.2,2I3,I2,I4)')kk,isthkk(kk),idhkk(kk),
628 *jmohkk(1,kk),jmohkk(2,kk),jdahkk(1,kk),jdahkk(2,kk),
629 *(phkk(ll,kk),ll=1,5),idres(kk),idxres(kk),nobam(kk),idbam(kk)
637 WRITE(6,*)
' Printout of important Parameters after DPMJET run.'
638 *,
' Please note for DPMJET input all numbers are floating point!'
639 WRITE(6,*)
'PROJPAR ',ip,ipz
640 WRITE(6,*)
'TARPAR ',it,itz
641 WRITE(6,*)
'MOMENTUM ',ppn
642 WRITE(6,*)
'ENERGY ',epn
643 WRITE(6,*)
'CMENERGY ',umo
644 WRITE(6,*)
'NOFINALE ',ifinal
645 WRITE(6,*)
'EVAPORAT ',ievap
646 WRITE(6,*)
'OUTLEVEL ',ipri,ipev,ippa,ipco,
init,iphkk
647 auauau=rd2out(iseed1,iseed2)
648 WRITE(6,*)
'RANDOMIZ ',iseed1,iseed2,
' Final RNDM (RM48) seeds'
649 WRITE(6,*)
'STRUCFUN ',istruf+100*istrut
650 WRITE(6,*)
'SAMPT ',isampt
651 WRITE(6,*)
'SELHARD ',0,iophrd, 0,dropjt,ptthr,ptthr2
652 WRITE(6,*)
'SIGMAPOM ',0,isig,ipim+10*icon,imax,mmax,
nmax
653 WRITE(6,*)
'PSHOWER ',ipshow
654 WRITE(6,*)
'CENTRAL ',icentr
655 WRITE(6,*)
'CMHISTO ',cmhis
656 WRITE(6,*)
'SEASU3 ',seasq
657 WRITE(6,*)
'RECOMBIN ',irecom
658 WRITE(6,*)
'SINGDIFF ',isingd
659 WRITE(6,*)
'TAUFOR ',taufor,ktauge,itauve
660 WRITE(6,*)
'POPCORN ',pdb
661 WRITE(6,*)
'POPCORCK ',ijpock,pdbck
662 WRITE(6,*)
'POPCORSE ',pdbse,pdbseu
663 WRITE(6,*)
'CASADIQU ',icasad,casaxx
664 WRITE(6,*)
'DIQUARKS ',idiqua,idiquu,amedd
665 WRITE(6,*)
'HADRONIZ ',ihadrz
666 WRITE(6,*)
'INTPT ',intpt
667 WRITE(6,*)
'PAULI ',lpauli
668 WRITE(6,*)
'FERMI ',fermp,fermod
669 WRITE(6,*)
'CRONINPT ',mkcron,cronco
670 WRITE(6,*)
'SEADISTR ',xseacu+0.95d0,unon,unom,unosea
671 WRITE(6,*)
'SEAQUARK ',seaqx,seaqxn
672 WRITE(6,*)
'SECINTER ',isecin
673 WRITE(6,*)
'XCUTS ',cvq,cdq,csea,ssmima
674 WRITE(6,*)
' Printout of important Parameters after DPMJET run.'
675 *,
' Please note for DPMJET input all numbers are floating point!'
681 auauau=rd2out(iseed1,iseed2)
682 WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
683 WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
684 WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
685 WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
697 WRITE(6,1100) irvv11,irvv12,irvv13,irvv14, irsv11,irsv12,irsv13,
698 + irsv14, irvs11,irvs12,irvs13,irvs14, irss11,irss12,irss13,irss14
699 1100
FORMAT (
' REJECTION COUNTERS FROM KKEVT',/, 5
x,
' V-V CHAINS',4i6/
700 +5
x,
' S-V CHAINS',4i6/ 5
x,
' V-S CHAINS',4i6/ 5
x,
' S-S CHAINS',4i6)
701 WRITE(6,
'(A,4I10)')
' POPCCK/SE/S3/S0 rejections ',
702 * irejck,irejse,irejs3,irejs0
703 WRITE(6,
'(A,4I10)')
' POPCCK/ASE/AS3/AS0 rejections ',
704 * irejsa,ireja3,ireja0
705 WRITE(6,
'(2A,8I6)')
' POPCCK ICK4,ICK6,IHAD4,IHAD6,ISE4,ISE6 ',
706 *
'ISE43,ISE63 ', ick4,ick6,ihad4,ihad6,ise4,ise6,ise43,ise63
707 WRITE(6,
'(2A,8I6)')
' POPCSAQ IHADA4,IHADA6,ISEA4,ISEA6 ',
708 *
'ISEA43,ISEA63 ', ihada4,ihada6,isea4,isea6,isea43,isea63
709 WRITE(6,*)
' NDVUU,NDVUS,NDVSS,NVDUU,NVDUS,NVDSS',
710 *
' NDSUU,NDSUS,NDSSS,NSDUU,NSDUS,NSDSS',
711 *
' NDZUU,NDZUS,NDZSS,NZDUU,NZDUS,NZDSS' ,
712 * ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
713 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
714 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
715 WRITE(6,*)
' NADVUU,NADVUS,NADVSS,NAVDUU,NAVDUS,NAVDSS',
716 *
' NADSUU,NADSUS,NADSSS,NASDUU,NASDUS,NASDSS',
717 *
' NADZUU,NADZUS,NADZSS,NAZDUU,NAZDUS,NAZDSS' ,
718 * nadvuu,nadvus,nadvss,navduu,navdus,navdss,
719 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
720 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
721 WRITE(6,*)
' NHSE1,NHSE2,NHSE3,NHASE1,NHASE2,NHASE3 ',
722 * nhse1,nhse2,nhse3,nhase1,nhase2,nhase3
729 WRITE(6,
'(A/7I8)')
' Diquark rejection IDIQRE(1-7),N,ss,su,ud',
730 & (idiqre(jj),jj=1,7)
731 WRITE(6,
'(A/7I8)')
' Diquark rejection IDIQRZ(1-7),N,ss,su,ud',
732 & (idiqrz(jj),jj=1,7)
733 WRITE(6,*)
' Diquark rej. IDVRE(1-3),ud,us,ss ',(idvre(jj),jj=1,3)
734 WRITE(6,*)
' Diquark rej. IVDRE(1-3),ud,us,ss ',(ivdre(jj),jj=1,3)
735 WRITE(6,*)
' Diquark rej. IDSRE(1-3),ud,us,ss ',(idsre(jj),jj=1,3)
736 WRITE(6,*)
' Diquark rej. ISDRE(1-3),ud,us,ss ',(isdre(jj),jj=1,3)
737 WRITE(6,*)
' Diquark rej. IDZRE(1-3),ud,us,ss ',(idzre(jj),jj=1,3)
738 WRITE(6,*)
' Diquark rej. IZDRE(1-3),ud,us,ss ',(izdre(jj),jj=1,3)
739 WRITE(6,*)
' NDZSU,NZDSU ',ndzsu,nzdsu
741 IF ((cmhis.EQ.1.d0).AND.(ioudif.EQ.1))
745 WRITE(6,*)
' Output of x-distribution survey',
746 *
'VXSP(II),VXST(II),VXSAP(II),VXSAT(II),',
747 *
'VXVP(II),VXVT(II),VXDP(II),VXDT(II)' ,
748 *nxsp,nxst,nxsap,nxsat,nxvp,nxvt,nxdp,nxdt
750 IF(nxsp.GE.1)vxsp(ii)=50.d0*vxsp(ii)/nxsp
751 IF(nxst.GE.1)vxst(ii)=50.d0*vxst(ii)/nxst
752 IF(nxsap.GE.1)vxsap(ii)=50.d0*vxsap(ii)/nxsap
753 IF(nxsat.GE.1)vxsat(ii)=50.d0*vxsat(ii)/nxsat
754 IF(nxvp.GE.1)vxvp(ii)=50.d0*vxvp(ii)/nxvp
755 IF(nxvt.GE.1)vxvt(ii)=50.d0*vxvt(ii)/nxvt
756 IF(nxdp.GE.1)vxdp(ii)=50.d0*vxdp(ii)/nxdp
757 IF(nxdt.GE.1)vxdt(ii)=50.d0*vxdt(ii)/nxdt
758 xxxxx=ii*0.02d0-0.01d0
765 fxvvv=(1.-xxxxx)**3/
sqrt(xxxxx)
766 fxddd=2.d0*xxxxx**3.0d0/
sqrt(1.d0-xxxxx)
767 vxsss(ii,1)=log10(vxsp(ii))
768 vxsss(ii,2)=log10(vxst(ii))
769 vxsss(ii,3)=log10(vxsap(ii))
770 vxsss(ii,4)=log10(vxsat(ii))
771 vxvvv(ii,1)=log10(vxvp(ii))
772 vxvvv(ii,2)=log10(vxvt(ii))
773 vxvvv(ii,3)=log10(vxdp(ii))
774 vxvvv(ii,4)=log10(vxdt(ii))
775 vxvvv(ii,5)=log10(fxvvv)
776 vxvvv(ii,6)=log10(fxddd)
777 axsp=axsp+0.02d0*vxsp(ii)*xxxxx
778 axst=axst+0.02d0*vxst(ii)*xxxxx
779 axsap=axsap+0.02d0*vxsap(ii)*xxxxx
780 axsat=axsat+0.02d0*vxsat(ii)*xxxxx
781 axvp=axvp+0.02d0*vxvp(ii)*xxxxx
782 axvt=axvt+0.02d0*vxvt(ii)*xxxxx
783 axdp=axdp+0.02d0*vxdp(ii)*xxxxx
784 axdt=axdt+0.02d0*vxdt(ii)*xxxxx
785 WRITE(6,*)vxsp(ii),vxst(ii),vxsap(ii),vxsat(ii),
786 * vxvp(ii),vxvt(ii),vxdp(ii),vxdt(ii)
789 *axsp,axst,axsap,axsat,axvp,axvt,axdp,axdt
790 CALL
plot(xxxx,vxsss,200,4,50,0.d0,0.02d0,-3.d0,0.05d0)
791 CALL
plot(xxxx,vxvvv,300,6,50,0.d0,0.02d0,-3.d0,0.05d0)
795 CALL
parpt(3,pt1,pt2,ipt,ncases)
799 fracxs=float(ncouxh)/(float(ncouxh)+float(ncouxt))
801 WRITE(6,*)
' Fraction of x-sect: ',fracxs,ncouxh,ncouxt
806 IF(ndsig.GE.1) dsigmc=dsigsu/ndsig
807 WRITE(6,*)
' Neutrino-nucleon cross section DSIGMC,NDSIG ',
808 & dsigmc,
' *10**(-38) cm**2 ',ndsig,
' evts'
841 IF(ishmal) CALL
shmak(3,nshmac,np,
nt,ip,it,umo,bimp)
842 IF(ishmal) CALL
shmak1(3,nshma2,np,
nt,ip,it,umo,bimp)
844 IF (ireso.EQ.1) CALL
distrp(3,ncases,ppn)
845 IF (cmhis.EQ.0.d0) CALL
distr(3,ncases,ppn,idummy)
846 IF (cmhis.EQ.1.d0) CALL
distrc(3,ncases,ppn,idummy)
847 IF (cmhis.EQ.2.d0) CALL
distco(3,ncases,ppn,idummy)
848 IF (ireso.EQ.1) CALL
disres(3,ncases,ppn)
851 IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN
852 CALL
plomb(5,pp,char,xfxfxf,itif,ijproj)
854 IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN
855 CALL
plombc(5,pp,char,xfxfxf,itif,ijproj)
866 SUBROUTINE dminit(NCASES,EPN,PPN,NCOUNT,IGLAUB)
867 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
979 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
1077 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
1078 +iibar(210),k1(210),k2(210)
1087 COMMON /paname/ btype(30)
1089 COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
1090 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
1092 COMMON /factmo/ ifacto
1094 COMMON /taufo/ taufor,ktauge,itauve,incmod
1096 COMMON /rptshm/ rproj,rtarg,bimpac
1098 COMMON /trafop/ gamp,bgamp,betp
1100 COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
1101 +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
1102 +irvs14, irvv11,irvv12,irvv13,irvv14
1104 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
1106 COMMON /dnun/ nn,np,
nt
1108 COMMON /nshmak/ nnshma,npshma,ntshma,nshmac,nshma2
1110 COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
1111 * bsite(0:1,200),nstatb,nsiteb
1113 COMMON /hadthr/ ehadth,inthad
1117 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
1118 COMMON /nuccc/ jt,jtz,jp,jpz,jjproj,jbproj,jjtarg,jbtarg
1120 COMMON /zentra/ icentr
1122 COMMON /cmhico/ cmhis
1124 COMMON /resona/ ireso
1126 LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
1128 COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
1129 +ipadis,ishmal,lpauli
1131 COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
1134 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
1135 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
1136 +prebin,taebin,fermod,etacou
1142 COMMON /projk/ iprojk
1144 parameter(lunber=14)
1147 COMMON /seaqxx/ seaqx,seaqxn
1148 COMMON /cronin/cronco,mkcron
1150 COMMON /seadiq/lseadi
1151 COMMON /final/ifinal
1152 COMMON /recom/irecom
1154 COMMON /neutyy/ neutyp,neudec
1155 COMMON /nstari/nstart
1156 common/popcor/pdb,ajsdef
1157 common/popcck/pdbck,pdbse,pdbseu,
1158 * ijpock,irejck,ick4,ihad4,ick6,ihad6
1159 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
1160 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
1161 *isea43,isea63,irejao
1163 COMMON /bufues/ bnnvv,bnnss,bnnsv,bnnvs,bnncc,
1164 * bnndv,bnnvd,bnnds,bnnsd,
1166 * bptvv,bptss,bptsv,bptvs,bptcc,bptdv,
1167 * bptvd,bptds,bptsd,
1169 * beevv,beess,beesv,beevs,beecc,beedv,
1170 * beevd,beeds,beesd,
1172 * ,bnndi,bptdi,beedi
1173 * ,bnnzd,bnndz,bptzd,bptdz,beezd,beedz
1174 COMMON /ncoucs/ bcouvv,bcouss,bcousv,bcouvs,
1175 * bcouzz,bcouhh,bcouds,bcousd,
1176 * bcoudz,bcouzd,bcoudi,
1177 * bcoudv,bcouvd,bcoucc
1178 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
1179 * anndv,annvd,annds,annsd,
1181 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
1183 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
1186 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
1187 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
1188 * acouzz,acouhh,acouds,acousd,
1189 * acoudz,acouzd,acoudi,
1190 * acoudv,acouvd,acoucc
1191 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
1192 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
1193 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
1194 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
1195 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
1196 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
1197 COMMON /hdjase/nhse1,nhse2,nhse3,nhase1,nhase2,nhase3
1199 COMMON /diffra/ isingd,idiftp,ioudif,iflagd
1200 COMMON /seasu3/seasq
1201 COMMON /ifragm/ifrag
1202 COMMON /fluctu/ifluct
1203 COMMON /diquax/amedd,idiqua,idiquu
1204 COMMON /vxsvd/vxsp(50),vxst(50),vxsap(50),vxsat(50),
1205 * vxvp(50),vxvt(50),vxdp(50),vxdt(50),
1206 * nxsp,nxst,nxsap,nxsat,nxvp,nxvt,nxdp,nxdt
1209 COMMON /nncms/ gamcm,bgcm,umo,pcm,eproj,pproj
1212 COMMON /xsecpt/ ptcut,sigs,dsigh
1213 COMMON /kglaub/jglaub
1214 COMMON /xsecnu/ecmuu,ecmoo,ngritt,nevtt
1258 CHARACTER*8 projty,targty
1261 COMMON /user1/titled,projty,targty
1262 COMMON /user2/cmener,sdfrac,ptlar,istruf,isingx,idubld
1286 COMMON /colle/nevhad,nvers,ihadrz,nfile
1298 common/collis/
s,ijprox,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
1307 common/booklt/btypex(30),nbook(30)
1325 COMMON /polmn0/pdifr,phard,psoft,alfah,betah,
1326 * sigtot,sigqel,sigel,sigine,sighin,sigd,sigdd
1333 COMMON /pomtyp/ipim,icon,isig,lmax,mmax,
nmax,difel,difnu
1340 COMMON /dropjj/dropjt,dropva
1341 COMMON /gluspl/nugluu,nsgluu
1344 COMMON /ptlarg/xsmax
1346 COMMON /ptsamp/ isampt
1347 COMMON /stars/istar2,istar3
1350 COMMON /strufu/istrum,istrut
1351 COMMON /cutofn/ncutox
1354 COMMON /harlun/ qlun,iharlu
1355 COMMON /pomtab/ipomta
1356 COMMON /sincha/isichaa
1358 COMMON /evappp/ievap
1360 parameter( frdiff = 0.2
d+00 )
1361 parameter( ethsea = 1.0
d+00 )
1363 LOGICAL ldiffr, linctv, levprt, lheavy, ldeexg, lgdhpr, lpreex,
1364 & lhlfix, lprfix, lparwv, lpower, lsngch, llvmod, lschdf
1365 COMMON / parevt / dpower, fsprd0, fshpfn, rn1gsc, rn2gsc,
1366 & ldiffr(39),lpower, linctv, levprt, lheavy,
1367 & ldeexg, lgdhpr, lpreex, lhlfix, lprfix, lparwv,
1368 & ilvmod, jlvmod, llvmod, lsngch, lschdf
1378 parameter(
mxpsst = 300 )
1379 parameter(
mxpsfb = 41000 )
1380 LOGICAL lfrmbk, lncmss
1381 COMMON / frbkcm / amufbk, eexfbk(
mxpsst), amfrbk(
mxpsst),
1383 & exmxfb, r0frbk, r0cfbk, c1cfbk, c2cfbk,
1388 & ifbcha(5,
mxpsfb), iposst, iposfb, ifbstf,
1389 & ifbfrb, nbufbk, lfrmbk, lncmss
1391 COMMON /inpflg/ iang,ifiss,ib0,igeom,istrag,keydk
1392 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
1393 COMMON /secint/isecin
1394 COMMON /nuclea/ pfermp(2),pfermn(2),fermdd,
1395 & ebindp(2),ebindn(2),epot(2,210),
1397 COMMON /ferfor/iferfo
1398 COMMON /casadi/casaxx,icasad
1399 COMMON /infore/ifrej
1409 CHARACTER*8 code,codewd,blank,sdum
1412 1
'TITLE ',
'PROJPAR ',
'TARPAR ',
'ENERGY ',
'HADRONIZ',
1413 2
'RANDOMIZ',
'FERMI ',
'EVENTAPE',
'START ',
'PARTEV ',
1414 3
'INTPT ',
'TECALBAM',
'RESONANC',
'VALVAL ',
'COMMENT ',
1415 4
'OUTLEVEL',
'LEPTOEVT',
'SEASEA ',
'PARTICLE',
'ALLPART ',
1416 5
'TAUFOR ',
'SEAVAL ',
'VALSEA ',
'MOMENTUM',
'PAULI ',
1417 6
'PROJKASK',
'CENTRAL ',
'SEADISTR',
'CMHISTO ',
'SIGTEST ',
1418 7
'XCUTS ',
'HADRIN ',
'FACTOMOM',
'COULOMB ',
'GLAUBERI',
1419 8
'EDENSITY',
'CMENERGY',
'INFOREJE',
'RECOMBIN',
'SINGDIFF',
1420 9
'NOFINALE',
'SEASU3 ',
'CRONINPT',
'POPCORN ',
'STOP ',
1421 9
'FLUCTUAT',
'DIQUARKS',
'HBOOKHIS',
'GLAUBERA',
'POMTABLE',
1422 9
'SINGLECH',
'HADRINTH',
'EVAPORAT',
'SEAQUARK',
'SECINTER',
1423 9
'POPCORCK',
'CASADIQU',
'POPCORSE',
'NEUTRINO',
'DIFFNUC ',
1424 9
'XSECNUC ',
' ',
' ',
' ',
' '/
1454 IF (ncount.EQ.1)
THEN
1462 1000
FORMAT(
'1 **************************************************',
1463 +
'**************************************************', //
1464 +
' DPMJET VERSION II.5 (Sept. 1999) ' /
1465 +
' DUAL PARTON MODEL FOR HADRON NUCLEUS COLLISIONS '/ /
1466 +
' AND NUCLEUS NUCLEUS COLLISIONS '/
1467 +
' INCLUDING A FORMATION TIME INTRANUCLEAR CASCADE'/
1468 4
' Minijets and DTUJET like multiple soft jets '/
1469 4
' Nuclear evaporation and residual target and '/
1470 4
' projectile nuclei '/
1471 +
' **************************************************',
1472 +
'**************************************************',//)
1574 IF (ihadrz.GE.2)
THEN
1616 ptthr=2.1+0.15*(log10(cmener/50.))**3
1618 ELSEIF(istrut.EQ.2)
THEN
1619 ptthr=2.5+0.12*(log10(cmener/50.))**3
1637 IF(ijproj.NE.0) nnpp=ijproj
1638 epn=
sqrt(ppn**2+aam(nnpp)**2)
1644 pproj =
sqrt((epn-amproj)*(epn+amproj))
1645 umo =
sqrt(amproj**2 + amtar**2 + 2.*amtar*eproj)
1647 gamcm = (eproj+amtar)/umo
1650 pcm=gamcm*pproj - bgcm*eproj
1652 ptthr=2.1+0.15*(log10(cmener/50.))**3
1654 ELSEIF(istrut.EQ.2)
THEN
1655 ptthr=2.5+0.12*(log10(cmener/50.))**3
1679 istruf=istruf-istrut*100
1776 IF(ieof.EQ.1) go to 40
1778 READ(5,1010)codewd,(
what(i),i=1,6),sdum
1779 WRITE(6,1020)codewd,(
what(i),i=1,6),sdum
1782 IF(codewd.EQ.
code(isw)) go to 30
1790 + 50 , 60 , 90 , 120 , 130 ,
1794 + 140 , 150 , 160 , 170 , 210 ,
1798 + 220 , 230 , 240 , 250 , 260 ,
1802 + 280 , 290 , 300 , 310 , 320 ,
1806 + 330 , 340 , 350 , 360 , 370 ,
1810 + 380 , 390 , 400 , 410 , 420 ,
1814 + 430 , 440 , 450 , 460 , 470 ,
1818 + 480 , 490 , 500 , 510 , 520 ,
1822 + 530 , 535 , 538, 539, 540 ,
1826 + 541 , 542 , 543, 544, 545,
1831 + 551 , 552 , 553 , 554 ,555 ,
1835 + 556 , 557 ,558 , 559 ,560,
1839 + 620 , 630 ,640 , 650 ,660,610),isw
1861 1010
FORMAT(a8,2
x,6e10.0,a8 )
1862 1020
FORMAT(
' *****NEXT CONTROL CARD ***** ',a10,6(1
x,g11.4), 2
x,a10)
1864 1030
FORMAT(/,
' UNKNOWN CODEWORD - CONTROL CARD IGNORED')
1865 1040
FORMAT(/,
' UNEXPECTED END OF INPUT - STOP ASSUMED.')
1866 1050
FORMAT(/,
' UNEXPECTED END OF INPUT - START ASSUMED.')
1906 IF (codewd.GT.
'-zzzzzzz')
1907 1
WRITE(6,91) codewd,(
what(i),i=1,6),sdum
1908 91
FORMAT(
' ---- control input card : ----'
1909 1 /1
x,a8,2
x,6(f10.3),a8)
1926 IF(codewd.EQ.
'STRUCFUN')
THEN
1932 istruf=istruf-istrut*100
1934 WRITE(6,*)
' ISTRUF,ISTRUT ',istruf,istrut
1948 ELSEIF(codewd.EQ.
'PSHOWER ')
THEN
1976 1070
FORMAT(//,5
x,a80,//)
1991 IF(sdum.EQ.blank)
THEN
1997 IF(ip.EQ.1) ijproj=1
1998 IF(ip.EQ.1) ijprox=1
2004 IF(ip.EQ.1) jjproj=1
2005 IF(ip.EQ.1) jjprox=1
2008 IF(sdum.EQ.btype(ii))
THEN
2011 ibproj=iibar(ijproj)
2016 jbproj=iibar(ijproj)
2022 WRITE(6,
'(A)')
' WRONG STRUCTURE OF PROJPAR CARD'
2039 IF(sdum.EQ.blank)
THEN
2046 IF(sdum.EQ.btype(ii))
THEN
2056 WRITE(6,
'(A)')
' WRONG STRUCTURE OF TARPAR CARD'
2072 IF(ijproj.NE.0) nnpp=ijproj
2073 ppn=
sqrt((epn-aam(nnpp))*(epn+aam(nnpp)))
2080 pproj =
sqrt((epn-amproj)*(epn+amproj))
2081 umo =
sqrt(amproj**2 + amtar**2 + 2.*amtar*eproj)
2084 ptthr=2.1+0.15*(log10(cmener/50.))**3
2086 ELSEIF(istrut.EQ.2)
THEN
2087 ptthr=2.5+0.12*(log10(cmener/50.))**3
2090 gamcm = (eproj+amtar)/umo
2093 pcm=gamcm*pproj - bgcm*eproj
2095 print 1033, eproj,pproj,
2096 +amproj,amtar,umo,gamcm,bgcm,pcm
2097 1033
FORMAT(
' CMS: ' ,
2098 +
' EPROJ,PPROJ,AMPROJ,AMTAR,UMO,GAMCM,BGCM,PCM'/8e22.13)
2119 IF (ihadrz.GE.2)
THEN
2138 auauau=rd2in(iseed1,iseed2)
2156 IF (
what(1).EQ.1.d0)
THEN
2163 IF(fermod.LT.0.0d0.OR.scafer.GT.2.0d0) scafer=1.0d0
2181 1080
FORMAT (
' THIS FILE CONTAINS EVENTS FROM KKEVT ')
2200 IF (ireso.EQ.1) CALL
distrp(1,ijproj,ppn)
2201 IF (cmhis.EQ.0.d0) CALL
distr(1,ijproj,ppn,idummy)
2202 IF (cmhis.EQ.1.d0) CALL
distrc(1,ijproj,ppn,idummy)
2203 IF (cmhis.EQ.2.d0) CALL
distco(1,ijproj,ppn,idummy)
2204 IF (ireso.EQ.1) CALL
disres(1,nhkkh1,ppn)
2205 IF (ipadis) CALL
distpa(1)
2206 IF (ioudif.EQ.1) CALL
diadif(1,0)
2207 CALL
shmak(1,nn,np,
nt,ip,it,umo,bimp)
2208 CALL
shmak1(1,nn,np,
nt,ip,it,umo,bimp)
2210 WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
2213 * ,form=
'UNFORMATTED')
2217 WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
2228 IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN
2229 CALL
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
2231 IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN
2232 CALL
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
2241 IF(lpauli .AND. (.NOT.fermp))
THEN
2243 +
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
2244 +
' IF FERMI ACTIVE',
' LPAULI CHANGED TO .FALSE.'
2251 IF(nevnts.LE.0) nevnts=1000
2255 IF(ncases.LE.0) ncases=100
2258 IF(iglaub.NE.1) iglaub=0
2261 IF(iglaub.EQ.1)
THEN
2264 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
2267 CALL
shmakf(ip,ipz,it,itz)
2276 IF(ipim.EQ.2)CALL
prblm2(cmener)
2345 ptthr=2.1+0.15*(log10(cmener/50.))**3
2347 ELSEIF(istrut.EQ.2)
THEN
2348 ptthr=2.5+0.12*(log10(cmener/50.))**3
2366 IF (
what(1).EQ.1.d0)
THEN
2384 IF (
what(1).EQ.1.d0)
THEN
2427 IF(
what(1).GT.0.5d0) ireso=1
2440 IF (
what(1).EQ.1.d0)
THEN
2462 270
WRITE(6,1120)
title
2515 OPEN(29,
file=
'lepto.evt',
2519 IF (ireso.EQ.1) CALL
distrp(1,ijproj,ppn)
2520 IF (cmhis.EQ.0.d0) CALL
distr(1,ijproj,ppn,idummy)
2521 IF (cmhis.EQ.1.d0) CALL
distrc(1,ijproj,ppn,idummy)
2522 IF (cmhis.EQ.2.d0) CALL
distco(1,ijproj,ppn,idummy)
2523 IF (ireso.EQ.1) CALL
disres(1,nhkkh1,ppn)
2524 IF (ipadis) CALL
distpa(1)
2525 IF (ioudif.EQ.1) CALL
diadif(1,0)
2527 WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
2530 * ,form=
'UNFORMATTED')
2534 WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
2538 WRITE(6,*)
' NEUTRINO: after INCINI call'
2542 IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN
2543 CALL
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
2545 IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN
2546 CALL
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
2555 IF(lpauli .AND. (.NOT.fermp))
THEN
2557 +
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
2558 +
' IF FERMI ACTIVE',
' LPAULI CHANGED TO .FALSE.'
2564 IF(nevnts.LE.0) nevnts=1000
2568 IF(ncases.LE.0) ncases=100
2574 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
2637 WRITE(6,*)
' NEUTRINO initialization finished'
2654 IF (
what(1).EQ.1.d0)
THEN
2684 IF (
what(1).EQ.1.d0)
THEN
2724 IF (
what(1).EQ.1.d0)
THEN
2741 IF (
what(1).EQ.1.d0)
THEN
2758 IF(ijproj.NE.0) nnpp=ijproj
2759 epn=
sqrt(ppn**2+aam(nnpp)**2)
2765 pproj =
sqrt((epn-amproj)*(epn+amproj))
2766 umo =
sqrt(amproj**2 + amtar**2 + 2.*amtar*eproj)
2769 ptthr=2.1+0.15*(log10(cmener/50.))**3
2771 ELSEIF(istrut.EQ.2)
THEN
2772 ptthr=2.5+0.12*(log10(cmener/50.))**3
2775 gamcm = (eproj+amtar)/umo
2778 pcm=gamcm*pproj - bgcm*eproj
2780 print 1033, eproj,pproj,
2781 +amproj,amtar,umo,gamcm,bgcm,pcm
2795 IF (
what(1).EQ.1.d0)
THEN
2855 IF(unon.LT.0.1d0) unon=2.0
2857 IF(unom.LT.0.1d0) unom=1.5
2859 IF(unosea.LT.0.1d0) unosea=2.0
2905 IF(cvq.LT.0.5d0) cvq=1.0
2907 IF(cdq.LT.1.0d0) cdq=2.0
2909 IF(csea.LT.0.1d0) csea =0.1
2911 IF(ssmima.LT.0.0d0) ssmima=0.14
2913 IF(
what(5).GT.2.0d0) vvmthr=
what(5)
2931 IF(inthad.LT.0 .OR. inthad.GT.2) inthad=0
2932 IF(inthad.EQ.1)
WRITE(6,
'(/5X,A/)')
2933 +
' FHAD: INELASTIC INTERACTION FORCED'
2934 IF(inthad.EQ.2)
WRITE(6,
'(/5X,A/)')
2935 +
' FHAD: ELASTIC INTERACTION FORCED'
2994 Write(47,473)ip,ipz,it,itz
2995 473
FORMAT(
' NUCLEUS ',4i10)
3008 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3011 WRITE(47,
'(4F10.5)') bmax,bstep,rproj,rtarg
3013 WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
3021 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3022 WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
3056 ptthr=2.1+0.15*(log10(cmener/50.))**3
3058 ELSEIF(istrut.EQ.2)
THEN
3059 ptthr=2.5+0.12*(log10(cmener/50.))**3
3066 IF(ijproj.NE.0) nnpp=ijproj
3067 epn=(cmener**2 + aam(nnpp)**2 - aam(1)**2)/(2.*aam(1))
3068 ppn=
sqrt((epn-aam(nnpp))*(epn+aam(nnpp)))
3069 ecm =
sqrt(aam(idp)**2+aam(1)**2+2.0d0*aam(1)*epn)
3074 pproj =
sqrt((epn-amproj)*(epn+amproj))
3075 eproj=
sqrt(pproj**2+amproj**2)
3079 gamcm = (eproj+amtar)/umo
3081 gamcm = (eproj+aam(1))/umo
3083 pcm=gamcm*pproj - bgcm*eproj
3084 print 1033, eproj,pproj,
3085 +amproj,amtar,umo,gamcm,bgcm,pcm
3117 IF (
what(1).EQ.1.d0) irecom=1
3118 IF (
what(1).NE.1.d0) irecom=0
3119 IF (
what(1).EQ.1.d0) lseadi=.true.
3153 IF (
what(1).EQ.1.d0) ifinal=1
3154 IF (
what(1).EQ.0.d0) ifinal=0
3218 IF(ifluct.EQ.1)CALL
fluini
3231 IF(
what(3).GT.0.d0)
THEN
3269 Write(47,1473)ip,ipz,it,itz
3270 1473
FORMAT(
' NUCLEUS ',4i10)
3279 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3282 WRITE(47,
'(4F10.5)') bmax,bstep,rproj,rtarg
3284 WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
3292 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3293 WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
3333 WRITE(6,
'(A,F10.2)')
' Threshold for HADRIN events = (GeV)',
3417 IF (
nint(
what(1)) .GE. 10000 )
THEN
3430 ELSE IF (
nint(whtsav) .NE. 0 )
THEN
3436 IF ( abs(
nint(
what(2))) .GE. 10 )
THEN
3440 ELSE IF (
nint(
what(2)) .NE. 0 )
THEN
3461 IF(
what(1).EQ.0.)
THEN
3488 IF (
nint(
what(2)) .LT. 0 ) ldeexg = .false.
3548 IF(
what(2).GE.0.1d0)
THEN
3589 OPEN(29,
file=
'qel.evt',
3593 IF (ireso.EQ.1) CALL
distrp(1,ijproj,ppn)
3594 IF (cmhis.EQ.0.d0) CALL
distr(1,ijproj,ppn,idummy)
3595 IF (cmhis.EQ.1.d0) CALL
distrc(1,ijproj,ppn,idummy)
3596 IF (cmhis.EQ.2.d0) CALL
distco(1,ijproj,ppn,idummy)
3597 IF (ireso.EQ.1) CALL
disres(1,nhkkh1,ppn)
3598 IF (ipadis) CALL
distpa(1)
3599 IF (ioudif.EQ.1) CALL
diadif(1,0)
3601 WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
3604 * ,form=
'UNFORMATTED')
3608 WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
3616 WRITE(6,*)
' NEUTRINO: after INCINI call'
3620 IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN
3621 CALL
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
3623 IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN
3624 CALL
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
3633 IF(lpauli .AND. (.NOT.fermp))
THEN
3635 +
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
3636 +
' IF FERMI ACTIVE',
' LPAULI CHANGED TO .FALSE.'
3643 IF(nevnts.LE.0) nevnts=1000
3647 IF(ncases.LE.0) ncases=100
3652 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3715 WRITE(6,*)
' NEUTRINO initialization finished'
3729 OPEN(29,
file=
'diffnuc.evt',
3733 IF (ireso.EQ.1) CALL
distrp(1,ijproj,ppn)
3734 IF (cmhis.EQ.0.d0) CALL
distr(1,ijproj,ppn,idummy)
3735 IF (cmhis.EQ.1.d0) CALL
distrc(1,ijproj,ppn,idummy)
3736 IF (cmhis.EQ.2.d0) CALL
distco(1,ijproj,ppn,idummy)
3737 IF (ireso.EQ.1) CALL
disres(1,nhkkh1,ppn)
3738 IF (ipadis) CALL
distpa(1)
3739 IF (ioudif.EQ.1) CALL
diadif(1,0)
3741 WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
3744 * ,form=
'UNFORMATTED')
3748 WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
3759 IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN
3760 CALL
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
3762 IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN
3763 CALL
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
3772 IF(lpauli .AND. (.NOT.fermp))
THEN
3774 +
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
3775 +
' IF FERMI ACTIVE',
' LPAULI CHANGED TO .FALSE.'
3782 IF(nevnts.LE.0) nevnts=1000
3786 IF(ncases.LE.0) ncases=100
3791 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3873 WRITE(6,*)
'call xsglau'
3874 CALL
xsglau(ip,it,ijproj,1)
3943 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3953 COMMON /paname/ btype(30)
3955 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
3956 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
3957 +prebin,taebin,fermod,etacou
3959 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
3961 LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
3963 COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
3964 +ipadis,ishmal,lpauli
3966 COMMON /hadthr/ ehadth,inthad
3968 COMMON /nshmak/ nnshma,npshma,ntshma,nshmac,nshma2
3970 COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
3971 +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
3972 +irvs14, irvv11,irvv12,irvv13,irvv14
3974 COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
3975 * bsite(0:1,200),nstatb,nsiteb
3978 COMMON /damp/ ca,ci,ga
3982 COMMON /diffra/ isingd,idiftp,ioudif,iflagd
3983 COMMON /nomije/ ptmije(10),nnmije(10)
3984 DATA ptmije /5.d0,7.d0,9.d0,11.d0,13.d0,15.d0,17.d0
3985 +,19.d0,21.d0,23.d0 /
3988 DATA irco1,irco2,irco3,irco4,irco5 /5*0/
3989 DATA irss11,irss12,irss13,irss14,irsv11,irsv12,irsv13,irsv14 /8*0/
3990 DATA irvs11,irvs12,irvs13,irvs14,irvv11,irvv12,irvv13,irvv14 /8*0/
3994 DATA prepot /210*0.0/
3995 DATA taepot /210*0.0/
3996 DATA taebin,prebin,fermod /2*0.0d0,0.6d0/
3998 DATA btype /
'PROTON ' ,
'APROTON ' ,
'ELECTRON' ,
'POSITRON' ,
3999 +
'NEUTRIE ' ,
'ANEUTRIE' ,
'PHOTON ' ,
'NEUTRON ' ,
'ANEUTRON' ,
4000 +
'MUON+ ' ,
'MUON- ' ,
'KAONLONG' ,
'PION+ ' ,
'PION- ' ,
4001 +
'KAON+ ' ,
'KAON- ' ,
'LAMBDA ' ,
'ALAMBDA ' ,
'KAONSHRT' ,
4002 +
'SIGMA- ' ,
'SIGMA+ ' ,
'SIGMAZER' ,
'PIZERO ' ,
'KAONZERO' ,
4003 +
'AKAONZER' ,
'RESERVED' ,
'BLANK ' ,
'BLANK ' ,
'BLANK ' ,
4006 DATA ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr /0, 0, 0, -1, 0,
4009 DATA intpt, fermp, ihadss,ihadsv,ihadvs,ihadvv, ihada /.true.,
4010 +.true., 4*.false., .true./
4011 DATA ipadis, ishmal, lpauli /.false., .false., .true./
4016 DATA nstatb, nsiteb /2000, 200/
4020 DATA isingd,idiftp,ioudif,iflagd /0,0,0,0/
4041 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4044 CHARACTER*8 projty,targty
4047 COMMON /user1/
title,projty,targty
4048 COMMON /user2/cmener,sdfrac,ptlar,istruf,isingd,idubld
4052 COMMON /colle/nevhad,nvers,ihadrz,nfile
4056 common/collis/
s,ijproj,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
4063 common/booklt/btype(30),nbook(30)
4068 COMMON /polmn0/pdifr,phard,psoft,alfah,betah,
4069 * sigtot,sigqel,sigel,sigine,sighin,sigd,sigdd
4073 COMMON /pomtyp/ipim,icon,isig,lmax,mmax,
nmax,difel,difnu
4078 COMMON /dropjj/dropjt,dropva
4079 COMMON /gluspl/nugluu,nsgluu
4080 COMMON /ptlarg/xsmax
4081 COMMON /ptsamp/ isampt
4082 COMMON /stars/istar2,istar3
4083 COMMON /strufu/istrum,istrut
4084 COMMON /popcor/pdb,ajsdef
4089 CHARACTER*8 codewd,sdum
4099 9
FORMAT(
' special code word was used ')
4177 IF(codewd.EQ.
'SIGMAPOM')
THEN
4196 IF (itest.EQ.1)CALL
pomdi
4210 ELSEIF(codewd.EQ.
'GLUSPLIT')
THEN
4233 ELSEIF(codewd.EQ.
'PARTEV ')
THEN
4236 IF (
what(2).EQ.0.d0)npev=30
4238 IF (
what(3).EQ.0.d0)nvers=1
4243 IF(ipim.EQ.2)CALL
prblm2(cmener)
4265 ELSEIF(codewd.EQ.
'SELHARD ')
THEN
4270 IF(
what(5).NE.0.d0)
THEN
4272 IF(cmener.LT.2000.0d0.AND.isig.EQ.3)ptthr=
what(5)
4273 IF (cmener.GE.2000.0d0.AND.isig.EQ.3)
4274 * ptthr=0.25*
log(cmener/2000.)+2.
4275 IF(ptthr2.LT.ptthr)ptthr2=ptthr
4277 ptthr=2.1+0.15*(log10(cmener/50.))**3
4279 ELSEIF(istrut.EQ.2)
THEN
4280 ptthr=2.5+0.12*(log10(cmener/50.))**3
4284 1244
FORMAT (
' THRESHOLD PT FOR HARD SCATTERING PTTHR=',f12.2)
4299 ELSEIF(codewd.EQ.
'XSLAPT ')
THEN
4314 ELSEIF(codewd.EQ.
'SAMPT ')
THEN
4317 IF( isampt.LT.0 .OR. isampt.GT.4 ) isampt=0
4355 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4358 common/booklt/btype(30),nbook(30)
4360 DATA btype /
'PROTON ' ,
'APROTON ' ,
'ELECTRON' ,
4361 1
'POSITRON' ,
'NEUTRIE ' ,
'ANEUTRIE' ,
4362 2
'PHOTON ' ,
'NEUTRON ' ,
'ANEUTRON' ,
4363 3
'MUON+ ' ,
'MUON- ' ,
'KAONLONG' ,
4364 4
'PION+ ' ,
'PION- ' ,
'KAON+ ' ,
4365 5
'KAON- ' ,
'LAMBDA ' ,
'ALAMBDA ' ,
4366 6
'KAONSHRT' ,
'SIGMA- ' ,
'SIGMA+ ' ,
4367 7
'SIGMAZER' ,
'PIZERO ' ,
'KAONZERO' ,
4368 9
'AKAONZER' ,
' ' ,
' ' ,
4372 DATA nbook / 2212 , -2212 , 11 ,
4374 2 22 , 2112 , -2112 ,
4376 4 211 , -211 , 321 ,
4377 5 -321 , 3122 , -3122 ,
4378 6 310 , 3114 , 3224 ,
4379 7 3214 , 111 , 311 ,
4396 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4398 parameter(
zero=0.d0,
one=1.d0)
4399 parameter( alfa=0.56268
d-01,
beta=0.17173
d+03 )
4400 parameter( acc = 0.0001d0 )
4401 COMMON /xsecpt/ ptcut,sigs,dsigh
4402 COMMON /sigma / sigsof,bs,zsof,sighar,
fill(7)
4403 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
4405 common/collis/
s,ijproj,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
4407 CHARACTER*8 projty,targty
4410 COMMON /user1/
title,projty,targty
4411 COMMON /user2/cmener,sdfrac,ptlar,istruf,isingd,idubld
4413 common/ptsamp/ isampt
4414 dimension pptt(50),dpptt(50)
4417 IF ( mode.EQ.0 )
THEN
4419 pptt(ii)=ii*ptcut/50.
4423 IF(ecm.LT.1000.)
THEN
4424 aacucu=0.85*(ecm-400.)/600.
4425 sigs=(1.-aacucu)*sigsof
4431 WRITE(6,
'(A,4E12.4)')
' SAMPPT:ECM,PTCUT,SIGS,DSIGH',
4432 * ecm,ptcut,sigs,dsigh
4433 IF( isampt.EQ.0 )
THEN
4434 c = dsigh/(2.*sigs*ptcut)
4436 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4437 * ,c,sigsof,sighar,rmin
4438 ELSEIF( isampt.EQ.1 )
THEN
4442 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4443 * ,c,sigsof,sighar,rmin
4444 ELSEIF( isampt.EQ.2 )
THEN
4446 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4447 * ,c,sigsof,sighar,rmin
4448 ELSEIF( isampt.EQ.3 )
THEN
4450 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4451 * ,c,sigsof,sighar,rmin
4452 ELSEIF( isampt.EQ.4)
THEN
4453 aaaa=ptcut**2*(sigsof+sighar)
4454 IF (aaaa.LE.0.00001d0)
THEN
4455 aaaa=abs(aaaa)+0.0002
4456 WRITE(6,5559)ptcut,sigsof,sighar
4457 5559
FORMAT(
' SAMPPT:PTCUT,SIDSOF.SIGHATD:',3e12.3)
4460 b = 0.5*
bsofpt(acc,c,ptcut)
4461 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4462 * ,c,sigsof,sighar,rmin
4465 rmin =
exp(b*ptcut**2)
4469 IF( ioutpa.GE.-1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4470 * ,c,sigsof,sighar,rmin
4471 9010
FORMAT(
' SAMPPT MODE,ISAMPT,PTCUT,SIGS,DSIGH,B,C,SIGSOF',
4475 ELSEIF ( mode.EQ.1 )
THEN
4476 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4477 * ,c,sigsof,sighar,rmin
4478 ptt =
log(1.0-
rndm(v)*(1.0-rmin))/(b+0.00001d0)
4480 iipt=
pt*50./ptcut+1.
4482 dpptt(iipt)=dpptt(iipt)+1./(
pt+0.000001d0)
4485 ELSEIF(mode.EQ.2)
THEN
4487 dpptt(ii)=log10(1.
e-8+dpptt(ii))
4489 IF(iouxev.GE.-1)
THEN
4491 203
FORMAT(
' PT DISTRIBUTION OF SOFT PARTONS AS SAMPLED IN BSOFPT')
4500 *
FUNCTION bsofpt(ACC,CC,PPTCUT)
4501 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4504 COMMON /bsoff1/c,ptcut
4505 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
4507 dimension
x(50),
y(50)
4515 IF(c.LT.1.
d-10)
THEN
4534 IF (.NOT.succes)
THEN
4535 IF (kkkk.EQ.0)go to 400
4538 IF(iouxev.GE.1)
WRITE(6,111)b1,b2
4549 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4551 COMMON /bsoff1/c,ptcut
4554 df=c*ptcut**2*aaa-aaa
4564 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4566 COMMON /bsoff1/c,ptcut
4569 IF(qqq.GT.-60.)
THEN
4572 bsof1=c*(aaa-1.)-b*aaa
4581 *
FUNCTION rtsafe(FUNCD,X1,X2,XACC)
4582 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4584 parameter(maxit=200,itepri=0)
4585 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
4586 CALL funcd(
x1,fl,df)
4588 CALL funcd(
x2,fh,df)
4590 IF(fl*fh.GE.0.) pause
'ROOT MUST BE BRACKETED'
4616 * .OR. abs(2.*
f).GT.abs(dxold*df) )
THEN
4629 IF(abs(
dx).LT.xacc)
RETURN
4639 pause
'RTSAFE EXCEEDING MAXIMUM ITERATIONS'
4641 9995
FORMAT(
' VR1,VR2:',2e12.5)
4642 9996
FORMAT(
' RTSAFE,XH,XL,DXOLD,F,DF IN LOOP 11 J=1,MAXIT')
4643 9997
FORMAT(3
x,6e10.3)
4644 9998
FORMAT(
' RTSAFE: RTSAFE,F,DF =',3e12.5)
4645 9999
FORMAT(
' RTSAFE: F,DF =',2e12.5)
4651 *
FUNCTION var(A,B,C,D)
4652 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4654 parameter( ambmax = 1.0
d+38, epsi = 1.2
d-38,
one=1.d0 )
4658 abl = log10( abl + epsi )
4661 ccl = log10( ccl + epsi )
4663 IF( rcheck .LE. 38.d0 )
THEN
4666 var = ambmax*siab*sicc -
d
4668 IF(
var .GT. 1.0
d+18 )
var = 1.0
e+18
4669 IF(
var .LT. -1.0
d+18 )
var = -1.0
e+18
4674 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4677 parameter(factor=1.6d0,ntry=50)
4679 IF(
x1.EQ.
x2)pause
'You have to guess an initial range'
4684 IF(
f1*
f2.LT.0.d0)
RETURN
4685 IF(abs(
f1).LT.abs(
f2))
THEN
4718 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4722 parameter(amuamu=0.93149432d0)
4725 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
4813 COMMON /delp/ delpx,delpy,delpz,delpe
4815 COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
4817 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
4819 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4824 IF(help.GT.5.d0)phelp=help-5.
4839 IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13)
THEN
4844 pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
4845 pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
4848 IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14)
THEN
4855 IF(isthkk(i).EQ.1)
THEN
4862 IF(isthkk(i).EQ.13.AND.jdahkk(2,i).EQ.0)
THEN
4867 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
4868 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
4870 IF(isthkk(i).EQ.14.AND.jdahkk(2,i).EQ.0)
THEN
4877 IF(isthkk(i).EQ.16)
THEN
4883 eext=eext + phkk(4,i) - phkk(4,imo)
4885 IF(isthkk(i).EQ.15)
THEN
4891 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
4892 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
4893 eext=eext + phkk(4,i) - phkk(4,imo)
4895 IF(isthkk(i).EQ.1)
THEN
4898 IF(isthkk(i).EQ.-1)
THEN
4899 eeem1=eeem1+phkk(4,i)
4901 IF(isthkk(i).EQ.1001)
THEN
4902 ee1001=ee1001+phkk(4,i)
4905 eee=eee1+eeem1+ee1001
4916 aip=aip+(ait*amuamu+1.
d-3*
energy(ait,aitz))/epnto
4922 IF(it.EQ.ip)tole=0.02
4925 IF(delle.GE.tole)irej=1
4928 IF(icheck.LE.100)
THEN
4929 WRITE(6,
'(A,I5,E10.3,5F10.4)')
4930 *
' IP,EPN,AEEE,AEEEE,AEEE1,AEEEM1,AEE101:',
4931 * ip,epn,aeee,aeeee,aeee1,aeeem1,aee101
4932 WRITE(6,
'(A,I5,E10.3,7E12.4)')
4933 *
' IP,EPN,EEE,EEEE,EEE1,EEEM1,EE1001,DELLE,ELLE:',
4934 * ip,epn,eee,eeee,eee1,eeem1,ee1001,delle,elle
4950 1000
FORMAT(
' CHECKF: PX,PY,PZ,PE,EEXT,EEXP',2f7.3,2e12.3,2f7.3
4951 * / 8
x,
' DELPX/Y/Z/E',4f7.3,i10,
' IORIG')
4957 IF(isthkk(i).EQ.11)
THEN
4958 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
4959 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
4960 + (vhkk(khkk,i),khkk=1,4)
4963 IF(isthkk(i).EQ.12)
THEN
4964 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
4965 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
4966 + (vhkk(khkk,i),khkk=1,4)
4969 IF(isthkk(i).EQ.1)
THEN
4970 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
4971 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
4972 + (vhkk(khkk,i),khkk=1,4)
4974 1010
FORMAT (i6,i4,5i6,9(1pe10.2))
4976 IF(isthkk(i).EQ.16)
THEN
4978 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
4979 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
4980 + (vhkk(khkk,i),khkk=1,4)
4993 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4997 parameter(amuamu=0.93149432d0)
5000 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
5088 COMMON /delp/ delpx,delpy,delpz,delpe
5090 COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
5092 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
5094 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
5099 IF(help.GT.5.d0)phelp=help-5.d0
5100 pthelp=12.d0+phelp*5.d0
5120 IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13)
THEN
5125 pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
5126 pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
5129 IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14)
THEN
5136 IF(isthkk(i).EQ.1)
THEN
5143 IF(isthkk(i).EQ.13.AND.jdahkk(2,i).EQ.0)
THEN
5148 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
5149 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
5151 IF(isthkk(i).EQ.14.AND.jdahkk(2,i).EQ.0)
THEN
5158 IF(isthkk(i).EQ.16)
THEN
5164 eext=eext + phkk(4,i) - phkk(4,imo)
5166 IF(isthkk(i).EQ.15)
THEN
5172 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
5173 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
5174 eext=eext + phkk(4,i) - phkk(4,imo)
5176 IF(isthkk(i).EQ.1)
THEN
5181 IF(isthkk(i).EQ.-1)
THEN
5182 eeem1=eeem1+phkk(4,i)
5186 IF(isthkk(i).EQ.1001)
THEN
5187 ee1001=ee1001+phkk(4,i)
5188 pz1001=pz1001+phkk(3,i)
5189 px1001=px1001+phkk(1,i)
5192 eee=eee1+eeem1+ee1001
5193 pzpz=pz1+pzm1+pz1001
5194 pxpx=px1+pxm1+px1001
5205 IF(isthkk(i).EQ.1001)
THEN
5206 phkk(3,i)=phkk(3,i)+delpz
5207 phkk(4,i)=
sqrt(phkk(1,i)**2+phkk(2,i)**2+phkk(3,i)**2
5209 ee1001=ee1001+phkk(4,i)
5212 eee=eee1+eeem1+ee1001
5220 bip=epn+(ait*amuamu+1.
d-3*
energy(ait,aitz))
5226 IF(delle.GE.tole)irej=1
5229 IF(icheck.LE.20)
THEN
5230 WRITE(6,
'(A,I5,E10.3,4F10.4)')
5231 *
' IP,EPN,PXPX,PX1,PXM1,PX1001:',
5232 * ip,epn,pxpx,px1,pxm1,px1001
5233 WRITE(6,
'(A,I5,E10.3,6F10.4)')
5234 *
' IP,PPN,PZPZ,PZ1,PZM1,PZ1001,BIP,BMI:',
5235 * ip,ppn,pzpz,pz1,pzm1,pz1001,bip,bmi
5236 WRITE(6,
'(A,I5,E10.3,5E12.4)')
5237 *
' IP,EPN,EEE,EEE1,EEEM1,EE1001,DELLE:',
5238 * ip,epn,eee,eee1,eeem1,ee1001,delle
5254 1000
FORMAT(
' CHECKF: PX,PY,PZ,PE,EEXT,EEXP',2f7.3,2e12.3,2f7.3
5255 * / 8
x,
' DELPX/Y/Z/E',4f7.3,i10,
' IORIG')
5261 IF(isthkk(i).EQ.11)
THEN
5262 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5263 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5264 + (vhkk(khkk,i),khkk=1,4)
5267 IF(isthkk(i).EQ.12)
THEN
5268 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5269 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5270 + (vhkk(khkk,i),khkk=1,4)
5273 IF(isthkk(i).EQ.1)
THEN
5274 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5275 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5276 + (vhkk(khkk,i),khkk=1,4)
5278 1010
FORMAT (i6,i4,5i6,9(1pe10.2))
5280 IF(isthkk(i).EQ.16)
THEN
5282 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5283 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5284 + (vhkk(khkk,i),khkk=1,4)
5296 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5300 parameter(amuamu=0.93149432d0)
5303 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
5390 COMMON /zentra/ icentr
5392 COMMON /delp/ delpx,delpy,delpz,delpe
5394 COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
5396 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
5398 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
5403 IF(help.GT.5.d0)phelp=help-5.
5418 IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13)
THEN
5423 pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
5424 pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
5427 IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14)
THEN
5434 IF(isthkk(i).EQ.1)
THEN
5441 IF(isthkk(i).EQ.13.AND.jdahkk(2,i).EQ.0)
THEN
5446 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
5447 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
5449 IF(isthkk(i).EQ.14.AND.jdahkk(2,i).EQ.0)
THEN
5456 IF(isthkk(i).EQ.16)
THEN
5462 eext=eext + phkk(4,i) - phkk(4,imo)
5464 IF(isthkk(i).EQ.15)
THEN
5470 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
5471 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
5472 eext=eext + phkk(4,i) - phkk(4,imo)
5474 IF(isthkk(i).EQ.1)
THEN
5504 IF(ip.EQ.it.AND.it.GT.1)tole=0.05d0*ip
5506 IF(epn.LE.5.d0)tole=3.d0*tole
5508 IF (abs(
px).GT.pthelp.OR. abs(
py).GT.pthelp.OR.
5509 * abs(
pz)/epn.GT.tole.
5510 + or. abs(pe)/epn.GT.tole)
THEN
5513 IF(icheck.LE.50.AND.irej.EQ.1)
THEN
5514 WRITE(6,1000)
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe,iorig
5516 1000
FORMAT(
' CHECKO: PX,PY,PZ,PE,EEXT,EEXP',2f7.3,2e12.3,2f7.3
5517 * / 8
x,
' DELPX/Y/Z/E',4f7.3,i10,
' IORIG')
5521 WRITE(6,1000)
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe,iorig
5525 IF(isthkk(i).EQ.11)
THEN
5526 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5527 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5528 + (vhkk(khkk,i),khkk=1,4)
5531 IF(isthkk(i).EQ.12)
THEN
5532 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5533 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5534 + (vhkk(khkk,i),khkk=1,4)
5537 IF(isthkk(i).EQ.1)
THEN
5538 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5539 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5540 + (vhkk(khkk,i),khkk=1,4)
5542 1010
FORMAT (i6,i4,5i6,9(1pe10.2))
5544 IF(isthkk(i).EQ.16)
THEN
5546 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5547 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5548 + (vhkk(khkk,i),khkk=1,4)
5555 WRITE(6,1000)
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe,iorig
5561 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5567 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
5655 COMMON /delp/ delpx,delpy,delpz,delpe
5657 COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
5659 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
5669 IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13)
THEN
5674 pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
5675 pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
5678 IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14)
THEN
5685 IF(isthkk(i).EQ.1)
THEN
5692 IF(isthkk(i).EQ.13.AND.jdahkk(1,i).EQ.0)
THEN
5697 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
5698 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
5700 IF(isthkk(i).EQ.14.AND.jdahkk(1,i).EQ.0)
THEN
5707 IF(isthkk(i).EQ.16)
THEN
5713 eext=eext + phkk(4,i) - phkk(4,imo)
5715 IF(isthkk(i).EQ.15)
THEN
5721 eext=eext + phkk(4,i) - phkk(4,imo)
5728 WRITE(6,1000)
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe
5729 1000
FORMAT(
' CHECKE: PX,PY,PZ,PE,EEXT,EEXP',6f7.3/ 8
x,
' DELPX/Y/Z/E',4
5731 WRITE(6,
'(8X,A,6F8.3)')
' TASUMA,TASUBI,TABI,TAMASU,TAMA,TAIMMA',
5732 +tasuma,tasubi,tabi,tamasu,tama,taimma
5734 IF (abs(
px).GT.0.004.OR. abs(
py).GT.0.004.OR. abs(
pz).GT.0.004.
5735 + or. abs(pe).GT.0.004)
THEN
5739 IF(isthkk(i).EQ.11)
THEN
5740 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5741 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5742 + (vhkk(khkk,i),khkk=1,4)
5745 IF(isthkk(i).EQ.12)
THEN
5746 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5747 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5748 + (vhkk(khkk,i),khkk=1,4)
5751 IF(isthkk(i).EQ.1)
THEN
5752 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5753 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5754 + (vhkk(khkk,i),khkk=1,4)
5756 1010
FORMAT (i6,i4,5i6,9(1pe10.2))
5758 IF(isthkk(i).EQ.16)
THEN
5760 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5761 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5762 + (vhkk(khkk,i),khkk=1,4)
5790 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5797 DATA a1,a2,a3,a4,a5 /0.01575, 0.0178, 0.000710, 0.0237, 0.034/
5800 IF(ia.LE.1.OR.
iz.EQ.0)
THEN
5805 ebind = a1*aa - a2*aa**0. 666667- a3*
iz*
iz*aa**(-0.333333) - a4
5807 IF (
mod(ia,2).EQ.1)
THEN
5809 ELSEIF (
mod(
iz,2).EQ.1)
THEN
5821 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5828 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
5926 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
5927 +iibar(210),k1(210),k2(210)
5930 COMMON /factmo/ ifacto
5932 COMMON /taufo/ taufor,ktauge,itauve,incmod
5934 COMMON /hadthr/ ehadth,inthad
5938 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
5939 COMMON /nuccc/ jt,jtz,jp,jpz,jjproj,jbproj,jjtarg,jbtarg
5941 COMMON /zentra/ icentr
5943 COMMON /cmhico/ cmhis
5945 COMMON /resona/ ireso
5947 COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
5954 COMMON /projk/ iprojk
5956 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
5957 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
5958 +prebin,taebin,fermod,etacou
5960 COMMON /recom/irecom
5970 ppn=
sqrt((epn-aam(ijproj))*(epn+aam(ijproj)))
5971 ibproj=iibar(ijproj)
5974 jbproj=iibar(ijproj)
6026 SUBROUTINE hadhad(EPN,PPN,NHKKH1,IHTAWW,ITTA,IREJFO)
6027 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6037 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
6125 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
6137 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
6138 +iibar(210),k1(210),k2(210)
6141 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
6143 parameter(maxfin=10)
6144 COMMON /dfinls/ itrh(maxfin),cxrh(maxfin),cyrh(maxfin), czrh
6145 +(maxfin),elrh(maxfin),plrh(maxfin),irh
6147 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
6148 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
6149 +prebin,taebin,fermod,etacou
6156 IF(ipri.GE.2)
WRITE(6,1001) ijproj,ijproj,ppn,epn,cccxp,cccyp,
6157 +ccczp,ihtaww,itta,ieline
6158 1001
FORMAT(
' HADHAD 1:',
6159 +
' IJPROJ,IJPROJ,PPN,EPN,CCCXP,CCCYP,CCCZP,IHTAWW,ITTA,IELINE'/ 2
6160 +i3,2e12.3,3f7.3,3i4)
6165 CALL
sihnin(ijproj,itta,ppn,sight)
6166 CALL
sihnel(ijproj,itta,ppn,sighte)
6167 sigtot=sight + sighte
6168 IF (sigtot*
rndm(bb).LE.sighte)ieline=1
6169 IF(ipri.GE.2)
WRITE(6,1000) ijproj,ijproj,ppn,epn,cccxp,cccyp,
6170 +ccczp,ihtaww,itta,ieline
6171 1000
FORMAT(
' HADHAD 2 nach si...:',
6172 +
' IJPROJ,IJPROJ,PPN,EPN,CCCXP,CCCYP,CCCZP,IHTAWW,ITTA,IELINE'/ 2
6173 +i3,2e12.3,3f7.3,3i4)
6177 IF(ipri.GE.2)
WRITE(6,1012) ijproj,ijproj,ppn,epn,cccxp,cccyp,
6178 +ccczp,ihtaww,itta,ieline
6179 1012
FORMAT(
' HADHAD 12 loop:',
6180 +
' IJPROJ,IJPROJ,PPN,EPN,CCCXP,CCCYP,CCCZP,IHTAWW,ITTA,IELINE'/ 2
6181 +i3,2e12.3,3f7.3,3i4)
6184 WRITE(6,
'(I3,5(1PE12.4),I5/3X,5(1PE12.4))') ii,elrh(ii),plrh
6185 + (ii),cxrh(ii),cyrh(ii),czrh(ii),itrh(ii), (phkk(jjj,nhkk),jjj
6190 CALL
fhad(ijproj,ijproj,ppn,epn,cccxp,cccyp,ccczp, ihtaww,itta,
6195 +
WRITE(6,
'(A)')
' exit from hadhad with irejfo=1 '
6201 IF (ihadha.LT.3)
THEN
6204 IF(itsec.EQ.1.AND.elrh(ii).LE.taefep+aam(itsec)) goto 12
6205 IF(itsec.EQ.8.AND.elrh(ii).LE.taefen+aam(itsec)) goto 12
6206 IF(iibar(itsec).NE.1.AND.elrh(ii)-aam(itsec)
6207 + .LE.taepot(itsec)) goto 12
6212 IF (ipri.GE.2)
WRITE (6,1010)irh,nhkkh1,ihtaww,itta
6213 1010
FORMAT (
' HADHAD IRH,NHKKH1,IHTAWW,ITTA = ',4i5)
6217 +
' HADHAD - PARTICLE TRANSFER FROM /FINLSP/ INTO /HKKEVT/',
6218 +
' II, ELRH, PLRH, CXRH, CYRH, CZRH / PHKK(1-5)'
6227 WRITE (6,
'(A,2I5)') .EQ.
' HADHAD:NHKKNMXHKK ',nhkk,
nmxhkk
6231 idhkk(nhkk)=
mpdgha(itsec)
6233 jmohkk(2,nhkk)=ihtaww
6236 phkk(1,nhkk)=plrh(ii)*cxrh(ii)
6237 phkk(2,nhkk)=plrh(ii)*cyrh(ii)
6238 phkk(3,nhkk)=plrh(ii)*czrh(ii)
6239 phkk(4,nhkk)=elrh(ii)
6240 IF(phkk(4,nhkk)-aam(itsec).LE.taepot(itsec).
6241 + and.iibar(itsec).EQ.1)
THEN
6246 phkk(5,nhkk)=aam(itrh(ii))
6249 WRITE(6,
'(I3,5(1PE12.4),I5/3X,5(1PE12.4),I5)')
6251 + (ii),cxrh(ii),cyrh(ii),czrh(ii),itrh(ii), (phkk(jjj,nhkk),jjj
6254 vhkk(1,nhkk)=vhkk(1,ihtaww)
6255 vhkk(2,nhkk)=vhkk(2,ihtaww)
6256 vhkk(3,nhkk)=vhkk(3,ihtaww)
6257 vhkk(4,nhkk)=vhkk(4,1)
6260 jdahkk(1,1)=nhkkh1+1
6262 jdahkk(1,ihtaww)=nhkkh1+1
6263 jdahkk(2,ihtaww)=nhkk
6267 +
WRITE(6,
'(A)')
' exit from hadhad with irejfo=0 '
6271 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6277 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
6377 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
6378 +iibar(210),k1(210),k2(210)
6381 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
6383 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
6385 COMMON /chabai/chargi,barnui
6386 COMMON /evappp/ievap
6401 IF (isthkk(i).EQ.13)
THEN
6403 IF(idhkk(i).EQ.2212)chaeve=chaeve+1.d0
6405 IF (isthkk(i).EQ.14)
THEN
6407 IF(idhkk(i).EQ.2212)chaeve=chaeve+1.d0
6410 DO 521 i=nhkkh1,nhkk
6411 IF (isthkk(i).EQ.1.OR.isthkk(i).EQ.15.OR.isthkk(i).EQ.16)
THEN
6414 IF (nrhkk.LE.0.OR.nrhkk.GT.410)
THEN
6415 WRITE(6,1389)nrhkk,i,idhkk(i),nhkkh1,nhkk
6416 1389
FORMAT (
' distr: NRHKK ERROR ',5i10)
6421 chaeve=chaeve+ichhkk
6425 ELSEIF(ievap.EQ.1)
THEN
6427 IF (isthkk(i).EQ.1)
THEN
6431 chaeve=chaeve+ichhkk
6437 IF (isthkk(i).EQ.-1)
THEN
6438 IF(idhkk(i).EQ.2112)
THEN
6442 IF(idhkk(i).EQ.2212)
THEN
6448 IF((idhkk(i).EQ.80000).AND.(isthkk(i).NE.1000))
THEN
6449 chaeve=chaeve+idxres(i)
6450 baeve=baeve+idres(i)
6455 IF(ievl.LE.10)
WRITE(6,
'(2A,4F10.2)')
' Event charge and B-number',
6456 *
'=',chaeve,baeve,chargi,barnui
6457 IF(chaeve-chargi.NE.0.d0.OR.baeve-barnui.NE.0.d0)
THEN
6459 IF(ievl.LE.1000)
WRITE(6,
'(2A,4F10.2)')
'Event charge and B-numb',
6460 *
'(violated) =',chaeve,baeve,chargi,barnui
6471 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6473 dimension
pt(50,10),ypt(50,10)
6479 pt(j,i)=j*dpt-dpt/2.
6486 IF(ipt1.GT.50)ipt1=50
6487 IF(ipt2.GT.50)ipt2=50
6488 ypt(ipt1,ipt)=ypt(ipt1,ipt)+1.
6489 ypt(ipt2,ipt)=ypt(ipt2,ipt)+1.
6490 ypt(ipt1,10)=ypt(ipt1,10)+1.
6491 ypt(ipt2,10)=ypt(ipt2,10)+1.
6496 ypt(j,i)=ypt(j,i)/nevt
6497 ypt(j,i)=log10(ypt(j,i)+1.
d-18)
6509 SUBROUTINE hkkfil(IST,ID,M1,M2,PX,PY,PZ,E,NHKKAU,KORMO,ICALL)
6511 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6513 parameter(lout=6,llook=9)
6514 parameter(tiny10=1.0
d-10,tiny4=1.0
d-3)
6520 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
6607 COMMON /nncms/ gacms,bgcms,umo,pcm,eproj,pproj
6608 COMMON /trafop/galab,bglab,blab
6609 COMMON /projk/ iprojk
6621 WRITE(lout,1000) nhkk
6622 1000
FORMAT(1
x,
'HKKFIL: NHKK exeeds NMXHKK = ',i7,
6623 &
'! program execution stopped..')
6626 IF (m1.LT.0) mo1 = nhkk+m1
6627 IF (m2.LT.0) mo2 = nhkk+m2
6630 IF(kormo.EQ.999)
THEN
6631 jmohkk(1,nhkk) = mo1
6632 jmohkk(2,nhkk) = mo2
6634 jmohkk(1,nhkk)=nhkkau+kormo-1
6640 IF (jdahkk(1,mo1).NE.0)
THEN
6641 jdahkk(2,mo1) = nhkk
6643 jdahkk(1,mo1) = nhkk
6645 jdahkk(1,mo1)=nhkkau
6648 IF (jdahkk(1,mo2).NE.0)
THEN
6649 jdahkk(2,mo2) = nhkk
6651 jdahkk(1,mo2) = nhkk
6653 jdahkk(1,mo2) = nhkkau
6659 phkk(5,nhkk) = phkk(4,nhkk)**2-phkk(1,nhkk)**2-
6660 & phkk(2,nhkk)**2-phkk(3,nhkk)**2
6661 IF ((phkk(5,nhkk).LT.0.0d0).AND.(abs(phkk(5,nhkk)).GT.tiny4))
6662 &
WRITE(lout,
'(1X,A,G10.3)')
6663 &
'HKKFIL: negative mass**2 ',phkk(5,nhkk)
6664 phkk(5,nhkk) =
sqrt(abs(phkk(5,nhkk)))
6665 IF (ist.EQ.88888.OR.ist.EQ.88887.OR.ist.EQ.88889)
THEN
6671 vhkk(i,nhkk) = vhkk(i,mo2)
6673 vhkk(4,nhkk) = vhkk(3,mo2)/blab-vhkk(3,mo1)/bglab
6677 vhkk(i,nhkk) = vhkk(i,mo1)
6678 IF (iprojk.EQ.1)
THEN
6679 whkk(i,nhkk) = whkk(i,mo1)
6692 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6698 common/
pyjets/
n,npad,k(4000,5),
p(4000,5),v(4000,5)
6699 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
6700 common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
6701 common/
pydat3/mdcy(500,3),mdme(4000,2),brat(4000),kfdp(4000,5)
6703 CHARACTER chap*16,chan*16,chad(5)*16
6707 WRITE(mstu(11),6800)
6711 IF(mstu(2).NE.0) kfmax=mstu(2)
6716 IF(kc.EQ.0) goto 220
6717 IF(mstu(14).EQ.0.AND.kf.GT.100.AND.kc.LE.100) goto 220
6718 IF(mstu(14).GT.0.AND.kf.GT.100.AND.max(
mod(kf/1000,10),
6719 &
mod(kf/100,10)).GT.mstu(14)) goto 220
6727 IF(kf.LE.100.AND.chap.EQ.
' '.AND.mdcy(kc,2).EQ.0) goto 220
6732 idc2=mdcy(kc,2)+mdcy(kc,3)-1
6733 WRITE(mstu(11),6900)kbam,
6734 & kf,kc,idc1,idc2,chap,chan,kchg(kc,1),kchg(kc,2),
6735 & kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
6737 & kf,kc,idc1,idc2,chap,chan,kchg(kc,1),kchg(kc,2),
6738 & kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
6742 IF(kf.GT.100.AND.kc.LE.100) goto 220
6743 DO 210 idc=mdcy(kc,2),mdcy(kc,2)+mdcy(kc,3)-1
6746 CALL
pyname(kfdp(idc,j),chad(j))
6748 kbamdp(j)=
mcihad(kfdp(idc,j))
6749 IF(kbamdp(j).EQ.26)kbamdp(j)=0
6751 WRITE(26,7001) idc,mdme(idc,1),mdme(idc,2),brat(idc),
6753 210
WRITE(mstu(11),7000) idc,mdme(idc,1),mdme(idc,2),brat(idc),
6756 IF(kabam.NE.410)
THEN
6757 WRITE(mstu(11),6900)kabam,
6758 & -kf,-kc,idc1,idc2,chan,chap,-kchg(kc,1),kchg(kc,2),
6759 & kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
6760 WRITE(26,6900)kabam,
6761 & -kf,-kc,idc1,idc2,chan,chap,-kchg(kc,1),kchg(kc,2),
6762 & kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
6763 DO 211 idc=mdcy(kc,2),mdcy(kc,2)+mdcy(kc,3)-1
6767 IF(kcdp.LE.0.OR.kcdp.GT.500)
THEN
6773 IF(kchg(kcdp,3).EQ.0)kfdpm=kfdp(idc,j)
6775 IF(kbamdp(j).EQ.26)kbamdp(j)=0
6777 CALL
pyname(kfdpm,chad(j))
6779 WRITE(26,7001) idc,mdme(idc,1),mdme(idc,2),brat(idc),
6781 211
WRITE(mstu(11),7000) idc,mdme(idc,1),mdme(idc,2),brat(idc),
6789 6800
FORMAT(///30
x,
'Particle/parton data table'//1
x,
'BAM',
6790 &1
x,
'ABAM',1
x,
'KF',1
x,
'KC',1
x,
'DCF',1
x,
'DCL',1
x,
6791 &
'particle',8
x,
'antiparticle',6
x,
'chg col anti',8
x,
'mass',7
x,
6792 &
'width',7
x,
'w-cut',5
x,
'lifetime',1
x,
'decay'/11
x,
'IDC',1
x,
'on/off',
6793 &1
x,
'ME',3
x,
'Br.rat.',4
x,
'decay products')
6794 6900
FORMAT(/1
x,i4,i6,i4,2i5,a16,a16,3i3,1
x,f12.5,2(1
x,f11.5),
6796 7000
FORMAT(10
x,i4,2
x,i3,2
x,i3,2
x,f8.5,4
x,5a16)
6797 7001
FORMAT(10
x,i4,2
x,i3,2
x,i3,2
x,f8.5,4
x,5i5)
REAL *8 function bsof1(B)
subroutine checkn(EPN, PPN, IREJ, IORIG)
subroutine shmak1(ICASE, NN, NNA, NNB, NA, NB, UMO, BIMP)
G4int nint(G4double number)
REAL *8 function bsofpt(ACC, CC, PPTCUT)
subroutine pyname(KF, CHAU)
DOUBLE PRECISION function rndm(RDUMMY)
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
subroutine plomb(I, PP, CHAR, XF, ITIF, IJPROJ)
subroutine disres(IOP, IJPROJ, PPN)
BasicVector3D< T > unit() const
void fill(G4double x, G4double weight=1.)
subroutine shmak(ICASE, NN, NNA, NNB, NA, NB, UMO, BIMP)
subroutine hadhad(EPN, PPN, NHKKH1, IHTAWW, ITTA, IREJFO)
subroutine shmaki(NA, NCA, NB, NCB, RPROJ, RTARG, PPN)
subroutine sihnin(IPROJ, ITAR, PO, SIIN)
subroutine distrc(IOP, NHKKH1, PO, IGENER)
DOUBLE PRECISION function ebind(IA, IZ)
G4int mod(G4int a, G4int b)
subroutine parpt(IFL, PT1, PT2, IPT, NEVT)
subroutine plot(X, Y, N, M, MM, XO, DX, YO, DY)
REAL *8 function var(A, B, C, D)
REAL *8 function rtsafe(FUNCD, X1, X2, XACC)
subroutine checkf(EPN, PPN, IREJ, IORIG)
subroutine fhad(IPRMOD, IPRO, PLAB, ELAB, CX, CY, CZ, ITHKK, ITTA, IELINE, IREJFH)
subroutine distco(IOP, IJPROJ, PPN, IDUMMY)
subroutine defaul(EPN, PPN)
subroutine shmakf(NA, NCA, NB, NCB)
subroutine samppt(MODE, PT)
subroutine diadif(IOP, NHKKH1)
const char * what(void) const
subroutine kkinc(EPN, NTMASS, NTCHAR, NPMASS, NPCHAR, IDP, KKMAT, IDT, NHKKH1, IREJ)
subroutine plombc(I, PP, CHAR, XF, ITIF, IJPROJ)
subroutine xsglau(NA, NB, IJPROJ, NTARG)
subroutine title(NA, NB, NCA, NCB)
static c2_log_p< float_type > & log()
make a *new object
subroutine checko(EPN, PPN, IREJ, IORIG)
subroutine dminit(NCASES, EPN, PPN, NCOUNT, IGLAUB)
subroutine dttest(CODEWD, WHAT, SDUM)
subroutine chebch(IREJ, NHKKH1)
G4double tau(G4double bet, G4double homega, G4double ef, G4double t)
static c2_sqrt_p< float_type > & sqrt()
make a *new object
subroutine checke(EPN, PPN)
subroutine hkkfil(IST, ID, M1, M2, PX, PY, PZ, E, NHKKAU, KORMO, ICALL)
subroutine distr(IOP, NHKKH1, PO, IGENER)
void print(const std::vector< T > &data)
subroutine defaux(EPN, PPN)
subroutine zbrac(FUNC, X1, X2, SUCCES)
subroutine bsofc1(B, F, DF)
subroutine distrp(IOP, NHKKH1, PO)
subroutine sihnel(IPROJ, ITAR, POO, SIEL)
static c2_exp_p< float_type > & exp()
make a *new object