| 1 |
REAL FUNCTION qlookraw() |
| 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, 6 |
| 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 |
enddo |
| 38 |
endif |
| 39 |
if (perror(l).eq.128.) ver(i,8)=ver(i,8)+1 |
| 40 |
if (perror(l).eq.129.) ver(i,9)=ver(i,9)+1 |
| 41 |
if (perror(l).eq.130.) ver(i,10)=ver(i,10)+1 |
| 42 |
if (perror(l).eq.131.) ver(i,11)=ver(i,11)+1 |
| 43 |
if (perror(l).eq.132.) ver(i,12)=ver(i,12)+1 |
| 44 |
if (perror(l).eq.133.) ver(i,13)=ver(i,13)+1 |
| 45 |
if (perror(l).eq.134.) ver(i,14)=ver(i,14)+1 |
| 46 |
if (perror(l).eq.135.) ver(i,15)=ver(i,15)+1 |
| 47 |
if (perror(l).eq.136.) ver(i,16)=ver(i,16)+1 |
| 48 |
if (perror(l).eq.137.) ver(i,17)=ver(i,17)+1 |
| 49 |
if (perror(l).eq.138.) ver(i,18)=ver(i,18)+1 |
| 50 |
if (perror(l).eq.139.) ver(i,19)=ver(i,19)+1 |
| 51 |
if (perror(l).eq.140.) ver(i,20)=ver(i,20)+1 |
| 52 |
if (perror(l).eq.141.) ver(i,21)=ver(i,21)+1 |
| 53 |
enddo |
| 54 |
do l = 1, 4 |
| 55 |
do i = 1, 11 |
| 56 |
shit(l,i) = 0 |
| 57 |
enddo |
| 58 |
enddo |
| 59 |
do l = 1, 2 |
| 60 |
do i = 1, 22 |
| 61 |
if (l.eq.1.and.mod(i,2).eq.0) then |
| 62 |
se = 1 |
| 63 |
pl = i / 2 |
| 64 |
endif |
| 65 |
if (l.eq.1.and.mod(i,2).ne.0) then |
| 66 |
se = 2 |
| 67 |
pl = 1 + (i - 1)/2 |
| 68 |
endif |
| 69 |
if (l.eq.2.and.mod(i,2).eq.0) then |
| 70 |
se = 3 |
| 71 |
pl = i / 2 |
| 72 |
endif |
| 73 |
if (l.eq.2.and.mod(i,2).ne.0) then |
| 74 |
se = 4 |
| 75 |
pl = 1 + (i - 1)/2 |
| 76 |
endif |
| 77 |
do k = 1, 96 |
| 78 |
if (dexy(l,i,k).gt.0.and.dexy(l,i,k).lt.10000.) then |
| 79 |
shit(se,pl) = shit(se,pl) + 1 |
| 80 |
c print *,'ou!',se ,pl |
| 81 |
ELSEIF (dexy(l,i,k).lt.0.) then |
| 82 |
ver(se,22) = ver(se,22) + 1 |
| 83 |
endif |
| 84 |
enddo |
| 85 |
enddo |
| 86 |
enddo |
| 87 |
shit3 = 0 |
| 88 |
do i = 100, 400, 100 |
| 89 |
shit2 = 0 |
| 90 |
do l = 1, 11 |
| 91 |
c print *,'i+l ',i+l |
| 92 |
c print *,'shit ',float(shit(i/100,l)) |
| 93 |
call hf1(i+l+1000,float(shit(i/100,l)),1.) |
| 94 |
shit2 = shit2 + shit(i/100,l) |
| 95 |
enddo |
| 96 |
call hf1(i+1050,float(shit2),1.) |
| 97 |
shit3 = shit3 + shit2 |
| 98 |
enddo |
| 99 |
call hf1(1500,float(shit3),1.) |
| 100 |
C |
| 101 |
END |
Parent Directory
| 
Revision Log