/[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.2 by mocchiut, Thu Jul 6 09:03:28 2006 UTC revision 1.4 by mocchiut, Sat Sep 9 05:17:20 2006 UTC
# Line 106  C     COPY THE ALFA VECTOR FROM AL_PP TO Line 106  C     COPY THE ALFA VECTOR FROM AL_PP TO
106           AL_P(i) = al_pp(i)           AL_P(i) = al_pp(i)
107        enddo        enddo
108        if (al_p(5).eq.0.) THEN        if (al_p(5).eq.0.) THEN
109           PRINT *,' TOF - WARNING F77: track with R = 0, discarded'  c         PRINT *,' TOF - WARNING F77: track with R = 0, discarded'
110           GOTO 969           GOTO 969
111        ENDIF        ENDIF
112  *     -------- *** tracking routine *** --------  *     -------- *** tracking routine *** --------
113        IFAIL = 0        IFAIL = 0
114        call DOTRACK(NPTOF,ZIN,XOUT,YOUT,AL_P,IFAIL)        call DOTRACK(NPTOF,ZIN,XOUT,YOUT,AL_P,IFAIL)
115        if(IFAIL.ne.0)then        if(IFAIL.ne.0)then
116           print *,' TOF - WARNING F77: tracking failed '  c         print *,' TOF - WARNING F77: tracking failed '
117           goto 969           goto 969
118        endif        endif
119  *     ------------------------------------------  *     ------------------------------------------
# Line 328  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 348  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 369  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 388  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 409  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 428  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
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.23