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

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

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by mocchiut, Sat Jun 17 12:14:56 2006 UTC revision 1.3 by mocchiut, Thu Aug 10 06:32:05 2006 UTC
# Line 32  c     define TOF Z-coordinates Line 32  c     define TOF Z-coordinates
32        INTEGER icount        INTEGER icount
33    
34        REAL beta_mean        REAL beta_mean
35          INTEGER j
36    
37        REAL theta12,theta13,theta23        REAL theta12,theta13,theta23
38  C--   DATA ZTOF/53.74,53.04,23.94,23.44,-23.49,-24.34/ !Sergio 9.05.2006  C--   DATA ZTOF/53.74,53.04,23.94,23.44,-23.49,-24.34/ !Sergio 9.05.2006
# Line 81  C     Line 81  C    
81           beta_a(i) = 100.           beta_a(i) = 100.
82        enddo        enddo
83    
84          do i=1,4
85             do j=1,12
86                adc_c(i,j) = 1000.
87             enddo
88          enddo
89    
90          do i=1,12
91             do j=1,4
92                tofmask(j,i) = 0
93             enddo
94          enddo
95    
96  C------ read  tracking routine  C------ read  tracking routine
97  *     igoodevent = igoodevent+1  *     igoodevent = igoodevent+1
98  *     assigned input  parameters for track routine  *     assigned input  parameters for track routine
# Line 229  C--------------     S32 ---------------- Line 241  C--------------     S32 ----------------
241        enddo        enddo
242        endif        endif
243    
244    C------------------------------------------------------------------
245    C---  calculate track position in paddle using timing difference
246    C------------------------------------------------------------------
247    
248          do i=1,3
249             xtofpos(i)=100.
250             ytofpos(i)=100.
251          enddo
252    C-----------------------------S1 --------------------------------
253          
254          IF (tof11_i.GT.none_find) THEN
255             ytofpos(1)  = ((tof11(1,tof11_i,itdc)-tof11(2,tof11_i,itdc))/2.
256         +        -y_coor_lin11(tof11_i,offset))/y_coor_lin11(tof11_i,slope)
257          endif
258    
259          IF (tof12_i.GT.none_find) THEN
260             xtofpos(1)  = ((tof12(1,tof12_i,itdc)-tof12(2,tof12_i,itdc))/2.
261         +        -x_coor_lin12(tof12_i,offset))/x_coor_lin12(tof12_i,slope)
262          endif
263          
264          
265    C-----------------------------S2 --------------------------------
266          
267          IF (tof21_i.GT.none_find) THEN
268             xtofpos(2) = ((tof21(1,tof21_i,itdc)-tof21(2,tof21_i,itdc))/2.
269         +        -x_coor_lin21(tof21_i,offset))/x_coor_lin21(tof21_i,slope)
270          endif
271          
272          IF (tof22_i.GT.none_find) THEN
273             ytofpos(2) = ((tof22(1,tof22_i,itdc)-tof22(2,tof22_i,itdc))/2.
274         +        -y_coor_lin22(tof22_i,offset))/y_coor_lin22(tof22_i,slope)
275          endif
276          
277          
278    C-----------------------------S3 --------------------------------
279          
280          IF (tof31_i.GT.none_find) THEN
281             ytofpos(3)  = ((tof31(1,tof31_i,itdc)-tof31(2,tof31_i,itdc))/2.
282         +        -y_coor_lin31(tof31_i,offset))/y_coor_lin31(tof31_i,slope)
283          endif
284          
285          IF (tof32_i.GT.none_find) THEN
286             xtofpos(3)  = ((tof32(1,tof32_i,itdc)-tof32(2,tof32_i,itdc))/2.
287         +        -x_coor_lin32(tof32_i,offset))/x_coor_lin32(tof32_i,slope)
288          endif
289          
290          
291          do i=1,3
292             if (abs(xtofpos(i)).gt.100.) then
293                xtofpos(i)=101.
294             endif
295             if (abs(ytofpos(i)).gt.100.) then
296                ytofpos(i)=101.
297             endif
298          enddo
299          
300  C----------------------------------------------------------------------  C----------------------------------------------------------------------
301  C--------------------Corrections on ADC-data -------------------------  C--------------------Corrections on ADC-data -------------------------
302  C---------------------zenith angle theta  ---------------------------  C---------------------zenith angle theta  ---------------------------
# Line 260  C-----------------------------S1 ------- Line 328  C-----------------------------S1 -------
328           i = tof11_i           i = tof11_i
329           xdummy=tof11(left,i,iadc)           xdummy=tof11(left,i,iadc)
330           tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta13)           tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta13)
331           if (tof11(left,i,iadc).lt.1000) then           if (tof11(left,i,iadc).lt.4095) then
332              xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))              xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))
333              xkorr0=adcx11(left,i,1)              xkorr0=adcx11(left,i,1)
334              adc_c(ch11a(i),hb11a(i))=tof11(left,i,iadc)/xkorr              adc_c(ch11a(i),hb11a(i))=tof11(left,i,iadc)/xkorr
335           endif           endif
336    
337           tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta13)           tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta13)
338           if (tof11(right,i,iadc).lt.1000) then           if (tof11(right,i,iadc).lt.4095) then
339              xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))              xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))
340              xkorr0=adcx11(right,i,1)              xkorr0=adcx11(right,i,1)
341              adc_c(ch11b(i),hb11b(i))=tof11(right,i,iadc)/xkorr              adc_c(ch11b(i),hb11b(i))=tof11(right,i,iadc)/xkorr
# Line 280  C-----------------------------S1 ------- Line 348  C-----------------------------S1 -------
348    
349           i = tof12_i           i = tof12_i
350           tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta13)           tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta13)
351           if (tof12(left,i,iadc).lt.1000) then           if (tof12(left,i,iadc).lt.4095) then
352              xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))              xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))
353              xkorr0=adcx12(left,i,1)              xkorr0=adcx12(left,i,1)
354              adc_c(ch12a(i),hb12a(i))=tof12(left,i,iadc)/xkorr              adc_c(ch12a(i),hb12a(i))=tof12(left,i,iadc)/xkorr
355           endif           endif
356    
357           tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta13)           tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta13)
358           if (tof12(right,i,iadc).lt.1000) then           if (tof12(right,i,iadc).lt.4095) then
359              xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))              xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))
360              xkorr0=adcx12(right,i,1)              xkorr0=adcx12(right,i,1)
361              adc_c(ch12b(i),hb12b(i))=tof12(right,i,iadc)/xkorr              adc_c(ch12b(i),hb12b(i))=tof12(right,i,iadc)/xkorr
# Line 301  C-----------------------------S2 ------- Line 369  C-----------------------------S2 -------
369    
370           i = tof21_i           i = tof21_i
371           tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta13)           tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta13)
372           if (tof21(left,i,iadc).lt.1000) then           if (tof21(left,i,iadc).lt.4095) then
373              xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))              xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))
374              xkorr0=adcx21(left,i,1)              xkorr0=adcx21(left,i,1)
375              adc_c(ch21a(i),hb21a(i))=tof21(left,i,iadc)/xkorr              adc_c(ch21a(i),hb21a(i))=tof21(left,i,iadc)/xkorr
376           endif           endif
377    
378           tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta13)           tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta13)
379           if (tof21(right,i,iadc).lt.1000) then           if (tof21(right,i,iadc).lt.4095) then
380              xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))              xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))
381              xkorr0=adcx21(right,i,1)              xkorr0=adcx21(right,i,1)
382              adc_c(ch21b(i),hb21b(i))=tof21(right,i,iadc)/xkorr              adc_c(ch21b(i),hb21b(i))=tof21(right,i,iadc)/xkorr
# Line 320  C-----------------------------S2 ------- Line 388  C-----------------------------S2 -------
388    
389           i = tof22_i           i = tof22_i
390           tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta13)           tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta13)
391           if (tof22(left,i,iadc).lt.1000) then           if (tof22(left,i,iadc).lt.4095) then
392              xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))              xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))
393              xkorr0=adcx22(left,i,1)              xkorr0=adcx22(left,i,1)
394              adc_c(ch22a(i),hb22a(i))=tof22(left,i,iadc)/xkorr              adc_c(ch22a(i),hb22a(i))=tof22(left,i,iadc)/xkorr
395           endif           endif
396    
397           tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta13)           tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta13)
398           if (tof22(right,i,iadc).lt.1000) then           if (tof22(right,i,iadc).lt.4095) then
399              xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))              xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))
400              xkorr0=adcx22(right,i,1)              xkorr0=adcx22(right,i,1)
401              adc_c(ch22b(i),hb22b(i))=tof22(right,i,iadc)/xkorr              adc_c(ch22b(i),hb22b(i))=tof22(right,i,iadc)/xkorr
# Line 341  C-----------------------------S3 ------- Line 409  C-----------------------------S3 -------
409    
410           i = tof31_i           i = tof31_i
411           tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta13)           tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta13)
412           if (tof31(left,i,iadc).lt.1000) then           if (tof31(left,i,iadc).lt.4095) then
413              xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))              xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))
414              xkorr0=adcx31(left,i,1)              xkorr0=adcx31(left,i,1)
415              adc_c(ch31a(i),hb31a(i))=tof31(left,i,iadc)/xkorr              adc_c(ch31a(i),hb31a(i))=tof31(left,i,iadc)/xkorr
416           endif           endif
417    
418           tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta13)           tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta13)
419           if (tof31(right,i,iadc).lt.1000) then           if (tof31(right,i,iadc).lt.4095) then
420              xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))              xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))
421              xkorr0=adcx31(right,i,1)              xkorr0=adcx31(right,i,1)
422              adc_c(ch31b(i),hb31b(i))=tof31(right,i,iadc)/xkorr              adc_c(ch31b(i),hb31b(i))=tof31(right,i,iadc)/xkorr
# Line 360  C-----------------------------S3 ------- Line 428  C-----------------------------S3 -------
428    
429           i = tof32_i           i = tof32_i
430           tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta13)           tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta13)
431           if (tof32(left,i,iadc).lt.1000) then           if (tof32(left,i,iadc).lt.4095) then
432              xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))              xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))
433              xkorr0=adcx32(left,i,1)              xkorr0=adcx32(left,i,1)
434              adc_c(ch32a(i),hb32a(i))=tof32(left,i,iadc)/xkorr              adc_c(ch32a(i),hb32a(i))=tof32(left,i,iadc)/xkorr
435           endif           endif
436    
437           tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta13)           tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta13)
438           if (tof32(right,i,iadc).lt.1000) then           if (tof32(right,i,iadc).lt.4095) then
439              xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))              xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))
440              xkorr0=adcx32(right,i,1)              xkorr0=adcx32(right,i,1)
441              adc_c(ch32b(i),hb32b(i))=tof32(right,i,iadc)/xkorr              adc_c(ch32b(i),hb32b(i))=tof32(right,i,iadc)/xkorr
# Line 395  C     S11 - S31 Line 463  C     S11 - S31
463           c2 = k_S11S31(2,ihelp)                   c2 = k_S11S31(2,ihelp)        
464           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
465       &   beta_a(1) = c2/(cos(theta13)*(ds-c1))       &   beta_a(1) = c2/(cos(theta13)*(ds-c1))
466    
467    C------- ToF Mask - S11 - S31
468    
469             tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
470         $        tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
471             tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
472         $        tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
473    
474             tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
475         $        tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
476             tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
477         $        tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
478    
479    C-------          
480    
481        ENDIF        ENDIF
482                
483  C     S11 - S32  C     S11 - S32
# Line 407  C     S11 - S32 Line 490  C     S11 - S32
490           c2 = k_S11S32(2,ihelp)                   c2 = k_S11S32(2,ihelp)        
491           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
492       &   beta_a(2) = c2/(cos(theta13)*(ds-c1))       &   beta_a(2) = c2/(cos(theta13)*(ds-c1))
493    
494    C------- ToF Mask - S11 - S32
495    
496             tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
497         $        tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
498             tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
499         $        tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
500    
501             tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
502         $        tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
503             tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
504         $        tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
505    
506    C-------          
507    
508        ENDIF        ENDIF
509    
510  C     S12 - S31  C     S12 - S31
# Line 419  C     S12 - S31 Line 517  C     S12 - S31
517           c2 = k_S12S31(2,ihelp)           c2 = k_S12S31(2,ihelp)
518           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
519       &   beta_a(3) = c2/(cos(theta13)*(ds-c1))       &   beta_a(3) = c2/(cos(theta13)*(ds-c1))
520    
521    C------- ToF Mask - S12 - S31
522    
523             tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
524         $        tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
525             tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
526         $        tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
527    
528             tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
529         $        tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
530             tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
531         $        tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
532    
533    C-------
534    
535        ENDIF        ENDIF
536    
537  C     S12 - S32  C     S12 - S32
# Line 431  C     S12 - S32 Line 544  C     S12 - S32
544           c2 = k_S12S32(2,ihelp)           c2 = k_S12S32(2,ihelp)
545           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
546       &   beta_a(4) = c2/(cos(theta13)*(ds-c1))       &   beta_a(4) = c2/(cos(theta13)*(ds-c1))
547    
548    C------- ToF Mask - S12 - S32
549    
550             tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
551         $        tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
552             tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
553         $        tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
554    
555             tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
556         $        tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
557             tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
558         $        tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
559    
560    C-------
561    
562        ENDIF        ENDIF
563    
564  C     S21 - S31  C     S21 - S31
# Line 443  C     S21 - S31 Line 571  C     S21 - S31
571           c2 = k_S21S31(2,ihelp)           c2 = k_S21S31(2,ihelp)
572           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
573       &   beta_a(5) = c2/(cos(theta23)*(ds-c1))       &   beta_a(5) = c2/(cos(theta23)*(ds-c1))
574    
575    C------- ToF Mask - S21 - S31
576    
577             tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
578         $        tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
579             tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
580         $        tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
581    
582             tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
583         $        tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
584             tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
585         $        tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
586    
587    C-------
588    
589        ENDIF        ENDIF
590    
591  C     S21 - S32  C     S21 - S32
# Line 455  C     S21 - S32 Line 598  C     S21 - S32
598           c2 = k_S21S32(2,ihelp)           c2 = k_S21S32(2,ihelp)
599           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
600       &   beta_a(6) = c2/(cos(theta23)*(ds-c1))       &   beta_a(6) = c2/(cos(theta23)*(ds-c1))
601    
602    C------- ToF Mask - S21 - S32
603    
604             tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
605         $        tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
606             tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
607         $        tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
608    
609             tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
610         $        tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
611             tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
612         $        tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
613    
614    C-------
615    
616        ENDIF        ENDIF
617    
618  C     S22 - S31  C     S22 - S31
# Line 467  C     S22 - S31 Line 625  C     S22 - S31
625           c2 = k_S22S31(2,ihelp)           c2 = k_S22S31(2,ihelp)
626           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
627       &   beta_a(7) = c2/(cos(theta13)*(ds-c1))       &   beta_a(7) = c2/(cos(theta13)*(ds-c1))
628    
629    C------- ToF Mask - S22 - S31
630    
631             tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
632         $        tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
633             tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
634         $        tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
635    
636             tofmask(ch31a(tof31_i),hb31a(tof31_i)) =
637         $        tofmask(ch31a(tof31_i),hb31a(tof31_i)) + 1
638             tofmask(ch31b(tof31_i),hb31b(tof31_i)) =
639         $        tofmask(ch31b(tof31_i),hb31b(tof31_i)) + 1
640    
641    C-------  
642    
643        ENDIF        ENDIF
644                
645  C     S22 - S32  C     S22 - S32
# Line 479  C     S22 - S32 Line 652  C     S22 - S32
652           c2 = k_S22S32(2,ihelp)           c2 = k_S22S32(2,ihelp)
653           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
654       &   beta_a(8) = c2/(cos(theta13)*(ds-c1))       &   beta_a(8) = c2/(cos(theta13)*(ds-c1))
655    
656    C------- ToF Mask - S22 - S32
657    
658             tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
659         $        tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
660             tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
661         $        tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
662    
663             tofmask(ch32a(tof32_i),hb32a(tof32_i)) =
664         $        tofmask(ch32a(tof32_i),hb32a(tof32_i)) + 1
665             tofmask(ch32b(tof32_i),hb32b(tof32_i)) =
666         $        tofmask(ch32b(tof32_i),hb32b(tof32_i)) + 1
667    
668    C-------  
669    
670        ENDIF        ENDIF
671    
672  C     S11 - S21  C     S11 - S21
# Line 491  C     S11 - S21 Line 679  C     S11 - S21
679           c2 = k_S11S21(2,ihelp)           c2 = k_S11S21(2,ihelp)
680           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
681       &   beta_a(9) = c2/(cos(theta12)*(ds-c1))       &   beta_a(9) = c2/(cos(theta12)*(ds-c1))
682    
683    C------- ToF Mask - S11 - S21
684    
685             tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
686         $        tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
687             tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
688         $        tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
689    
690             tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
691         $        tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
692             tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
693         $        tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
694    
695    C-------  
696    
697        ENDIF        ENDIF
698                
699  C     S11 - S22  C     S11 - S22
# Line 503  C     S11 - S22 Line 706  C     S11 - S22
706           c2 = k_S11S22(2,ihelp)           c2 = k_S11S22(2,ihelp)
707           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
708       &   beta_a(10) = c2/(cos(theta12)*(ds-c1))       &   beta_a(10) = c2/(cos(theta12)*(ds-c1))
709    
710    C------- ToF Mask - S11 - S22
711    
712             tofmask(ch11a(tof11_i),hb11a(tof11_i)) =
713         $        tofmask(ch11a(tof11_i),hb11a(tof11_i)) + 1
714             tofmask(ch11b(tof11_i),hb11b(tof11_i)) =
715         $        tofmask(ch11b(tof11_i),hb11b(tof11_i)) + 1
716    
717             tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
718         $        tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
719             tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
720         $        tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
721    
722    C-------
723    
724        ENDIF        ENDIF
725    
726  C     S12 - S21  C     S12 - S21
# Line 515  C     S12 - S21 Line 733  C     S12 - S21
733           c2 = k_S12S21(2,ihelp)           c2 = k_S12S21(2,ihelp)
734           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
735       &   beta_a(11) = c2/(cos(theta12)*(ds-c1))       &   beta_a(11) = c2/(cos(theta12)*(ds-c1))
736    
737    C------- ToF Mask - S12 - S21
738    
739             tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
740         $        tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
741             tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
742         $        tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
743    
744             tofmask(ch21a(tof21_i),hb21a(tof21_i)) =
745         $        tofmask(ch21a(tof21_i),hb21a(tof21_i)) + 1
746             tofmask(ch21b(tof21_i),hb21b(tof21_i)) =
747         $        tofmask(ch21b(tof21_i),hb21b(tof21_i)) + 1
748    
749    C-------  
750    
751        ENDIF        ENDIF
752    
753  C     S12 - S22  C     S12 - S22
# Line 527  C     S12 - S22 Line 760  C     S12 - S22
760           c2 = k_S12S22(2,ihelp)           c2 = k_S12S22(2,ihelp)
761           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))           if ((xhelp1.lt.8000.).and.(xhelp2.lt.8000))
762       &   beta_a(12) = c2/(cos(theta12)*(ds-c1))       &   beta_a(12) = c2/(cos(theta12)*(ds-c1))
763    
764    C------- ToF Mask - S12 - S22
765    
766             tofmask(ch12a(tof12_i),hb12a(tof12_i)) =
767         $        tofmask(ch12a(tof12_i),hb12a(tof12_i)) + 1
768             tofmask(ch12b(tof12_i),hb12b(tof12_i)) =
769         $        tofmask(ch12b(tof12_i),hb12b(tof12_i)) + 1
770    
771             tofmask(ch22a(tof22_i),hb22a(tof22_i)) =
772         $        tofmask(ch22a(tof22_i),hb22a(tof22_i)) + 1
773             tofmask(ch22b(tof22_i),hb22b(tof22_i)) =
774         $        tofmask(ch22b(tof22_i),hb22b(tof22_i)) + 1
775    
776    C-------  
777    
778        ENDIF        ENDIF
779                
780  C-------  C-------
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23