/[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.5 by mocchiut, Tue Sep 12 13:58:24 2006 UTC revision 1.7 by pam-de, Wed Dec 6 07:21:32 2006 UTC
# Line 47  c     define TOF Z-coordinates Line 47  c     define TOF Z-coordinates
47                
48        INTEGER j        INTEGER j
49    
50        REAL theta12,theta13,theta23  c      REAL theta12,theta13,theta23
51          REAL theta13
52  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
53        REAL tofarm12        REAL tofarm12
54        PARAMETER (tofarm12 = 29.70) ! from 53.39 to 23.69        PARAMETER (tofarm12 = 29.70) ! from 53.39 to 23.69
# Line 359  C-----------------------------S1 ------- Line 360  C-----------------------------S1 -------
360           i = tof11_i           i = tof11_i
361           xdummy=tof11(left,i,iadc)           xdummy=tof11(left,i,iadc)
362           tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta13)           tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta13)
363           if (tof11(left,i,iadc).lt.1000) then           if (tof11(left,i,iadc).lt.4095) then
364              xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))              xkorr=adcx11(left,i,1)*exp(-yhelp/adcx11(left,i,2))
365              xkorr0=adcx11(left,i,1)              xkorr0=adcx11(left,i,1)
366              adc_c(ch11a(i),hb11a(i))=tof11(left,i,iadc)/xkorr              adc_c(ch11a(i),hb11a(i))=tof11(left,i,iadc)/xkorr
367           endif           endif
368    
369           tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta13)           tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta13)
370           if (tof11(right,i,iadc).lt.1000) then           if (tof11(right,i,iadc).lt.4095) then
371              xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))              xkorr=adcx11(right,i,1)*exp(yhelp/adcx11(right,i,2))
372              xkorr0=adcx11(right,i,1)              xkorr0=adcx11(right,i,1)
373              adc_c(ch11b(i),hb11b(i))=tof11(right,i,iadc)/xkorr              adc_c(ch11b(i),hb11b(i))=tof11(right,i,iadc)/xkorr
# Line 379  C-----------------------------S1 ------- Line 380  C-----------------------------S1 -------
380    
381           i = tof12_i           i = tof12_i
382           tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta13)           tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta13)
383           if (tof12(left,i,iadc).lt.1000) then           if (tof12(left,i,iadc).lt.4095) then
384              xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))              xkorr=adcx12(left,i,1)*exp(-xhelp/adcx12(left,i,2))
385              xkorr0=adcx12(left,i,1)              xkorr0=adcx12(left,i,1)
386              adc_c(ch12a(i),hb12a(i))=tof12(left,i,iadc)/xkorr              adc_c(ch12a(i),hb12a(i))=tof12(left,i,iadc)/xkorr
387           endif           endif
388    
389           tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta13)           tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta13)
390           if (tof12(right,i,iadc).lt.1000) then           if (tof12(right,i,iadc).lt.4095) then
391              xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))              xkorr=adcx12(right,i,1)*exp(xhelp/adcx12(right,i,2))
392              xkorr0=adcx12(right,i,1)              xkorr0=adcx12(right,i,1)
393              adc_c(ch12b(i),hb12b(i))=tof12(right,i,iadc)/xkorr              adc_c(ch12b(i),hb12b(i))=tof12(right,i,iadc)/xkorr
# Line 400  C-----------------------------S2 ------- Line 401  C-----------------------------S2 -------
401    
402           i = tof21_i           i = tof21_i
403           tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta13)           tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta13)
404           if (tof21(left,i,iadc).lt.1000) then           if (tof21(left,i,iadc).lt.4095) then
405              xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))              xkorr=adcx21(left,i,1)*exp(-xhelp/adcx21(left,i,2))
406              xkorr0=adcx21(left,i,1)              xkorr0=adcx21(left,i,1)
407              adc_c(ch21a(i),hb21a(i))=tof21(left,i,iadc)/xkorr              adc_c(ch21a(i),hb21a(i))=tof21(left,i,iadc)/xkorr
408           endif           endif
409    
410           tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta13)           tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta13)
411           if (tof21(right,i,iadc).lt.1000) then           if (tof21(right,i,iadc).lt.4095) then
412              xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))              xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))
413              xkorr0=adcx21(right,i,1)              xkorr0=adcx21(right,i,1)
414              adc_c(ch21b(i),hb21b(i))=tof21(right,i,iadc)/xkorr              adc_c(ch21b(i),hb21b(i))=tof21(right,i,iadc)/xkorr
# Line 419  C-----------------------------S2 ------- Line 420  C-----------------------------S2 -------
420    
421           i = tof22_i           i = tof22_i
422           tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta13)           tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta13)
423           if (tof22(left,i,iadc).lt.1000) then           if (tof22(left,i,iadc).lt.4095) then
424              xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))              xkorr=adcx22(left,i,1)*exp(-yhelp/adcx22(left,i,2))
425              xkorr0=adcx22(left,i,1)              xkorr0=adcx22(left,i,1)
426              adc_c(ch22a(i),hb22a(i))=tof22(left,i,iadc)/xkorr              adc_c(ch22a(i),hb22a(i))=tof22(left,i,iadc)/xkorr
427           endif           endif
428    
429           tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta13)           tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta13)
430           if (tof22(right,i,iadc).lt.1000) then           if (tof22(right,i,iadc).lt.4095) then
431              xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))              xkorr=adcx22(right,i,1)*exp(yhelp/adcx22(right,i,2))
432              xkorr0=adcx22(right,i,1)              xkorr0=adcx22(right,i,1)
433              adc_c(ch22b(i),hb22b(i))=tof22(right,i,iadc)/xkorr              adc_c(ch22b(i),hb22b(i))=tof22(right,i,iadc)/xkorr
# Line 440  C-----------------------------S3 ------- Line 441  C-----------------------------S3 -------
441    
442           i = tof31_i           i = tof31_i
443           tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta13)           tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta13)
444           if (tof31(left,i,iadc).lt.1000) then           if (tof31(left,i,iadc).lt.4095) then
445              xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))              xkorr=adcx31(left,i,1)*exp(-yhelp/adcx31(left,i,2))
446              xkorr0=adcx31(left,i,1)              xkorr0=adcx31(left,i,1)
447              adc_c(ch31a(i),hb31a(i))=tof31(left,i,iadc)/xkorr              adc_c(ch31a(i),hb31a(i))=tof31(left,i,iadc)/xkorr
448           endif           endif
449    
450           tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta13)           tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta13)
451           if (tof31(right,i,iadc).lt.1000) then           if (tof31(right,i,iadc).lt.4095) then
452              xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))              xkorr=adcx31(right,i,1)*exp(yhelp/adcx31(right,i,2))
453              xkorr0=adcx31(right,i,1)              xkorr0=adcx31(right,i,1)
454              adc_c(ch31b(i),hb31b(i))=tof31(right,i,iadc)/xkorr              adc_c(ch31b(i),hb31b(i))=tof31(right,i,iadc)/xkorr
# Line 459  C-----------------------------S3 ------- Line 460  C-----------------------------S3 -------
460    
461           i = tof32_i           i = tof32_i
462           tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta13)           tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta13)
463           if (tof32(left,i,iadc).lt.1000) then           if (tof32(left,i,iadc).lt.4095) then
464              xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))              xkorr=adcx32(left,i,1)*exp(-xhelp/adcx32(left,i,2))
465              xkorr0=adcx32(left,i,1)              xkorr0=adcx32(left,i,1)
466              adc_c(ch32a(i),hb32a(i))=tof32(left,i,iadc)/xkorr              adc_c(ch32a(i),hb32a(i))=tof32(left,i,iadc)/xkorr
467           endif           endif
468    
469           tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta13)           tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta13)
470           if (tof32(right,i,iadc).lt.1000) then           if (tof32(right,i,iadc).lt.4095) then
471              xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))              xkorr=adcx32(right,i,1)*exp(xhelp/adcx32(right,i,2))
472              xkorr0=adcx32(right,i,1)              xkorr0=adcx32(right,i,1)
473              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.5  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.23