IAP GITLAB

Commit fc2d1f70 authored by Colin Baus's avatar Colin Baus

sergey said it was ok to change it to 208. so qgsjetII is fine, however,...

sergey said it was ok to change it to 208. so qgsjetII is fine, however, epos-con.f was also changed

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@3194 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 8f725328
......@@ -310,7 +310,7 @@ CRMCoptions::GetOutputFileName() const
case 1 : rootFileName << "p"; break;
case -1 : rootFileName << "antip"; break;
case 12 : rootFileName << "C"; break;
case 207 : rootFileName << "Pb"; break;
case 208 : rootFileName << "Pb"; break;
default:
// update this
rootFileName << "pdg" ;
......@@ -318,10 +318,12 @@ CRMCoptions::GetOutputFileName() const
}
switch (fTargetId) {
case 120 : rootFileName << "pi"; break;
case -120 : rootFileName << "antipi"; break;
case 1 : rootFileName << "_p"; break;
case -1 : rootFileName << "_antip"; break;
case 12 : rootFileName << "_C"; break;
case 207 : rootFileName << "_Pb"; break;
case 208 : rootFileName << "_Pb"; break;
default:
// update this
rootFileName << "_" << "pdg" ;
......
......@@ -138,7 +138,7 @@ c define vec to boost from cm. to cms frame
boostvec3=dsinh(dble(ycm2det))
boostvec4=dcosh(dble(ycm2det))
boostvec5=1d0
write(*,*)nevhep,nhep,boostvec3,boostvec4,ycm2det
c write(*,*)nevhep,nhep,boostvec3,boostvec4,ycm2det
do i=1,nhep
c boost output to cms frame
call utlob2(-1,boostvec1,boostvec2,boostvec3
......
......@@ -775,7 +775,7 @@ c-----------------------------------------------------------------------
common/nucl3/phi,bimp
common /cncl/xproj(mamx),yproj(mamx),zproj(mamx)
*,xtarg(mamx),ytarg(mamx),ztarg(mamx)
parameter(iapmax=207)
parameter(iapmax=208)
double precision bqgs2,bmaxqgs2,bmaxnex2,bminnex2,xan,xbn
common /qgsIInex1/xan(iapmax,3),xbn(iapmax,3)
*,bqgs2,bmaxqgs2,bmaxnex2,bminnex2
......
......@@ -68,7 +68,7 @@ c-----------------------------------------------------------------------------
integer debug
character*7 ty
character*2 tyq
parameter(iapmax=207)
parameter(iapmax=208)
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
common /qgarr3/ rmin,emax,eev
......@@ -255,7 +255,7 @@ c-----------------------------------------------------------------------------
integer debug
character *7 ty
logical lcalc
parameter(iapmax=207)
parameter(iapmax=208)
dimension mij(40,40,4),nij(40,40,4),cs1(40,40,160)
*,evs(40,100,3,2),ixemax(40,3,2),gz0(5),gz1(3)
*,qfan0(11,14),qfan2(11,11,3),fann(14)
......@@ -2006,7 +2006,7 @@ c iat - target mass number
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
common /qgarr4/ ey0(3)
......@@ -4321,7 +4321,7 @@ c icz - hadron class
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
common /qgarr6/ pi,bm,amws
common /qgarr15/ fp(3),rq(2,3),cd(2,3),gsoft(3)
common /qgarr16/ cc(2,3),iddp(iapmax),iddt(iapmax)
......@@ -4372,7 +4372,7 @@ c icz - hadron class
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
common /qgarr6/ pi,bm,amws
common /qgarr15/ fp(3),rq(2,3),cd(2,3),gsoft(3)
common /qgarr16/ cc(2,3),iddp(iapmax),iddt(iapmax)
......@@ -4425,7 +4425,7 @@ c jj=1 - q
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
common /qgarr6/ pi,bm,amws
common /qgarr15/ fp(3),rq(2,3),cd(2,3),gsoft(3)
common /qgarr16/ cc(2,3),iddp(iapmax),iddt(iapmax)
......@@ -4569,7 +4569,7 @@ c icz - hadron class
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
common /qgarr6/ pi,bm,amws
common /qgarr15/ fp(3),rq(2,3),cd(2,3),gsoft(3)
common /qgarr16/ cc(2,3),iddp(iapmax),iddt(iapmax)
......@@ -4637,7 +4637,7 @@ c icz - hadron class
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
common /qgarr6/ pi,bm,amws
common /qgarr15/ fp(3),rq(2,3),cd(2,3),gsoft(3)
common /qgarr16/ cc(2,3),iddp(iapmax),iddt(iapmax)
......@@ -4929,7 +4929,7 @@ c hadron-hadron and hadron-nucleus cross sections calculation
c----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension gz(5),wt1(3),wt2(3)
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
......@@ -5125,7 +5125,7 @@ c (used for total and diffractive hadron-nucleus cross-sections calculation)
c----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension xb(iapmax,3),vabs(2)
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
......@@ -5177,7 +5177,7 @@ c vvxtl = 1 - exp[-sum_{j<J} chi_targ(j)]
c------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension vpac(iapmax),vtac(iapmax)
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
......@@ -5356,7 +5356,7 @@ c to it-th targ. nucleon
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207,npbmax=1000,npnmax=900,npmax=900,legmax=900)
parameter(iapmax=208,npbmax=1000,npnmax=900,npmax=900,legmax=900)
dimension xas(iapmax,3),vabs(2),vabsi(2,iapmax),wdifi(iapmax)
*,vpac(iapmax),vtac(iapmax),xpomip(npmax),xpomim(npmax)
*,vvxim(npmax),bpomim(npmax),xpompi(legmax),xpomti(legmax)
......@@ -6130,7 +6130,7 @@ c iret=1 - reject configuration
c-----------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207,npbmax=1000,npnmax=900,npmax=900
parameter(iapmax=208,npbmax=1000,npnmax=900,npmax=900
*,levmax=20,legmax=900)
dimension vpac(iapmax),vtac(iapmax)
*,vpac0(iapmax),vtac0(iapmax),vpact(iapmax),vtact(iapmax)
......@@ -10890,7 +10890,7 @@ c ncolb - number of inel.-wounded targ. nucleons
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207,npbmax=1000,npnmax=900,npmax=900
parameter(iapmax=208,npbmax=1000,npnmax=900,npmax=900
*,legmax=900,njmax=50000)
dimension wppr0(iapmax),wmtg0(iapmax),wppr1(iapmax),wmtg1(iapmax)
*,wppr2(iapmax),wmtg2(iapmax),izp(iapmax),izt(iapmax)
......@@ -11544,7 +11544,7 @@ c imod = 1 - propose values according to x^delf * (1 - sum_i x_i)^ahl
c-------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207,npbmax=1000,npnmax=900,npmax=900,legmax=900)
parameter(iapmax=208,npbmax=1000,npnmax=900,npmax=900,legmax=900)
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
common /qgarr6/ pi,bm,amws
......@@ -11688,7 +11688,7 @@ c imod = 1 - propose values according to x^delf * (1 - sum_i x_i)^ahl
c-------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207,npbmax=1000,npnmax=900,npmax=900,legmax=900)
parameter(iapmax=208,npbmax=1000,npnmax=900,npmax=900,legmax=900)
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
common /qgarr6/ pi,bm,amws
......@@ -11790,7 +11790,7 @@ c---------------------------------------------------------------------------
*ldau1(30,49),lpar1(30,50),
*qv2(30,50),zv2(30,50),qm2(30,50),iqv2(30,50),
*ldau2(30,49),lpar2(30,50)
parameter(iapmax=207,npbmax=1000,npnmax=900,npmax=900,legmax=900)
parameter(iapmax=208,npbmax=1000,npnmax=900,npmax=900,legmax=900)
common /qgarr2/ scm,wp0,wm0
common /qgarr6/ pi,bm,amws
common /qgarr8/ wwm,be(4),dc(5),deta,almpt,ptdif,ptndi
......@@ -13355,7 +13355,7 @@ c ia - number of nucleons to be considered
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension xa(iapmax,3)
common /qgarr5/ rnuc(2),wsnuc(2),wbnuc(2),anorm
*,cr1(2),cr2(2),cr3(2)
......@@ -14879,7 +14879,7 @@ c integrands for hadron-hadron and hadron-nucleus cross-sections calculation
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension gz(3),gz0(5)
common /qgarr1/ ia(2),icz,icp
common /qgarr16/ cc(2,3),iddp(iapmax),iddt(iapmax)
......@@ -14928,7 +14928,7 @@ c multifragmentation
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension sa(iapmax,3)
common /qgarr13/ nsf,iaf(iapmax)
common /qgarr43/ moniou
......@@ -14999,7 +14999,7 @@ c xa - array for spectator nucleons positions
c ns - total number of spectators
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
parameter(iapmax=207)
parameter(iapmax=208)
dimension xa(iapmax,3)
integer debug
common /qgarr1/ ia(2),icz,icp
......@@ -16209,7 +16209,7 @@ c (target diffraction is not treated explicitely and contributes to
c gdd, gqel, gcoh).
c-------------------------------------------------------------------------------
implicit double precision (a-h,o-z)
parameter(iapmax=207)
parameter(iapmax=208)
dimension wabs(28),wdd(28),wqel(28),wcoh(28)
*,wprod(28),b0(28),ai(28),xa(iapmax,3),xb(iapmax,3)
common /qgarr1/ ia(2),icz,icp
......@@ -16299,7 +16299,7 @@ c-------------------------------------------------------------------------------
subroutine qggcr(b,gabs,gdd,gqel,gcoh,xa,xb,ia)
c-------------------------------------------------------------------------------
implicit double precision (a-h,o-z)
parameter(iapmax=207)
parameter(iapmax=208)
dimension xa(iapmax,3),xb(iapmax,3),vabs(2)
gabs=1.
......
......@@ -7,7 +7,7 @@ c-----------------------------------------------------------------------
c-----------------------------------------------------------------------
c Primary initialization for QGSJET-II
c-----------------------------------------------------------------------
parameter(iapmax=207)
parameter(iapmax=208)
include 'epos.inc'
integer debug
common /qgdebug/ debug
......@@ -44,7 +44,7 @@ c-----------------------------------------------------------------------
c-----------------------------------------------------------------------
c Initialization for each type of event (for given proj, targ and egy)
c-----------------------------------------------------------------------
parameter(iapmax=207)
parameter(iapmax=208)
include 'epos.inc'
common/geom/rmproj,rmtarg,bmax,bkmx
c QGSJET-II Common
......@@ -54,7 +54,7 @@ c QGSJET-II Common
if(matarg.gt.iapmax.or.maproj.gt.iapmax)
& call utstop('Nucleus too big for QGSJET-II (Mmax=207) !&')
& call utstop('Nucleus too big for QGSJET-II (Mmax=208) !&')
call iclass(idproj,iclpro)
call iclass(idtarg,icltar)
icp=idtrafo('nxs','qgs',idproj)
......@@ -81,7 +81,7 @@ c-----------------------------------------------------------------------
c-----------------------------------------------------------------------
c call qgsjet-II to simulate interaction
c-----------------------------------------------------------------------
parameter(iapmax=207,nptmax=95000)
parameter(iapmax=208,nptmax=95000)
include 'epos.inc'
include 'epos.incems'
common/geom/rmproj,rmtarg,bmax,bkmx
......
......@@ -57,7 +57,7 @@ c-----------------------------------------------------------------------------
integer debug
character*7 ty
character*2 tyq
parameter(iapmax=207)
parameter(iapmax=208)
common /qgarr3/ rmin,emax,eev
common /qgarr6/ pi,bm,amws
common /qgarr7/ xa(iapmax,3),xb(iapmax,3),b
......@@ -225,7 +225,7 @@ c-----------------------------------------------------------------------------
integer debug
character *7 ty
logical lcalc
parameter(iapmax=207)
parameter(iapmax=208)
dimension mij(40,40,4),nij(40,40,4,2),cs1(40,40,160,2)
*,evs(40,100,3,2),ixemax(40,3,2),qfan0(51,11,6),fann(6)
*,gz0(5),gz1(3)
......@@ -1489,7 +1489,7 @@ c iat - target mass number
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
common /qgarr4/ ey0(3)
......@@ -2482,7 +2482,7 @@ c gz(i), i=1,5 - contributions of different cuts (elastic, inel., and diffr.)
c----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension gz(5),wt1(3),wt2(3)
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
......@@ -2878,7 +2878,7 @@ c it - targ. index
c------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
ctp dimension vpau0(iapmax),vtau0(iapmax),vpau(iapmax),vtau(iapmax)
dimension vpau0(iapmax),vtau0(iapmax)
*,vpac(iapmax),vtac(iapmax)
......@@ -4309,7 +4309,7 @@ c x,y - projectile nucleon coordinates
c----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension xb(iapmax,3),vabs(2)
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
......@@ -4413,7 +4413,7 @@ c to it-th targ. nucleon
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207,npbmax=1000,npnmax=1000,npmax=5000
parameter(iapmax=208,npbmax=1000,npnmax=1000,npmax=5000
*,legmax=900) !so030206
dimension xas(iapmax,3),vabs(2),vabsi(2,iapmax),wdifi(iapmax)
*,vpac(iapmax),vtac(iapmax),xpomip(npmax),xpomim(npmax)
......@@ -5016,7 +5016,7 @@ c iret=1 - reject configuration
c-----------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207,npmax=5000,levmax=20,legmax=900) !so030206
parameter(iapmax=208,npmax=5000,levmax=20,legmax=900) !so030206
dimension vpac(iapmax),vtac(iapmax)
*,vpac0(iapmax),vtac0(iapmax),vpact0(iapmax),vtact0(iapmax)
*,xpomip(npmax),xpomim(npmax),xpompi(legmax),xpomti(legmax)
......@@ -6415,7 +6415,7 @@ c ncolb - number of inel.-wounded targ. nucleons
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207,npbmax=1000,npnmax=1000,npmax=5000
parameter(iapmax=208,npbmax=1000,npnmax=1000,npmax=5000
*,legmax=900,njmax=50000) !so030206
dimension wppr0(iapmax),wmtg0(iapmax),wppr1(iapmax),wmtg1(iapmax)
*,wppr2(iapmax),wmtg2(iapmax),izp(iapmax),izt(iapmax)
......@@ -7095,7 +7095,7 @@ c imod = 1 - propose values according to x^(delf-alpd) * (1 - sum_i x_i)^ahl
c-------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207,npbmax=1000,npnmax=1000,npmax=5000
parameter(iapmax=208,npbmax=1000,npnmax=1000,npmax=5000
*,legmax=900) !so030206
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
......@@ -7237,7 +7237,7 @@ c imod = 1 - propose values according to x^(delf-alpd) * (1 - sum_i x_i)^ahl
c-------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207,npbmax=1000,npnmax=1000,npmax=5000
parameter(iapmax=208,npbmax=1000,npnmax=1000,npmax=5000
*,legmax=900) !so030206
common /qgarr1/ ia(2),icz,icp
common /qgarr2/ scm,wp0,wm0
......@@ -8741,7 +8741,7 @@ c jj=2 - target
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension xa(iapmax,3)
common /qgarr5/ rnuc(2),wsnuc(2),wbnuc(2),anorm
*,cr1(2),cr2(2),cr3(2)
......@@ -10277,7 +10277,7 @@ c qgfau - integrands for hadron-hadron and hadron-nucleus cross-sections calc.
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension gz(3),gz0(5)
common /qgarr1/ ia(2),icz,icp
common /qgarr16/ cc(2,3),iddp(iapmax),iddt(iapmax)
......@@ -10324,7 +10324,7 @@ c qgfrag - multifragmentation (search for connected nucleon clasters)
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension sa(iapmax,3)
common /qgarr13/ nsf,iaf(iapmax)
common /qgarr43/ moniou
......@@ -10394,7 +10394,7 @@ c nsf - number of secondary fragments;
c iaf(i) - mass of the i-th fragment
c-----------------------------------------------------------------------------
implicit double precision (a-h,o-z)
parameter(iapmax=207)
parameter(iapmax=208)
dimension xa(iapmax,3)
integer debug
common /qgarr1/ ia(2),icz,icp
......@@ -11512,7 +11512,7 @@ c gdd, gqel, gcoh).
c-------------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension wabs(28),wdd(28),wqel(28),wcoh(28)
*,wprod(28),b0(28),ai(28),xa(iapmax,3),xb(iapmax,3)
common /qgarr1/ ia(2),icz,icp
......@@ -11610,7 +11610,7 @@ c b - impact parameter
c-------------------------------------------------------------------------------
implicit double precision (a-h,o-z)
integer debug
parameter(iapmax=207)
parameter(iapmax=208)
dimension xa(iapmax,3),xb(iapmax,3),vabs(2)
common /qgarr43/ moniou
common /qgdebug/ debug
......
......@@ -7,7 +7,7 @@ c-----------------------------------------------------------------------
c-----------------------------------------------------------------------
c Primary initialization for QGSJET-II
c-----------------------------------------------------------------------
parameter(iapmax=207)
parameter(iapmax=208)
include 'epos.inc'
integer debug
common /qgdebug/ debug
......@@ -42,7 +42,7 @@ c-----------------------------------------------------------------------
c-----------------------------------------------------------------------
c Initialization for each type of event (for given proj, targ and egy)
c-----------------------------------------------------------------------
parameter(iapmax=207)
parameter(iapmax=208)
include 'epos.inc'
common/geom/rmproj,rmtarg,bmax,bkmx
c QGSJET-II Common
......@@ -52,7 +52,7 @@ c QGSJET-II Common
if(matarg.gt.iapmax.or.maproj.gt.iapmax)
& call utstop('Nucleus too big for QGSJET-II (Mmax=207) !&')
& call utstop('Nucleus too big for QGSJET-II (Mmax=208) !&')
call iclass(idproj,iclpro)
call iclass(idtarg,icltar)
icp=idtrafo('nxs','qgs',idproj)
......@@ -77,7 +77,7 @@ c-----------------------------------------------------------------------
c-----------------------------------------------------------------------
c call qgsjet-II to simulate interaction
c-----------------------------------------------------------------------
parameter(iapmax=207,nptmax=95000)
parameter(iapmax=208,nptmax=95000)
include 'epos.inc'
include 'epos.incems'
common/geom/rmproj,rmtarg,bmax,bkmx
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment