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 /inxdpm/intdpm
193 COMMON /nstari/nstart
194 COMMON /ncshxx/ncouxh,ncouxt
196 COMMON /nucros/dsigsu,dsigmc,ndsig
197 COMMON /nncms/ gamcm,bgcm,umo,pcm,eproj,pproj
198 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
199 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
200 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
201 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
202 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
203 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
204 COMMON /hdjase/nhse1,nhse2,nhse3,nhase1,nhase2,nhase3
205 COMMON /casadi/casaxx,icasad
206 COMMON /vxsvd/vxsp(50),vxst(50),vxsap(50),vxsat(50),
207 * vxvp(50),vxvt(50),vxdp(50),vxdt(50),
208 * nxsp,nxst,nxsap,nxsat,nxvp,nxvt,nxdp,nxdt
209 dimension vxsss(50,6),vxvvv(50,6),xxxx(50,6)
210 dimension xb(200),bimpp(200)
215 CHARACTER*8 projty,targty
216 COMMON /user1/titled,projty,targty
217 COMMON /user2/cmener,sdfrac,ptlar,istruf,isingx,idubld
218 COMMON /strufu/istrum,istrut
219 COMMON /ptsamp/ isampt
220 common/collis/
s,ijprox,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
221 COMMON /dropjj/dropjt,dropva
222 COMMON /pomtyp/ipim,icon,isig,lmax,mmax,
nmax,difel,difnu
224 COMMON /zentra/ icentr
228 COMMON /taufo/ taufor,ktauge,itauve,incmod
229 common/popcor/pdb,ajsdef
230 COMMON /diquax/amedd,idiqua,idiquu
231 COMMON /colle/nevhad,nvers,ihadrz,nfile
232 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
233 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
234 +prebin,taebin,fermod,etacou
235 COMMON /cronin/cronco,mkcron
236 COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
238 COMMON /secint/isecin
268 OPEN(47,
file=
'GLAUBTAR.DAT',
270 OPEN(37,
file=
'GLAUBCROSSPB.DAT',
308 CALL rluxgo(lux_level,iseed,0,0)
388 CALL
parpt(1,pt1,pt2,ipt,nevt)
391 CALL
dminit(ncases,multe,epnn,ppnn,ncount,iglaub)
411 IF( iglaub.EQ.1)
THEN
419 WRITE(6,*)
' Printout of important Parameters before DPMJET run.'
420 *,
' Please note for DPMJET input all numbers are floating point!'
421 WRITE(6,*)
'PROJPAR ',ip,ipz
422 WRITE(6,*)
'TARPAR ',it,itz
423 WRITE(6,*)
'MOMENTUM ',ppn
424 WRITE(6,*)
'ENERGY ',epn
425 WRITE(6,*)
'CMENERGY ',umo
426 WRITE(6,*)
'NOFINALE ',ifinal
427 WRITE(6,*)
'EVAPORAT ',ievap
428 WRITE(6,*)
'OUTLEVEL ',ipri,ipev,ippa,ipco,
init,iphkk
429 auauau=rd2out(iseed1,iseed2)
430 WRITE(6,*)
'RANDOMIZ ',iseed1,iseed2,
' Initial RNDM (RM48) seeds'
431 WRITE(6,*)
'STRUCFUN ',istruf+100*istrut
432 WRITE(6,*)
'SAMPT ',isampt
433 WRITE(6,*)
'SELHARD ',0,iophrd, 0,dropjt,ptthr,ptthr2
434 WRITE(6,*)
'SIGMAPOM ',0,isig,ipim+10*icon,imax,mmax,
nmax
435 WRITE(6,*)
'PSHOWER ',ipshow
436 WRITE(6,*)
'CENTRAL ',icentr
437 WRITE(6,*)
'CMHISTO ',cmhis
438 WRITE(6,*)
'SEASU3 ',seasq
439 WRITE(6,*)
'RECOMBIN ',irecom
440 WRITE(6,*)
'SINGDIFF ',isingd
441 WRITE(6,*)
'TAUFOR ',taufor,ktauge,itauve
442 WRITE(6,*)
'POPCORN ',pdb
443 WRITE(6,*)
'POPCORCK ',ijpock,pdbck
444 WRITE(6,*)
'POPCORSE ',pdbse,pdbseu
445 WRITE(6,*)
'CASADIQU ',icasad,casaxx
446 WRITE(6,*)
'DIQUARKS ',idiqua,idiquu,amedd
447 WRITE(6,*)
'HADRONIZ ',ihadrz
448 WRITE(6,*)
'INTPT ',intpt
449 WRITE(6,*)
'PAULI ',lpauli
450 WRITE(6,*)
'FERMI ',fermp,fermod
451 WRITE(6,*)
'CRONINPT ',mkcron,cronco
452 WRITE(6,*)
'SEADISTR ',xseacu+0.95d0,unon,unom,unosea
453 WRITE(6,*)
'SEAQUARK ',seaqx,seaqxn
454 WRITE(6,*)
'SECINTER ',isecin
455 WRITE(6,*)
'XCUTS ',cvq,cdq,csea,ssmima
456 WRITE(6,*)
'START ',ncases
457 WRITE(6,*)
' Printout of important Parameters before DPMJET run.'
458 *,
' Please note for DPMJET input all numbers are floating point!'
464 ndone=(iiii-1)*ncaset
466 WRITE(6,
'(A,4I5)')
' KKINC call ',iit,iitz,iip,iipz
468 1111
FORMAT(
' NDONE= ',i10)
475 ELSEIF(multe.EQ.1)
THEN
476 epn=0.1d0*epnn+
rndm(v)*(1.9d0*epnn)
478 IF(ijproj.NE.0) nnpp=ijproj
479 ppn=
sqrt((epn-aam(nnpp))*(epn+aam(nnpp)))
486 pproj =
sqrt((epn-amproj)*(epn+amproj))
487 umo =
sqrt(amproj**2 + amtar**2 + 2.*amtar*eproj)
490 ptthr=2.1+0.15*(log10(cmener/50.))**3
492 ELSEIF(istrut.EQ.2)
THEN
493 ptthr=2.5+0.12*(log10(cmener/50.))**3
496 gamcm = (eproj+amtar)/umo
499 pcm=gamcm*pproj - bgcm*eproj
503 1033
FORMAT(
' CMS: ' ,
504 +
' EPROJ,PPROJ,AMPROJ,AMTAR,UMO,GAMCM,BGCM,PCM'/8e22.13)
507 numev = i+(iiii-1)*ncaset
508 IF ((i.EQ.486).OR.(i.EQ.803).OR.(i.EQ.1368).OR.
509 & (i.EQ.1465).OR.(i.EQ.1693).OR.(i.EQ.1808))
THEN
580 CALL
kkinc(epn,iit,iitz,iip,iipz,iiproj,kkmat,
581 * iitarg,nhkkh1,irej)
582 ELSEIF(intdpm.EQ.1)
THEN
600 CALL
dpmevt(elabt,iiipro,iiip,iiipz,iiit,iiitz,kkmat,nhkkh1)
603 IF(irej.EQ.1)go to 765
692 WRITE(6,*)
' Printout of important Parameters after DPMJET run.'
693 *,
' Please note for DPMJET input all numbers are floating point!'
694 WRITE(6,*)
'PROJPAR ',ip,ipz
695 WRITE(6,*)
'TARPAR ',it,itz
696 WRITE(6,*)
'MOMENTUM ',ppn
697 WRITE(6,*)
'ENERGY ',epn
698 WRITE(6,*)
'CMENERGY ',umo
699 WRITE(6,*)
'NOFINALE ',ifinal
700 WRITE(6,*)
'EVAPORAT ',ievap
701 WRITE(6,*)
'OUTLEVEL ',ipri,ipev,ippa,ipco,
init,iphkk
702 auauau=rd2out(iseed1,iseed2)
703 WRITE(6,*)
'RANDOMIZ ',iseed1,iseed2,
' Final RNDM (RM48) seeds'
704 WRITE(6,*)
'STRUCFUN ',istruf+100*istrut
705 WRITE(6,*)
'SAMPT ',isampt
706 WRITE(6,*)
'SELHARD ',0,iophrd, 0,dropjt,ptthr,ptthr2
707 WRITE(6,*)
'SIGMAPOM ',0,isig,ipim+10*icon,imax,mmax,
nmax
708 WRITE(6,*)
'PSHOWER ',ipshow
709 WRITE(6,*)
'CENTRAL ',icentr
710 WRITE(6,*)
'CMHISTO ',cmhis
711 WRITE(6,*)
'SEASU3 ',seasq
712 WRITE(6,*)
'RECOMBIN ',irecom
713 WRITE(6,*)
'SINGDIFF ',isingd
714 WRITE(6,*)
'TAUFOR ',taufor,ktauge,itauve
715 WRITE(6,*)
'POPCORN ',pdb
716 WRITE(6,*)
'POPCORCK ',ijpock,pdbck
717 WRITE(6,*)
'POPCORSE ',pdbse,pdbseu
718 WRITE(6,*)
'CASADIQU ',icasad,casaxx
719 WRITE(6,*)
'DIQUARKS ',idiqua,idiquu,amedd
720 WRITE(6,*)
'HADRONIZ ',ihadrz
721 WRITE(6,*)
'INTPT ',intpt
722 WRITE(6,*)
'PAULI ',lpauli
723 WRITE(6,*)
'FERMI ',fermp,fermod
724 WRITE(6,*)
'CRONINPT ',mkcron,cronco
725 WRITE(6,*)
'SEADISTR ',xseacu+0.95d0,unon,unom,unosea
726 WRITE(6,*)
'SEAQUARK ',seaqx,seaqxn
727 WRITE(6,*)
'SECINTER ',isecin
728 WRITE(6,*)
'XCUTS ',cvq,cdq,csea,ssmima
729 WRITE(6,*)
' Printout of important Parameters after DPMJET run.'
730 *,
' Please note for DPMJET input all numbers are floating point!'
736 auauau=rd2out(iseed1,iseed2)
737 WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
738 WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
739 WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
740 WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
752 WRITE(6,1100) irvv11,irvv12,irvv13,irvv14, irsv11,irsv12,irsv13,
753 + irsv14, irvs11,irvs12,irvs13,irvs14, irss11,irss12,irss13,irss14
754 1100
FORMAT (
' REJECTION COUNTERS FROM KKEVT',/, 5
x,
' V-V CHAINS',4i6/
755 +5
x,
' S-V CHAINS',4i6/ 5
x,
' V-S CHAINS',4i6/ 5
x,
' S-S CHAINS',4i6)
756 WRITE(6,
'(A,4I10)')
' POPCCK/SE/S3/S0 rejections ',
757 * irejck,irejse,irejs3,irejs0
758 WRITE(6,
'(A,4I10)')
' POPCCK/ASE/AS3/AS0 rejections ',
759 * irejsa,ireja3,ireja0
760 WRITE(6,
'(2A,8I6)')
' POPCCK ICK4,ICK6,IHAD4,IHAD6,ISE4,ISE6 ',
761 *
'ISE43,ISE63 ', ick4,ick6,ihad4,ihad6,ise4,ise6,ise43,ise63
762 WRITE(6,
'(2A,8I6)')
' POPCSAQ IHADA4,IHADA6,ISEA4,ISEA6 ',
763 *
'ISEA43,ISEA63 ', ihada4,ihada6,isea4,isea6,isea43,isea63
764 WRITE(6,*)
' NDVUU,NDVUS,NDVSS,NVDUU,NVDUS,NVDSS',
765 *
' NDSUU,NDSUS,NDSSS,NSDUU,NSDUS,NSDSS',
766 *
' NDZUU,NDZUS,NDZSS,NZDUU,NZDUS,NZDSS' ,
767 * ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
768 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
769 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
770 WRITE(6,*)
' NADVUU,NADVUS,NADVSS,NAVDUU,NAVDUS,NAVDSS',
771 *
' NADSUU,NADSUS,NADSSS,NASDUU,NASDUS,NASDSS',
772 *
' NADZUU,NADZUS,NADZSS,NAZDUU,NAZDUS,NAZDSS' ,
773 * nadvuu,nadvus,nadvss,navduu,navdus,navdss,
774 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
775 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
776 WRITE(6,*)
' NHSE1,NHSE2,NHSE3,NHASE1,NHASE2,NHASE3 ',
777 * nhse1,nhse2,nhse3,nhase1,nhase2,nhase3
784 WRITE(6,
'(A/7I8)')
' Diquark rejection IDIQRE(1-7),N,ss,su,ud',
785 & (idiqre(jj),jj=1,7)
786 WRITE(6,
'(A/7I8)')
' Diquark rejection IDIQRZ(1-7),N,ss,su,ud',
787 & (idiqrz(jj),jj=1,7)
788 WRITE(6,*)
' Diquark rej. IDVRE(1-3),ud,us,ss ',(idvre(jj),jj=1,3)
789 WRITE(6,*)
' Diquark rej. IVDRE(1-3),ud,us,ss ',(ivdre(jj),jj=1,3)
790 WRITE(6,*)
' Diquark rej. IDSRE(1-3),ud,us,ss ',(idsre(jj),jj=1,3)
791 WRITE(6,*)
' Diquark rej. ISDRE(1-3),ud,us,ss ',(isdre(jj),jj=1,3)
792 WRITE(6,*)
' Diquark rej. IDZRE(1-3),ud,us,ss ',(idzre(jj),jj=1,3)
793 WRITE(6,*)
' Diquark rej. IZDRE(1-3),ud,us,ss ',(izdre(jj),jj=1,3)
794 WRITE(6,*)
' NDZSU,NZDSU ',ndzsu,nzdsu
796 IF ((cmhis.EQ.1.d0).AND.(ioudif.EQ.1))
800 WRITE(6,*)
' Output of x-distribution survey',
801 *
'VXSP(II),VXST(II),VXSAP(II),VXSAT(II),',
802 *
'VXVP(II),VXVT(II),VXDP(II),VXDT(II)' ,
803 *nxsp,nxst,nxsap,nxsat,nxvp,nxvt,nxdp,nxdt
805 IF(nxsp.GE.1)vxsp(ii)=50.d0*vxsp(ii)/nxsp
806 IF(nxst.GE.1)vxst(ii)=50.d0*vxst(ii)/nxst
807 IF(nxsap.GE.1)vxsap(ii)=50.d0*vxsap(ii)/nxsap
808 IF(nxsat.GE.1)vxsat(ii)=50.d0*vxsat(ii)/nxsat
809 IF(nxvp.GE.1)vxvp(ii)=50.d0*vxvp(ii)/nxvp
810 IF(nxvt.GE.1)vxvt(ii)=50.d0*vxvt(ii)/nxvt
811 IF(nxdp.GE.1)vxdp(ii)=50.d0*vxdp(ii)/nxdp
812 IF(nxdt.GE.1)vxdt(ii)=50.d0*vxdt(ii)/nxdt
813 xxxxx=ii*0.02d0-0.01d0
820 fxvvv=(1.-xxxxx)**3/
sqrt(xxxxx)
821 fxddd=2.d0*xxxxx**3.0d0/
sqrt(1.d0-xxxxx)
822 vxsss(ii,1)=log10(vxsp(ii))
823 vxsss(ii,2)=log10(vxst(ii))
824 vxsss(ii,3)=log10(vxsap(ii))
825 vxsss(ii,4)=log10(vxsat(ii))
826 vxvvv(ii,1)=log10(vxvp(ii))
827 vxvvv(ii,2)=log10(vxvt(ii))
828 vxvvv(ii,3)=log10(vxdp(ii))
829 vxvvv(ii,4)=log10(vxdt(ii))
830 vxvvv(ii,5)=log10(fxvvv)
831 vxvvv(ii,6)=log10(fxddd)
832 axsp=axsp+0.02d0*vxsp(ii)*xxxxx
833 axst=axst+0.02d0*vxst(ii)*xxxxx
834 axsap=axsap+0.02d0*vxsap(ii)*xxxxx
835 axsat=axsat+0.02d0*vxsat(ii)*xxxxx
836 axvp=axvp+0.02d0*vxvp(ii)*xxxxx
837 axvt=axvt+0.02d0*vxvt(ii)*xxxxx
838 axdp=axdp+0.02d0*vxdp(ii)*xxxxx
839 axdt=axdt+0.02d0*vxdt(ii)*xxxxx
840 WRITE(6,*)vxsp(ii),vxst(ii),vxsap(ii),vxsat(ii),
841 * vxvp(ii),vxvt(ii),vxdp(ii),vxdt(ii)
844 *axsp,axst,axsap,axsat,axvp,axvt,axdp,axdt
845 CALL
plot(xxxx,vxsss,200,4,50,0.d0,0.02d0,-3.d0,0.05d0)
846 CALL
plot(xxxx,vxvvv,300,6,50,0.d0,0.02d0,-3.d0,0.05d0)
850 CALL
parpt(3,pt1,pt2,ipt,ncases)
854 fracxs=float(ncouxh)/(float(ncouxh)+float(ncouxt))
856 WRITE(6,*)
' Fraction of x-sect: ',fracxs,ncouxh,ncouxt
861 IF(ndsig.GE.1) dsigmc=dsigsu/ndsig
862 WRITE(6,*)
' Neutrino-nucleon cross section DSIGMC,NDSIG ',
863 & dsigmc,
' *10**(-38) cm**2 ',ndsig,
' evts'
896 IF(ishmal) CALL
shmak(3,nshmac,np,
nt,ip,it,umo,bimp)
897 IF(ishmal) CALL
shmak1(3,nshma2,np,
nt,ip,it,umo,bimp)
899 IF (ireso.EQ.1) CALL
distrp(3,ncases,ppn)
900 IF (cmhis.EQ.0.d0) CALL
distr(3,ncases,ppn,idummy)
901 IF (cmhis.EQ.1.d0) CALL
distrc(3,ncases,ppn,idummy)
902 IF (cmhis.EQ.2.d0) CALL
distco(3,ncases,ppn,idummy)
903 IF (ireso.EQ.1) CALL
disres(3,ncases,ppn)
906 IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN
907 CALL
plomb(5,pp,char,xfxfxf,itif,ijproj)
909 IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN
910 CALL
plombc(5,pp,char,xfxfxf,itif,ijproj)
921 SUBROUTINE dminit(NCASES,MULTE,EPN,PPN,NCOUNT,IGLAUB)
922 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1034 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
1132 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
1133 +iibar(210),k1(210),k2(210)
1142 COMMON /paname/ btype(30)
1144 COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
1145 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
1147 COMMON /factmo/ ifacto
1149 COMMON /taufo/ taufor,ktauge,itauve,incmod
1151 COMMON /rptshm/ rproj,rtarg,bimpac
1153 COMMON /trafop/ gamp,bgamp,betp
1155 COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
1156 +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
1157 +irvs14, irvv11,irvv12,irvv13,irvv14
1159 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
1161 COMMON /dnun/ nn,np,
nt
1163 COMMON /nshmak/ nnshma,npshma,ntshma,nshmac,nshma2
1165 COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
1166 * bsite(0:1,200),nstatb,nsiteb
1168 COMMON /hadthr/ ehadth,inthad
1172 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
1173 COMMON /nuccc/ jt,jtz,jp,jpz,jjproj,jbproj,jjtarg,jbtarg
1175 COMMON /zentra/ icentr
1177 COMMON /cmhico/ cmhis
1179 COMMON /resona/ ireso
1181 LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
1183 COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
1184 +ipadis,ishmal,lpauli
1186 COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
1189 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
1190 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
1191 +prebin,taebin,fermod,etacou
1197 COMMON /projk/ iprojk
1199 parameter(lunber=14)
1202 COMMON /seaqxx/ seaqx,seaqxn
1203 COMMON /cronin/cronco,mkcron
1205 COMMON /seadiq/lseadi
1206 COMMON /final/ifinal
1207 COMMON /recom/irecom
1209 COMMON /neutyy/ neutyp,neudec
1210 COMMON /nstari/nstart
1211 common/popcor/pdb,ajsdef
1212 common/popcck/pdbck,pdbse,pdbseu,
1213 * ijpock,irejck,ick4,ihad4,ick6,ihad6
1214 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
1215 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
1216 *isea43,isea63,irejao
1218 COMMON /bufues/ bnnvv,bnnss,bnnsv,bnnvs,bnncc,
1219 * bnndv,bnnvd,bnnds,bnnsd,
1221 * bptvv,bptss,bptsv,bptvs,bptcc,bptdv,
1222 * bptvd,bptds,bptsd,
1224 * beevv,beess,beesv,beevs,beecc,beedv,
1225 * beevd,beeds,beesd,
1227 * ,bnndi,bptdi,beedi
1228 * ,bnnzd,bnndz,bptzd,bptdz,beezd,beedz
1229 COMMON /ncoucs/ bcouvv,bcouss,bcousv,bcouvs,
1230 * bcouzz,bcouhh,bcouds,bcousd,
1231 * bcoudz,bcouzd,bcoudi,
1232 * bcoudv,bcouvd,bcoucc
1233 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
1234 * anndv,annvd,annds,annsd,
1236 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
1238 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
1241 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
1242 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
1243 * acouzz,acouhh,acouds,acousd,
1244 * acoudz,acouzd,acoudi,
1245 * acoudv,acouvd,acoucc
1246 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
1247 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
1248 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
1249 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
1250 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
1251 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
1252 COMMON /hdjase/nhse1,nhse2,nhse3,nhase1,nhase2,nhase3
1254 COMMON /diffra/ isingd,idiftp,ioudif,iflagd
1255 COMMON /seasu3/seasq
1256 COMMON /ifragm/ifrag
1257 COMMON /fluctu/ifluct
1258 COMMON /diquax/amedd,idiqua,idiquu
1259 COMMON /inxdpm/intdpm
1260 COMMON /vxsvd/vxsp(50),vxst(50),vxsap(50),vxsat(50),
1261 * vxvp(50),vxvt(50),vxdp(50),vxdt(50),
1262 * nxsp,nxst,nxsap,nxsat,nxvp,nxvt,nxdp,nxdt
1265 COMMON /nncms/ gamcm,bgcm,umo,pcm,eproj,pproj
1268 COMMON /xsecpt/ ptcut,sigs,dsigh
1269 COMMON /kglaub/jglaub
1270 COMMON /xsecnu/ecmuu,ecmoo,ngritt,nevtt
1314 CHARACTER*8 projty,targty
1317 COMMON /user1/titled,projty,targty
1318 COMMON /user2/cmener,sdfrac,ptlar,istruf,isingx,idubld
1342 COMMON /colle/nevhad,nvers,ihadrz,nfile
1354 common/collis/
s,ijprox,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
1363 common/booklt/btypex(30),nbook(30)
1381 COMMON /polmn0/pdifr,phard,psoft,alfah,betah,
1382 * sigtot,sigqel,sigel,sigine,sighin,sigd,sigdd
1389 COMMON /pomtyp/ipim,icon,isig,lmax,mmax,
nmax,difel,difnu
1396 COMMON /dropjj/dropjt,dropva
1397 COMMON /gluspl/nugluu,nsgluu
1400 COMMON /ptlarg/xsmax
1402 COMMON /ptsamp/ isampt
1403 COMMON /stars/istar2,istar3
1406 COMMON /strufu/istrum,istrut
1407 COMMON /cutofn/ncutox
1410 COMMON /harlun/ qlun,iharlu
1411 COMMON /pomtab/ipomta
1412 COMMON /sincha/isichaa
1414 COMMON /evappp/ievap
1416 parameter( frdiff = 0.2
d+00 )
1417 parameter( ethsea = 1.0
d+00 )
1419 LOGICAL ldiffr, linctv, levprt, lheavy, ldeexg, lgdhpr, lpreex,
1420 & lhlfix, lprfix, lparwv, lpower, lsngch, llvmod, lschdf
1421 COMMON / parevt / dpower, fsprd0, fshpfn, rn1gsc, rn2gsc,
1422 & ldiffr(39),lpower, linctv, levprt, lheavy,
1423 & ldeexg, lgdhpr, lpreex, lhlfix, lprfix, lparwv,
1424 & ilvmod, jlvmod, llvmod, lsngch, lschdf
1434 parameter(
mxpsst = 300 )
1435 parameter(
mxpsfb = 41000 )
1436 LOGICAL lfrmbk, lncmss
1437 COMMON / frbkcm / amufbk, eexfbk(
mxpsst), amfrbk(
mxpsst),
1439 & exmxfb, r0frbk, r0cfbk, c1cfbk, c2cfbk,
1444 & ifbcha(5,
mxpsfb), iposst, iposfb, ifbstf,
1445 & ifbfrb, nbufbk, lfrmbk, lncmss
1447 COMMON /inpflg/ iang,ifiss,ib0,igeom,istrag,keydk
1448 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
1449 COMMON /secint/isecin
1450 COMMON /nuclea/ pfermp(2),pfermn(2),fermdd,
1451 & ebindp(2),ebindn(2),epot(2,210),
1453 COMMON /ferfor/iferfo
1454 COMMON /casadi/casaxx,icasad
1455 COMMON /infore/ifrej
1465 CHARACTER*8 code,codewd,blank,sdum
1468 1
'TITLE ',
'PROJPAR ',
'TARPAR ',
'ENERGY ',
'HADRONIZ',
1469 2
'RANDOMIZ',
'FERMI ',
'EVENTAPE',
'START ',
'PARTEV ',
1470 3
'INTPT ',
'TECALBAM',
'RESONANC',
'VALVAL ',
'COMMENT ',
1471 4
'OUTLEVEL',
'LEPTOEVT',
'SEASEA ',
'PARTICLE',
'ALLPART ',
1472 5
'TAUFOR ',
'SEAVAL ',
'VALSEA ',
'MOMENTUM',
'PAULI ',
1473 6
'PROJKASK',
'CENTRAL ',
'SEADISTR',
'CMHISTO ',
'SIGTEST ',
1474 7
'XCUTS ',
'HADRIN ',
'FACTOMOM',
'COULOMB ',
'GLAUBERI',
1475 8
'EDENSITY',
'CMENERGY',
'INFOREJE',
'RECOMBIN',
'SINGDIFF',
1476 9
'NOFINALE',
'SEASU3 ',
'CRONINPT',
'POPCORN ',
'STOP ',
1477 9
'FLUCTUAT',
'DIQUARKS',
'HBOOKHIS',
'GLAUBERA',
'POMTABLE',
1478 9
'SINGLECH',
'HADRINTH',
'EVAPORAT',
'SEAQUARK',
'SECINTER',
1479 9
'POPCORCK',
'CASADIQU',
'POPCORSE',
'NEUTRINO',
'DIFFNUC ',
1480 9
'XSECNUC ',
'INTERDPM',
' ',
' ',
' '/
1510 IF (ncount.EQ.1)
THEN
1518 1000
FORMAT(
'1 **************************************************',
1519 +
'**************************************************', //
1520 +
' DPMJET VERSION II.5 (Sept. 1999) ' /
1521 +
' DUAL PARTON MODEL FOR HADRON NUCLEUS COLLISIONS '/ /
1522 +
' AND NUCLEUS NUCLEUS COLLISIONS '/
1523 +
' INCLUDING A FORMATION TIME INTRANUCLEAR CASCADE'/
1524 4
' Minijets and DTUJET like multiple soft jets '/
1525 4
' Nuclear evaporation and residual target and '/
1526 4
' projectile nuclei '/
1527 +
' **************************************************',
1528 +
'**************************************************',//)
1633 IF (ihadrz.GE.2)
THEN
1675 ptthr=2.1+0.15*(log10(cmener/50.))**3
1677 ELSEIF(istrut.EQ.2)
THEN
1678 ptthr=2.5+0.12*(log10(cmener/50.))**3
1696 IF(ijproj.NE.0) nnpp=ijproj
1697 epn=
sqrt(ppn**2+aam(nnpp)**2)
1703 pproj =
sqrt((epn-amproj)*(epn+amproj))
1704 umo =
sqrt(amproj**2 + amtar**2 + 2.*amtar*eproj)
1706 gamcm = (eproj+amtar)/umo
1709 pcm=gamcm*pproj - bgcm*eproj
1711 ptthr=2.1+0.15*(log10(cmener/50.))**3
1713 ELSEIF(istrut.EQ.2)
THEN
1714 ptthr=2.5+0.12*(log10(cmener/50.))**3
1738 istruf=istruf-istrut*100
1835 IF(ieof.EQ.1) go to 40
1837 READ(5,1010)codewd,(
what(i),i=1,6),sdum
1838 WRITE(6,1020)codewd,(
what(i),i=1,6),sdum
1841 IF(codewd.EQ.
code(isw)) go to 30
1849 + 50 , 60 , 90 , 120 , 130 ,
1853 + 140 , 150 , 160 , 170 , 210 ,
1857 + 220 , 230 , 240 , 250 , 260 ,
1861 + 280 , 290 , 300 , 310 , 320 ,
1865 + 330 , 340 , 350 , 360 , 370 ,
1869 + 380 , 390 , 400 , 410 , 420 ,
1873 + 430 , 440 , 450 , 460 , 470 ,
1877 + 480 , 490 , 500 , 510 , 520 ,
1881 + 530 , 535 , 538, 539, 540 ,
1885 + 541 , 542 , 543, 544, 545,
1890 + 551 , 552 , 553 , 554 ,555 ,
1894 + 556 , 557 ,558 , 559 ,560,
1898 + 620 , 630 ,640 , 650 ,660,610),isw
1920 1010
FORMAT(a8,2
x,6e10.0,a8 )
1921 1020
FORMAT(
' *****NEXT CONTROL CARD ***** ',a10,6(1
x,g11.4), 2
x,a10)
1923 1030
FORMAT(/,
' UNKNOWN CODEWORD - CONTROL CARD IGNORED')
1924 1040
FORMAT(/,
' UNEXPECTED END OF INPUT - STOP ASSUMED.')
1925 1050
FORMAT(/,
' UNEXPECTED END OF INPUT - START ASSUMED.')
1965 IF (codewd.GT.
'-zzzzzzz')
1966 1
WRITE(6,91) codewd,(
what(i),i=1,6),sdum
1967 91
FORMAT(
' ---- control input card : ----'
1968 1 /1
x,a8,2
x,6(f10.3),a8)
1985 IF(codewd.EQ.
'STRUCFUN')
THEN
1991 istruf=istruf-istrut*100
1993 WRITE(6,*)
' ISTRUF,ISTRUT ',istruf,istrut
2007 ELSEIF(codewd.EQ.
'PSHOWER ')
THEN
2035 1070
FORMAT(//,5
x,a80,//)
2050 IF(sdum.EQ.blank)
THEN
2056 IF(ip.EQ.1) ijproj=1
2057 IF(ip.EQ.1) ijprox=1
2063 IF(ip.EQ.1) jjproj=1
2064 IF(ip.EQ.1) jjprox=1
2067 IF(sdum.EQ.btype(ii))
THEN
2070 ibproj=iibar(ijproj)
2075 jbproj=iibar(ijproj)
2081 WRITE(6,
'(A)')
' WRONG STRUCTURE OF PROJPAR CARD'
2098 IF(sdum.EQ.blank)
THEN
2105 IF(sdum.EQ.btype(ii))
THEN
2115 WRITE(6,
'(A)')
' WRONG STRUCTURE OF TARPAR CARD'
2131 IF(ijproj.NE.0) nnpp=ijproj
2132 ppn=
sqrt((epn-aam(nnpp))*(epn+aam(nnpp)))
2139 pproj =
sqrt((epn-amproj)*(epn+amproj))
2140 umo =
sqrt(amproj**2 + amtar**2 + 2.*amtar*eproj)
2143 ptthr=2.1+0.15*(log10(cmener/50.))**3
2145 ELSEIF(istrut.EQ.2)
THEN
2146 ptthr=2.5+0.12*(log10(cmener/50.))**3
2149 gamcm = (eproj+amtar)/umo
2152 pcm=gamcm*pproj - bgcm*eproj
2154 print 1033, eproj,pproj,
2155 +amproj,amtar,umo,gamcm,bgcm,pcm
2156 1033
FORMAT(
' CMS: ' ,
2157 +
' EPROJ,PPROJ,AMPROJ,AMTAR,UMO,GAMCM,BGCM,PCM'/8e22.13)
2178 IF (ihadrz.GE.2)
THEN
2197 auauau=rd2in(iseed1,iseed2)
2215 IF (
what(1).EQ.1.d0)
THEN
2222 IF(fermod.LT.0.0d0.OR.scafer.GT.2.0d0) scafer=1.0d0
2240 1080
FORMAT (
' THIS FILE CONTAINS EVENTS FROM KKEVT ')
2260 IF (ireso.EQ.1) CALL
distrp(1,ijproj,ppn)
2261 IF (cmhis.EQ.0.d0) CALL
distr(1,ijproj,ppn,idummy)
2262 IF (cmhis.EQ.1.d0) CALL
distrc(1,ijproj,ppn,idummy)
2263 IF (cmhis.EQ.2.d0) CALL
distco(1,ijproj,ppn,idummy)
2264 IF (ireso.EQ.1) CALL
disres(1,nhkkh1,ppn)
2265 IF (ipadis) CALL
distpa(1)
2266 IF (ioudif.EQ.1) CALL
diadif(1,0)
2267 CALL
shmak(1,nn,np,
nt,ip,it,umo,bimp)
2268 CALL
shmak1(1,nn,np,
nt,ip,it,umo,bimp)
2270 WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
2273 * ,form=
'UNFORMATTED')
2277 WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
2288 IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN
2289 CALL
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
2291 IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN
2292 CALL
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
2301 IF(lpauli .AND. (.NOT.fermp))
THEN
2303 +
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
2304 +
' IF FERMI ACTIVE',
' LPAULI CHANGED TO .FALSE.'
2311 IF(nevnts.LE.0) nevnts=1000
2315 IF(ncases.LE.0) ncases=100
2320 IF(iglaub.NE.1) iglaub=0
2323 IF(iglaub.EQ.1)
THEN
2326 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
2329 CALL
shmakf(ip,ipz,it,itz)
2338 IF(ipim.EQ.2)CALL
prblm2(cmener)
2407 ptthr=2.1+0.15*(log10(cmener/50.))**3
2409 ELSEIF(istrut.EQ.2)
THEN
2410 ptthr=2.5+0.12*(log10(cmener/50.))**3
2428 IF (
what(1).EQ.1.d0)
THEN
2446 IF (
what(1).EQ.1.d0)
THEN
2489 IF(
what(1).GT.0.5d0) ireso=1
2502 IF (
what(1).EQ.1.d0)
THEN
2524 270
WRITE(6,1120)
title
2577 OPEN(29,
file=
'lepto.evt',
2581 IF (ireso.EQ.1) CALL
distrp(1,ijproj,ppn)
2582 IF (cmhis.EQ.0.d0) CALL
distr(1,ijproj,ppn,idummy)
2583 IF (cmhis.EQ.1.d0) CALL
distrc(1,ijproj,ppn,idummy)
2584 IF (cmhis.EQ.2.d0) CALL
distco(1,ijproj,ppn,idummy)
2585 IF (ireso.EQ.1) CALL
disres(1,nhkkh1,ppn)
2586 IF (ipadis) CALL
distpa(1)
2587 IF (ioudif.EQ.1) CALL
diadif(1,0)
2589 WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
2592 * ,form=
'UNFORMATTED')
2596 WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
2600 WRITE(6,*)
' NEUTRINO: after INCINI call'
2604 IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN
2605 CALL
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
2607 IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN
2608 CALL
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
2617 IF(lpauli .AND. (.NOT.fermp))
THEN
2619 +
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
2620 +
' IF FERMI ACTIVE',
' LPAULI CHANGED TO .FALSE.'
2626 IF(nevnts.LE.0) nevnts=1000
2630 IF(ncases.LE.0) ncases=100
2636 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
2699 WRITE(6,*)
' NEUTRINO initialization finished'
2716 IF (
what(1).EQ.1.d0)
THEN
2746 IF (
what(1).EQ.1.d0)
THEN
2786 IF (
what(1).EQ.1.d0)
THEN
2803 IF (
what(1).EQ.1.d0)
THEN
2820 IF(ijproj.NE.0) nnpp=ijproj
2821 epn=
sqrt(ppn**2+aam(nnpp)**2)
2827 pproj =
sqrt((epn-amproj)*(epn+amproj))
2828 umo =
sqrt(amproj**2 + amtar**2 + 2.*amtar*eproj)
2831 ptthr=2.1+0.15*(log10(cmener/50.))**3
2833 ELSEIF(istrut.EQ.2)
THEN
2834 ptthr=2.5+0.12*(log10(cmener/50.))**3
2837 gamcm = (eproj+amtar)/umo
2840 pcm=gamcm*pproj - bgcm*eproj
2842 print 1033, eproj,pproj,
2843 +amproj,amtar,umo,gamcm,bgcm,pcm
2857 IF (
what(1).EQ.1.d0)
THEN
2917 IF(unon.LT.0.1d0) unon=2.0
2919 IF(unom.LT.0.1d0) unom=1.5
2921 IF(unosea.LT.0.1d0) unosea=2.0
2967 IF(cvq.LT.0.5d0) cvq=1.0
2969 IF(cdq.LT.1.0d0) cdq=2.0
2971 IF(csea.LT.0.1d0) csea =0.1
2973 IF(ssmima.LT.0.0d0) ssmima=0.14
2975 IF(
what(5).GT.2.0d0) vvmthr=
what(5)
2993 IF(inthad.LT.0 .OR. inthad.GT.2) inthad=0
2994 IF(inthad.EQ.1)
WRITE(6,
'(/5X,A/)')
2995 +
' FHAD: INELASTIC INTERACTION FORCED'
2996 IF(inthad.EQ.2)
WRITE(6,
'(/5X,A/)')
2997 +
' FHAD: ELASTIC INTERACTION FORCED'
3056 Write(47,473)ip,ipz,it,itz
3057 473
FORMAT(
' NUCLEUS ',4i10)
3070 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3073 WRITE(47,
'(4F10.5)') bmax,bstep,rproj,rtarg
3075 WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
3083 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3084 WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
3118 ptthr=2.1+0.15*(log10(cmener/50.))**3
3120 ELSEIF(istrut.EQ.2)
THEN
3121 ptthr=2.5+0.12*(log10(cmener/50.))**3
3128 IF(ijproj.NE.0) nnpp=ijproj
3129 epn=(cmener**2 + aam(nnpp)**2 - aam(1)**2)/(2.*aam(1))
3130 ppn=
sqrt((epn-aam(nnpp))*(epn+aam(nnpp)))
3131 ecm =
sqrt(aam(idp)**2+aam(1)**2+2.0d0*aam(1)*epn)
3136 pproj =
sqrt((epn-amproj)*(epn+amproj))
3137 eproj=
sqrt(pproj**2+amproj**2)
3141 gamcm = (eproj+amtar)/umo
3143 gamcm = (eproj+aam(1))/umo
3145 pcm=gamcm*pproj - bgcm*eproj
3146 print 1033, eproj,pproj,
3147 +amproj,amtar,umo,gamcm,bgcm,pcm
3179 IF (
what(1).EQ.1.d0) irecom=1
3180 IF (
what(1).NE.1.d0) irecom=0
3181 IF (
what(1).EQ.1.d0) lseadi=.true.
3215 IF (
what(1).EQ.1.d0) ifinal=1
3216 IF (
what(1).EQ.0.d0) ifinal=0
3280 IF(ifluct.EQ.1)CALL
fluini
3293 IF(
what(3).GT.0.d0)
THEN
3331 Write(47,1473)ip,ipz,it,itz
3332 1473
FORMAT(
' NUCLEUS ',4i10)
3341 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3344 WRITE(47,
'(4F10.5)') bmax,bstep,rproj,rtarg
3346 WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
3354 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3355 WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
3395 WRITE(6,
'(A,F10.2)')
' Threshold for HADRIN events = (GeV)',
3479 IF (
nint(
what(1)) .GE. 10000 )
THEN
3492 ELSE IF (
nint(whtsav) .NE. 0 )
THEN
3498 IF ( abs(
nint(
what(2))) .GE. 10 )
THEN
3502 ELSE IF (
nint(
what(2)) .NE. 0 )
THEN
3523 IF(
what(1).EQ.0.)
THEN
3550 IF (
nint(
what(2)) .LT. 0 ) ldeexg = .false.
3610 IF(
what(2).GE.0.1d0)
THEN
3651 OPEN(29,
file=
'qel.evt',
3655 IF (ireso.EQ.1) CALL
distrp(1,ijproj,ppn)
3656 IF (cmhis.EQ.0.d0) CALL
distr(1,ijproj,ppn,idummy)
3657 IF (cmhis.EQ.1.d0) CALL
distrc(1,ijproj,ppn,idummy)
3658 IF (cmhis.EQ.2.d0) CALL
distco(1,ijproj,ppn,idummy)
3659 IF (ireso.EQ.1) CALL
disres(1,nhkkh1,ppn)
3660 IF (ipadis) CALL
distpa(1)
3661 IF (ioudif.EQ.1) CALL
diadif(1,0)
3663 WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
3666 * ,form=
'UNFORMATTED')
3670 WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
3678 WRITE(6,*)
' NEUTRINO: after INCINI call'
3682 IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN
3683 CALL
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
3685 IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN
3686 CALL
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
3695 IF(lpauli .AND. (.NOT.fermp))
THEN
3697 +
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
3698 +
' IF FERMI ACTIVE',
' LPAULI CHANGED TO .FALSE.'
3705 IF(nevnts.LE.0) nevnts=1000
3709 IF(ncases.LE.0) ncases=100
3714 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3777 WRITE(6,*)
' NEUTRINO initialization finished'
3791 OPEN(29,
file=
'diffnuc.evt',
3795 IF (ireso.EQ.1) CALL
distrp(1,ijproj,ppn)
3796 IF (cmhis.EQ.0.d0) CALL
distr(1,ijproj,ppn,idummy)
3797 IF (cmhis.EQ.1.d0) CALL
distrc(1,ijproj,ppn,idummy)
3798 IF (cmhis.EQ.2.d0) CALL
distco(1,ijproj,ppn,idummy)
3799 IF (ireso.EQ.1) CALL
disres(1,nhkkh1,ppn)
3800 IF (ipadis) CALL
distpa(1)
3801 IF (ioudif.EQ.1) CALL
diadif(1,0)
3803 WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
3806 * ,form=
'UNFORMATTED')
3810 WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
3821 IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN
3822 CALL
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
3824 IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN
3825 CALL
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
3834 IF(lpauli .AND. (.NOT.fermp))
THEN
3836 +
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
3837 +
' IF FERMI ACTIVE',
' LPAULI CHANGED TO .FALSE.'
3844 IF(nevnts.LE.0) nevnts=1000
3848 IF(ncases.LE.0) ncases=100
3853 CALL
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
3935 WRITE(6,*)
'call xsglau'
3936 CALL
xsglau(ip,it,ijproj,1)
4008 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4018 COMMON /paname/ btype(30)
4020 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
4021 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
4022 +prebin,taebin,fermod,etacou
4024 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4026 LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
4028 COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
4029 +ipadis,ishmal,lpauli
4031 COMMON /hadthr/ ehadth,inthad
4033 COMMON /nshmak/ nnshma,npshma,ntshma,nshmac,nshma2
4035 COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
4036 +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
4037 +irvs14, irvv11,irvv12,irvv13,irvv14
4039 COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
4040 * bsite(0:1,200),nstatb,nsiteb
4043 COMMON /damp/ ca,ci,ga
4047 COMMON /diffra/ isingd,idiftp,ioudif,iflagd
4048 COMMON /nomije/ ptmije(10),nnmije(10)
4049 DATA ptmije /5.d0,7.d0,9.d0,11.d0,13.d0,15.d0,17.d0
4050 +,19.d0,21.d0,23.d0 /
4053 DATA irco1,irco2,irco3,irco4,irco5 /5*0/
4054 DATA irss11,irss12,irss13,irss14,irsv11,irsv12,irsv13,irsv14 /8*0/
4055 DATA irvs11,irvs12,irvs13,irvs14,irvv11,irvv12,irvv13,irvv14 /8*0/
4059 DATA prepot /210*0.0/
4060 DATA taepot /210*0.0/
4061 DATA taebin,prebin,fermod /2*0.0d0,0.6d0/
4063 DATA btype /
'PROTON ' ,
'APROTON ' ,
'ELECTRON' ,
'POSITRON' ,
4064 +
'NEUTRIE ' ,
'ANEUTRIE' ,
'PHOTON ' ,
'NEUTRON ' ,
'ANEUTRON' ,
4065 +
'MUON+ ' ,
'MUON- ' ,
'KAONLONG' ,
'PION+ ' ,
'PION- ' ,
4066 +
'KAON+ ' ,
'KAON- ' ,
'LAMBDA ' ,
'ALAMBDA ' ,
'KAONSHRT' ,
4067 +
'SIGMA- ' ,
'SIGMA+ ' ,
'SIGMAZER' ,
'PIZERO ' ,
'KAONZERO' ,
4068 +
'AKAONZER' ,
'RESERVED' ,
'BLANK ' ,
'BLANK ' ,
'BLANK ' ,
4071 DATA ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr /0, 0, 0, -1, 0,
4074 DATA intpt, fermp, ihadss,ihadsv,ihadvs,ihadvv, ihada /.true.,
4075 +.true., 4*.false., .true./
4076 DATA ipadis, ishmal, lpauli /.false., .false., .true./
4081 DATA nstatb, nsiteb /2000, 200/
4085 DATA isingd,idiftp,ioudif,iflagd /0,0,0,0/
4106 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4109 CHARACTER*8 projty,targty
4112 COMMON /user1/
title,projty,targty
4113 COMMON /user2/cmener,sdfrac,ptlar,istruf,isingd,idubld
4117 COMMON /colle/nevhad,nvers,ihadrz,nfile
4121 common/collis/
s,ijproj,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
4128 common/booklt/btype(30),nbook(30)
4133 COMMON /polmn0/pdifr,phard,psoft,alfah,betah,
4134 * sigtot,sigqel,sigel,sigine,sighin,sigd,sigdd
4138 COMMON /pomtyp/ipim,icon,isig,lmax,mmax,
nmax,difel,difnu
4143 COMMON /dropjj/dropjt,dropva
4144 COMMON /gluspl/nugluu,nsgluu
4145 COMMON /ptlarg/xsmax
4146 COMMON /ptsamp/ isampt
4147 COMMON /stars/istar2,istar3
4148 COMMON /strufu/istrum,istrut
4149 COMMON /popcor/pdb,ajsdef
4154 CHARACTER*8 codewd,sdum
4164 9
FORMAT(
' special code word was used ')
4242 IF(codewd.EQ.
'SIGMAPOM')
THEN
4261 IF (itest.EQ.1)CALL
pomdi
4275 ELSEIF(codewd.EQ.
'GLUSPLIT')
THEN
4298 ELSEIF(codewd.EQ.
'PARTEV ')
THEN
4301 IF (
what(2).EQ.0.d0)npev=30
4303 IF (
what(3).EQ.0.d0)nvers=1
4308 IF(ipim.EQ.2)CALL
prblm2(cmener)
4330 ELSEIF(codewd.EQ.
'SELHARD ')
THEN
4335 IF(
what(5).NE.0.d0)
THEN
4337 IF(cmener.LT.2000.0d0.AND.isig.EQ.3)ptthr=
what(5)
4338 IF (cmener.GE.2000.0d0.AND.isig.EQ.3)
4339 * ptthr=0.25*
log(cmener/2000.)+2.
4340 IF(ptthr2.LT.ptthr)ptthr2=ptthr
4342 ptthr=2.1+0.15*(log10(cmener/50.))**3
4344 ELSEIF(istrut.EQ.2)
THEN
4345 ptthr=2.5+0.12*(log10(cmener/50.))**3
4349 1244
FORMAT (
' THRESHOLD PT FOR HARD SCATTERING PTTHR=',f12.2)
4364 ELSEIF(codewd.EQ.
'XSLAPT ')
THEN
4379 ELSEIF(codewd.EQ.
'SAMPT ')
THEN
4382 IF( isampt.LT.0 .OR. isampt.GT.4 ) isampt=0
4420 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4423 common/booklt/btype(30),nbook(30)
4425 DATA btype /
'PROTON ' ,
'APROTON ' ,
'ELECTRON' ,
4426 1
'POSITRON' ,
'NEUTRIE ' ,
'ANEUTRIE' ,
4427 2
'PHOTON ' ,
'NEUTRON ' ,
'ANEUTRON' ,
4428 3
'MUON+ ' ,
'MUON- ' ,
'KAONLONG' ,
4429 4
'PION+ ' ,
'PION- ' ,
'KAON+ ' ,
4430 5
'KAON- ' ,
'LAMBDA ' ,
'ALAMBDA ' ,
4431 6
'KAONSHRT' ,
'SIGMA- ' ,
'SIGMA+ ' ,
4432 7
'SIGMAZER' ,
'PIZERO ' ,
'KAONZERO' ,
4433 9
'AKAONZER' ,
' ' ,
' ' ,
4437 DATA nbook / 2212 , -2212 , 11 ,
4439 2 22 , 2112 , -2112 ,
4441 4 211 , -211 , 321 ,
4442 5 -321 , 3122 , -3122 ,
4443 6 310 , 3114 , 3224 ,
4444 7 3214 , 111 , 311 ,
4461 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4463 parameter(
zero=0.d0,
one=1.d0)
4464 parameter( alfa=0.56268
d-01,
beta=0.17173
d+03 )
4465 parameter( acc = 0.0001d0 )
4466 COMMON /xsecpt/ ptcut,sigs,dsigh
4467 COMMON /sigma / sigsof,bs,zsof,sighar,
fill(7)
4468 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
4470 common/collis/
s,ijproj,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
4472 CHARACTER*8 projty,targty
4475 COMMON /user1/
title,projty,targty
4476 COMMON /user2/cmener,sdfrac,ptlar,istruf,isingd,idubld
4478 common/ptsamp/ isampt
4479 dimension pptt(50),dpptt(50)
4483 ptthr=2.5+0.12*(log10(cmener/50.))**3
4486 IF ( mode.EQ.0 )
THEN
4488 pptt(ii)=ii*ptcut/50.
4492 IF(ecm.LT.1000.)
THEN
4493 aacucu=0.85*(ecm-400.)/600.
4494 sigs=(1.-aacucu)*sigsof
4506 5559
FORMAT(
' SAMPPT:PTCUT,SIDSOF.SIGHATD,ISAMPT:',3e12.3,i5)
4508 IF( isampt.EQ.0 )
THEN
4509 c = dsigh/(2.*sigs*ptcut)
4511 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4512 * ,c,sigsof,sighar,rmin
4513 ELSEIF( isampt.EQ.1 )
THEN
4517 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4518 * ,c,sigsof,sighar,rmin
4519 ELSEIF( isampt.EQ.2 )
THEN
4521 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4522 * ,c,sigsof,sighar,rmin
4523 ELSEIF( isampt.EQ.3 )
THEN
4525 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4526 * ,c,sigsof,sighar,rmin
4527 ELSEIF( isampt.EQ.4)
THEN
4528 aaaa=ptcut**2*(sigsof+sighar)
4529 IF (aaaa.LE.0.00001d0)
THEN
4530 aaaa=abs(aaaa)+0.0002
4535 b = 0.5*
bsofpt(acc,c,ptcut)
4536 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4537 * ,c,sigsof,sighar,rmin
4541 rmin =
exp(b*ptcut**2)
4545 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4546 * ,c,sigsof,sighar,rmin
4547 9010
FORMAT(
' SAMPPT MODE,ISAMPT,PTCUT,SIGS,DSIGH,B,C,SIGSOF',
4551 ELSEIF ( mode.EQ.1 )
THEN
4552 IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,b
4553 * ,c,sigsof,sighar,rmin
4554 ptt =
log(1.0-
rndm(v)*(1.0-rmin))/(b+0.00001d0)
4556 iipt=
pt*50./ptcut+1.
4558 dpptt(iipt)=dpptt(iipt)+1./(
pt+0.000001d0)
4561 ELSEIF(mode.EQ.2)
THEN
4563 dpptt(ii)=log10(1.
e-8+dpptt(ii))
4565 IF(iouxev.GE.-1)
THEN
4567 203
FORMAT(
' PT DISTRIBUTION OF SOFT PARTONS AS SAMPLED IN BSOFPT')
4576 *
FUNCTION bsofpt(ACC,CC,PPTCUT)
4577 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4580 COMMON /bsoff1/c,ptcut
4581 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
4583 dimension
x(50),
y(50)
4591 IF(c.LT.1.
d-10)
THEN
4610 IF (.NOT.succes)
THEN
4611 IF (kkkk.EQ.0)go to 400
4614 IF(iouxev.GE.1)
WRITE(6,111)b1,b2
4625 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4627 COMMON /bsoff1/c,ptcut
4630 df=c*ptcut**2*aaa-aaa
4640 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4642 COMMON /bsoff1/c,ptcut
4645 IF(qqq.GT.-60.)
THEN
4648 bsof1=c*(aaa-1.)-b*aaa
4657 *
FUNCTION rtsafe(FUNCD,X1,X2,XACC)
4658 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4660 parameter(maxit=200,itepri=0)
4661 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
4662 CALL funcd(
x1,fl,df)
4664 CALL funcd(
x2,fh,df)
4666 IF(fl*fh.GE.0.) pause
'ROOT MUST BE BRACKETED'
4692 * .OR. abs(2.*
f).GT.abs(dxold*df) )
THEN
4705 IF(abs(
dx).LT.xacc)
RETURN
4715 pause
'RTSAFE EXCEEDING MAXIMUM ITERATIONS'
4717 9995
FORMAT(
' VR1,VR2:',2e12.5)
4718 9996
FORMAT(
' RTSAFE,XH,XL,DXOLD,F,DF IN LOOP 11 J=1,MAXIT')
4719 9997
FORMAT(3
x,6e10.3)
4720 9998
FORMAT(
' RTSAFE: RTSAFE,F,DF =',3e12.5)
4721 9999
FORMAT(
' RTSAFE: F,DF =',2e12.5)
4727 *
FUNCTION var(A,B,C,D)
4728 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4730 parameter( ambmax = 1.0
d+38, epsi = 1.2
d-38,
one=1.d0 )
4734 abl = log10( abl + epsi )
4737 ccl = log10( ccl + epsi )
4739 IF( rcheck .LE. 38.d0 )
THEN
4742 var = ambmax*siab*sicc -
d
4744 IF(
var .GT. 1.0
d+18 )
var = 1.0
e+18
4745 IF(
var .LT. -1.0
d+18 )
var = -1.0
e+18
4750 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
4753 parameter(factor=1.6d0,ntry=50)
4755 IF(
x1.EQ.
x2)pause
'You have to guess an initial range'
4760 IF(
f1*
f2.LT.0.d0)
RETURN
4761 IF(abs(
f1).LT.abs(
f2))
THEN
4794 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4798 parameter(amuamu=0.93149432d0)
4801 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
4889 COMMON /delp/ delpx,delpy,delpz,delpe
4891 COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
4893 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
4895 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4900 IF(help.GT.5.d0)phelp=help-5.
4915 IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13)
THEN
4920 pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
4921 pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
4924 IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14)
THEN
4931 IF(isthkk(i).EQ.1)
THEN
4938 IF(isthkk(i).EQ.13.AND.jdahkk(2,i).EQ.0)
THEN
4943 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
4944 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
4946 IF(isthkk(i).EQ.14.AND.jdahkk(2,i).EQ.0)
THEN
4953 IF(isthkk(i).EQ.16)
THEN
4959 eext=eext + phkk(4,i) - phkk(4,imo)
4961 IF(isthkk(i).EQ.15)
THEN
4967 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
4968 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
4969 eext=eext + phkk(4,i) - phkk(4,imo)
4971 IF(isthkk(i).EQ.1)
THEN
4974 IF(isthkk(i).EQ.-1)
THEN
4975 eeem1=eeem1+phkk(4,i)
4977 IF(isthkk(i).EQ.1001)
THEN
4978 ee1001=ee1001+phkk(4,i)
4981 eee=eee1+eeem1+ee1001
4992 aip=aip+(ait*amuamu+1.
d-3*
energy(ait,aitz))/epnto
4998 IF(it.EQ.ip)tole=0.02
5001 IF(delle.GE.tole)irej=1
5004 IF(icheck.LE.100)
THEN
5005 WRITE(6,
'(A,I5,E10.3,5F10.4)')
5006 *
' IP,EPN,AEEE,AEEEE,AEEE1,AEEEM1,AEE101:',
5007 * ip,epn,aeee,aeeee,aeee1,aeeem1,aee101
5008 WRITE(6,
'(A,I5,E10.3,7E12.4)')
5009 *
' IP,EPN,EEE,EEEE,EEE1,EEEM1,EE1001,DELLE,ELLE:',
5010 * ip,epn,eee,eeee,eee1,eeem1,ee1001,delle,elle
5026 1000
FORMAT(
' CHECKF: PX,PY,PZ,PE,EEXT,EEXP',2f7.3,2e12.3,2f7.3
5027 * / 8
x,
' DELPX/Y/Z/E',4f7.3,i10,
' IORIG')
5033 IF(isthkk(i).EQ.11)
THEN
5034 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5035 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5036 + (vhkk(khkk,i),khkk=1,4)
5039 IF(isthkk(i).EQ.12)
THEN
5040 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5041 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5042 + (vhkk(khkk,i),khkk=1,4)
5045 IF(isthkk(i).EQ.1)
THEN
5046 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5047 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5048 + (vhkk(khkk,i),khkk=1,4)
5050 1010
FORMAT (i6,i4,5i6,9(1pe10.2))
5052 IF(isthkk(i).EQ.16)
THEN
5054 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5055 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5056 + (vhkk(khkk,i),khkk=1,4)
5069 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5073 parameter(amuamu=0.93149432d0)
5076 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
5164 COMMON /delp/ delpx,delpy,delpz,delpe
5166 COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
5168 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
5170 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
5175 IF(help.GT.5.d0)phelp=help-5.d0
5176 pthelp=12.d0+phelp*5.d0
5196 IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13)
THEN
5201 pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
5202 pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
5205 IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14)
THEN
5212 IF(isthkk(i).EQ.1)
THEN
5219 IF(isthkk(i).EQ.13.AND.jdahkk(2,i).EQ.0)
THEN
5224 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
5225 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
5227 IF(isthkk(i).EQ.14.AND.jdahkk(2,i).EQ.0)
THEN
5234 IF(isthkk(i).EQ.16)
THEN
5240 eext=eext + phkk(4,i) - phkk(4,imo)
5242 IF(isthkk(i).EQ.15)
THEN
5248 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
5249 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
5250 eext=eext + phkk(4,i) - phkk(4,imo)
5252 IF(isthkk(i).EQ.1)
THEN
5257 IF(isthkk(i).EQ.-1)
THEN
5258 eeem1=eeem1+phkk(4,i)
5262 IF(isthkk(i).EQ.1001)
THEN
5263 ee1001=ee1001+phkk(4,i)
5264 pz1001=pz1001+phkk(3,i)
5265 px1001=px1001+phkk(1,i)
5268 eee=eee1+eeem1+ee1001
5269 pzpz=pz1+pzm1+pz1001
5270 pxpx=px1+pxm1+px1001
5281 IF(isthkk(i).EQ.1001)
THEN
5282 phkk(3,i)=phkk(3,i)+delpz
5283 phkk(4,i)=
sqrt(phkk(1,i)**2+phkk(2,i)**2+phkk(3,i)**2
5285 ee1001=ee1001+phkk(4,i)
5288 eee=eee1+eeem1+ee1001
5296 bip=epn+(ait*amuamu+1.
d-3*
energy(ait,aitz))
5302 IF(delle.GE.tole)irej=1
5305 IF(icheck.LE.20)
THEN
5306 WRITE(6,
'(A,I5,E10.3,4F10.4)')
5307 *
' IP,EPN,PXPX,PX1,PXM1,PX1001:',
5308 * ip,epn,pxpx,px1,pxm1,px1001
5309 WRITE(6,
'(A,I5,E10.3,6F10.4)')
5310 *
' IP,PPN,PZPZ,PZ1,PZM1,PZ1001,BIP,BMI:',
5311 * ip,ppn,pzpz,pz1,pzm1,pz1001,bip,bmi
5312 WRITE(6,
'(A,I5,E10.3,5E12.4)')
5313 *
' IP,EPN,EEE,EEE1,EEEM1,EE1001,DELLE:',
5314 * ip,epn,eee,eee1,eeem1,ee1001,delle
5330 1000
FORMAT(
' CHECKF: PX,PY,PZ,PE,EEXT,EEXP',2f7.3,2e12.3,2f7.3
5331 * / 8
x,
' DELPX/Y/Z/E',4f7.3,i10,
' IORIG')
5337 IF(isthkk(i).EQ.11)
THEN
5338 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5339 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5340 + (vhkk(khkk,i),khkk=1,4)
5343 IF(isthkk(i).EQ.12)
THEN
5344 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5345 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5346 + (vhkk(khkk,i),khkk=1,4)
5349 IF(isthkk(i).EQ.1)
THEN
5350 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5351 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5352 + (vhkk(khkk,i),khkk=1,4)
5354 1010
FORMAT (i6,i4,5i6,9(1pe10.2))
5356 IF(isthkk(i).EQ.16)
THEN
5358 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5359 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5360 + (vhkk(khkk,i),khkk=1,4)
5372 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5376 parameter(amuamu=0.93149432d0)
5379 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
5466 COMMON /zentra/ icentr
5468 COMMON /delp/ delpx,delpy,delpz,delpe
5470 COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
5472 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
5474 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
5479 IF(help.GT.5.d0)phelp=help-5.
5494 IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13)
THEN
5499 pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
5500 pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
5503 IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14)
THEN
5510 IF(isthkk(i).EQ.1)
THEN
5517 IF(isthkk(i).EQ.13.AND.jdahkk(2,i).EQ.0)
THEN
5522 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
5523 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
5525 IF(isthkk(i).EQ.14.AND.jdahkk(2,i).EQ.0)
THEN
5532 IF(isthkk(i).EQ.16)
THEN
5538 eext=eext + phkk(4,i) - phkk(4,imo)
5540 IF(isthkk(i).EQ.15)
THEN
5546 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
5547 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
5548 eext=eext + phkk(4,i) - phkk(4,imo)
5550 IF(isthkk(i).EQ.1)
THEN
5580 IF(ip.EQ.it.AND.it.GT.1)tole=0.05d0*ip
5582 IF(epn.LE.5.d0)tole=3.d0*tole
5584 IF (abs(
px).GT.pthelp.OR. abs(
py).GT.pthelp.OR.
5585 * abs(
pz)/epn.GT.tole.
5586 + or. abs(pe)/epn.GT.tole)
THEN
5589 IF(icheck.LE.50.AND.irej.EQ.1)
THEN
5590 WRITE(6,1000)
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe,iorig
5592 1000
FORMAT(
' CHECKO: PX,PY,PZ,PE,EEXT,EEXP',2f7.3,2e12.3,2f7.3
5593 * / 8
x,
' DELPX/Y/Z/E',4f7.3,i10,
' IORIG')
5597 WRITE(6,1000)
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe,iorig
5601 IF(isthkk(i).EQ.11)
THEN
5602 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5603 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5604 + (vhkk(khkk,i),khkk=1,4)
5607 IF(isthkk(i).EQ.12)
THEN
5608 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5609 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5610 + (vhkk(khkk,i),khkk=1,4)
5613 IF(isthkk(i).EQ.1)
THEN
5614 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5615 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5616 + (vhkk(khkk,i),khkk=1,4)
5618 1010
FORMAT (i6,i4,5i6,9(1pe10.2))
5620 IF(isthkk(i).EQ.16)
THEN
5622 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5623 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5624 + (vhkk(khkk,i),khkk=1,4)
5631 WRITE(6,1000)
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe,iorig
5637 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5643 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
5731 COMMON /delp/ delpx,delpy,delpz,delpe
5733 COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
5735 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
5745 IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13)
THEN
5750 pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
5751 pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
5754 IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14)
THEN
5761 IF(isthkk(i).EQ.1)
THEN
5768 IF(isthkk(i).EQ.13.AND.jdahkk(1,i).EQ.0)
THEN
5773 pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
5774 pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
5776 IF(isthkk(i).EQ.14.AND.jdahkk(1,i).EQ.0)
THEN
5783 IF(isthkk(i).EQ.16)
THEN
5789 eext=eext + phkk(4,i) - phkk(4,imo)
5791 IF(isthkk(i).EQ.15)
THEN
5797 eext=eext + phkk(4,i) - phkk(4,imo)
5804 WRITE(6,1000)
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe
5805 1000
FORMAT(
' CHECKE: PX,PY,PZ,PE,EEXT,EEXP',6f7.3/ 8
x,
' DELPX/Y/Z/E',4
5807 WRITE(6,
'(8X,A,6F8.3)')
' TASUMA,TASUBI,TABI,TAMASU,TAMA,TAIMMA',
5808 +tasuma,tasubi,tabi,tamasu,tama,taimma
5810 IF (abs(
px).GT.0.004.OR. abs(
py).GT.0.004.OR. abs(
pz).GT.0.004.
5811 + or. abs(pe).GT.0.004)
THEN
5815 IF(isthkk(i).EQ.11)
THEN
5816 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5817 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5818 + (vhkk(khkk,i),khkk=1,4)
5821 IF(isthkk(i).EQ.12)
THEN
5822 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5823 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5824 + (vhkk(khkk,i),khkk=1,4)
5827 IF(isthkk(i).EQ.1)
THEN
5828 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5829 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5830 + (vhkk(khkk,i),khkk=1,4)
5832 1010
FORMAT (i6,i4,5i6,9(1pe10.2))
5834 IF(isthkk(i).EQ.16)
THEN
5836 WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
5837 + (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
5838 + (vhkk(khkk,i),khkk=1,4)
5866 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5873 DATA a1,a2,a3,a4,a5 /0.01575, 0.0178, 0.000710, 0.0237, 0.034/
5876 IF(ia.LE.1.OR.
iz.EQ.0)
THEN
5881 ebind = a1*aa - a2*aa**0. 666667- a3*
iz*
iz*aa**(-0.333333) - a4
5883 IF (
mod(ia,2).EQ.1)
THEN
5885 ELSEIF (
mod(
iz,2).EQ.1)
THEN
5897 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5904 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
6002 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
6003 +iibar(210),k1(210),k2(210)
6006 COMMON /factmo/ ifacto
6008 COMMON /taufo/ taufor,ktauge,itauve,incmod
6010 COMMON /hadthr/ ehadth,inthad
6014 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
6015 COMMON /nuccc/ jt,jtz,jp,jpz,jjproj,jbproj,jjtarg,jbtarg
6017 COMMON /zentra/ icentr
6019 COMMON /cmhico/ cmhis
6021 COMMON /resona/ ireso
6023 COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
6030 COMMON /projk/ iprojk
6032 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
6033 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
6034 +prebin,taebin,fermod,etacou
6036 COMMON /recom/irecom
6046 ppn=
sqrt((epn-aam(ijproj))*(epn+aam(ijproj)))
6047 ibproj=iibar(ijproj)
6050 jbproj=iibar(ijproj)
6102 SUBROUTINE hadhad(EPN,PPN,NHKKH1,IHTAWW,ITTA,IREJFO)
6103 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6113 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
6201 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
6213 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
6214 +iibar(210),k1(210),k2(210)
6217 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
6219 parameter(maxfin=10)
6220 COMMON /dfinls/ itrh(maxfin),cxrh(maxfin),cyrh(maxfin), czrh
6221 +(maxfin),elrh(maxfin),plrh(maxfin),irh
6223 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
6224 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
6225 +prebin,taebin,fermod,etacou
6232 IF(ipri.GE.2)
WRITE(6,1001) ijproj,ijproj,ppn,epn,cccxp,cccyp,
6233 +ccczp,ihtaww,itta,ieline
6234 1001
FORMAT(
' HADHAD 1:',
6235 +
' IJPROJ,IJPROJ,PPN,EPN,CCCXP,CCCYP,CCCZP,IHTAWW,ITTA,IELINE'/ 2
6236 +i3,2e12.3,3f7.3,3i4)
6241 CALL
sihnin(ijproj,itta,ppn,sight)
6242 CALL
sihnel(ijproj,itta,ppn,sighte)
6243 sigtot=sight + sighte
6244 IF (sigtot*
rndm(bb).LE.sighte)ieline=1
6245 IF(ipri.GE.2)
WRITE(6,1000) ijproj,ijproj,ppn,epn,cccxp,cccyp,
6246 +ccczp,ihtaww,itta,ieline
6247 1000
FORMAT(
' HADHAD 2 nach si...:',
6248 +
' IJPROJ,IJPROJ,PPN,EPN,CCCXP,CCCYP,CCCZP,IHTAWW,ITTA,IELINE'/ 2
6249 +i3,2e12.3,3f7.3,3i4)
6253 IF(ipri.GE.2)
WRITE(6,1012) ijproj,ijproj,ppn,epn,cccxp,cccyp,
6254 +ccczp,ihtaww,itta,ieline
6255 1012
FORMAT(
' HADHAD 12 loop:',
6256 +
' IJPROJ,IJPROJ,PPN,EPN,CCCXP,CCCYP,CCCZP,IHTAWW,ITTA,IELINE'/ 2
6257 +i3,2e12.3,3f7.3,3i4)
6260 WRITE(6,
'(I3,5(1PE12.4),I5/3X,5(1PE12.4))') ii,elrh(ii),plrh
6261 + (ii),cxrh(ii),cyrh(ii),czrh(ii),itrh(ii), (phkk(jjj,nhkk),jjj
6266 CALL
fhad(ijproj,ijproj,ppn,epn,cccxp,cccyp,ccczp, ihtaww,itta,
6271 +
WRITE(6,
'(A)')
' exit from hadhad with irejfo=1 '
6277 IF (ihadha.LT.3)
THEN
6280 IF(itsec.EQ.1.AND.elrh(ii).LE.taefep+aam(itsec)) goto 12
6281 IF(itsec.EQ.8.AND.elrh(ii).LE.taefen+aam(itsec)) goto 12
6282 IF(iibar(itsec).NE.1.AND.elrh(ii)-aam(itsec)
6283 + .LE.taepot(itsec)) goto 12
6288 IF (ipri.GE.2)
WRITE (6,1010)irh,nhkkh1,ihtaww,itta
6289 1010
FORMAT (
' HADHAD IRH,NHKKH1,IHTAWW,ITTA = ',4i5)
6293 +
' HADHAD - PARTICLE TRANSFER FROM /FINLSP/ INTO /HKKEVT/',
6294 +
' II, ELRH, PLRH, CXRH, CYRH, CZRH / PHKK(1-5)'
6303 WRITE (6,
'(A,2I5)') .EQ.
' HADHAD:NHKKNMXHKK ',nhkk,
nmxhkk
6307 idhkk(nhkk)=
mpdgha(itsec)
6309 jmohkk(2,nhkk)=ihtaww
6312 phkk(1,nhkk)=plrh(ii)*cxrh(ii)
6313 phkk(2,nhkk)=plrh(ii)*cyrh(ii)
6314 phkk(3,nhkk)=plrh(ii)*czrh(ii)
6315 phkk(4,nhkk)=elrh(ii)
6316 IF(phkk(4,nhkk)-aam(itsec).LE.taepot(itsec).
6317 + and.iibar(itsec).EQ.1)
THEN
6322 phkk(5,nhkk)=aam(itrh(ii))
6325 WRITE(6,
'(I3,5(1PE12.4),I5/3X,5(1PE12.4),I5)')
6327 + (ii),cxrh(ii),cyrh(ii),czrh(ii),itrh(ii), (phkk(jjj,nhkk),jjj
6330 vhkk(1,nhkk)=vhkk(1,ihtaww)
6331 vhkk(2,nhkk)=vhkk(2,ihtaww)
6332 vhkk(3,nhkk)=vhkk(3,ihtaww)
6333 vhkk(4,nhkk)=vhkk(4,1)
6336 jdahkk(1,1)=nhkkh1+1
6338 jdahkk(1,ihtaww)=nhkkh1+1
6339 jdahkk(2,ihtaww)=nhkk
6343 +
WRITE(6,
'(A)')
' exit from hadhad with irejfo=0 '
6347 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6353 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
6453 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
6454 +iibar(210),k1(210),k2(210)
6457 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
6459 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
6461 COMMON /chabai/chargi,barnui
6462 COMMON /evappp/ievap
6477 IF (isthkk(i).EQ.13)
THEN
6479 IF(idhkk(i).EQ.2212)chaeve=chaeve+1.d0
6481 IF (isthkk(i).EQ.14)
THEN
6483 IF(idhkk(i).EQ.2212)chaeve=chaeve+1.d0
6486 DO 521 i=nhkkh1,nhkk
6487 IF (isthkk(i).EQ.1.OR.isthkk(i).EQ.15.OR.isthkk(i).EQ.16)
THEN
6490 IF (nrhkk.LE.0.OR.nrhkk.GT.410)
THEN
6491 WRITE(6,1389)nrhkk,i,idhkk(i),nhkkh1,nhkk
6492 1389
FORMAT (
' distr: NRHKK ERROR ',5i10)
6497 chaeve=chaeve+ichhkk
6501 ELSEIF(ievap.EQ.1)
THEN
6503 IF (isthkk(i).EQ.1)
THEN
6507 chaeve=chaeve+ichhkk
6513 IF (isthkk(i).EQ.-1)
THEN
6514 IF(idhkk(i).EQ.2112)
THEN
6518 IF(idhkk(i).EQ.2212)
THEN
6524 IF((idhkk(i).EQ.80000).AND.(isthkk(i).NE.1000))
THEN
6525 chaeve=chaeve+idxres(i)
6526 baeve=baeve+idres(i)
6531 IF(ievl.LE.10)
WRITE(6,
'(2A,4F10.2)')
' Event charge and B-number',
6532 *
'=',chaeve,baeve,chargi,barnui
6533 IF(chaeve-chargi.NE.0.d0.OR.baeve-barnui.NE.0.d0)
THEN
6535 IF(ievl.LE.1000)
WRITE(6,
'(2A,4F10.2)')
'Event charge and B-numb',
6536 *
'(violated) =',chaeve,baeve,chargi,barnui
6547 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6549 dimension
pt(50,10),ypt(50,10)
6555 pt(j,i)=j*dpt-dpt/2.
6562 IF(ipt1.GT.50)ipt1=50
6563 IF(ipt2.GT.50)ipt2=50
6564 ypt(ipt1,ipt)=ypt(ipt1,ipt)+1.
6565 ypt(ipt2,ipt)=ypt(ipt2,ipt)+1.
6566 ypt(ipt1,10)=ypt(ipt1,10)+1.
6567 ypt(ipt2,10)=ypt(ipt2,10)+1.
6572 ypt(j,i)=ypt(j,i)/nevt
6573 ypt(j,i)=log10(ypt(j,i)+1.
d-18)
6585 SUBROUTINE hkkfil(IST,ID,M1,M2,PX,PY,PZ,E,NHKKAU,KORMO,ICALL)
6587 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6589 parameter(lout=6,llook=9)
6590 parameter(tiny10=1.0
d-10,tiny4=1.0
d-3)
6596 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
6683 COMMON /nncms/ gacms,bgcms,umo,pcm,eproj,pproj
6684 COMMON /trafop/galab,bglab,blab
6685 COMMON /projk/ iprojk
6698 WRITE(lout,1000) nhkk
6699 1000
FORMAT(1
x,
'HKKFIL: NHKK exeeds NMXHKK = ',i7,
6700 &
'! program execution stopped..')
6703 IF (m1.LT.0) mo1 = nhkk+m1
6704 IF (m2.LT.0) mo2 = nhkk+m2
6707 IF(kormo.EQ.999)
THEN
6708 jmohkk(1,nhkk) = mo1
6709 jmohkk(2,nhkk) = mo2
6711 jmohkk(1,nhkk)=nhkkau+kormo-1
6714 IF(nhkk.LE.jmohkk(1,nhkk))
THEN
6716 WRITE(6,*)
' HKKFIL(IST,ID,M1,M2,PX,PY,PZ,E,NHKKAU,KORMO)',
6717 * nhkk,ist,id,m1,m2,
px,
py,
pz,
e,nhkkau,kormo,icall,jmohkk(1,nhkk)
6722 IF (jdahkk(1,mo1).NE.0)
THEN
6723 jdahkk(2,mo1) = nhkk
6725 jdahkk(1,mo1) = nhkk
6727 jdahkk(1,mo1)=nhkkau
6730 IF (jdahkk(1,mo2).NE.0)
THEN
6731 jdahkk(2,mo2) = nhkk
6733 jdahkk(1,mo2) = nhkk
6735 jdahkk(1,mo2) = nhkkau
6741 phkk(5,nhkk) = phkk(4,nhkk)**2-phkk(1,nhkk)**2-
6742 & phkk(2,nhkk)**2-phkk(3,nhkk)**2
6743 IF ((phkk(5,nhkk).LT.0.0d0).AND.(abs(phkk(5,nhkk)).GT.tiny4))
6744 &
WRITE(lout,
'(1X,A,G10.3)')
6745 &
'HKKFIL: negative mass**2 ',phkk(5,nhkk)
6746 phkk(5,nhkk) =
sqrt(abs(phkk(5,nhkk)))
6747 IF (ist.EQ.88888.OR.ist.EQ.88887.OR.ist.EQ.88889)
THEN
6753 vhkk(i,nhkk) = vhkk(i,mo2)
6755 vhkk(4,nhkk) = vhkk(3,mo2)/blab-vhkk(3,mo1)/bglab
6759 vhkk(i,nhkk) = vhkk(i,mo1)
6760 IF (iprojk.EQ.1)
THEN
6761 whkk(i,nhkk) = whkk(i,mo1)
6774 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6780 common/
pyjets/
n,npad,k(4000,5),
p(4000,5),v(4000,5)
6781 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
6782 common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
6783 common/
pydat3/mdcy(500,3),mdme(4000,2),brat(4000),kfdp(4000,5)
6785 CHARACTER chap*16,chan*16,chad(5)*16
6789 WRITE(mstu(11),6800)
6793 IF(mstu(2).NE.0) kfmax=mstu(2)
6798 IF(kc.EQ.0) goto 220
6799 IF(mstu(14).EQ.0.AND.kf.GT.100.AND.kc.LE.100) goto 220
6800 IF(mstu(14).GT.0.AND.kf.GT.100.AND.max(
mod(kf/1000,10),
6801 &
mod(kf/100,10)).GT.mstu(14)) goto 220
6809 IF(kf.LE.100.AND.chap.EQ.
' '.AND.mdcy(kc,2).EQ.0) goto 220
6814 idc2=mdcy(kc,2)+mdcy(kc,3)-1
6815 WRITE(mstu(11),6900)kbam,
6816 & kf,kc,idc1,idc2,chap,chan,kchg(kc,1),kchg(kc,2),
6817 & kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
6819 & kf,kc,idc1,idc2,chap,chan,kchg(kc,1),kchg(kc,2),
6820 & kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
6824 IF(kf.GT.100.AND.kc.LE.100) goto 220
6825 DO 210 idc=mdcy(kc,2),mdcy(kc,2)+mdcy(kc,3)-1
6828 CALL
pyname(kfdp(idc,j),chad(j))
6830 kbamdp(j)=
mcihad(kfdp(idc,j))
6831 IF(kbamdp(j).EQ.26)kbamdp(j)=0
6833 WRITE(26,7001) idc,mdme(idc,1),mdme(idc,2),brat(idc),
6835 210
WRITE(mstu(11),7000) idc,mdme(idc,1),mdme(idc,2),brat(idc),
6838 IF(kabam.NE.410)
THEN
6839 WRITE(mstu(11),6900)kabam,
6840 & -kf,-kc,idc1,idc2,chan,chap,-kchg(kc,1),kchg(kc,2),
6841 & kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
6842 WRITE(26,6900)kabam,
6843 & -kf,-kc,idc1,idc2,chan,chap,-kchg(kc,1),kchg(kc,2),
6844 & kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
6845 DO 211 idc=mdcy(kc,2),mdcy(kc,2)+mdcy(kc,3)-1
6849 IF(kcdp.LE.0.OR.kcdp.GT.500)
THEN
6855 IF(kchg(kcdp,3).EQ.0)kfdpm=kfdp(idc,j)
6857 IF(kbamdp(j).EQ.26)kbamdp(j)=0
6859 CALL
pyname(kfdpm,chad(j))
6861 WRITE(26,7001) idc,mdme(idc,1),mdme(idc,2),brat(idc),
6863 211
WRITE(mstu(11),7000) idc,mdme(idc,1),mdme(idc,2),brat(idc),
6871 6800
FORMAT(///30
x,
'Particle/parton data table'//1
x,
'BAM',
6872 &1
x,
'ABAM',1
x,
'KF',1
x,
'KC',1
x,
'DCF',1
x,
'DCL',1
x,
6873 &
'particle',8
x,
'antiparticle',6
x,
'chg col anti',8
x,
'mass',7
x,
6874 &
'width',7
x,
'w-cut',5
x,
'lifetime',1
x,
'decay'/11
x,
'IDC',1
x,
'on/off',
6875 &1
x,
'ME',3
x,
'Br.rat.',4
x,
'decay products')
6876 6900
FORMAT(/1
x,i4,i6,i4,2i5,a16,a16,3i3,1
x,f12.5,2(1
x,f11.5),
6878 7000
FORMAT(10
x,i4,2
x,i3,2
x,i3,2
x,f8.5,4
x,5a16)
6879 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 dpmevt(ELABT, IIPROJ, IIP, IIPZ, IIT, IITZ, KKMAT, NHKKH1)
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