root2paw/src/optrklev1.for

C------------------------------------------------
      FUNCTION OPTRKLEV1(name,trnev)
C------------------------------------------------
C     Written by Emiliano Mocchiutti
C     v. 0.2
C       
      INCLUDE 'COMMON.TXT'
C
      CHARACTER*100 name
C
      INTEGER lrec, optrklev1,trnev
      INTEGER istat
      PARAMETER (lrec=8190)
      REAL hmemor(9000000)
      integer Iquest(100)    
      COMMON /pawc/hmemor
      Common /QUEST/ Iquest
      CALL HLIMIT(9000000)
      Iquest(10) = 256000
C
      optrklev1 = 0
      trnev = 0
C
C     CHANGE here (hbname and hbnt if you want another ntuple id number
C
      call HROPEN(79,'LEVEL1',name,'QP',4096,istat)
      if(istat.ne.0) goto 19      
      call HRIN(21,9999,20)       
      call access_level1
      call HNOENT(41,trnev)

 19   continue

C
c
      RETURN
      END


      SUBROUTINE access_level1
*
      include 'COMMON.TXT'

      CALL HBNAME(41,' ',0,'$CLEAR')
*
      CALL HBNAME(41,'EVENT',GOOD1,'$SET')
*
      CALL HBNAME(41,'CPU',PKT_TYPE1,'$SET')
*
      CALL HBNAME(41,'CLUSTER',NCLSTR1,'$SET')
*
      CALL HBNAME(41,'CNOISE',CNEV,'$SET')
*

*
      return
      END
1	C------------------------------------------------
2	FUNCTION OPTRKLEV1(name,trnev)
3	C------------------------------------------------
4	C Written by Emiliano Mocchiutti
5	C v. 0.2
6	C
7	INCLUDE 'COMMON.TXT'
8	C
9	CHARACTER*100 name
10	C
11	INTEGER lrec, optrklev1,trnev
12	INTEGER istat
13	PARAMETER (lrec=8190)
14	REAL hmemor(9000000)
15	integer Iquest(100)
16	COMMON /pawc/hmemor
17	Common /QUEST/ Iquest
18	CALL HLIMIT(9000000)
19	Iquest(10) = 256000
20	C
21	optrklev1 = 0
22	trnev = 0
23	C
24	C CHANGE here (hbname and hbnt if you want another ntuple id number
25	C
26	call HROPEN(79,'LEVEL1',name,'QP',4096,istat)
27	if(istat.ne.0) goto 19
28	call HRIN(21,9999,20)
29	call access_level1
30	call HNOENT(41,trnev)
31
32	19 continue
33
34	C
35	c
36	RETURN
37	END
38
39
40	SUBROUTINE access_level1
41	*
42	include 'COMMON.TXT'
43
44	CALL HBNAME(41,' ',0,'$CLEAR')
45	*
46	CALL HBNAME(41,'EVENT',GOOD1,'$SET')
47	*
48	CALL HBNAME(41,'CPU',PKT_TYPE1,'$SET')
49	*
50	CALL HBNAME(41,'CLUSTER',NCLSTR1,'$SET')
51	*
52	CALL HBNAME(41,'CNOISE',CNEV,'$SET')
53	*
54
55	*
56	return
57	END