/[PAMELA software]/chewbacca/PamOffLineSW/forroutines/calorimeter/calunpack.for
ViewVC logotype

Diff of /chewbacca/PamOffLineSW/forroutines/calorimeter/calunpack.for

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

revision 1.2 by mocchiut, Thu Dec 4 13:53:15 2008 UTC revision 1.4 by mocchiut, Tue Dec 30 14:21:19 2008 UTC
# Line 172  C Line 172  C
172  C          C        
173  C Begin !  C Begin !
174  C  C
175    c      dumpo = iev
176        start = 320        start = 320
177        firsttime = 1        firsttime = 1
178        SOGLIA0 = 70        SOGLIA0 = 70
179          sic = 0
180          sicb = 0
181   2    continue   2    continue
182  C  C
183  C     input length must be > 0, if not go out with error code 142  C     input length must be > 0, if not go out with error code 142
# Line 193  C Line 196  C
196  C  C
197  C     DEBUG: PRINT OUT THE INPUT VECTOR  C     DEBUG: PRINT OUT THE INPUT VECTOR
198  C  C
199        if (iev.eq.dumpo) then  c      if (iev.eq.dumpo) then
200           do l=1,lung  c         do l=1,lung
201              write(*,17)l,vecta(l)  c            write(*,17)l,vecta(l)
202           enddo  c         enddo
203        endif  c      endif
204  C      dumpo = iev  C      dumpo = iev
205  C  C
206  C     DETERMINE LENGTH IN WORDS FROM LENGTH IN BYTES  C     DETERMINE LENGTH IN WORDS FROM LENGTH IN BYTES
# Line 295  C Line 298  C
298  c  c
299   32   continue   32   continue
300  C  C
301          if ( ic .lt. 1 ) then
302             if (dumpo.eq.iev) print *,' AGH IC = ',IC      
303             ic = 1
304          endif
305          if ( icsave .lt. 0 ) then
306             if (dumpo.eq.iev) print *,' AGH ICSAVE = ',ICSAVE      
307             icsave = 0
308          endif
309    C
310        ke = 0        ke = 0
311        chis = chi        chis = chi
312        icold = ic        icold = ic
# Line 835  c     contr = contr + 1             Line 847  c     contr = contr + 1            
847              headcor = 1              headcor = 1
848              ichc = ic - 1              ichc = ic - 1
849              if (iev.eq.dumpo)              if (iev.eq.dumpo)
850       &           print *,'crc is wrong ',ic,       &           print *,' A crc is wrong ',ic,
851       &           ' search section ',contr,' coco = ',coco       &           ' search section ',contr,' coco = ',coco
852              goto 32                  goto 32    
853           else           else
# Line 1126  c     Line 1138  c    
1138           ichc = sic - 1           ichc = sic - 1
1139           icb = sicb           icb = sicb
1140           if (iev.eq.dumpo)           if (iev.eq.dumpo)
1141       &        print *,'crc is wrong ',sic,       &        print *,' B crc is wrong ',sic,
1142       &        ' search section ',contr,' coco = ',coco       &        ' search section ',contr,' coco = ',coco
1143           goto 32               goto 32    
1144  c  c
# Line 1253  c Line 1265  c
1265                 enddo                 enddo
1266              enddo                enddo  
1267           endif           endif
1268           do l=1,lung  c         do l=1,lung
1269              write(*,17)l,vecta(l)  c            write(*,17)l,vecta(l)
1270           enddo  c         enddo
1271        endif        endif
1272        iev = iev + 1              iev = iev + 1      
1273        RETURN        RETURN
# Line 1352  C Line 1364  C
1364        i = inf        i = inf
1365  c  c
1366   10   continue   10   continue
1367        if (i.gt.sup) then        if (i.gt.sup.or.i.gt.120000) then
1368           RETURN           RETURN
1369        endif        endif
1370  C  C
# Line 1392  c Line 1404  c
1404           ipl = int(st/6) + 1           ipl = int(st/6) + 1
1405           ipr = st - (ipl - 1) * 6 + 1           ipr = st - (ipl - 1) * 6 + 1
1406           i = i + 1           i = i + 1
1407           if (i.gt.sup) RETURN           if (i.gt.sup.or.i.gt.120000) RETURN
1408           if ( ipl.ge.1.and.ipl.le.11.and.ipr.ge.1.and.ipr.le.6)           if ( ipl.ge.1.and.ipl.le.11.and.ipr.ge.1.and.ipr.le.6)
1409       +        basse(ipl,ipr) = vect(i)                   +        basse(ipl,ipr) = vect(i)            
1410  c          c        
1411   20      continue   20      continue
1412           if (i.gt.sup) RETURN           if (i.gt.sup.or.i.gt.120000) RETURN
1413  C  C
1414           i = i + 1           i = i + 1
1415           if (i.gt.sup) RETURN           if (i.gt.sup.or.i.gt.120000) RETURN
1416           if ((vect(i)+1).gt.16.or.(vect(i)+1).lt.1) then           if ((vect(i)+1).gt.16.or.(vect(i)+1).lt.1) then
1417              goto 10              goto 10
1418           endif           endif
1419           ist = vect(i) + 1 + 16 * (ipr - 1)           ist = vect(i) + 1 + 16 * (ipr - 1)
1420           i = i + 1           i = i + 1
1421           if (i.gt.sup) RETURN           if (i.gt.sup.or.i.gt.120000) RETURN
1422           if ( ipl.ge.1.and.ipl.le.11.and.ist.ge.1.and.ist.le.96)           if ( ipl.ge.1.and.ipl.le.11.and.ist.ge.1.and.ist.le.96)
1423       +        dedx(ipl,ist) = vect(i)       +        dedx(ipl,ist) = vect(i)
1424           goto 20           goto 20
# Line 1418  C Line 1430  C
1430           if ( ipl.ge.1.and.ipl.le.11 ) then           if ( ipl.ge.1.and.ipl.le.11 ) then
1431              do j = 1,16              do j = 1,16
1432                 i = i + 1                 i = i + 1
1433                 if (i.gt.sup) RETURN                 if (i.gt.sup.or.i.gt.120000) RETURN
1434                 ist = j + 16 * (ipr - 1)                 ist = j + 16 * (ipr - 1)
1435                 if (ist.ge.1.and.ist.le.96) dedx(ipl,ist) = vect(i)                 if (ist.ge.1.and.ist.le.96) dedx(ipl,ist) = vect(i)
1436              enddo              enddo
1437           endif           endif
1438           i = i + 1           i = i + 1
1439           if (i.gt.sup) RETURN           if (i.gt.sup.or.i.gt.120000) RETURN
1440           goto 10           goto 10
1441  C          C        
1442        endif        endif
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.23