1 gaskelld 1.1 subroutine NtupleInit(filename)
2 implicit none
3 save
4
5 include 'hbook.inc'
6 include 'simulate.inc'
7
8 character*80 filename,directory
9 character*16 NtupleTag(80),name,title
10 integer*4 m,io,recl,bank,id,status
11 parameter(recl = 1024)
12 parameter(bank = 8000)
13 parameter(io = 29)
14 parameter(name = 'SimcNtuple')
15 parameter(title = 'SIMTUPLE')
16
17 NtupleID = defaultID
18 id = NtupleID
19 NtupleIO = io
20 NtupleName = name
21
22 gaskelld 1.1 call HCDIR(directory,'R') !CERNLIB read current directory
23 call HROPEN(io,name,filename,'N',recl,status) !CERNLIB
24 !directory set to "//TUPLE"
25 if (status.ne.0)then
26 write(6,*) 'HROPEN error: istat=',status
27 stop
28 endif
29
30 m = 0
31 m = m+1
32 NtupleTag(m) = 'hsdelta' ! 1
33 m = m+1
34 NtupleTag(m) = 'hsyptar' ! 2
35 m = m+1
36 NtupleTag(m) = 'hsxptar' ! 3
37 m = m+1
38 NtupleTag(m) = 'hsytar' ! 4
39 m = m+1
40 NtupleTag(m) = 'hsxfp' ! 5
41 m = m+1
42 NtupleTag(m) = 'hsxpfp' ! 6
43 gaskelld 1.1 m = m+1
44 NtupleTag(m) = 'hsyfp' ! 7
45 m = m+1
46 NtupleTag(m) = 'hsypfp' ! 8
47 m = m+1
48 NtupleTag(m) = 'hsdeltai' ! 9
49 m = m+1
50 NtupleTag(m) = 'hsyptari' ! 10
51 m = m+1
52 NtupleTag(m) = 'hsxptari' ! 11
53 m = m+1
54 NtupleTag(m) = 'hsytari' ! 12
55 m = m+1
56 NtupleTag(m) = 'ssdelta' ! 13
57 m = m+1
58 NtupleTag(m) = 'ssyptar' ! 14
59 m = m+1
60 NtupleTag(m) = 'ssxptar' ! 15
61 m = m+1
62 NtupleTag(m) = 'ssytar' ! 16
63 m = m+1
64 gaskelld 1.1 NtupleTag(m) = 'ssxfp' ! 17
65 m = m+1
66 NtupleTag(m) = 'ssxpfp' ! 18
67 m = m+1
68 NtupleTag(m) = 'ssyfp' ! 19
69 m = m+1
70 NtupleTag(m) = 'ssypfp' ! 20
71 m = m+1
72 NtupleTag(m) = 'ssdeltai' ! 21
73 m = m+1
74 NtupleTag(m) = 'ssyptari' ! 22
75 m = m+1
76 NtupleTag(m) = 'ssxptari' ! 23
77 m = m+1
78 NtupleTag(m) = 'ssytari' ! 24
79 m = m+1
80 NtupleTag(m) = 'q' ! 25
81 m = m+1
82 NtupleTag(m) = 'nu' ! 26
83 m = m+1
84 NtupleTag(m) = 'Q2' ! 27
85 gaskelld 1.1 m = m+1
86 NtupleTag(m) = 'W' ! 28
87 m = m+1
88 NtupleTag(m) = 'epsilon' ! 29
89 m = m+1
90 NtupleTag(m) = 'Em' ! 30
91 m = m+1
92 NtupleTag(m) = 'Pm' ! 31
93 m = m+1
94 NtupleTag(m) = 'thetapq' ! 32
95 m = m+1
96 NtupleTag(m) = 'phipq' ! 33
97
98 if (doing_pion .or. doing_kaon .or. doing_delta) then
99 m = m+1
100 NtupleTag(m) = 'missmass' ! 34
101 m = m+1
102 NtupleTag(m) = 'mmnuc' ! 35
103 m = m+1
104 NtupleTag(m) = 'phad' ! 36
105 m = m+1
106 gaskelld 1.1 NtupleTag(m) = 't' ! 37
107 m = m+1
108 NtupleTag(m) = 'pmpar' ! 38
109 m = m+1
110 NtupleTag(m) = 'pmper' ! 39
111 m = m+1
112 NtupleTag(m) = 'pmoop' ! 40
113 m = m+1
114 NtupleTag(m) = 'fry' ! 41 !+y is up.
115 m = m+1
116 NtupleTag(m) = 'radphot' ! 42
117 m = m+1
118 NtupleTag(m) = 'pfermi' ! 43
119 m = m+1
120 NtupleTag(m) = 'siglab' ! 44
121 m = m+1
122 NtupleTag(m) = 'sigcm' ! 45
123 m = m+1
124 NtupleTag(m) = 'Weight' ! 46
125 m = m+1
126 NtupleTag(m) = 'jacobian' ! 47
127 gaskelld 1.1 m = m+1
128 NtupleTag(m) = 'decdist' ! 48
129 m = m+1
130 NtupleTag(m) = 'Mhadron' ! 49
131 m = m+1
132 NtupleTag(m) = 'pdotqhat' ! 50
133 m = m+1
134 NtupleTag(m) = 'Q2i' ! 51
135 m = m+1
136 NtupleTag(m) = 'Wi' ! 52
137 m = m+1
138 NtupleTag(m) = 'ti' ! 53
139 m = m+1
140 NtupleTag(m) = 'thetapqi' ! 54
141 m = m+1
142 NtupleTag(m) = 'phipqi' ! 55
143 m = m+1
144 NtupleTag(m) = 'thetacm' ! 56
145 m = m+1
146 NtupleTag(m) = 'phicm' ! 57
147 m = m+1
148 gaskelld 1.1 NtupleTag(m) = 'phadcm' ! 58
149 if (doing_kaon) then
150 m = m+1
151 NtupleTag(m) = 'saghai' ! 59
152 m = m+1
153 NtupleTag(m) = 'factor' ! 60
154 endif
155 else if (doing_semi.or.doing_rho) then
156 m = m+1
157 NtupleTag(m) = 'missmass' ! 34 <- Wprime for semi-inclusive folks
158 m = m+1
159 NtupleTag(m) = 'ppi' ! 35
160 m = m+1
161 NtupleTag(m) = 't' ! 36
162 m = m+1
163 NtupleTag(m) = 'fry' ! 37 !+y is up.
164 m = m+1
165 NtupleTag(m) = 'radphot' ! 38
166 m = m+1
167 NtupleTag(m) = 'siglab' ! 39
168 m = m+1
169 gaskelld 1.1 NtupleTag(m) = 'sigcent' ! 40
170 m = m+1
171 NtupleTag(m) = 'Weight' ! 41
172 m = m+1
173 NtupleTag(m) = 'decdist' ! 42
174 m = m+1
175 NtupleTag(m) = 'Mhadron' ! 43
176 m = m+1
177 NtupleTag(m) = 'z' ! 44
178 m = m+1
179 NtupleTag(m) = 'zi' ! 45
180 m = m+1
181 NtupleTag(m) = 'pt2' ! 46
182 m = m+1
183 NtupleTag(m) = 'pt2i' ! 47
184 m = m+1
185 NtupleTag(m) = 'xbj' ! 48
186 m = m+1
187 NtupleTag(m) = 'xbji' ! 49
188 m = m+1
189 NtupleTag(m) = 'thqi' ! 50
190 gaskelld 1.1 m = m+1
191 NtupleTag(m) = 'sighad' ! 51
192 m = m+1
193 NtupleTag(m) = 'jacobian' ! 52
194 m = m+1
195 NtupleTag(m) = 'centjac' ! 53
196 m = m+1
197 NtupleTag(m) = 'pfermi' ! 54
198 m = m+1
199 NtupleTag(m) = 'xfermi' ! 55
200
201 if(doing_rho) then
202 m = m+1
203 NtupleTag(m) = 'Mrho' ! 56
204 m = m+1
205 NtupleTag(m) = 'Thrho' ! 57
206 endif
207
208 else if (doing_hyd_elast .or. doing_deuterium .or. doing_heavy) then
209 m = m+1
210 NtupleTag(m) = 'corrsing' ! 34
211 gaskelld 1.1 m = m+1
212 NtupleTag(m) = 'Pmx' ! 35 !for Heepcheck
213 m = m+1
214 NtupleTag(m) = 'Pmy' ! 36 !for Heepcheck
215 m = m+1
216 NtupleTag(m) = 'Pmz' ! 37 !for Heepcheck
217 m = m+1
218 NtupleTag(m) = 'PmPar' ! 38
219 m = m+1
220 NtupleTag(m) = 'PmPer' ! 39
221 m = m+1
222 NtupleTag(m) = 'PmOop' ! 40
223 m = m+1
224 NtupleTag(m) = 'fry' ! 41 !+y is up.
225 m = m+1
226 NtupleTag(m) = 'radphot' ! 42
227 m = m+1
228 NtupleTag(m) = 'sigcc' ! 43
229 m = m+1
230 NtupleTag(m) = 'Weight' ! 44
231 endif
232 gaskelld 1.1
233 ! else !used to be the if (doing_phsp) option.
234 ! m=m+1
235 ! NtupleTag(m)='gd'
236 ! m=m+1
237 ! NtupleTag(m)='gt'
238 ! m=m+1
239 ! NtupleTag(m)='gp'
240 ! m=m+1
241 ! NtupleTag(m)='gy'
242 ! m=m+1
243 ! NtupleTag(m)='rd'
244 ! m=m+1
245 ! NtupleTag(m)='rt'
246 ! m=m+1
247 ! NtupleTag(m)='rp'
248 ! m=m+1
249 ! NtupleTag(m)='ry'
250 ! m=m+1
251 ! NtupleTag(m)='w'
252 ! endif
253 gaskelld 1.1
254 NtupleSize = m
255
256 call HBOOKN(id,title,NtupleSize,name,bank,NtupleTag) !create Ntuple
257
258 call HCDIR(NtupleDirectory,'R') !record Ntuple directory
259
260 call HCDIR(directory,' ') !reset CERNLIB directory
261
262 return
263 end
264 r=r=r= A 5s?^Y?^Y? r=r= P
A @ T?r=r=
265 9o?9o?9o? 3 ?== /usr/local/cern/ = &=&=&= /" $R 2H? H?2H? H? dTg 2d?== /w/cache | q?== /u/apps | q?== /u/site |
266 L?== /u/scratch | D?== /w/work | 2d?== /u/mirror | K?== /w/mss | =I==I= |
|