/[PAMELA software]/calo/unpacking/qlook.for
ViewVC logotype

Contents of /calo/unpacking/qlook.for

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1.1.1 - (show annotations) (download) (vendor branch)
Mon Dec 5 16:23:21 2005 UTC (19 years ago) by mocchiut
Branch: MAIN, unpacking
CVS Tags: start, v1r00, HEAD
Changes since 1.1: +0 -0 lines
Error occurred while calculating annotation data.
Imported sources

1 REAL FUNCTION qlook()
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 ver(4,22)
17 integer l, bit, bi, i, k, se
18 INTEGER st4, shit(4,11), shit2, shit3
19 *
20 *-- Enter user code here
21 *
22 do l=1,4
23 i = l
24 st4 = 0
25 c st4 = IAND(stwerr(l) ! ,'00FF'x)
26 st4 = stwerr(l) ! ,'00FF'x)
27 if (st4.ne.0) then
28 do bit=0, 17
29 bi = ibits(st4,bit,1)
30 if (bit.eq.0.and.bi.ne.0) ver(i,7)=ver(i,7)+1
31 if (bit.eq.1.and.bi.ne.0) ver(i,6)=ver(i,6)+1
32 if (bit.eq.2.and.bi.ne.0) ver(i,5)=ver(i,5)+1
33 if (bit.eq.3.and.bi.ne.0) ver(i,4)=ver(i,4)+1
34 if (bit.eq.4.and.bi.ne.0) ver(i,3)=ver(i,3)+1
35 if (bit.eq.5.and.bi.ne.0) ver(i,2)=ver(i,2)+1
36 if (bit.eq.6.and.bi.ne.0) ver(i,1)=ver(i,1)+1
37 if (bit.eq.16.and.bi.ne.0) ver(i,17)=ver(i,17)+1
38 if (bit.eq.17.and.bi.ne.0) ver(i,18)=ver(i,18)+1
39 enddo
40 endif
41 if (perror(l).eq.128.) ver(i,8)=ver(i,8)+1
42 if (perror(l).eq.129.) ver(i,9)=ver(i,9)+1
43 if (perror(l).eq.130.) ver(i,10)=ver(i,10)+1
44 if (perror(l).eq.131.) ver(i,11)=ver(i,11)+1
45 if (perror(l).eq.132.) ver(i,12)=ver(i,12)+1
46 if (perror(l).eq.133.) ver(i,13)=ver(i,13)+1
47 if (perror(l).eq.134.) ver(i,14)=ver(i,14)+1
48 if (perror(l).eq.135.) ver(i,15)=ver(i,15)+1
49 if (perror(l).eq.136.) ver(i,16)=ver(i,16)+1
50 c if (perror(l).eq.137.) ver(i,17)=ver(i,17)+1
51 c if (perror(l).eq.138.) ver(i,18)=ver(i,18)+1
52 if (perror(l).eq.139.) ver(i,19)=ver(i,19)+1
53 if (perror(l).eq.140.) ver(i,20)=ver(i,20)+1
54 if (perror(l).eq.141.) ver(i,21)=ver(i,21)+1
55 enddo
56 do l = 1, 4
57 do i = 1, 11
58 shit(l,i) = 0
59 enddo
60 enddo
61 do l = 1, 2
62 do i = 1, 22
63 if (l.eq.1.and.mod(i,2).eq.0) then
64 se = 1
65 pl = i / 2
66 endif
67 if (l.eq.1.and.mod(i,2).ne.0) then
68 se = 2
69 pl = 1 + (i - 1)/2
70 endif
71 if (l.eq.2.and.mod(i,2).eq.0) then
72 se = 3
73 pl = i / 2
74 endif
75 if (l.eq.2.and.mod(i,2).ne.0) then
76 se = 4
77 pl = 1 + (i - 1)/2
78 endif
79 do k = 1, 96
80 if (dexyc(l,i,k).gt.0.) then
81 shit(se,pl) = shit(se,pl) + 1
82 ELSEIF (dexyc(l,i,k).lt.0.) then
83 ver(se,22) = ver(se,22) + 1
84 endif
85 enddo
86 enddo
87 enddo
88 shit3 = 0
89 do i = 100, 400, 100
90 shit2 = 0
91 do l = 1, 11
92 c print *,'i+l ',i+l
93 c print *,'shit ',float(shit(i/100,l))
94 call hf1(i+l+1000,float(shit(i/100,l)),1.)
95 shit2 = shit2 + shit(i/100,l)
96 enddo
97 call hf1(i+1050,float(shit2),1.)
98 shit3 = shit3 + shit2
99 enddo
100 call hf1(1500,float(shit3),1.)
101 C
102 END

  ViewVC Help
Powered by ViewVC 1.1.23