CCC Common FACFAC for Factorials
integer factorial_max
- parameter (factorial_max = 10000) ! max # multiplicity per event;
+ parameter (factorial_max = 80000) ! max # multiplicity per event;
CCC ! for any specific particle ID;
CCC ! also used for log(n!).
Common/FACFAC/ FACLOG(factorial_max)
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')