September 16, 2002

PAW version 2.13/08

The PAW version 2.13/08 is part of the 2002 release of CERNLIB. The source files, binaries and libraries are available from the CERN Program Library home page.

The PAW support is done by:

Olivier Couet.

The PAW related packages versions are:

Package	Version	Date
PAW	2.13/08	19/09/2002
KUIP	2.08/03	23/05/2001
HBOOK	4.27/03	12/12/2001
COMIS	1.23/01	19/06/2002
HPLOT	5.24/02	12/11/2002
HIGZ	1.28/07	22/05/2002
SIGMA	1.10/03	10/02/1994

The following subjects can be found in this page:

PAW

Improvements
Help improvements
Bugs fixed

HBOOK

Improvements
Bugs fixed

HIGZ

Improvements
Bugs fixed

HPLOT

Bugs fixed

COMIS

Bugs fixed

PAW

Improvements

New routines callable via COMIS: DBESJ0, DBESJ1, DBESY0, DBESY1, DBESI0, DBESI1, DBESK0, DBESK1, DEBSI0, DEBSI1, DEBSK0, DEBSK1.

Help improvements

In the command 1DHISTO and 2DHISTO the way to define axis titles together with the histogram title is now documented.

Bugs fixed

In case of alphanumeric labels on 1D or 2D histograms, the alphabetical sorting done by the internal query processor procedure smap_sort was wrong. With the ntuple generated by the following example ...

      COMMON /PAWC/  SPACE(50000)

      CHARACTER*40   CHFILE
      CHARACTER*4    TYPE
      COMMON /CEXAM/ TYPE

      CALL HLIMIT ( 50000 )

      CHFILE = 'test.cwn'
      CALL HROPEN ( 11, 'MYFILE', CHFILE, 'N', 1024, ISTAT )
      CALL HBNT   ( 10, 'TEST NTUPLE', ' ' )
      CALL HBNAMC ( 10, 'TYPE', TYPE, 'TYPE:C' )

      CALL FILL('BB',1)
      CALL FILL('CC',2)
      CALL FILL('EE',3)
      CALL FILL('FF',4)
      CALL FILL('GG',5)
      CALL FILL('HH',6)
      CALL FILL('IT',7)
      CALL FILL('JJ',8)
      CALL FILL('LL',9)
      CALL FILL('MM',8)
      CALL FILL('NN',7)
      CALL FILL('OO',6)
      CALL FILL('SE',5)
      CALL FILL('ZZ',4)
      CALL FILL('AA',3)

      CALL HROUT  ( 10, ICYCLE, ' ' )
      CALL HREND  ( 'MYFILE' )
      END

      SUBROUTINE FILL (C,N)
      CHARACTER*(*) C
      CHARACTER*4    TYPE
      COMMON /CEXAM/ TYPE
      TYPE = C
      DO I=1,N
         CALL HFNT(10)
      ENDDO
      END

... the entries in the SE channel were dropped into the underflow channel when one use the command:

PAW >  nt/plot 10.type

The smap_sort sorting is now replaced by the sorting done by the HBOOK routine HLABEL.

The option S in the NT/PLOT command, used the limit of the histogram created by the previous NT/PLOT command.
Example:
```
 nt/plot 10.age age<40
 null 0 100 0 100
 nt/plot 10.age    
 
```
The last NT/PLOT command was also done with the cut age<40.

The following macro produced wrong results:

