/[PAMELA software]/DarthVader/ToFLevel2/src/toftrk.for
ViewVC logotype

Annotation of /DarthVader/ToFLevel2/src/toftrk.for

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.6 - (hide annotations) (download)
Thu Nov 9 17:05:49 2006 UTC (18 years ago) by mocchiut
Branch: MAIN
CVS Tags: v2r01
Changes since 1.5: +2 -1 lines 
Added unpackError in all classes, ToF .rz bug fixed, install_GL_PARAM time screw bug fixed
1 mocchiut 1.1 *****************************************************************************
2     INTEGER FUNCTION TOFTRK()
3 mocchiut 1.5
4     C****************************************************************************
5     C 31-08-06 WM
6     C Changed to use DOTRACK2
7     C Beta calculation: now the flightpath (instead of cos(theta)) is used
8     C Beta calculation: all 4 TDV measurements must be < 4095 (in the old
9     C routine it was (t1+t2)<8000
10     C
11     C****************************************************************************
12 mocchiut 1.1 IMPLICIT NONE
13     C
14     include 'input_tof.txt'
15     include 'output_tof.txt'
16     include 'tofcomm.txt'
17     C
18    
19     c =======================================
20     c variables for tracking routine
21     c =======================================
22     integer NPOINT_MAX
23     parameter(NPOINT_MAX=100)
24    
25     c define TOF Z-coordinates
26     integer NPTOF
27     parameter (NPTOF=6)
28     DOUBLE PRECISION ZTOF(NPTOF)
29     DATA ZTOF/53.74,53.04,23.94,23.44,-23.49,-24.34/ !Sergio 9.05.2006
30    
31     integer itof
32    
33     DOUBLE PRECISION al_p(5),
34 mocchiut 1.5 & xout(NPOINT_MAX),yout(NPOINT_MAX),zin(NPTOF),
35     & THXOUT(NPOINT_MAX),THYOUT(NPOINT_MAX),TLOUT(NPOINT_MAX)
36    
37 mocchiut 1.1
38     INTEGER IFAIL
39     REAL dx,dy,dr,ds
40 mocchiut 1.5 REAL t1,t2,t3,t4
41 mocchiut 1.1 REAL yhelp,xdummy,xkorr0,xhelp,xhelp1,xhelp2
42     REAL c1,c2,sw,sxw,w_i
43 mocchiut 1.5 REAL dist,dl,F
44     INTEGER icount,ievent
45 mocchiut 1.1
46     REAL beta_mean
47 mocchiut 1.5
48 mocchiut 1.2 INTEGER j
49 mocchiut 1.1
50 mocchiut 1.6 c REAL theta12,theta13,theta23
51     REAL theta13
52 mocchiut 1.1 C-- DATA ZTOF/53.74,53.04,23.94,23.44,-23.49,-24.34/ !Sergio 9.05.2006
53     REAL tofarm12
54     PARAMETER (tofarm12 = 29.70) ! from 53.39 to 23.69
55     REAL tofarm23
56     PARAMETER (tofarm23 = 47.61) ! from 23.69 to -23.92
57     REAL tofarm13
58     PARAMETER (tofarm13 = 77.31) ! from 53.39 to -23.92
59    
60    
61    
62    
63     INTEGER ihelp
64     REAL xkorr
65    
66     REAL yl,yh,xl,xh
67     C
68     REAL hmemor(9000000)
69     INTEGER Iquest(100)
70     C
71 mocchiut 1.5 DATA ievent / 0 /
72    
73 mocchiut 1.1 COMMON / pawcd / hmemor
74     save / pawcd /
75     C
76     Common / QUESTd / Iquest
77     save / questd /
78     C
79     C Begin !
80     C
81     TOFTRK = 0
82    
83     *******************************************************************
84    
85 mocchiut 1.5 ievent = ievent +1
86    
87    
88 mocchiut 1.1 offset = 1
89     slope = 2
90     left = 1
91     right = 2
92     none_ev = 0
93     none_find = 0
94     tdc_ev = 1
95     adc_ev = 1
96     itdc = 1
97     iadc = 2
98    
99     do i=1,13
100     beta_a(i) = 100.
101     enddo
102    
103 mocchiut 1.2 do i=1,4
104     do j=1,12
105     adc_c(i,j) = 1000.
106     enddo
107     enddo
108    
109     do i=1,12
110     do j=1,4
111     tofmask(j,i) = 0
112     enddo
113     enddo
114    
115 mocchiut 1.1 C------ read tracking routine
116     * igoodevent = igoodevent+1
117     * assigned input parameters for track routine
118     * 1) Z-coordinates where the trajectory is evaluated
119     do itof=1,NPTOF
120     ZIN(itof) = ZTOF(itof)
121     enddo
122     * 2) track status vector
123     C COPY THE ALFA VECTOR FROM AL_PP TO AL_P FOR THE TRACK "T"
124     do i=1,5
125     AL_P(i) = al_pp(i)
126     enddo
127 mocchiut 1.5
128     c write(*,*) AL_P
129    
130 mocchiut 1.1 if (al_p(5).eq.0.) THEN
131 mocchiut 1.5 PRINT *,' TOF - WARNING F77: track with R = 0, discarded'
132 mocchiut 1.1 GOTO 969
133     ENDIF
134     * -------- *** tracking routine *** --------
135     IFAIL = 0
136 mocchiut 1.5 C call DOTRACK(NPTOF,ZIN,XOUT,YOUT,AL_P,IFAIL)
137     call DOTRACK2(NPTOF,ZIN,XOUT,YOUT,THXOUT,THYOUT,TLOUT,AL_P,IFAIL)
138    
139    
140     C write(*,*) (TLOUT(i),i=1,6)
141    
142 mocchiut 1.1 if(IFAIL.ne.0)then
143 mocchiut 1.5 print *,' TOF - WARNING F77: tracking failed '
144 mocchiut 1.1 goto 969
145     endif
146     * ------------------------------------------
147    
148     969 continue
149    
150     C-------------------------------------------------------------
151     C------- check which paddle penetrated the track -----------
152     C-------------------------------------------------------------
153     c middle y (or x) position of the upper and middle ToF-Paddle
154     c DATA tof11_x/ -17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85/
155     c DATA tof12_y/ -13.75,-8.25,-2.75,2.75,8.25,13.75/
156     c DATA tof21_y/ -3.75,3.75/
157     c DATA tof22_x/ -4.5,4.5/
158     c DATA tof31_x/ -6.0,0.,6.0/
159     c DATA tof32_y/ -5.0,0.0,5.0/
160     c
161     c S11 8 paddles 33.0 x 5.1 cm
162     c S12 6 paddles 40.8 x 5.5 cm
163     c S21 2 paddles 18.0 x 7.5 cm
164     c S22 2 paddles 15.0 x 9.0 cm
165     c S31 3 paddles 15.0 x 6.0 cm
166     c S32 3 paddles 18.0 x 5.0 cm
167    
168 mocchiut 1.5 c write(*,*) xout(1),xout(2),xout(3),xout(4),xout(5),xout(6)
169     c write(*,*) yout(1),yout(2),yout(3),yout(4),yout(5),yout(6)
170 mocchiut 1.1
171     C-------------- S11 --------------------------------------
172    
173     tof11_i = none_find
174    
175     yl = -33.0/2.
176     yh = 33.0/2.
177    
178     if ((yout(1).gt.yl).and.(yout(1).lt.yh)) then
179     do i=1,8
180     xl = tof11_x(i) - 5.1/2.
181     xh = tof11_x(i) + 5.1/2.
182     if ((xout(1).gt.xl).and.(xout(1).le.xh)) then
183     tof11_i=i
184     endif
185     enddo
186     endif
187    
188     C-------------- S12 --------------------------------------
189    
190     tof12_i = none_find
191    
192     xl = -40.8/2.
193     xh = 40.8/2.
194    
195     if ((xout(2).gt.xl).and.(xout(2).lt.xh)) then
196     do i=1,6
197     yl = tof12_y(i) - 5.5/2.
198     yh = tof12_y(i) + 5.5/2.
199     if ((yout(2).gt.yl).and.(yout(2).le.yh)) then
200     tof12_i=i
201     endif
202     enddo
203     endif
204    
205     C-------------- S21 --------------------------------------
206    
207     tof21_i = none_find
208    
209     xl = -18./2.
210     xh = 18./2.
211    
212     if ((xout(3).gt.xl).and.(xout(3).lt.xh)) then
213     do i=1,2
214     yl = tof21_y(i) - 7.5/2.
215     yh = tof21_y(i) + 7.5/2.
216     if ((yout(3).gt.yl).and.(yout(3).le.yh)) then
217     tof21_i=i
218     endif
219     enddo
220     endif
221    
222     C-------------- S22 --------------------------------------
223    
224     tof22_i = none_find
225    
226     yl = -15./2.
227     yh = 15./2.
228    
229     if ((yout(4).gt.yl).and.(yout(4).lt.yh)) then
230     do i=1,2
231     xl = tof22_x(i) - 9.0/2.
232     xh = tof22_x(i) + 9.0/2.
233     if ((xout(4).gt.xl).and.(xout(4).le.xh)) then
234     tof22_i=i
235     endif
236     enddo
237     endif
238    
239     C-------------- S31 --------------------------------------
240    
241     tof31_i = none_find
242    
243     yl = -15.0/2.
244     yh = 15.0/2.
245    
246     if ((yout(5).gt.yl).and.(yout(5).lt.yh)) then
247     do i=1,3
248     xl = tof31_x(i) - 6.0/2.
249     xh = tof31_x(i) + 6.0/2.
250     if ((xout(5).gt.xl).and.(xout(5).le.xh)) then
251     tof31_i=i
252     endif
253     enddo
254     endif
255    
256     C-------------- S32 --------------------------------------
257    
258     tof32_i = none_find
259    
260     xl = -18.0/2.
261     xh = 18.0/2.
262    
263     if ((xout(6).gt.xl).and.(xout(6).lt.xh)) then
264     do i=1,3
265     yl = tof32_y(i) - 5.0/2.
266     yh = tof32_y(i) + 5.0/2.
267     if ((yout(6).gt.yl).and.(yout(6).le.yh)) then
268     tof32_i=i
269     endif
270     enddo
271     endif
272    
273 mocchiut 1.5
274     C write(*,*) tof11_i,tof12_i,tof21_i,tof22_i,tof31_i,tof32_i
275    
276 mocchiut 1.2 C------------------------------------------------------------------
277     C--- calculate track position in paddle using timing difference
278     C------------------------------------------------------------------
279    
280     do i=1,3
281     xtofpos(i)=100.
282     ytofpos(i)=100.
283     enddo
284     C-----------------------------S1 --------------------------------
285    
286     IF (tof11_i.GT.none_find) THEN
287     ytofpos(1) = ((tof11(1,tof11_i,itdc)-tof11(2,tof11_i,itdc))/2.
288     + -y_coor_lin11(tof11_i,offset))/y_coor_lin11(tof11_i,slope)
289     endif
290    
291     IF (tof12_i.GT.none_find) THEN
292     xtofpos(1) = ((tof12(1,tof12_i,itdc)-tof12(2,tof12_i,itdc))/2.
293     + -x_coor_lin12(tof12_i,offset))/x_coor_lin12(tof12_i,slope)
294     endif
295    
296    
297     C-----------------------------S2 --------------------------------
298    
299     IF (tof21_i.GT.none_find) THEN
300     xtofpos(2) = ((tof21(1,tof21_i,itdc)-tof21(2,tof21_i,itdc))/2.
301     + -x_coor_lin21(tof21_i,offset))/x_coor_lin21(tof21_i,slope)
302     endif
303    
304     IF (tof22_i.GT.none_find) THEN
305     ytofpos(2) = ((tof22(1,tof22_i,itdc)-tof22(2,tof22_i,itdc))/2.
306     + -y_coor_lin22(tof22_i,offset))/y_coor_lin22(tof22_i,slope)
307     endif
308    
309    
310     C-----------------------------S3 --------------------------------
311    
312     IF (tof31_i.GT.none_find) THEN
313     ytofpos(3) = ((tof31(1,tof31_i,itdc)-tof31(2,tof31_i,itdc))/2.
314     + -y_coor_lin31(tof31_i,offset))/y_coor_lin31(tof31_i,slope)
315     endif
316    
317     IF (tof32_i.GT.none_find) THEN
318     xtofpos(3) = ((tof32(1,tof32_i,itdc)-tof32(2,tof32_i,itdc))/2.
319     + -x_coor_lin32(tof32_i,offset))/x_coor_lin32(tof32_i,slope)
320     endif
321    
322    
323     do i=1,3
324     if (abs(xtofpos(i)).gt.100.) then
325     xtofpos(i)=101.
326     endif
327     if (abs(ytofpos(i)).gt.100.) then
328     ytofpos(i)=101.
329     endif
330     enddo
331    
332 mocchiut 1.1 C----------------------------------------------------------------------
333     C--------------------Corrections on ADC-data -------------------------
334     C---------------------zenith angle theta ---------------------------
335     C----------------------------------------------------------------------
336    
337    
338     dx=0.
339     dy=0.
340     dr=0.
341     theta13 = 0.
342    
343     if (xout(1).lt.100.) then
344     dx = xout(1)-xout(6)
345     dy = yout(1)-yout(6)
346     dr = sqrt(dx*dx+dy*dy)
347     theta13 = atan(dr/tofarm13)
348     endif
349    
350    
351     C----------------------------------------------------------------------
352     C------------------angle and ADC(x) correction
353     C----------------------------------------------------------------------
354     C-----------------------------S1 --------------------------------
355    
356     yhelp=yout(1)
357    
358     IF (tof11_i.GT.none_find.AND.yhelp.lt.100) THEN
359    
360     i = tof11_i
361     xdummy=tof11(left,i,iadc)
362     tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta13)
363 mocchiut 1.5 if (tof11(left,i,iadc).lt.1000) then
364 mocchiut 1.1 xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))
365     xkorr0=adcx11(left,i,1)
366     adc_c(ch11a(i),hb11a(i))=tof11(left,i,iadc)/xkorr
367     endif
368    
369     tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta13)
370 mocchiut 1.5 if (tof11(right,i,iadc).lt.1000) then
371 mocchiut 1.1 xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))
372     xkorr0=adcx11(right,i,1)
373     adc_c(ch11b(i),hb11b(i))=tof11(right,i,iadc)/xkorr
374     endif
375     ENDIF
376    
377    
378     xhelp=xout(2)
379     IF (tof12_i.GT.none_find.AND.xhelp.lt.100) THEN
380    
381     i = tof12_i
382     tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta13)
383 mocchiut 1.5 if (tof12(left,i,iadc).lt.1000) then
384 mocchiut 1.1 xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))
385     xkorr0=adcx12(left,i,1)
386     adc_c(ch12a(i),hb12a(i))=tof12(left,i,iadc)/xkorr
387     endif
388    
389     tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta13)
390 mocchiut 1.5 if (tof12(right,i,iadc).lt.1000) then
391 mocchiut 1.1 xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))
392     xkorr0=adcx12(right,i,1)
393     adc_c(ch12b(i),hb12b(i))=tof12(right,i,iadc)/xkorr
394     endif
395     ENDIF
396    
397     C-----------------------------S2 --------------------------------
398    
399     xhelp=xout(3)
400     IF (tof21_i.GT.none_find.AND.xhelp.lt.100) THEN
401    
402     i = tof21_i
403     tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta13)
404 mocchiut 1.5 if (tof21(left,i,iadc).lt.1000) then
405 mocchiut 1.1 xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))
406     xkorr0=adcx21(left,i,1)
407     adc_c(ch21a(i),hb21a(i))=tof21(left,i,iadc)/xkorr
408     endif
409    
410     tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta13)
411 mocchiut 1.5 if (tof21(right,i,iadc).lt.1000) then
412 mocchiut 1.1 xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))
413     xkorr0=adcx21(right,i,1)
414     adc_c(ch21b(i),hb21b(i))=tof21(right,i,iadc)/xkorr
415     endif
416     ENDIF
417    
418     yhelp=yout(4)
419     IF (tof22_i.GT.none_find.AND.yhelp.lt.100) THEN
420    
421     i = tof22_i
422     tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta13)
423 mocchiut 1.5 if (tof22(left,i,iadc).lt.1000) then
424 mocchiut 1.1 xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))
425     xkorr0=adcx22(left,i,1)
426     adc_c(ch22a(i),hb22a(i))=tof22(left,i,iadc)/xkorr
427     endif
428    
429     tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta13)
430 mocchiut 1.5 if (tof22(right,i,iadc).lt.1000) then
431 mocchiut 1.1 xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))
432     xkorr0=adcx22(right,i,1)
433     adc_c(ch22b(i),hb22b(i))=tof22(right,i,iadc)/xkorr
434     endif
435     ENDIF
436    
437     C-----------------------------S3 --------------------------------
438    
439     yhelp=yout(5)
440     IF (tof31_i.GT.none_find.AND.yhelp.lt.100) THEN
441    
442     i = tof31_i
443     tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta13)
444 mocchiut 1.5 if (tof31(left,i,iadc).lt.1000) then
445 mocchiut 1.1 xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))
446     xkorr0=adcx31(left,i,1)
447     adc_c(ch31a(i),hb31a(i))=tof31(left,i,iadc)/xkorr
448     endif
449    
450     tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta13)
451 mocchiut 1.5 if (tof31(right,i,iadc).lt.1000) then
452 mocchiut 1.1 xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))
453     xkorr0=adcx31(right,i,1)
454     adc_c(ch31b(i),hb31b(i))=tof31(right,i,iadc)/xkorr
455     endif
456     ENDIF
457    
458     xhelp=xout(6)
459     IF (tof32_i.GT.none_find.AND.xhelp.lt.100) THEN
460    
461     i = tof32_i
462     tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta13)
463 mocchiut 1.5 if (tof32(left,i,iadc).lt.1000) then
464 mocchiut 1.1 xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))
465     xkorr0=adcx32(left,i,1)
466     adc_c(ch32a(i),hb32a(i))=tof32(left,i,iadc)/xkorr
467     endif
468    
469     tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta13)
470 mocchiut 1.5 if (tof32(right,i,iadc).lt.1000) then
471 mocchiut 1.1 xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))
472     xkorr0=adcx32(right,i,1)
473     adc_c(ch32b(i),hb32b(i))=tof32(right,i,iadc)/xkorr
474     endif
475     ENDIF
476    
477     C-----------------------------------------------------------------------
478     C----------------------calculate Beta ------------------------------
479     C-----------------------------------------------------------------------
480     C-------------------difference of sums ---------------------------
481     C
482     C DS = (t1+t2) - t3+t4)
483     C DS = c1 + c2/beta*cos(theta)
484     C c2 = 2d/c gives c2 = 2d/(c*TDCresolution) TDC=50ps/channel
485     C => c2 = ca.60 for 0.45 m c2 = ca.109 for 0.81 m
486     C since TDC resolution varies slightly c2 has to be calibrated
487 mocchiut 1.5 C instead of cos(theta) use factor F:
488     C F = pathlength/d
489     C => beta = c2*F/(DS-c1))
490    
491     dist = ZTOF(1) - ZTOF(5)
492     dl = 0.
493     DO I=1,5
494     dl = dl + TLOUT(i)
495     ENDDO
496     F = dl/dist
497 mocchiut 1.1
498     C S11 - S31
499     IF (tof11_i.GT.none_find.AND.tof31_i.GT.none_find) THEN
500 mocchiut 1.5 t1 = tof11(1,tof11_i,itdc)
501     t2 = tof11(2,tof11_i,itdc)
502     t3 = tof31(1,tof31_i,itdc)
503     t4 = tof31(2,tof31_i,itdc)
504     IF ((t1.lt.4095).and.(t2.lt.4095).and.
505     & (t3.lt.4095).and.(t4.lt.4095)) THEN
506 mocchiut 1.1 xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
507     xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
508     ds = xhelp1-xhelp2
509     ihelp=(tof11_i-1)*3+tof31_i
510     c1 = k_S11S31(1,ihelp)
511     c2 = k_S11S31(2,ihelp)
512 mocchiut 1.5 beta_a(1) = c2*F/(ds-c1)
513     C write(*,*) 'S11-S31 ',xhelp1,xhelp2, beta_a(1)
514 mocchiut 1.2 C------- ToF Mask - S11 - S31
515    
516     tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
517     $ tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
518     tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
519     $ tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
520    
521     tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
522     $ tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
523     tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
524     $ tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
525    
526 mocchiut 1.5 ENDIF
527     ENDIF
528 mocchiut 1.2
529 mocchiut 1.5 dist = ZTOF(1) - ZTOF(6)
530     dl = 0.
531     DO I=1,6
532     dl = dl + TLOUT(i)
533     ENDDO
534     F = dl/dist
535 mocchiut 1.1
536     C S11 - S32
537     IF (tof11_i.GT.none_find.AND.tof32_i.GT.none_find) THEN
538 mocchiut 1.5 t1 = tof11(1,tof11_i,itdc)
539     t2 = tof11(2,tof11_i,itdc)
540     t3 = tof32(1,tof32_i,itdc)
541     t4 = tof32(2,tof32_i,itdc)
542     IF ((t1.lt.4095).and.(t2.lt.4095).and.
543     & (t3.lt.4095).and.(t4.lt.4095)) THEN
544 mocchiut 1.1 xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
545     xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
546     ds = xhelp1-xhelp2
547     ihelp=(tof11_i-1)*3+tof32_i
548     c1 = k_S11S32(1,ihelp)
549     c2 = k_S11S32(2,ihelp)
550 mocchiut 1.5 beta_a(2) = c2*F/(ds-c1)
551     C write(*,*) 'S11-S32 ',xhelp1,xhelp2, beta_a(2)
552 mocchiut 1.2
553     C------- ToF Mask - S11 - S32
554    
555     tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
556     $ tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
557     tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
558     $ tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
559    
560     tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
561     $ tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
562     tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
563     $ tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
564    
565     C-------
566    
567 mocchiut 1.1 ENDIF
568 mocchiut 1.5 ENDIF
569 mocchiut 1.1
570     C S12 - S31
571 mocchiut 1.5 dist = ZTOF(2) - ZTOF(5)
572     dl = 0.
573     DO I=2,5
574     dl = dl + TLOUT(i)
575     ENDDO
576     F = dl/dist
577    
578 mocchiut 1.1 IF (tof12_i.GT.none_find.AND.tof31_i.GT.none_find) THEN
579 mocchiut 1.5 t1 = tof12(1,tof12_i,itdc)
580     t2 = tof12(2,tof12_i,itdc)
581     t3 = tof31(1,tof31_i,itdc)
582     t4 = tof31(2,tof31_i,itdc)
583     IF ((t1.lt.4095).and.(t2.lt.4095).and.
584     & (t3.lt.4095).and.(t4.lt.4095)) THEN
585 mocchiut 1.1 xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
586     xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
587     ds = xhelp1-xhelp2
588     ihelp=(tof12_i-1)*3+tof31_i
589     c1 = k_S12S31(1,ihelp)
590     c2 = k_S12S31(2,ihelp)
591 mocchiut 1.5 beta_a(3) = c2*F/(ds-c1)
592     C write(*,*) 'S12-S31 ',xhelp1,xhelp2, beta_a(3)
593 mocchiut 1.2
594     C------- ToF Mask - S12 - S31
595    
596     tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
597     $ tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
598     tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
599     $ tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
600    
601     tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
602     $ tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
603     tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
604     $ tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
605    
606     C-------
607    
608 mocchiut 1.1 ENDIF
609 mocchiut 1.5 ENDIF
610 mocchiut 1.1
611     C S12 - S32
612 mocchiut 1.5
613     dist = ZTOF(2) - ZTOF(6)
614     dl = 0.
615     DO I=2,6
616     dl = dl + TLOUT(i)
617     ENDDO
618     F = dl/dist
619    
620 mocchiut 1.1 IF (tof12_i.GT.none_find.AND.tof32_i.GT.none_find) THEN
621 mocchiut 1.5 t1 = tof12(1,tof12_i,itdc)
622     t2 = tof12(2,tof12_i,itdc)
623     t3 = tof32(1,tof32_i,itdc)
624     t4 = tof32(2,tof32_i,itdc)
625     IF ((t1.lt.4095).and.(t2.lt.4095).and.
626     & (t3.lt.4095).and.(t4.lt.4095)) THEN
627 mocchiut 1.1 xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
628     xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
629     ds = xhelp1-xhelp2
630     ihelp=(tof12_i-1)*3+tof32_i
631     c1 = k_S12S32(1,ihelp)
632     c2 = k_S12S32(2,ihelp)
633 mocchiut 1.5 beta_a(4) = c2*F/(ds-c1)
634     C write(*,*) 'S12-S32 ',xhelp1,xhelp2, beta_a(4)
635 mocchiut 1.2
636     C------- ToF Mask - S12 - S32
637    
638     tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
639     $ tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
640     tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
641     $ tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
642    
643     tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
644     $ tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
645     tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
646     $ tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
647    
648     C-------
649    
650 mocchiut 1.1 ENDIF
651 mocchiut 1.5 ENDIF
652 mocchiut 1.1
653     C S21 - S31
654 mocchiut 1.5
655     dist = ZTOF(3) - ZTOF(5)
656     dl = 0.
657     DO I=3,5
658     dl = dl + TLOUT(i)
659     ENDDO
660     F = dl/dist
661    
662 mocchiut 1.1 IF (tof21_i.GT.none_find.AND.tof31_i.GT.none_find) THEN
663 mocchiut 1.5 t1 = tof21(1,tof21_i,itdc)
664     t2 = tof21(2,tof21_i,itdc)
665     t3 = tof31(1,tof31_i,itdc)
666     t4 = tof31(2,tof31_i,itdc)
667     IF ((t1.lt.4095).and.(t2.lt.4095).and.
668     & (t3.lt.4095).and.(t4.lt.4095)) THEN
669 mocchiut 1.1 xhelp1 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
670     xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
671     ds = xhelp1-xhelp2
672     ihelp=(tof21_i-1)*3+tof31_i
673     c1 = k_S21S31(1,ihelp)
674     c2 = k_S21S31(2,ihelp)
675 mocchiut 1.5 beta_a(5) = c2*F/(ds-c1)
676 mocchiut 1.2
677     C------- ToF Mask - S21 - S31
678    
679     tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
680     $ tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
681     tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
682     $ tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
683    
684     tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
685     $ tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
686     tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
687     $ tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
688    
689     C-------
690    
691 mocchiut 1.1 ENDIF
692 mocchiut 1.5 ENDIF
693 mocchiut 1.1
694     C S21 - S32
695 mocchiut 1.5
696     dist = ZTOF(3) - ZTOF(6)
697     dl = 0.
698     DO I=3,6
699     dl = dl + TLOUT(i)
700     ENDDO
701     F = dl/dist
702    
703 mocchiut 1.1 IF (tof21_i.GT.none_find.AND.tof32_i.GT.none_find) THEN
704 mocchiut 1.5 t1 = tof21(1,tof21_i,itdc)
705     t2 = tof21(2,tof21_i,itdc)
706     t3 = tof32(1,tof32_i,itdc)
707     t4 = tof32(2,tof32_i,itdc)
708     IF ((t1.lt.4095).and.(t2.lt.4095).and.
709     & (t3.lt.4095).and.(t4.lt.4095)) THEN
710 mocchiut 1.1 xhelp1 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
711     xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
712     ds = xhelp1-xhelp2
713     ihelp=(tof21_i-1)*3+tof32_i
714     c1 = k_S21S32(1,ihelp)
715     c2 = k_S21S32(2,ihelp)
716 mocchiut 1.5 beta_a(6) = c2*F/(ds-c1)
717 mocchiut 1.2
718     C------- ToF Mask - S21 - S32
719    
720     tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
721     $ tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
722     tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
723     $ tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
724    
725     tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
726     $ tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
727     tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
728     $ tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
729    
730     C-------
731    
732 mocchiut 1.1 ENDIF
733 mocchiut 1.5 ENDIF
734 mocchiut 1.1
735     C S22 - S31
736 mocchiut 1.5
737     dist = ZTOF(4) - ZTOF(5)
738     dl = 0.
739     DO I=4,5
740     dl = dl + TLOUT(i)
741     ENDDO
742     F = dl/dist
743    
744 mocchiut 1.1 IF (tof22_i.GT.none_find.AND.tof31_i.GT.none_find) THEN
745 mocchiut 1.5 t1 = tof22(1,tof22_i,itdc)
746     t2 = tof22(2,tof22_i,itdc)
747     t3 = tof31(1,tof31_i,itdc)
748     t4 = tof31(2,tof31_i,itdc)
749     IF ((t1.lt.4095).and.(t2.lt.4095).and.
750     & (t3.lt.4095).and.(t4.lt.4095)) THEN
751 mocchiut 1.1 xhelp1 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
752     xhelp2 = tof31(1,tof31_i,itdc)+tof31(2,tof31_i,itdc)
753     ds = xhelp1-xhelp2
754     ihelp=(tof22_i-1)*3+tof31_i
755     c1 = k_S22S31(1,ihelp)
756     c2 = k_S22S31(2,ihelp)
757 mocchiut 1.5 beta_a(7) = c2*F/(ds-c1)
758 mocchiut 1.2
759     C------- ToF Mask - S22 - S31
760    
761     tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
762     $ tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
763     tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
764     $ tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
765    
766     tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
767     $ tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
768     tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
769     $ tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
770    
771     C-------
772    
773 mocchiut 1.1 ENDIF
774 mocchiut 1.5 ENDIF
775 mocchiut 1.1
776     C S22 - S32
777 mocchiut 1.5
778     dist = ZTOF(4) - ZTOF(6)
779     dl = 0.
780     DO I=4,6
781     dl = dl + TLOUT(i)
782     ENDDO
783     F = dl/dist
784    
785 mocchiut 1.1 IF (tof22_i.GT.none_find.AND.tof32_i.GT.none_find) THEN
786 mocchiut 1.5 t1 = tof22(1,tof22_i,itdc)
787     t2 = tof22(2,tof22_i,itdc)
788     t3 = tof32(1,tof32_i,itdc)
789     t4 = tof32(2,tof32_i,itdc)
790     IF ((t1.lt.4095).and.(t2.lt.4095).and.
791     & (t3.lt.4095).and.(t4.lt.4095)) THEN
792 mocchiut 1.1 xhelp1 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
793     xhelp2 = tof32(1,tof32_i,itdc)+tof32(2,tof32_i,itdc)
794     ds = xhelp1-xhelp2
795     ihelp=(tof22_i-1)*3+tof32_i
796     c1 = k_S22S32(1,ihelp)
797     c2 = k_S22S32(2,ihelp)
798 mocchiut 1.5 beta_a(8) = c2*F/(ds-c1)
799 mocchiut 1.2
800     C------- ToF Mask - S22 - S32
801    
802     tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
803     $ tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
804     tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
805     $ tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
806    
807     tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
808     $ tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
809     tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
810     $ tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
811    
812     C-------
813    
814 mocchiut 1.1 ENDIF
815 mocchiut 1.5 ENDIF
816 mocchiut 1.1
817     C S11 - S21
818 mocchiut 1.5
819     dist = ZTOF(1) - ZTOF(3)
820     dl = 0.
821     DO I=1,3
822     dl = dl + TLOUT(i)
823     ENDDO
824     F = dl/dist
825    
826 mocchiut 1.1 IF (tof11_i.GT.none_find.AND.tof21_i.GT.none_find) THEN
827 mocchiut 1.5 t1 = tof11(1,tof11_i,itdc)
828     t2 = tof11(2,tof11_i,itdc)
829     t3 = tof21(1,tof21_i,itdc)
830     t4 = tof21(2,tof21_i,itdc)
831     IF ((t1.lt.4095).and.(t2.lt.4095).and.
832     & (t3.lt.4095).and.(t4.lt.4095)) THEN
833 mocchiut 1.1 xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
834     xhelp2 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
835     ds = xhelp1-xhelp2
836     ihelp=(tof11_i-1)*2+tof21_i
837     c1 = k_S11S21(1,ihelp)
838     c2 = k_S11S21(2,ihelp)
839 mocchiut 1.5 beta_a(9) = c2*F/(ds-c1)
840 mocchiut 1.2
841     C------- ToF Mask - S11 - S21
842    
843     tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
844     $ tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
845     tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
846     $ tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
847    
848     tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
849     $ tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
850     tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
851     $ tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
852    
853     C-------
854    
855 mocchiut 1.1 ENDIF
856 mocchiut 1.5 ENDIF
857 mocchiut 1.1
858     C S11 - S22
859 mocchiut 1.5
860     dist = ZTOF(1) - ZTOF(4)
861     dl = 0.
862     DO I=1,4
863     dl = dl + TLOUT(i)
864     ENDDO
865     F = dl/dist
866    
867 mocchiut 1.1 IF (tof11_i.GT.none_find.AND.tof22_i.GT.none_find) THEN
868 mocchiut 1.5 t1 = tof11(1,tof11_i,itdc)
869     t2 = tof11(2,tof11_i,itdc)
870     t3 = tof22(1,tof22_i,itdc)
871     t4 = tof22(2,tof22_i,itdc)
872     IF ((t1.lt.4095).and.(t2.lt.4095).and.
873     & (t3.lt.4095).and.(t4.lt.4095)) THEN
874 mocchiut 1.1 xhelp1 = tof11(1,tof11_i,itdc)+tof11(2,tof11_i,itdc)
875     xhelp2 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
876     ds = xhelp1-xhelp2
877     ihelp=(tof11_i-1)*2+tof22_i
878     c1 = k_S11S22(1,ihelp)
879     c2 = k_S11S22(2,ihelp)
880 mocchiut 1.5 beta_a(10) = c2*F/(ds-c1)
881 mocchiut 1.2
882     C------- ToF Mask - S11 - S22
883    
884     tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
885     $ tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
886     tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
887     $ tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
888    
889     tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
890     $ tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
891     tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
892     $ tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
893    
894     C-------
895    
896 mocchiut 1.1 ENDIF
897 mocchiut 1.5 ENDIF
898 mocchiut 1.1
899     C S12 - S21
900 mocchiut 1.5
901     dist = ZTOF(2) - ZTOF(3)
902     dl = 0.
903     DO I=2,3
904     dl = dl + TLOUT(i)
905     ENDDO
906     F = dl/dist
907    
908 mocchiut 1.1 IF (tof12_i.GT.none_find.AND.tof21_i.GT.none_find) THEN
909 mocchiut 1.5 t1 = tof12(1,tof12_i,itdc)
910     t2 = tof12(2,tof12_i,itdc)
911     t3 = tof21(1,tof21_i,itdc)
912     t4 = tof21(2,tof21_i,itdc)
913     IF ((t1.lt.4095).and.(t2.lt.4095).and.
914     & (t3.lt.4095).and.(t4.lt.4095)) THEN
915 mocchiut 1.1 xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
916     xhelp2 = tof21(1,tof21_i,itdc)+tof21(2,tof21_i,itdc)
917     ds = xhelp1-xhelp2
918     ihelp=(tof12_i-1)*2+tof21_i
919     c1 = k_S12S21(1,ihelp)
920     c2 = k_S12S21(2,ihelp)
921 mocchiut 1.5 beta_a(11) = c2*F/(ds-c1)
922 mocchiut 1.2
923     C------- ToF Mask - S12 - S21
924    
925     tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
926     $ tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
927     tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
928     $ tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
929    
930     tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
931     $ tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
932     tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
933     $ tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
934    
935     C-------
936    
937 mocchiut 1.1 ENDIF
938 mocchiut 1.5 ENDIF
939 mocchiut 1.1
940     C S12 - S22
941 mocchiut 1.5
942     dist = ZTOF(2) - ZTOF(4)
943     dl = 0.
944     DO I=2,4
945     dl = dl + TLOUT(i)
946     ENDDO
947     F = dl/dist
948    
949 mocchiut 1.1 IF (tof12_i.GT.none_find.AND.tof22_i.GT.none_find) THEN
950 mocchiut 1.5 t1 = tof12(1,tof12_i,itdc)
951     t2 = tof12(2,tof12_i,itdc)
952     t3 = tof22(1,tof22_i,itdc)
953     t4 = tof22(2,tof22_i,itdc)
954     IF ((t1.lt.4095).and.(t2.lt.4095).and.
955     & (t3.lt.4095).and.(t4.lt.4095)) THEN
956 mocchiut 1.1 xhelp1 = tof12(1,tof12_i,itdc)+tof12(2,tof12_i,itdc)
957     xhelp2 = tof22(1,tof22_i,itdc)+tof22(2,tof22_i,itdc)
958     ds = xhelp1-xhelp2
959     ihelp=(tof12_i-1)*2+tof22_i
960     c1 = k_S12S22(1,ihelp)
961     c2 = k_S12S22(2,ihelp)
962 mocchiut 1.5 beta_a(12) = c2*F/(ds-c1)
963 mocchiut 1.2
964     C------- ToF Mask - S12 - S22
965    
966     tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
967     $ tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
968     tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
969     $ tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
970    
971     tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
972     $ tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
973     tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
974     $ tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
975    
976     C-------
977    
978 mocchiut 1.1 ENDIF
979 mocchiut 1.5 ENDIF
980 mocchiut 1.1
981     C-------
982    
983     icount=0
984     sw=0.
985     sxw=0.
986     beta_mean=100.
987    
988     do i=1,12
989     if ((beta_a(i).gt.-1.5).and.(beta_a(i).lt.1.5)) then
990     icount= icount+1
991     if (i.le.4) w_i=1./(0.13**2.)
992     if ((i.ge.5).and.(i.le.8)) w_i=1./(0.16**2.)
993     if (i.ge.9) w_i=1./(0.25**2.) ! to be checked
994     sxw=sxw + beta_a(i)*w_i
995     sw =sw + w_i
996     endif
997     enddo
998    
999     if (icount.gt.0) beta_mean=sxw/sw
1000     beta_a(13) = beta_mean
1001    
1002 mocchiut 1.5 C write(*,*)'TOFTRK ',ievent,beta_a(1),beta_a(2),beta_a(3),beta_a(4)
1003    
1004 mocchiut 1.1 RETURN
1005     END
1006    
1007    
  ViewVC Help
Powered by ViewVC 1.1.23