+#ifdef __APPLE__
+#define STOP CALL EXIT !
+#define stop CALL EXIT !
+#endif
+ subroutine ah
+ STOP
+ end
C Program Mult_Gen
SUBROUTINE multgen
implicit none
real*4 MultFac_xfunc_save(nmax_integ)
CCC Open I/O Files:
- open(unit=4,type='old',access='sequential',name='mult_gen.in')
+ open(unit=4,status='old',access='sequential',file='mult_gen.in')
C--> Commented for AliRoot direct COMMON block access
C open(unit=7,type='new',access='sequential',name='mult_gen.out')
C open(unit=8,type='new',access='sequential',name='mult_gen.log')
integer reac_plane_cntrl,iphi
integer n_integ_pts
+ external ran
+
do i = 1,factorial_max
do j = 1,4
pout(pid,j,i) = 0.0
pz = (y/abs(y))*pzmag
end if
pmag = sqrt(pt*pt + pz*pz)
- if(pt.ne.0.0) then
+ if( (pt.ne.0.0) .and. (pmag .ne.pz) .and. (pmag.ne.(-pz)) ) then
+
pseudorapidity = 0.5*log((pmag+pz)/(pmag-pz))
- else if (pt.eq.0.0) then
+ else
+CCC if (pt.eq.0.0) then
pseudorapidity = 86.0
C--> write(8,10)
10 Format(10x,'Function pseudorapidity called with pt=0')