opt zfl
do i=1,1000
   v/del *
   sigma a=array(200,200#1)
   sigma b=array(200,1#200)
   sigma c=b*2.
   null 0. 200. 0. 400.
   hplot/symb  b a  200 20 0.21
   hplot/symb  b c  200 24 0.21
enddo

A possible workaround was to not use opt zfl, but opt zfl1 instead.

When a vector of length 1 was created, for example as the result of a VEC/READ command, the V/PLOT command didn't recognise it as a vector. Commands such as VEC/HFILL or VEC/PLOT also returned an "unknown vector" error.
PAW on Linux is now compiled with the version 2.95 of the egcs compiler. That version runs on machine having the version 2.91 (typically RH6.x machines). With 2.95 the division of complex numbers is done with inline code whereas it was done using a function in version 2.91. Due to this difference one got undefined symbols (c_div and z_div) when a COMIS function doing a division of two complex numbers, was called in compile mode (extension .f77) from a PAW module created with egcs 2.95 on a machine having egcs 2.91.

HBOOK

Improvements

New option E in HRFILE. Used with the option Q, this option gives the possibility to create HBOOK files with up to 2**32 records. This uses of the option N of RZMAKE. This option (E) is also available via HROPEN. For example:
```
      IQUEST(10)=256000
      CALL HROPEN(1,'EXAMPLE','example.rz','NQE',LREC,ISTAT)
```

Bugs fixed

The number of entries of a non existing histogram returned by HNOENT was not 0.

The following program didn't work. It could not access the shared memory area named 'ASA' created by an other process.

      Program shm_read
      Common /PAWC/ A(100000)
      Call Hlimit(100000)
      Call Hlimap(0,'ASA') 
      End

A possible workaround was to put the shared memory name in a character variable, ie:

      Program shm_read
      Common /PAWC/ A(100000)
      Character*8 C
      Call Hlimit(100000)
      C = 'ASA' 
      Call Hlimap(0,C)
      End

When HLABEL was called in a loop, the program crashed after a while with a ZEBRA error. The following little test program demonstrated the problem.

      program test
      PARAMETER (nword_paw=50000)
      COMMON/pawc/rpaw(nword_paw)
      CHARACTER*4 clab(1)
      clab(1) = 'A...'
      CALL hlimit(nword_paw)
      DO i=1,20000
         CALL hbook1(1,'Test',100.,0.,1.,0.)
         CALL hlabel(1,1,clab,'NX')
         CALL hdelet(1)
      enddo
      end

This was due to a misplaced call to MZFORM in HLABEL.

HIGZ

Improvements

Several improvements have been done in the 3D algorithms drawing routines (Lego, surface etc ...). This do not change the functionalities but make the code more robust.

Bugs fixed

When a HIST/PLOT command (in PAW) was performed with the option S with markers or lines, if the second bin of the histogram was outside the plot frame, the first bin was clipped too even if it is inside. The following PAW macro illustrated the problem:
```
  v/cr a(5) r 1 0 1 0.5 1
  set mtyp 20
  1d 100 ' ' 5 0. 5.
  put/cont 100 a
  null 0. 5. .1 1.5
  h/pl 100 s
  h/pl 100 sp
  h/pl 100 sl
```
The clipping was done in IGHIST and was causing this problem. This code was not necessary since the clipping is correctly handled elsewhere in HIGZ for all the basic primitives.

The "work" buffers were too small in the HIGZ routine IGRAPH and the following macro didn't not work for more that 101 points.

   vec/del *
   n = [1]
   vec/cr time([n])
   vec/cr etime([n]) r [n]*0.0
   vec/cr y([n]) r
   vec/cr ey([n]) r [n]*0.7
   do i = 1, [n]
     vec/input time([i]) [i]
     vec/input y([i]) $eval(3.0+0.01*[i])
   enddo
   null -10 160 0 10
   hplot/errors time y etime ey [n] 20 0.1 3

HPLOT

Bugs fixed

The fix done is the previous release (Some rounding problem on IEEE machines) introduced problems in case of negative division number. This is now fixed.
Bug fixed in HPLAER (asymetric errors drawing). Due some changes in parameters passing to this routine, the errors along the Y axis were symetric.

COMIS

Bugs fixed

Due to some bug in the DecOSF f77 compiler, the following COMIS program:

        subroutine decbug
        z = 2.0
        y = x**z
        print *,y
        end

produced the following error message:

  PAW > call  decbug.f
  forrtl: error (65): floating invalid
   NaN
  PAW >

-PIC option removed from the list of default c compiler options on LINUX.