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

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

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

revision 1.2 by pam-de, Fri Oct 20 12:22:18 2006 UTC revision 1.8 by mocchiut, Wed Nov 23 21:19:38 2011 UTC
# Line 1  Line 1 
1  C------------------------------------------------  C------------------------------------------------
2        FUNCTION RDTOFCAL(calname)        REAL FUNCTION RDTOFCAL(calname,cnlen)
3  C------------------------------------------------  C------------------------------------------------
4    
5        include 'input_tof.txt'        include 'input_tof.txt'
6    
7  c----- HBOOK  c----- HBOOK
8        INTEGER HMEM        INTEGER HMEM
9    c      REAL RDTOFCAL
10        parameter (NWPAWC=8500000)        parameter (NWPAWC=8500000)
11        common/PAWC/HMEM(NWPAWC)        common/PAWC/HMEM(NWPAWC)
12    
# Line 15  c     local variables Line 16  c     local variables
16  c------------------------------------------------------------------------  c------------------------------------------------------------------------
17        character*250 calname        character*250 calname
18        character*250 calname1        character*250 calname1
19        character ch  c      character ch
20        character*2 ch2  c      character*2 ch2
21          integer cnlen
22  c------------------------------------------------------------------  c------------------------------------------------------------------
23    
24        COMMON/QUEST/IQUEST(100)        COMMON/QUEST/IQUEST(100)
# Line 27  C Line 29  C
29    
30        INTEGER lrec        INTEGER lrec
31        INTEGER istat, ierr, icycle        INTEGER istat, ierr, icycle
32    
33    c      REAL x_coor_lin11(8,2),y_coor_lin12(6,2)
34        PARAMETER (lrec=4096)        PARAMETER (lrec=4096)
35    
36  C-------------------------------------------------------------------------  C-------------------------------------------------------------------------
# Line 35  c     HBOOK initialization Line 39  c     HBOOK initialization
39  c  c
40  c------------------------------------------------------------------------  c------------------------------------------------------------------------
41    
42          RDTOFCAL = 0.
43        call HLIMIT(NWPAWC)        call HLIMIT(NWPAWC)
44    
45  c------------------------------------------------------------------------  c------------------------------------------------------------------------
46  c---------------    open ToF calib rz file  c---------------    open ToF calib rz file
47  c------------------------------------------------------------------------  c------------------------------------------------------------------------
48    
49         do ii=1,250  c       do ii=1,250
50  c         ch=calname(ii:ii)  c         ch=calname(ii:ii)
          ch2=calname(ii:ii+1)  
 c         write(*,*) ii,ch  
51  c         if (ch.eq.'.') goto 100  c         if (ch.eq.'.') goto 100
52           if (ch2.eq.'rz') goto 100  c       enddo  
53         enddo    c100   continue
54  100   continue  c      calname1 = calname(1:ii)//'rz'
55    
56    c      do ii=1,250
57    cc         ch=calname(ii:ii)
58    c         ch2=calname(ii:ii+1)
59    cc         write(*,*) ii,ch
60    cc         if (ch.eq.'.') goto 100
61    c         if (ch2.eq.'rz') goto 100
62    c       enddo
63    c100   continue
64  c      calname1 = calname(1:ii)//'rz'  c      calname1 = calname(1:ii)//'rz'
65         calname1 = calname(1:ii+1)  c       calname1 = calname(1:ii+1)
66    
67          calname1 = calname(1:cnlen)
68    
69    
70    
 c      print*,'__________ opening TOF calib rz file __________'  
 c      print*,calname  
 c      print*,calname1  
71    
72         CALL HROPEN(59,'TOF K1A',calname1,'QP',4096,istat)         CALL HROPEN(59,'TOFK1131',calname1,'QP',4096,istat)
73    
74          if (istat.ne.0) then  ! check if HROPEN was OK          if (istat.ne.0) then  ! check if HROPEN was OK
75                write(*,*) 'Can''t open correct ToF calibration File !!!'                write(*,*) 'Can''t open correct ToF calibration File !!!'
76          stop          stop
77          endif          endif
78    
 c        print*,' reading TOF CALIB n-tuple...'  
   
