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

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

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

revision 1.9 by pamelats, Thu Nov 27 13:47:54 2008 UTC revision 1.13 by mocchiut, Mon Nov 23 09:50:50 2009 UTC
# Line 31  C             call comes from "tofl2com" Line 31  C             call comes from "tofl2com"
31  C  mar-08 WM: Bug found in dEdx if check_charge>1  C  mar-08 WM: Bug found in dEdx if check_charge>1
32  C  oct-08 WM: Calculation of zenith angle debugged, sometimes strange values  C  oct-08 WM: Calculation of zenith angle debugged, sometimes strange values
33  C             were possible  C             were possible
34    C  nov-09 WM: the dEdx part ("adctof_c") moved to the new dEdx routine from Napoli
35  C******************************************************************************  C******************************************************************************
36    
37        INTEGER FUNCTION TOFL2COM()        INTEGER FUNCTION TOFL2COM()
# Line 107  c     second index : 1... number of padd Line 108  c     second index : 1... number of padd
108  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
109        REAL tofarm12        REAL tofarm12
110        PARAMETER (tofarm12 = 29.70) ! from 53.39 to 23.69        PARAMETER (tofarm12 = 29.70) ! from 53.39 to 23.69
111        REAL tofarm23                REAL tofarm23        
112        PARAMETER (tofarm23 = 47.61)  ! from 23.69 to -23.92        PARAMETER (tofarm23 = 47.61)  ! from 23.69 to -23.92
113        REAL tofarm13        REAL tofarm13
114        PARAMETER (tofarm13 = 77.31)  ! from 53.39 to -23.92        PARAMETER (tofarm13 = 77.31)  ! from 53.39 to -23.92
# Line 636  C-----------------------------S3 ------- Line 637  C-----------------------------S3 -------
637  C----------------------------------------------------------------------  C----------------------------------------------------------------------
638  C---------------------  zenith angle theta  ---------------------------  C---------------------  zenith angle theta  ---------------------------
639  C----------------------------------------------------------------------  C----------------------------------------------------------------------
 C-----------------------------  old  ----------------------------------  
 c      dx=0.  
 c      dy=0.  
 c      dr=0.  
 c      theta13 = 0.  
 c  
 c         IF ((tof12_i.GT.none_find).AND.(tof32_i.GT.none_find))  
 c     &        dx  = xtofpos(1) - xtofpos(3)  
 c         IF ((tof11_i.GT.none_find).AND.(tof31_i.GT.none_find))  
 c     &        dy  = ytofpos(1) - ytofpos(3)  
 c         dr = sqrt(dx*dx+dy*dy)  
 c         theta13 = atan(dr/tofarm13)  
 c  
 c  
 C-----------------------------  new  ------------------------------  
