/[PAMELA software]/yoda/techmodel/forroutines/calorimeter/calpulse.for
ViewVC logotype

Diff of /yoda/techmodel/forroutines/calorimeter/calpulse.for

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

revision 2.0 by kusanagi, Tue Sep 21 20:51:10 2004 UTC revision 2.1 by kusanagi, Sun Oct 17 12:28:43 2004 UTC
# Line 1  Line 1 
1  C  C
2  C  Written by Mirko Boezio and Emiliano Mocchiutti  C  Written by Mirko Boezio and Emiliano Mocchiutti
3  C  C
4  C     * Version: 2.18.2 *  C     * Version: 2.18.3 *
5  C  C
6  C     (see calunpack for changelog)  C     (see calunpack for changelog)
7  C  C
# Line 66  C Line 66  C
66    
67        save / calib /        save / calib /
68    
69        COMMON /VARIE/ error, CONTR, e2, dump        COMMON /VARIE/ dump, CONTR
70        SAVE /VARIE/        SAVE /VARIE/
71            
72    
# Line 88  c Line 88  c
88        endif        endif
89  c  c
90        if (lunga.gt.60000.and.dump.gt.0) then        if (lunga.gt.60000.and.dump.gt.0) then
91  c         print *,'Calorimeter WARNING: more than 30000 words!'           print *,'Calorimeter WARNING: more than 30000 words!'
92           lunga = 60000           lunga = 60000
93        endif        endif
94  c  c
# Line 167  c     if so go out of this loop and go o Line 167  c     if so go out of this loop and go o
167                 m = ic                 m = ic
168                 call riempi(m,lunga,lleng,lung,pari,vect,vecta)                 call riempi(m,lunga,lleng,lung,pari,vect,vecta)
169                 icb = 1                 icb = 1
170                 if (st2b.ne.0) then  c               if (st2b.ne.0) then
171                    E2(contr) = vect(icb)                    E2(contr) = vect(icb)
172                 endif  c               endif
173                 goto 9                 goto 9
174              else              else
175  c    if not, is it one of the next sections? did we miss a section?  c    if not, is it one of the next sections? did we miss a section?
# Line 194  c Line 194  c
194                 m = ic                 m = ic
195                 call riempi(m,lunga,lleng,lung,pari,vect,vecta)                 call riempi(m,lunga,lleng,lung,pari,vect,vecta)
196                 icb = 1                 icb = 1
197                 if (st2b.ne.0) then  c               if (st2b.ne.0) then
198                    E2(contr) = vect(icb)                    E2(contr) = vect(icb)
199                 endif  c               endif
200                 goto 9                 goto 9
201              else              else
202                 if (st1b.eq.XE.or.st1b.eq.XO) then                 if (st1b.eq.XE.or.st1b.eq.XO) then
# Line 218  c Line 218  c
218                 m = ic                 m = ic
219                 call riempi(m,lunga,lleng,lung,pari,vect,vecta)                 call riempi(m,lunga,lleng,lung,pari,vect,vecta)
220                 icb = 1                 icb = 1
221                 if (st2b.ne.0) then  c               if (st2b.ne.0) then
222                    E2(contr) = vect(icb)                    E2(contr) = vect(icb)
223                 endif  c               endif
224                 goto 9                 goto 9
225              else              else
226                 if (st1b.eq.XO) then                 if (st1b.eq.XO) then
# Line 242  c Line 242  c
242                 m = ic                 m = ic
243                 call riempi(m,lunga,lleng,lung,pari,vect,vecta)                 call riempi(m,lunga,lleng,lung,pari,vect,vecta)
244                 icb = 1                 icb = 1
245                 if (st2b.ne.0) then  c               if (st2b.ne.0) then
246                    E2(contr) = vect(icb)                    E2(contr) = vect(icb)
247                 endif  c               endif
248              else              else
249  c     we should never arrive here (in case we run out of vector if section  c     we should never arrive here (in case we run out of vector if section
250  c     four is missing!)... however here it is in case of bugs!  c     four is missing!)... however here it is in case of bugs!
# Line 364  C Line 364  C
364              else              else
365                 bi = 0                 bi = 0
366              endif              endif
367              if (bi.eq.1) stwerr(l) = ibset(stwerr(l),bit)              if (l.ne.contr) bi = 0
368                if (bi.eq.1) then
369                   stwerr(l) = ibset(stwerr(l),bit)
370                else
371                   stwerr(l) = ibclr(stwerr(l),bit)
372                endif
373           enddo           enddo
374           perror(l) = float(error(l))           if (l.eq.contr) then
375                perror(l) = float(error(l))
376             else
377                perror(l) = 0.
378             endif
379        enddo        enddo
380  C  C
381        if (me.eq.0) iev2 = iev2 + 1        if (me.eq.0) iev2 = iev2 + 1
382        RETURN        RETURN
383        END        END
384    
385    
Legend:
Removed from v.2.0  
changed lines
  Added in v.2.1
  ViewVC Help
Powered by ViewVC 1.1.23