c3po/F77/ntuples_functions.f

C-----------------------------------------------------------------------
C     Subroutines to book and fill ntuples for Level2 variables
C-----------------------------------------------------------------------
      
C-----------------------------------------------------------------------
C     Tracker ntupla
C-----------------------------------------------------------------------

      function booktrackerntupla()
 
      include './common/common_level2.f'


      call HBNT(trkntp,'TRACKER',' ')
      
      call HBNAME(trkntp,'tracker',good2,'good2:I::[0,1]
     $     ,crc(12):I::[0,2],nclsx:I::[0,100],planex(nclsx):I::[0,6]
     $     ,xs(2,nclsx):R,sgnlxs(nclsx):R,nclsy:I::[0,100]
     $     ,planey(nclsy):I::[0,6],ys(2,nclsy):R,sgnlys(nclsy):R
     $     ,ntrk:I::[0,10],image(ntrk):I::[0,10],xm_nt(6,ntrk):R
     $     ,ym_nt(6,ntrk):R,zm_nt(6,ntrk):R,resx_nt(6,ntrk):R
     $     ,resy_nt(6,ntrk):R,al_nt(5,ntrk):R,coval(5,5,ntrk):R
     $     ,chi2_nt(ntrk):R,xgood_nt(6,ntrk):I::[0,1]
     $     ,ygood_nt(6,ntrk):I::[0,1],xv_nt(6,ntrk):R
     $     ,yv_nt(6,ntrk):R,zv_nt(6,ntrk):R,axv_nt(6,ntrk):R
     $     ,ayv_nt(6,ntrk):R,dedx_x(6,ntrk):R
     $     ,dedx_y(6,ntrk):R,bdl(ntrk):R')

      return
      end

C-----------------------------------------------------------------------
      function filltrackerntupla()

      include './common/common_level2.f'

      call HFNT(trkntp)

      return
      end
C-----------------------------------------------------------------------
      function closetrackerntupla()

      include './common/common_level2.f'
      
      call HROUT(trkntp,ICYCLE,' ')

      return
      end
C-----------------------------------------------------------------------
C-----------------------------------------------------------------------
C     Anticounter ntupla
C-----------------------------------------------------------------------

      function bookanticntupla()
 
      include './common/common_level2.f'

      call HBNT(acntp,'ANTICOUNTER',' ')
      
      call HBNAME(acntp,'antic',status,'status(2):I::[0,65535]
     $     ,hitmap(2):I::[0,65535],hitstatus(2):I::[0,65535]
     $     ,trigger(2):I::[0,65535]')

      return
      end

C-----------------------------------------------------------------------
      function fillanticntupla()

      include './common/common_level2.f'

      call HFNT(acntp)

      return
      end
C-----------------------------------------------------------------------
      function closeanticntupla()

      include './common/common_level2.f'
      
      call HROUT(acntp,ICYCLE,' ')

      return
      end
C-----------------------------------------------------------------------
C-----------------------------------------------------------------------
C     ToF ntupla
C-----------------------------------------------------------------------

      function booktofntupla()
 
      include './common/common_level2.f'

      call HBNT(tofntp,'TOF',' ')
      
      call HBNAME(tofntp,'tof',tof_j_flag,'tof_j_flag(6):I
     $     ,ntoftrk:I::[0,10],toftrkseqno(ntoftrk):I
     $     ,npmttdc(ntoftrk):I::[0,48],pmttdc(48,ntoftrk):I::[0,48]
     $     ,beta(13,ntoftrk):R,npmtadc(ntoftrk):I::[0,48]
     $     ,pmtadc(48,ntoftrk):I::[0,48]
     $     ,dedx(48,ntoftrk):R,xtofpos(3,ntoftrk):R
     $     ,ytofpos(3,ntoftrk):R,npmt:I::[0,48],pmt_id(npmt):I::[0,48]
     $     ,adc(npmt):R,tdc_tw(npmt):R')

      return
      end

C-----------------------------------------------------------------------
      function filltofntupla()

      include './common/common_level2.f'
    
      call HFNT(tofntp)

      return
      end
C-----------------------------------------------------------------------
      function closetofntupla()

      include './common/common_level2.f'
      
      call HROUT(tofntp,ICYCLE,' ')

      return
      end
C-----------------------------------------------------------------------
C-----------------------------------------------------------------------
C    Trigger ntupla
C-----------------------------------------------------------------------

      function booktriggerntupla()
 
      include './common/common_level2.f'

      call HBNT(trigntp,'TRIGGER',' ')
      
      call HBNAME(trigntp,'trigger',evcount,'evcount:I
     $     ,pmtpl(3):I,trigrate(6):I,dltime(2):I
     $     ,s4calcount(2):I,pmtcount1(24):I
     $     ,pmtcount2(24):I,patternbusy(3):I
     $     ,patterntrig(6):I,trigconf:I')

      return
      end

C-----------------------------------------------------------------------
      function filltriggerntupla()

      include './common/common_level2.f'

      call HFNT(trigntp)

      return
      end
C-----------------------------------------------------------------------
      function closetriggerntupla()

      include './common/common_level2.f'
      
      call HROUT(trigntp,ICYCLE,' ')

      return
      end
C-----------------------------------------------------------------------

C-----------------------------------------------------------------------
C    S4 ntupla
C-----------------------------------------------------------------------

      function books4ntupla()
 
      include './common/common_level2.f'

      call HBNT(s4ntp,'S4',' ')
      
      call HBNAME(s4ntp,'s4',s4calibrated,'s4calibrated:R
     $     ,s4adc:R')

      return
      end

C-----------------------------------------------------------------------
      function fills4ntupla()

      include './common/common_level2.f'

      call HFNT(s4ntp)

      return
      end
C-----------------------------------------------------------------------
      function closes4ntupla()

      include './common/common_level2.f'
      
      call HROUT(s4ntp,ICYCLE,' ')

      return
      end
C-----------------------------------------------------------------------

C-----------------------------------------------------------------------
C    ND ntupla
C-----------------------------------------------------------------------

      function bookndntupla()
 
      include './common/common_level2.f'

      call HBNT(ndntp,'ND',' ')
      
      call HBNAME(ndntp,'nd',upperback,'upperback:R
     $     ,bottomback:R,trigphysics:R')

      return
      end

C-----------------------------------------------------------------------
      function fillndntupla()

      include './common/common_level2.f'

      call HFNT(ndntp)

      return
      end
C-----------------------------------------------------------------------
      function closendntupla()

      include './common/common_level2.f'
      
      call HROUT(ndntp,ICYCLE,' ')

      return
      end
C-----------------------------------------------------------------------

C-----------------------------------------------------------------------
C    OrbitalInfo ntupla
C-----------------------------------------------------------------------

      function bookorbitalinfontupla()
 
      include './common/common_level2.f'

      call HBNT(orbntp,'ORBINFO',' ')
      
      call HBNAME(orbntp,'orbinfo',abstime,'abstime:I*8
     $     ,obt:I*8,pkt_num:I*8')

      return
      end

C-----------------------------------------------------------------------
      function fillorbitalinfontupla()

      include './common/common_level2.f'

      call HFNT(orbntp)

      return
      end
C-----------------------------------------------------------------------
      function closeorbitalinfontupla()

      include './common/common_level2.f'
      
      call HROUT(orbntp,ICYCLE,' ')

      return
      end
C-----------------------------------------------------------------------

C-----------------------------------------------------------------------
C    RunInfo and SoftInfo ntupla
C-----------------------------------------------------------------------

      function bookruninfontupla()
 
      include './common/common_level2.f'

      call HBNT(runntp,'RUNINFO',' ')
      
      call HBNAME(runntp,'glrun',id,'id:I*8
     $     ,id_reg_run:I*8,id_reg_run_l2:I*8   
     $     ,runheader_time:I*8,runtrailer_time:I*8
     $     ,ev_reg_phys_from:I,ev_reg_phys_to :I   
     $     ,ev_reg_runheader:I,ev_reg_runtraile:I 
     $     ,trk_calib_used:I,eff_wrk_schedule:I
     $     ,prh_var_trg_mode_a:I,prh_var_trg_mode_b:I
     $     ,acq_build_info:I,acq_var_info:I')

      call HBNAMC(runntp,'softinfo',trig_ver,'trig_ver:C*8
     $     ,tof_ver:C*8,ac_ver:C*8,trk_ver:C*8,cal_ver:C*8
     $     ,s4_ver:C*8,nd_ver:C*8,orb_ver:C*8')
      return
      end

C-----------------------------------------------------------------------
      function fillruninfontupla()

      include './common/common_level2.f'

      call HFNT(runntp)

      return
      end
C-----------------------------------------------------------------------
      function closeruninfontupla()

      include './common/common_level2.f'
      
      call HROUT(runntp,ICYCLE,' ')

      return
      end
C-----------------------------------------------------------------------
C-----------------------------------------------------------------------
C     Calorimeter ntupla
C-----------------------------------------------------------------------

      function bookcalontupla()
 
      include './common/common_level2.f'

      call HBNT(calontp,'CALO',' ')
      
      call HBNAME(calontp,'calo',calntrk,'calntrk:I::[0,11]
     $     ,caltrkseqno(calntrk):I::[-1,11],ncore(calntrk):I::[0,500000]
     $     ,noint(calntrk):I::[0,2000],ncyl(calntrk):I::[0,2000]
     $     ,nlast(calntrk):I::[0,2000],npre(calntrk):I::[0,2000]
     $     ,npresh(calntrk):I::[0,2000],ntr(calntrk):I::[0,2000]
     $     ,planetot(calntrk):I::[0,45],nlow(calntrk):I::[0,4225]
     $     ,qcore(calntrk):R,qcyl(calntrk):R,qlast(calntrk):R
     $     ,qpre(calntrk):R,qpresh(calntrk):R,qtr(calntrk):R
     $     ,qtrack(calntrk):R,qtrackx(calntrk):R,qtracky(calntrk):R
     $     ,dxtrack(calntrk):R,dytrack(calntrk):R,qmean(calntrk):R
     $     ,qlow(calntrk):R,dX0l(calntrk):R,tbar(2,22,calntrk):R
     $     ,good:I::[0,1],perr(4):I::[0,500],swerr(4):I::[0,65535]
     $     ,calcrc(4):I::[0,1],selftrigger:I::[0,1],nstrip:I::[0,4225]
     $     ,nx22:I::[0,97],planemax(2):I::[0,45],qtot:R,qx22:R,qmax:R
     $     ,qq(4):R,npcfit(2):I::[0,45],impx:R,impy:R,tanx:R,tany:R
     $     ,varcfit(2):R,elen:R,selen:R')

      return
      end

C-----------------------------------------------------------------------
      function fillcalontupla()

      include './common/common_level2.f'
    
      call HFNT(calontp)

      return
      end
C-----------------------------------------------------------------------
      function closecalontupla()

      include './common/common_level2.f'
      
      call HROUT(calontp,ICYCLE,' ')

      return
      end
C-----------------------------------------------------------------------
1	pam-fi	1.1	C-----------------------------------------------------------------------
2			C Subroutines to book and fill ntuples for Level2 variables
3			C-----------------------------------------------------------------------
4
5			C-----------------------------------------------------------------------
6			C Tracker ntupla
7			C-----------------------------------------------------------------------
8
9			function booktrackerntupla()
10
11			include './common/common_level2.f'
12
13
14			call HBNT(trkntp,'TRACKER',' ')
15
16			call HBNAME(trkntp,'tracker',good2,'good2:I::[0,1]
17			$ ,crc(12):I::[0,2],nclsx:I::[0,100],planex(nclsx):I::[0,6]
18			$ ,xs(2,nclsx):R,sgnlxs(nclsx):R,nclsy:I::[0,100]
19			$ ,planey(nclsy):I::[0,6],ys(2,nclsy):R,sgnlys(nclsy):R
20			$ ,ntrk:I::[0,10],image(ntrk):I::[0,10],xm_nt(6,ntrk):R
21			$ ,ym_nt(6,ntrk):R,zm_nt(6,ntrk):R,resx_nt(6,ntrk):R
22			$ ,resy_nt(6,ntrk):R,al_nt(5,ntrk):R,coval(5,5,ntrk):R
23			$ ,chi2_nt(ntrk):R,xgood_nt(6,ntrk):I::[0,1]
24			$ ,ygood_nt(6,ntrk):I::[0,1],xv_nt(6,ntrk):R
25			$ ,yv_nt(6,ntrk):R,zv_nt(6,ntrk):R,axv_nt(6,ntrk):R
26			$ ,ayv_nt(6,ntrk):R,dedx_x(6,ntrk):R
27			$ ,dedx_y(6,ntrk):R,bdl(ntrk):R')
28
29			return
30			end
31
32			C-----------------------------------------------------------------------
33			function filltrackerntupla()
34
35			include './common/common_level2.f'
36
37			call HFNT(trkntp)
38
39			return
40			end
41			C-----------------------------------------------------------------------
42			function closetrackerntupla()
43
44			include './common/common_level2.f'
45
46			call HROUT(trkntp,ICYCLE,' ')
47
48			return
49			end
50			C-----------------------------------------------------------------------
51			C-----------------------------------------------------------------------
52			C Anticounter ntupla
53			C-----------------------------------------------------------------------
54
55			function bookanticntupla()
56
57			include './common/common_level2.f'
58
59			call HBNT(acntp,'ANTICOUNTER',' ')
60
61			call HBNAME(acntp,'antic',status,'status(2):I::[0,65535]
62			$ ,hitmap(2):I::[0,65535],hitstatus(2):I::[0,65535]
63			$ ,trigger(2):I::[0,65535]')
64
65			return
66			end
67
68			C-----------------------------------------------------------------------
69			function fillanticntupla()
70
71			include './common/common_level2.f'
72
73			call HFNT(acntp)
74
75			return
76			end
77			C-----------------------------------------------------------------------
78			function closeanticntupla()
79
80			include './common/common_level2.f'
81
82			call HROUT(acntp,ICYCLE,' ')
83
84			return
85			end
86			C-----------------------------------------------------------------------
87			C-----------------------------------------------------------------------
88			C ToF ntupla
89			C-----------------------------------------------------------------------
90
91			function booktofntupla()
92
93			include './common/common_level2.f'
94
95			call HBNT(tofntp,'TOF',' ')
96
97			call HBNAME(tofntp,'tof',tof_j_flag,'tof_j_flag(6):I
98			$ ,ntoftrk:I::[0,10],toftrkseqno(ntoftrk):I
99	pam-fi	1.3	$ ,npmttdc(ntoftrk):I::[0,48],pmttdc(48,ntoftrk):I::[0,48]
100			$ ,beta(13,ntoftrk):R,npmtadc(ntoftrk):I::[0,48]
101			$ ,pmtadc(48,ntoftrk):I::[0,48]
102			$ ,dedx(48,ntoftrk):R,xtofpos(3,ntoftrk):R
103			$ ,ytofpos(3,ntoftrk):R,npmt:I::[0,48],pmt_id(npmt):I::[0,48]
104	pam-fi	1.1	$ ,adc(npmt):R,tdc_tw(npmt):R')
105
106			return
107			end
108
109			C-----------------------------------------------------------------------
110			function filltofntupla()
111
112			include './common/common_level2.f'
113
114			call HFNT(tofntp)
115
116			return
117			end
118			C-----------------------------------------------------------------------
119			function closetofntupla()
120
121			include './common/common_level2.f'
122
123			call HROUT(tofntp,ICYCLE,' ')
124
125			return
126			end
127			C-----------------------------------------------------------------------
128			C-----------------------------------------------------------------------
129			C Trigger ntupla
130			C-----------------------------------------------------------------------
131
132			function booktriggerntupla()
133
134			include './common/common_level2.f'
135
136			call HBNT(trigntp,'TRIGGER',' ')
137
138			call HBNAME(trigntp,'trigger',evcount,'evcount:I
139			$ ,pmtpl(3):I,trigrate(6):I,dltime(2):I
140			$ ,s4calcount(2):I,pmtcount1(24):I
141			$ ,pmtcount2(24):I,patternbusy(3):I
142			$ ,patterntrig(6):I,trigconf:I')
143
144			return
145			end
146
147			C-----------------------------------------------------------------------
148			function filltriggerntupla()
149
150			include './common/common_level2.f'
151
152			call HFNT(trigntp)
153
154			return
155			end
156			C-----------------------------------------------------------------------
157			function closetriggerntupla()
158
159			include './common/common_level2.f'
160
161			call HROUT(trigntp,ICYCLE,' ')
162
163			return
164			end
165			C-----------------------------------------------------------------------
166
167			C-----------------------------------------------------------------------
168			C S4 ntupla
169			C-----------------------------------------------------------------------
170
171			function books4ntupla()
172
173			include './common/common_level2.f'
174
175			call HBNT(s4ntp,'S4',' ')
176
177			call HBNAME(s4ntp,'s4',s4calibrated,'s4calibrated:R
178			$ ,s4adc:R')
179
180			return
181			end
182
183			C-----------------------------------------------------------------------
184			function fills4ntupla()
185
186			include './common/common_level2.f'
187
188			call HFNT(s4ntp)
189
190			return
191			end
192			C-----------------------------------------------------------------------
193			function closes4ntupla()
194
195			include './common/common_level2.f'
196
197			call HROUT(s4ntp,ICYCLE,' ')
198
199			return
200			end
201			C-----------------------------------------------------------------------
202
203			C-----------------------------------------------------------------------
204			C ND ntupla
205			C-----------------------------------------------------------------------
206
207			function bookndntupla()
208
209			include './common/common_level2.f'
210
211			call HBNT(ndntp,'ND',' ')
212
213			call HBNAME(ndntp,'nd',upperback,'upperback:R
214			$ ,bottomback:R,trigphysics:R')
215
216			return
217			end
218
219			C-----------------------------------------------------------------------
220			function fillndntupla()
221
222			include './common/common_level2.f'
223
224			call HFNT(ndntp)
225
226			return
227			end
228			C-----------------------------------------------------------------------
229			function closendntupla()
230
231			include './common/common_level2.f'
232
233			call HROUT(ndntp,ICYCLE,' ')
234
235			return
236			end
237			C-----------------------------------------------------------------------
238
239			C-----------------------------------------------------------------------
240			C OrbitalInfo ntupla
241			C-----------------------------------------------------------------------
242
243			function bookorbitalinfontupla()
244
245			include './common/common_level2.f'
246
247			call HBNT(orbntp,'ORBINFO',' ')
248
249			call HBNAME(orbntp,'orbinfo',abstime,'abstime:I*8
250			$ ,obt:I8,pkt_num:I8')
251
252			return
253			end
254
255			C-----------------------------------------------------------------------
256			function fillorbitalinfontupla()
257
258			include './common/common_level2.f'
259
260			call HFNT(orbntp)
261
262			return
263			end
264			C-----------------------------------------------------------------------
265			function closeorbitalinfontupla()
266
267			include './common/common_level2.f'
268
269			call HROUT(orbntp,ICYCLE,' ')
270
271			return
272			end
273			C-----------------------------------------------------------------------
274
275			C-----------------------------------------------------------------------
276			C RunInfo and SoftInfo ntupla
277			C-----------------------------------------------------------------------
278
279			function bookruninfontupla()
280
281			include './common/common_level2.f'
282
283			call HBNT(runntp,'RUNINFO',' ')
284
285			call HBNAME(runntp,'glrun',id,'id:I*8
286			$ ,id_reg_run:I8,id_reg_run_l2:I8
287			$ ,runheader_time:I8,runtrailer_time:I8
288			$ ,ev_reg_phys_from:I,ev_reg_phys_to :I
289			$ ,ev_reg_runheader:I,ev_reg_runtraile:I
290			$ ,trk_calib_used:I,eff_wrk_schedule:I
291			$ ,prh_var_trg_mode_a:I,prh_var_trg_mode_b:I
292			$ ,acq_build_info:I,acq_var_info:I')
293
294			call HBNAMC(runntp,'softinfo',trig_ver,'trig_ver:C*8
295			$ ,tof_ver:C8,ac_ver:C8,trk_ver:C8,cal_ver:C8
296			$ ,s4_ver:C8,nd_ver:C8,orb_ver:C*8')
297			return
298			end
299
300			C-----------------------------------------------------------------------
301			function fillruninfontupla()
302
303			include './common/common_level2.f'
304
305			call HFNT(runntp)
306
307			return
308			end
309			C-----------------------------------------------------------------------
310			function closeruninfontupla()
311
312			include './common/common_level2.f'
313
314			call HROUT(runntp,ICYCLE,' ')
315
316			return
317			end
318			C-----------------------------------------------------------------------
319			C-----------------------------------------------------------------------
320			C Calorimeter ntupla
321			C-----------------------------------------------------------------------
322
323			function bookcalontupla()
324
325			include './common/common_level2.f'
326
327			call HBNT(calontp,'CALO',' ')
328
329	pam-fi	1.4	call HBNAME(calontp,'calo',calntrk,'calntrk:I::[0,11]
330	pam-fi	1.2	$ ,caltrkseqno(calntrk):I::[-1,11],ncore(calntrk):I::[0,500000]
331	pam-fi	1.1	$ ,noint(calntrk):I::[0,2000],ncyl(calntrk):I::[0,2000]
332			$ ,nlast(calntrk):I::[0,2000],npre(calntrk):I::[0,2000]
333			$ ,npresh(calntrk):I::[0,2000],ntr(calntrk):I::[0,2000]
334			$ ,planetot(calntrk):I::[0,45],nlow(calntrk):I::[0,4225]
335			$ ,qcore(calntrk):R,qcyl(calntrk):R,qlast(calntrk):R
336			$ ,qpre(calntrk):R,qpresh(calntrk):R,qtr(calntrk):R
337			$ ,qtrack(calntrk):R,qtrackx(calntrk):R,qtracky(calntrk):R
338			$ ,dxtrack(calntrk):R,dytrack(calntrk):R,qmean(calntrk):R
339	pam-fi	1.4	$ ,qlow(calntrk):R,dX0l(calntrk):R,tbar(2,22,calntrk):R
340			$ ,good:I::[0,1],perr(4):I::[0,500],swerr(4):I::[0,65535]
341			$ ,calcrc(4):I::[0,1],selftrigger:I::[0,1],nstrip:I::[0,4225]
342			$ ,nx22:I::[0,97],planemax(2):I::[0,45],qtot:R,qx22:R,qmax:R
343			$ ,qq(4):R,npcfit(2):I::[0,45],impx:R,impy:R,tanx:R,tany:R
344	pam-fi	1.1	$ ,varcfit(2):R,elen:R,selen:R')
345
346			return
347			end
348
349			C-----------------------------------------------------------------------
350			function fillcalontupla()
351
352			include './common/common_level2.f'
353
354			call HFNT(calontp)
355
356			return
357			end
358			C-----------------------------------------------------------------------
359			function closecalontupla()
360
361			include './common/common_level2.f'
362
363			call HROUT(calontp,ICYCLE,' ')
364
365			return
366			end
367			C-----------------------------------------------------------------------