640    
641         xhelp1=0.         xhelp1=0.
642         if (tof11_i.GT.none_find) xhelp1=tof11_x(tof11_i)         if (tof11_i.GT.none_find) xhelp1=tof11_x(tof11_i)
# Line 958  C---- Line 944  C----
944           tdc_c(ch32b(i),hb32b(i))=tof32(right,i,itdc)           tdc_c(ch32b(i),hb32b(i))=tof32(right,i,itdc)
945                                               ENDIF                                               ENDIF
946        ENDDO        ENDDO
947                  
948    
949  C---------------------------------------------------------------  C---------------------------------------------------------------
950  C--- calculate track position in paddle using timing difference  C--- calculate track position in paddle using timing difference
# Line 1070  C-------------------  improved calculati Line 1056  C-------------------  improved calculati
1056    
1057    
1058  C------------------------------------------------------------------  C------------------------------------------------------------------
1059  C----------------------------------------------------------------------  C------------------------------------------------------------------
1060  C------------------angle and ADC(x) correction  C-------angle and ADC(x) correction: moved to new dEdx routine
1061  C----------------------------------------------------------------------  C------------------------------------------------------------------
1062  C-----------------------------S1 --------------------------------  C------------------------------------------------------------------
 c middle y (or x) position of the upper and middle ToF-Paddle  
 c       DATA tof11_x/ -17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85/  
 c       DATA tof12_y/ -13.75,-8.25,-2.75,2.75,8.25,13.75/  
 c       DATA tof21_y/  3.75,-3.75/     ! paddles in different order  
 c       DATA tof22_x/ -4.5,4.5/  
 c       DATA tof31_x/ -6.0,0.,6.0/  
 c       DATA tof32_y/ -5.0,0.0,5.0/  
   
       yhelp=0.  
 c      yhelp=100.  
       if (tof12_i.GT.none_find) yhelp=tof12_y(tof12_i)  
       if (ytofpos(1).lt.100)  yhelp=ytofpos(1)  
   
       IF (tof11_i.GT.none_find.AND.abs(yhelp).lt.100) THEN  
   
          i = tof11_i  
          if (tof11(left,i,iadc).lt.3786) then  
 c          if (adc(ch11a(i),hb11a(i)).lt.4095) then  
             tof11(left,i,iadc) = tof11(left,i,iadc)*cos(theta13)  
             xkorr = atten(left,11,i,yhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch11a(i),hb11a(i))=tof11(left,i,iadc)/xkorr  
          endif  
   
          if (tof11(right,i,iadc).lt.3786) then  
 c          if (adc(ch11b(i),hb11b(i)).lt.4095) then  
             tof11(right,i,iadc) = tof11(right,i,iadc)*cos(theta13)  
             xkorr = atten(right,11,i,yhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch11b(i),hb11b(i))=tof11(right,i,iadc)/xkorr  
          endif  
       ENDIF  
   
       xhelp=0.  
 c      xhelp=100.  
       if (tof11_i.GT.none_find) xhelp=tof11_x(tof11_i)  
       if (xtofpos(1).lt.100)  xhelp=xtofpos(1)  
   
       IF (tof12_i.GT.none_find.AND.abs(xhelp).lt.100) THEN  
   
          i = tof12_i  
          if (tof12(left,i,iadc).lt.3786) then  
 c          if (adc(ch12a(i),hb12a(i)).lt.4095) then  
             tof12(left,i,iadc) = tof12(left,i,iadc)*cos(theta13)  
             xkorr = atten(left,12,i,xhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch12a(i),hb12a(i))=tof12(left,i,iadc)/xkorr  
          endif  
   
          if (tof12(right,i,iadc).lt.3786) then  
 c          if (adc(ch12b(i),hb12b(i)).lt.4095) then  
             tof12(right,i,iadc) = tof12(right,i,iadc)*cos(theta13)  
             xkorr = atten(right,12,i,xhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch12b(i),hb12b(i))=tof12(right,i,iadc)/xkorr  
          endif  
       ENDIF  
   
 C-----------------------------S2 --------------------------------  
   
       xhelp=0.  
 c      xhelp=100.  
       if (tof22_i.GT.none_find) xhelp=tof22_x(tof22_i)  
       if (xtofpos(2).lt.100)  xhelp=xtofpos(2)  
   
       IF (tof21_i.GT.none_find.AND.abs(xhelp).lt.100) THEN  
   
          i = tof21_i  
          if (tof21(left,i,iadc).lt.3786) then  
 c          if (adc(ch21a(i),hb21a(i)).lt.4095) then  
             tof21(left,i,iadc) = tof21(left,i,iadc)*cos(theta13)  
             xkorr = atten(left,21,i,xhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch21a(i),hb21a(i))=tof21(left,i,iadc)/xkorr  
          endif  
   
          if (tof21(right,i,iadc).lt.3786) then  
 c          if (adc(ch21b(i),hb21b(i)).lt.4095) then  
             tof21(right,i,iadc) = tof21(right,i,iadc)*cos(theta13)  
             xkorr=adcx21(right,i,1)*exp(xhelp/adcx21(right,i,2))  
             xkorr = atten(right,21,i,xhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch21b(i),hb21b(i))=tof21(right,i,iadc)/xkorr  
          endif  
       ENDIF  
   
   
       yhelp=0.  
 c      yhelp=100.  
       if (tof21_i.GT.none_find) yhelp=tof21_y(tof21_i)  
       if (ytofpos(2).lt.100)  yhelp=ytofpos(2)  
   
       IF (tof22_i.GT.none_find.AND.abs(yhelp).lt.100) THEN  
   
          i = tof22_i  
          if (tof22(left,i,iadc).lt.3786) then  
 c          if (adc(ch22a(i),hb22a(i)).lt.4095) then  
             tof22(left,i,iadc) = tof22(left,i,iadc)*cos(theta13)  
             xkorr = atten(left,22,i,yhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch22a(i),hb22a(i))=tof22(left,i,iadc)/xkorr  
          endif  
   
          if (tof22(right,i,iadc).lt.3786) then  
 c          if (adc(ch22b(i),hb22b(i)).lt.4095) then  
             tof22(right,i,iadc) = tof22(right,i,iadc)*cos(theta13)  
             xkorr = atten(right,22,i,yhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch22b(i),hb22b(i))=tof22(right,i,iadc)/xkorr  
          endif  
       ENDIF  
   
 C-----------------------------S3 --------------------------------  
   
       yhelp=0.  
 c      yhelp=100.  
       if (tof32_i.GT.none_find) yhelp=tof32_y(tof32_i)  
       if (ytofpos(3).lt.100)  yhelp=ytofpos(3)  
   
       IF (tof31_i.GT.none_find.AND.abs(yhelp).lt.100) THEN  
   
          i = tof31_i  
          if (tof31(left,i,iadc).lt.3786) then  
 c          if (adc(ch31a(i),hb31a(i)).lt.4095) then  
             tof31(left,i,iadc) = tof31(left,i,iadc)*cos(theta13)  
             xkorr = atten(left,31,i,yhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch31a(i),hb31a(i))=tof31(left,i,iadc)/xkorr  
          endif  
   
          if (tof31(right,i,iadc).lt.3786) then  
 c          if (adc(ch31b(i),hb31b(i)).lt.4095) then  
             tof31(right,i,iadc) = tof31(right,i,iadc)*cos(theta13)  
             xkorr = atten(right,31,i,yhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch31b(i),hb31b(i))=tof31(right,i,iadc)/xkorr  
          endif  
       ENDIF  
   
       xhelp=0.  
 c      xhelp=100.  
       if (tof31_i.GT.none_find) xhelp=tof31_x(tof31_i)  
       if (xtofpos(3).lt.100)  xhelp=xtofpos(3)  
   
       IF (tof32_i.GT.none_find.AND.abs(xhelp).lt.100) THEN  
   
          i = tof32_i  
          if (tof32(left,i,iadc).lt.3786) then  
 c          if (adc(ch32a(i),hb32a(i)).lt.4095) then  
             tof32(left,i,iadc) = tof32(left,i,iadc)*cos(theta13)  
             xkorr = atten(left,32,i,xhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch32a(i),hb32a(i))=tof32(left,i,iadc)/xkorr  
          endif  
   
          if (tof32(right,i,iadc).lt.3786) then  
 c          if (adc(ch32b(i),hb32b(i)).lt.4095) then  
             tof32(right,i,iadc) = tof32(right,i,iadc)*cos(theta13)  
             xkorr = atten(right,32,i,xhelp)  
             xkorr=xkorr/hepratio  
             adctof_c(ch32b(i),hb32b(i))=tof32(right,i,iadc)/xkorr  
          endif  
       ENDIF  
