calo/unpacking/calibc.for

      REAL FUNCTION calibc()
*********************************************************
*                                                       *
* This file was generated by HUWFUN.                    *
*                                                       *
*********************************************************
*
*     Ntuple Id:      2    
*     Ntuple Title:   Pamela data
*     Creation:       13/07/2004 11.39.53
*
*********************************************************
*
      INCLUDE 'comcheck.txt'
*
      vector pede(2,22,96)
      vector matmip(2,22,96)
      real bas, valore
      integer m, inf, sup
*
*--   Enter user code here
*     
      do i=1,2
         do j=1,22
c
c     commenta questo per la mip totale globale
c
c            if (i.eq.2.and.j.eq.1) then
               pre = 0
               do k=1,96
                  if (mod(k+15,16).eq.0) pre = pre + 1
c     print *,k, pre
                  if (dexyc(i,j,k).gt.0) then
                     calbase = 0.
                     
                     if (base(i,j,pre).eq.0.) then
                        inf = (pre - 1) * 16 + 1
                        sup = pre * 16 
C                     print *,'i s ',inf,sup,pre
                        do m = inf, sup
                           valore = valore + dexyc(i,j,m) 
                        enddo
                        bas = valore / 16.
c                     bas = base(i,j,pre)
                     else
                        bas = base(i,j,pre)
                     endif
                     
                    calbase=(dexyc(i,j,k)-bas-pede(i,j,k))/matmip(i,j,k)
                     call hf1(1600,calbase,1.)
                  endif
                  call hf1(1700,dexyc(i,j,k),1.)
                  call hf1(1800,dexy(i,j,k),1.)
                  call hf1(1900,base(i,j,pre),1.)
c               if (cal_rms(i,j,k).ne.0.)
c     &         print *,'rms ',cal_rms(i,j,k),' i=',i,' j=',j,' k=',k
               enddo
c
c     e commenta anche qua ovviamente!
c
c            endif
         enddo
      enddo      
c
      calib = 1.
*
      END
1	mocchiut	1.1	REAL FUNCTION calibc()
2			*********************************************************
3			* *
4			* This file was generated by HUWFUN. *
5			* *
6			*********************************************************
7			*
8			* Ntuple Id: 2
9			* Ntuple Title: Pamela data
10			* Creation: 13/07/2004 11.39.53
11			*
12			*********************************************************
13			*
14			INCLUDE 'comcheck.txt'
15			*
16			vector pede(2,22,96)
17			vector matmip(2,22,96)
18			real bas, valore
19			integer m, inf, sup
20			*
21			*-- Enter user code here
22			*
23			do i=1,2
24			do j=1,22
25			c
26			c commenta questo per la mip totale globale
27			c
28			c if (i.eq.2.and.j.eq.1) then
29			pre = 0
30			do k=1,96
31			if (mod(k+15,16).eq.0) pre = pre + 1
32			c print *,k, pre
33			if (dexyc(i,j,k).gt.0) then
34			calbase = 0.
35
36			if (base(i,j,pre).eq.0.) then
37			inf = (pre - 1) * 16 + 1
38			sup = pre * 16
39			C print *,'i s ',inf,sup,pre
40			do m = inf, sup
41			valore = valore + dexyc(i,j,m)
42			enddo
43			bas = valore / 16.
44			c bas = base(i,j,pre)
45			else
46			bas = base(i,j,pre)
47			endif
48
49			calbase=(dexyc(i,j,k)-bas-pede(i,j,k))/matmip(i,j,k)
50			call hf1(1600,calbase,1.)
51			endif
52			call hf1(1700,dexyc(i,j,k),1.)
53			call hf1(1800,dexy(i,j,k),1.)
54			call hf1(1900,base(i,j,pre),1.)
55			c if (cal_rms(i,j,k).ne.0.)
56			c & print *,'rms ',cal_rms(i,j,k),' i=',i,' j=',j,' k=',k
57			enddo
58			c
59			c e commenta anche qua ovviamente!
60			c
61			c endif
62			enddo
63			enddo
64			c
65			calib = 1.
66			*
67			END