1 subroutine NumberModel(cmodel,model)
4 if(cmodel(1:n).eq.'qgsjet' ) model=2
5 if(cmodel(1:n).eq.'gheisha') model=3
6 c if(cmodel(1:n).eq.'gheisha' )stop'wrong choice!!!!' ! model=4
7 if(cmodel(1:n).eq.'pythia' )stop'wrong choice!!!!' ! model=4
8 if(cmodel(1:n).eq.'hijing' )stop'wrong choice!!!!' ! model=5
9 if(cmodel(1:n).eq.'sibyll' ) model=6
10 if(cmodel(1:n).eq.'IIqgsjet' ) model=7
13 subroutine IniModel(model)
14 if(model.eq.2)call IniQGSjet
15 if(model.eq.3)call IniGheisha
16 cxxx if(model.eq.4)call IniPythia
17 cxxx if(model.eq.5)call IniHijing
18 if(model.eq.6)call IniSibyll
19 if(model.eq.7)call IniQGSJetII
22 subroutine IniEvtModel
24 if(model.eq.2)call IniEvtQGS
25 if(model.eq.3)call IniEvtGhe
26 cxxx if(model.eq.4)then
28 cxxx if(engy.lt.egymin)engy=egymin
34 c if(engy.lt.egymin)engy=egymin
38 if(model.eq.6)call IniEvtSib
39 if(model.eq.7)call IniEvtQGSII
42 subroutine emsaaaModel(model,iret)
44 if(idtargin.eq.0)call IniEvtQGS
47 if(model.eq.3)call emsghe(iret)
48 cxxx if(model.eq.4)call emspyt(iret)
49 c if(model.eq.5)call emshij(iret)
50 if(model.eq.6)call emssib(iret)
52 if(idtargin.eq.0)call IniEvtQGSII
57 function crseModel(model,ekin,maproj,matarg,idtarg)
58 if(model.eq.2)crseModel=qgscrse(ekin,maproj,matarg,idtarg)
59 cxxx if(model.eq.4)??????????????????
60 c if(model.eq.5)crseModel=hijcrse(ekin,maproj,matarg,idtarg)
61 if(model.eq.6)crseModel=sibcrse(ekin,maproj,matarg,idtarg)
62 if(model.eq.7)crseModel=qgsIIcrse(ekin,maproj,matarg,idtarg)
66 subroutine m2XXFZ( a,b)
71 subroutine m3SIGMA(ek,idpro,idtar,latar,matar,sigi,sige)
72 call ghecrse(ek,idpro,idtar,latar,matar,sigi,sige)
75 subroutine m6SIGMA(icl,engy,stot,sela,sine,sdifr,slela,Rho)
83 call SIB_SIGMA_HP(L,engy,stot,sela,sine,sdifr,slela,Rho)
87 subroutine m7SIGMA(stot,scut,sine,slela)
88 double precision GzZ0(5)
89 common /qgarr1/ ia(2),icz,icp
92 call qgfz(0.d0,gzz0,0,0)
93 scut=real(gzz0(2))/2. !cut pomerons cross-section
94 stot=real(gzz0(1)) !tot cross-section
95 sine=real(gzz0(2)+gzz0(3)+gzz0(4))/2. !inelastic cross section