1063    
1064  C--------------------------------------------------------------------  C--------------------------------------------------------------------
1065  C----------------------calculate Beta  ------------------------------  C----------------------calculate Beta  ------------------------------
# Line 1609  C         if ((betatof_a(i).gt.-1.5).and Line 1432  C         if ((betatof_a(i).gt.-1.5).and
1432  C            icount= icount+1  C            icount= icount+1
1433  C            if (i.le.4) w_i=1./(0.13**2.)  C            if (i.le.4) w_i=1./(0.13**2.)
1434  C            if ((i.ge.5).and.(i.le.8)) w_i=1./(0.16**2.)  C            if ((i.ge.5).and.(i.le.8)) w_i=1./(0.16**2.)
1435  C           if (i.ge.9) w_i=1./(0.25**2.)     ! to be checked  C         if (i.ge.9) w_i=1./(0.25**2.)     ! to be checked
1436  C            sxw=sxw + betatof_a(i)*w_i  C            sxw=sxw + betatof_a(i)*w_i
1437  C            sw =sw + w_i  C            sw =sw + w_i
1438  C         endif  C         endif
# Line 1639  c      write(*,*) ytofpos Line 1462  c      write(*,*) ytofpos
1462  c      write(*,*) xtr_tof  c      write(*,*) xtr_tof
1463  c      write(*,*) ytr_tof  c      write(*,*) ytr_tof
1464                
1465   100  continue  c 100  continue
1466            continue
1467    
1468  C  C
1469        RETURN        RETURN
# Line 1783  C k2 constant is taken to be the standar Line 1607  C k2 constant is taken to be the standar
1607          REAL  a1,a2          REAL  a1,a2
1608          INTEGER jj          INTEGER jj
1609    
1610    c get rid of warnings EMILIANO
1611            i = 0
1612            slope = 0
1613            offset = 0
1614            none_find = 0
1615            none_ev = 0
1616            adc_ev = 0
1617            tdc_ev = 0
1618            iadc = 0
1619            itdc = 0
1620            right = 0
1621            left = 0
1622            tof12_y(1) = tof12_y(1)
1623            tof11_x(1) = tof11_x(1)
1624            tof21_y(1) = tof21_y(1)
1625            tof22_x(1) = tof22_x(1)
1626            tof32_y(1) = tof32_y(1)
1627            tof31_x(1) = tof31_x(1)
1628    c get rid of warnings
1629    
1630  C-----------------------------------------------------------  C-----------------------------------------------------------
1631  C--- get data  C--- get data
1632  C-----------------------------------------------------------  C-----------------------------------------------------------
# Line 2166  C*************************************** Line 2010  C***************************************
2010          DATA itop /1,1,2,2,3,3,4,4,1,1,2,2/          DATA itop /1,1,2,2,3,3,4,4,1,1,2,2/
2011          DATA ibot /5,6,5,6,5,6,5,6,3,4,3,4/          DATA ibot /5,6,5,6,5,6,5,6,3,4,3,4/
2012    
2013    
2014    c get rid of warnings EMILIANO
2015            slope = 0
2016            offset = 0
2017            none_find = 0
2018            none_ev = 0
2019            adc_ev = 0
2020            tdc_ev = 0
2021            iadc = 0
2022            itdc = 0
2023            right = 0
2024            left = 0
2025            tof12_y(1) = tof12_y(1)
2026            tof11_x(1) = tof11_x(1)
2027            tof21_y(1) = tof21_y(1)
2028            tof22_x(1) = tof22_x(1)
2029            tof32_y(1) = tof32_y(1)
2030            tof31_x(1) = tof31_x(1)
2031    c get rid of warnings
2032    
2033  C====================================================================  C====================================================================
2034    
2035          tof11_i = hitvec(1)          tof11_i = hitvec(1)
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.23