| 1 |
subroutine ptlevel1() |
| 2 |
|
| 3 |
c LOGICAL CHAIN |
| 4 |
c CHARACTER*128 CFILE |
| 5 |
c INTEGER IDNEVT,NCHEVT,ICHEVT |
| 6 |
c REAL OBS(13) |
| 7 |
* |
| 8 |
c COMMON /PAWIDN/ IDNEVT,OBS |
| 9 |
c COMMON /PAWCHN/ CHAIN, NCHEVT, ICHEVT |
| 10 |
c COMMON /PAWCHC/ CFILE |
| 11 |
* |
| 12 |
*-- Ntuple Variable Declarations |
| 13 |
* |
| 14 |
REAL DEDX(500),CLSIGNAL(8500),CNEV(12,24) |
| 15 |
INTEGER NEV1,NCLSTR1,VIEW(500),LADDER(500),MAXS(500),MULT(500) |
| 16 |
& ,INDSTART(500),INDMAX(500),TOTCLLENGTH |
| 17 |
LOGICAL GOOD1 |
| 18 |
* |
| 19 |
integer ncl(2) |
| 20 |
real clp(3,500) |
| 21 |
real cls(2,500) |
| 22 |
real clz(2,500) |
| 23 |
real sigma_matrix(12,3072) |
| 24 |
c |
| 25 |
Integer pkt_type1; |
| 26 |
Integer pkt_num1; |
| 27 |
Integer obt1; |
| 28 |
Integer which_calib1; |
| 29 |
|
| 30 |
c |
| 31 |
c parameter(NWPAWC=8) |
| 32 |
c common/pawc/hmem(NWPAWC) |
| 33 |
c save/pawc/ |
| 34 |
c common/quest/iquest(1) |
| 35 |
c save/quest/ |
| 36 |
c |
| 37 |
c |
| 38 |
|
| 39 |
COMMON /PWCR4/ GOOD1,NEV1, |
| 40 |
& pkt_type1,pkt_num1,obt1,which_calib1, |
| 41 |
& NCLSTR1,VIEW,LADDER,MAXS,MULT,DEDX |
| 42 |
& ,INDSTART,INDMAX,TOTCLLENGTH,CLSIGNAL,CNEV, |
| 43 |
& ncl,clp,cls,clz,sigma_matrix |
| 44 |
SAVE /PWCR4/ |
| 45 |
|
| 46 |
* |
| 47 |
|
| 48 |
*-- Enter user code here |
| 49 |
* |
| 50 |
c integer ,nstripsview(12) |
| 51 |
integer ind,ic,cllength,lows,highs,nstrips,is |
| 52 |
c integer sflag(12,3072) !to avoid putting more than one time a |
| 53 |
c strip in a histogram when a cluster is more than one time contained in |
| 54 |
c clsignal array (usually these are the double peaked clusters) |
| 55 |
real sign,noise,clinclcut |
| 56 |
|
| 57 |
real zlevel(12) |
| 58 |
data zlevel/22.28,22.28,13.38,13.38,4.48,4.48, |
| 59 |
+ -4.42,-4.42,-13.32,-13.32,-22.23,-22.23/ |
| 60 |
|
| 61 |
parameter(incut=5) |
| 62 |
parameter(signcut=50.) |
| 63 |
|
| 64 |
c parameter(traslys=2151) |
| 65 |
parameter(traslys=2115) !INVERSIONE????????????? |
| 66 |
|
| 67 |
c call hlimit(nwpawc) |
| 68 |
|
| 69 |
|
| 70 |
* *********** |
| 71 |
|
| 72 |
c PRINT *,'TRACKER: evento ',nev1 |
| 73 |
c event_num_trk(1) = nev1 |
| 74 |
c event_num_trk(1) = mod(nev1,1000000) |
| 75 |
c |
| 76 |
c do i = 1, 12 |
| 77 |
c do j = 1, 3072 |
| 78 |
c if ( sigma_matrix(i,j).ne.0. ) |
| 79 |
c & print *,'sigma(',i,',',j,') = ',sigma_matrix(i,j) |
| 80 |
c enddo |
| 81 |
c enddo |
| 82 |
c print *,' nclstr = ',nclstr1 |
| 83 |
c do i = 1, nclstr1 |
| 84 |
c print *,' view(',i,') ',view(i) |
| 85 |
c print *,' ladder(',i,') ',ladder(i) |
| 86 |
c enddo |
| 87 |
c |
| 88 |
|
| 89 |
ind=0 |
| 90 |
|
| 91 |
nclx = 0 |
| 92 |
ncly = 0 |
| 93 |
|
| 94 |
c print*,'acidenti a me ',event_num_trk(1),nclstr1 |
| 95 |
|
| 96 |
do ic=1,nclstr1 |
| 97 |
|
| 98 |
if(ic.ne.nclstr1) then |
| 99 |
cllength=indstart(ic+1)-indstart(ic) |
| 100 |
else |
| 101 |
cllength=(totCLlength+1)-indstart(ic) |
| 102 |
endif |
| 103 |
|
| 104 |
lows=maxs(ic)-(indmax(ic)-indstart(ic)) |
| 105 |
highs=lows+cllength-1 |
| 106 |
|
| 107 |
sign=0. |
| 108 |
noise=0. |
| 109 |
nstrips=0 |
| 110 |
posx = 0. |
| 111 |
posy1 = 0. |
| 112 |
posy2 = 0. |
| 113 |
posz = 0. |
| 114 |
pposx = 0. |
| 115 |
pposy1 = 0. |
| 116 |
pposy2 = 0. |
| 117 |
pposz = 0. |
| 118 |
|
| 119 |
C print *,'qua ',ic |
| 120 |
do is=lows,highs |
| 121 |
ind=ind+1 |
| 122 |
C print *,'qui view(ic) ',view(ic),' is ',is |
| 123 |
C print *,'qui clsignal(ind) ',clsignal(ind),' ind ',is |
| 124 |
clinclcut=incut*sigma_matrix(view(ic),is) |
| 125 |
if (clsignal(ind).gt.clinclcut) then |
| 126 |
C print *,'qua ci arrivo' |
| 127 |
sign=sign+clsignal(ind) |
| 128 |
noise=noise+sigma_matrix(view(ic),is) |
| 129 |
nstrips=nstrips+1 |
| 130 |
posz = zlevel(view(ic)) |
| 131 |
C print *,'qua anche sign ',sign,' noise ',noise |
| 132 |
if (mod(view(ic),2).eq.1) then !!! Y Y Y Y |
| 133 |
if (view(ic).ne.11) then |
| 134 |
str=float(mod(is-1,1024)+1) !Y |
| 135 |
posy1 = -6.9541 + (str-1)*0.00665 |
| 136 |
posy2 = 0.1 + (str-1)*0.00665 |
| 137 |
else !??????????????????INVERSIONE?????????? |
| 138 |
str=float(mod(is-1,1024)+1) !Y11 |
| 139 |
c str=traslys-str |
| 140 |
str=1024-str + 67 |
| 141 |
posy1 = -6.9541 + (str-1)*0.00665 !-0.5 |
| 142 |
posy2 = 0.1 + (str-1)*0.00665 !-0.5 |
| 143 |
endif |
| 144 |
pposy1 = pposy1 + posy1*clsignal(ind) |
| 145 |
pposy2 = pposy2 + posy2*clsignal(ind) |
| 146 |
else |
| 147 |
str=float(is) !X |
| 148 |
srt = mod(is-1,1024)+1 |
| 149 |
lad = int((is-1)/1024) + 1 |
| 150 |
if(lad.eq.1) posx=-7.8798+(srt-1)*0.005 |
| 151 |
if(lad.eq.2) posx=-2.5575+(srt-1)*0.005 |
| 152 |
if(lad.eq.3) posx= 2.7648+(srt-1)*0.005 |
| 153 |
pposx = pposx + posx*clsignal(ind) |
| 154 |
endif |
| 155 |
endif |
| 156 |
enddo |
| 157 |
if(mod(view(ic),2).eq.1)then |
| 158 |
C print *,'qui ci sono e pposy1 ',pposy1,' pposy2 ',pposy2,' sign ',sign |
| 159 |
ncly = ncly + 1 |
| 160 |
clp(2,ncly) = pposy1/sign |
| 161 |
clp(3,ncly) = pposy2/sign |
| 162 |
clz(2,ncly) = posz |
| 163 |
cls(2,ncly) = sign |
| 164 |
else |
| 165 |
nclx = nclx + 1 |
| 166 |
clp(1,nclx) = pposx/sign |
| 167 |
clz(1,nclx) = posz |
| 168 |
cls(1,nclx) = sign |
| 169 |
c call hfill(99,posx,posz,sign) |
| 170 |
endif |
| 171 |
enddo |
| 172 |
c do i = 1, 19 |
| 173 |
c print *,'cls(2,',i,') = ',cls(2,i) |
| 174 |
c print *,'cls(1,',i,') = ',cls(1,i) |
| 175 |
c enddo |
| 176 |
ncl(1) = nclx |
| 177 |
ncl(2) = ncly |
| 178 |
c project_trk_level1 = 1 |
| 179 |
C print *,'****' |
| 180 |
return |
| 181 |
end |
Parent Directory
| 
Revision Log