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(10,48):I::[0,48]
     $     ,beta(10,13):R,npmtadc(ntoftrk):I::[0,48]
     $     ,pmtadc(10,48):I::[0,48]
     $     ,dedx(10,48):R,xtofpos(10,3):R
     $     ,ytofpos(10,3):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,51]
     $     ,caltrkseqno(calntrk):I::[0,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,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	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	$ ,npmttdc(ntoftrk):I::[0,48],pmttdc(10,48):I::[0,48]
100	$ ,beta(10,13):R,npmtadc(ntoftrk):I::[0,48]
101	$ ,pmtadc(10,48):I::[0,48]
102	$ ,dedx(10,48):R,xtofpos(10,3):R
103	$ ,ytofpos(10,3):R,npmt:I::[0,48],pmt_id(npmt):I::[0,48]
104	$ ,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	call HBNAME(calontp,'calo',calntrk,'calntrk:I::[0,51]
330	$ ,caltrkseqno(calntrk):I::[0,11],ncore(calntrk):I::[0,500000]
331	$ ,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	$ ,qlow(calntrk):R,dX0l(calntrk):R,good:I::[0,1],perr(4)
340	$ :I::[0,500],swerr(4):I::[0,65535],calcrc(4):I::[0,1]
341	$ ,selftrigger:I::[0,1],nstrip:I::[0,4225],nx22:I::[0,97]
342	$ ,planemax(2):I::[0,45],qtot:R,qx22:R,qmax:R,qq(4):R
343	$ ,npcfit(2):I::[0,45],impx:R,impy:R,tanx:R,tany:R
344	$ ,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-----------------------------------------------------------------------