79          call HRIN(ntp_tofcalib,9999,0)          call HRIN(ntp_tofcalib,9999,0)
80    
81          call HBNAME(ntp_tofcalib,' ',0,'$CLEAR')          call HBNAME(ntp_tofcalib,' ',0,'$CLEAR')
82    
83          call HBNAME(ntp_tofcalib,'TOFKA',k_s11s31,'$SET')          call HBNAME(ntp_tofcalib,'TOFK1131',k_s11s31,'$SET')
84          call HBNAME(ntp_tofcalib,'TOFKB',k_s12s32,'$SET')          call HBNAME(ntp_tofcalib,'TOFK1132',k_s11s32,'$SET')
85          call HBNAME(ntp_tofcalib,'TOFKC',k_s21s31,'$SET')          call HBNAME(ntp_tofcalib,'TOFK1231',k_s12s31,'$SET')
86          call HBNAME(ntp_tofcalib,'TOFKD',k_s22s32,'$SET')          call HBNAME(ntp_tofcalib,'TOFK1232',k_s12s32,'$SET')
87    
88            call HBNAME(ntp_tofcalib,'TOFK2131',k_s21s31,'$SET')
89            call HBNAME(ntp_tofcalib,'TOFK2132',k_s21s32,'$SET')
90            call HBNAME(ntp_tofcalib,'TOFK2231',k_s22s31,'$SET')
91            call HBNAME(ntp_tofcalib,'TOFK2232',k_s22s32,'$SET')
92    
93            call HBNAME(ntp_tofcalib,'TOFK1121',k_s11s21,'$SET')
94            call HBNAME(ntp_tofcalib,'TOFK1122',k_s11s22,'$SET')
95            call HBNAME(ntp_tofcalib,'TOFK1221',k_s12s21,'$SET')
96            call HBNAME(ntp_tofcalib,'TOFK1222',k_s12s22,'$SET')
97    
98    C--------------
99    
100          call HBNAME(ntp_tofcalib,'TOFLIN11',y_coor_lin11,'$SET')          call HBNAME(ntp_tofcalib,'TOFLIN11',y_coor_lin11,'$SET')
101          call HBNAME(ntp_tofcalib,'TOFLIN12',x_coor_lin12,'$SET')          call HBNAME(ntp_tofcalib,'TOFLIN12',x_coor_lin12,'$SET')
102    c        call HBNAME(ntp_tofcalib,'TOFLIN11',x_coor_lin11,'$SET')
103    c        call HBNAME(ntp_tofcalib,'TOFLIN12',y_coor_lin12,'$SET')
104          call HBNAME(ntp_tofcalib,'TOFLIN21',x_coor_lin21,'$SET')          call HBNAME(ntp_tofcalib,'TOFLIN21',x_coor_lin21,'$SET')
105          call HBNAME(ntp_tofcalib,'TOFLIN22',y_coor_lin22,'$SET')          call HBNAME(ntp_tofcalib,'TOFLIN22',y_coor_lin22,'$SET')
106          call HBNAME(ntp_tofcalib,'TOFLIN31',y_coor_lin31,'$SET')          call HBNAME(ntp_tofcalib,'TOFLIN31',y_coor_lin31,'$SET')
# Line 98  c        print*,' reading TOF CALIB n-tu Line 123  c        print*,' reading TOF CALIB n-tu
123  c        call HPRNTU(ntp_tofcalib)  c        call HPRNTU(ntp_tofcalib)
124    
125          call HNOENT(ntp_tofcalib,iemax_cal)          call HNOENT(ntp_tofcalib,iemax_cal)
 c        write(*,*) 'Number of Events CALIB ',iemax_cal  
126    
127          do iev=1,iemax_cal          do iev=1,iemax_cal
128             call HGNT(ntp_tofcalib,iev,ierr) !reads an event             call HGNT(ntp_tofcalib,iev,ierr) !reads an event
129  c           write(34,*)iev,k1_s11s31,y_coor_lin11,tw11  c           print *,' ierr ',ierr
130          enddo          enddo
131    
132          call hrout(ntp_tofcalib,icycle,' ')          call hrout(ntp_tofcalib,icycle,' ')
133          call hrend('TOF K1A')          call hrend('TOFK1131')
134          close(0)          close(0)
135    
136           return          return
137          end          end
138    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.23