6 |
INCLUDE 'INTEST.TXT' |
INCLUDE 'INTEST.TXT' |
7 |
C |
C |
8 |
DOUBLE PRECISION al_p(5), |
DOUBLE PRECISION al_p(5), |
9 |
& xout(npla),yout(npla),zin(npla) |
& xout(nplav),yout(nplav),zin(nplav) |
10 |
C |
C |
11 |
REAL PIANO(22), VARFIT(2) |
REAL PIANO(NPLAV), VARFIT(2) |
12 |
REAL TX, TY, SHIFT |
REAL TX, TY, SHIFT |
13 |
REAL BAR(2,NPLA), DISTY |
REAL BAR(2,NPLAV), DISTY |
14 |
REAL DISTX, Y(NPLA), YY(NPLA) |
REAL DISTX, Y(NPLAV), YY(NPLAV) |
15 |
REAL RIG, PPLANEMAX, RMASS |
REAL RIG, PPLANEMAX, RMASS |
16 |
REAL RNSS, QTOTT, RQT, MX, MY |
REAL RNSS, QTOTT, RQT, MX, MY |
17 |
REAL CHECK, ENER, CX, CY |
REAL CHECK, ENER, CX, CY |
21 |
REAL ax,bx,eax,ebx,chi2x |
REAL ax,bx,eax,ebx,chi2x |
22 |
REAL ay,by,eay,eby,chi2y |
REAL ay,by,eay,eby,chi2y |
23 |
REAL parzen3, TMISD |
REAL parzen3, TMISD |
24 |
INTEGER Nfitx,Nfity |
INTEGER Nfitx,Nfity, MNPLA |
25 |
C |
C |
26 |
INTEGER INDEX, NTOT(2), NPIANI, GTR |
INTEGER INDEX, NTOT(2), NPIANI, GTR |
27 |
INTEGER j, m, i, IWPL(2), timpx, timpy, T, nn |
INTEGER j, m, i, IWPL(2), timpx, timpy, T, nn |
28 |
INTEGER IPLANE, NNX, NNY, INFX, INFY, ISUPX, ISUPY |
INTEGER IPLANE, NNX, NNY, INFX, INFY, ISUPX, ISUPY |
29 |
INTEGER IBAR(2,NPLA), NPFIT(2), CHTRACK,IWPLU |
INTEGER IBAR(2,NPLAV), NPFIT(2), CHTRACK,IWPLU |
30 |
INTEGER Iquest(100), ICONTROL5, nin, IFAIL |
INTEGER Iquest(100), ICONTROL5, nin, IFAIL |
31 |
C |
C |
32 |
PARAMETER (X01PL=0.74) |
PARAMETER (X01PL=0.74) |
66 |
C |
C |
67 |
C Begin ! |
C Begin ! |
68 |
C |
C |
69 |
|
c print *,' sono qui' |
70 |
CALOL2TR = 0; |
CALOL2TR = 0; |
71 |
NCORE = 0. |
NCORE = 0. |
72 |
QCORE = 0. |
QCORE = 0. |
92 |
QMEAN = 0. |
QMEAN = 0. |
93 |
C SELFTRIGGER = 0 |
C SELFTRIGGER = 0 |
94 |
C |
C |
95 |
C BEGIN WITH THE FISRT TRACK IF WE HAVE A TRACK FROM TRACKER |
C BEGIN WITH THE FIRST TRACK IF WE HAVE A TRACK FROM TRACKER |
96 |
C |
C |
97 |
T = 1 |
T = 1 |
98 |
C |
C |
99 |
10 CONTINUE |
c 10 CONTINUE |
100 |
|
CONTINUE |
101 |
C |
C |
102 |
IF (GOOD2.EQ.1) THEN |
IF (GOOD2.EQ.1) THEN |
103 |
C |
C |
104 |
CHTRACK = 0 |
CHTRACK = 0 |
105 |
C |
C |
106 |
CALL VZERO(IWPL,2) |
CALL VZERO(IWPL,2) |
107 |
CALL VZERO(BAR,2*NPLA) |
CALL VZERO(BAR,2*NPLAV) |
108 |
CALL VZERO(IBAR,2*NPLA) |
CALL VZERO(IBAR,2*NPLAV) |
109 |
CALL VZERO(TBAR,2*NPLA) |
CALL VZERO(TBAR,2*NPLAV) |
110 |
CALL VZERO(TIBAR,2*NPLA) |
CALL VZERO(TIBAR,2*NPLAV) |
111 |
CALL VZERO(Y,NPLA) |
CALL VZERO(Y,NPLAV) |
112 |
CALL VZERO(YY,NPLA) |
CALL VZERO(YY,NPLAV) |
113 |
CALL VZERO(XOUT,NPLA) |
CALL VZERO(XOUT,NPLAV) |
114 |
CALL VZERO(YOUT,NPLA) |
CALL VZERO(YOUT,NPLAV) |
115 |
do m = 1, 5 |
do m = 1, 5 |
116 |
al_p(m) = al_pp(t,m) |
al_p(m) = al_pp(t,m) |
117 |
|
c print *,' al_p(',m,') = ',al_p(m) |
118 |
enddo |
enddo |
119 |
if (al_p(5).eq.0.) THEN |
if (al_p(5).eq.0.) THEN |
120 |
PRINT *,' CALORIMETER - WARNING F77: track with R = 0, discarded' |
PRINT *,' CALORIMETER - WARNING F77: track with R = 0, discarded' |
127 |
YOUT(I) = 0. |
YOUT(I) = 0. |
128 |
IF (MOD(M,2).EQ.0) THEN |
IF (MOD(M,2).EQ.0) THEN |
129 |
DISTX = PIANO(I) + ZALIG |
DISTX = PIANO(I) + ZALIG |
130 |
|
c print *,'T Y PLANE I= ',I,' Z = ',DISTX |
131 |
ELSE |
ELSE |
132 |
DISTX = PIANO(I) - 5.81 + ZALIG |
DISTX = PIANO(I) - 5.81 + ZALIG |
133 |
|
c print *,'T X PLANE I= ',I,' Z = ',DISTX |
134 |
ENDIF |
ENDIF |
135 |
ZIN(I) = distx / 10. |
ZIN(I) = distx / 10. |
136 |
|
c print *,' ZIN(',I,') = ',ZIN(I) |
137 |
TBAR(M,I) = 0. |
TBAR(M,I) = 0. |
138 |
TIBAR(M,I) = 0 |
TIBAR(M,I) = 0 |
139 |
enddo |
enddo |
148 |
TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4)) |
TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4)) |
149 |
DO I = 1, NPLA |
DO I = 1, NPLA |
150 |
NN = 0 |
NN = 0 |
151 |
C????? IF (M.EQ.2) NN = 1 |
c IF (M.EQ.2) NN = 1 |
152 |
IF (MOD(I,2).EQ.NN) THEN |
IF (MOD(I,2).EQ.NN) THEN |
153 |
SHIFT = +0.5 |
IF (REVERSE.EQ.0) THEN |
154 |
|
SHIFT = -0.5 |
155 |
|
ELSE |
156 |
|
SHIFT = +0.5 |
157 |
|
ENDIF |
158 |
ELSE |
ELSE |
159 |
SHIFT = -0.5 |
IF (REVERSE.EQ.0) THEN |
160 |
|
SHIFT = +0.5 |
161 |
|
ELSE |
162 |
|
SHIFT = -0.5 |
163 |
|
ENDIF |
164 |
ENDIF |
ENDIF |
165 |
C |
C |
166 |
C CHECK IF XOUT OR YOUT ARE NaN |
C CHECK IF XOUT OR YOUT ARE NaN |
183 |
Y(I) = CX |
Y(I) = CX |
184 |
BAR(M,I) = Y(I) |
BAR(M,I) = Y(I) |
185 |
TBAR(M,I) = (Y(I) - XALIG)/10. |
TBAR(M,I) = (Y(I) - XALIG)/10. |
186 |
IF (I.EQ.22) MX=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22)) |
IF (I.EQ.NPLA) MX=ABS(Y(1)-Y(NPLA))/ |
187 |
|
& ABS(ZIN(1)-ZIN(NPLA)) |
188 |
ELSE |
ELSE |
189 |
YY(I) = CY |
YY(I) = CY |
190 |
BAR(M,I) = YY(I) |
BAR(M,I) = YY(I) |
191 |
TBAR(M,I) = (-YALIG + YY(I))/10. |
TBAR(M,I) = (-YALIG + YY(I))/10. |
192 |
IF (I.EQ.22) MY=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22)) |
IF (I.EQ.NPLA) MY=ABS(Y(1)-Y(NPLA))/ |
193 |
|
& ABS(ZIN(1)-ZIN(NPLA)) |
194 |
ENDIF |
ENDIF |
195 |
CALL LASTRISCIA(BAR(M,I),IBAR(M,I)) |
CALL LASTRISCIA(BAR(M,I),IBAR(M,I)) |
196 |
tibar(M,I) = ibar(m,i) |
tibar(M,I) = ibar(m,i) |
251 |
CALL VZERO(IWPL,2) |
CALL VZERO(IWPL,2) |
252 |
CALL VZERO(VARCFIT,2) |
CALL VZERO(VARCFIT,2) |
253 |
CALL VZERO(NPCFIT,2) |
CALL VZERO(NPCFIT,2) |
254 |
CALL VZERO(TBAR,2*NPLA) |
CALL VZERO(TBAR,2*NPLAV) |
255 |
CALL VZERO(TIBAR,2*NPLA) |
CALL VZERO(TIBAR,2*NPLAV) |
256 |
CALL VZERO(BAR,2*NPLA) |
CALL VZERO(BAR,2*NPLAV) |
257 |
CALL VZERO(IBAR,2*NPLA) |
CALL VZERO(IBAR,2*NPLAV) |
258 |
CALL VZERO(Y,NPLA) |
CALL VZERO(Y,NPLAV) |
259 |
CALL VZERO(YY,NPLA) |
CALL VZERO(YY,NPLAV) |
260 |
CALL VZERO(XOUT,NPLA) |
CALL VZERO(XOUT,NPLAV) |
261 |
CALL VZERO(YOUT,NPLA) |
CALL VZERO(YOUT,NPLAV) |
262 |
C |
C |
263 |
CALL SELFTRIG() |
CALL SELFTRIG() |
264 |
ELEN = PARZEN3 |
ELEN = PARZEN3 |
269 |
C |
C |
270 |
DO M = 1,2 |
DO M = 1,2 |
271 |
C |
C |
272 |
|
c print *,' ax ',ax,' ay ',ay |
273 |
|
c print *,' bx ',bx,' by ',by |
274 |
IF (NPCFIT(M).GE.2) THEN |
IF (NPCFIT(M).GE.2) THEN |
275 |
IF (M.EQ.1) THEN |
IF (M.EQ.1) THEN |
276 |
VARCFIT(1) = CHI2X |
VARCFIT(1) = CHI2X |
277 |
IMPX = 10. * ( AX + 12.1 ) |
IMPX = AX + BX * (ZALIG/10.) ! PAMELA REF |
278 |
TANX = BX |
TANX = BX |
279 |
ELSE |
ELSE |
280 |
VARCFIT(2) = CHI2Y |
VARCFIT(2) = CHI2Y |
281 |
IMPY = 10. * ( AY + 12.1 ) |
IMPY = AY + BY * (ZALIG/10.) ! PAMELA REF |
282 |
TANY = BY |
TANY = BY |
283 |
ENDIF |
ENDIF |
284 |
C |
C |
285 |
DO I = 1,NPLA |
DO I = 1,NPLA |
286 |
NN = 0 |
NN = 0 |
287 |
C????? IF (M.EQ.2) NN = 1 |
c IF (M.EQ.2) NN = 1 |
288 |
IF (MOD(I,2).EQ.NN) THEN |
IF (MOD(I,2).EQ.NN) THEN |
289 |
SHIFT = +0.5 |
IF (REVERSE.EQ.0) THEN |
290 |
|
SHIFT = -0.5 |
291 |
|
ELSE |
292 |
|
SHIFT = +0.5 |
293 |
|
ENDIF |
294 |
ELSE |
ELSE |
295 |
SHIFT = -0.5 |
IF (REVERSE.EQ.0) THEN |
296 |
|
SHIFT = +0.5 |
297 |
|
ELSE |
298 |
|
SHIFT = -0.5 |
299 |
|
ENDIF |
300 |
ENDIF |
ENDIF |
301 |
C |
C |
302 |
IF (M.EQ.1) THEN |
IF (M.EQ.1) THEN |
303 |
DISTX = PIANO(I) - 5.81 |
DISTX = PIANO(I) - 5.81 |
304 |
Y(I) = DISTX * TANX + IMPX |
Y(I) = (DISTX * TANX) + IMPX*10. !- XALIG |
305 |
c CBAR(M,I) = Y(I) |
c CBAR(M,I) = Y(I) |
306 |
BAR(M,I) = Y(I) |
BAR(M,I) = Y(I) + XALIG |
307 |
CBAR(M,I) = (Y(I) - XALIG)/10. |
CBAR(M,I) = Y(I) / 10. |
308 |
IF (I.EQ.22) MX=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22)) |
IF (I.EQ.NPLA) MX=ABS(Y(1)-Y(NPLA))/ |
309 |
|
& ABS(ZIN(1)-ZIN(NPLA)) |
310 |
C |
C |
311 |
ELSE |
ELSE |
312 |
DISTY = PIANO(I) |
DISTY = PIANO(I) |
313 |
YY(I) = DISTY * TANY + IMPY |
YY(I) = (DISTY * TANY) + IMPY*10. !- YALIG |
314 |
|
c print *,' I ',i,' YY ',YY(I) |
315 |
c CBAR(M,I) = YY(I) |
c CBAR(M,I) = YY(I) |
316 |
BAR(M,I) = YY(I) |
BAR(M,I) = YY(I) + YALIG |
317 |
CBAR(M,I) = (-YALIG + YY(I))/10. |
CBAR(M,I) = YY(I) / 10. |
318 |
IF (I.EQ.22) MY=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22)) |
IF (I.EQ.NPLA) MY=ABS(Y(1)-Y(NPLA))/ |
319 |
|
& ABS(ZIN(1)-ZIN(NPLA)) |
320 |
C |
C |
321 |
ENDIF |
ENDIF |
322 |
CALL LASTRISCIA(BAR(M,I),IBAR(M,I)) |
CALL LASTRISCIA(BAR(M,I),IBAR(M,I)) |
399 |
RNSS = 0. |
RNSS = 0. |
400 |
QTOTT = 0. |
QTOTT = 0. |
401 |
C |
C |
402 |
PPLANEMAX = 1.01*(LOG(ABS(RIG)/0.0081)-1.) |
PPLANEMAX = 1.01*(LOG(ABS(RIG)/0.0081)-1.) / 0.74 |
403 |
C |
C |
404 |
IPLANE = INT(ANINT(PPLANEMAX)) + 5 |
IPLANE = INT(ANINT(PPLANEMAX)) + 5 |
405 |
C |
C |
406 |
IF (IPLANE.GT.NPLA) IPLANE=NPLA |
IF (IPLANE.GT.NPLA) IPLANE=NPLA |
407 |
IF (IPLANE.LT.1) IPLANE = 1 |
IF (IPLANE.LT.1) IPLANE = 1 |
408 |
|
c print *,' calcolo...' |
409 |
C |
C |
410 |
C CALCULATE QLOW AND NLOW |
C CALCULATE QLOW AND NLOW |
411 |
C |
C |
429 |
C |
C |
430 |
DO J = 1,IPLANE |
DO J = 1,IPLANE |
431 |
NNX = IBAR(1,J) |
NNX = IBAR(1,J) |
432 |
|
RNSS = 0. ! BACO!! |
433 |
|
QTOTT = 0. ! BACO!! |
434 |
IF (NNX.NE.-1) THEN |
IF (NNX.NE.-1) THEN |
435 |
IF (NNX.LT.9) NNX = 9 |
IF (NNX.LT.9) NNX = 9 |
436 |
IF (NNX.GT.88) NNX = 88 |
IF (NNX.GT.88) NNX = 88 |
481 |
INFX = NNX - 8 |
INFX = NNX - 8 |
482 |
ISUPX = NNX + 8 |
ISUPX = NNX + 8 |
483 |
DO I = INFX,ISUPX |
DO I = INFX,ISUPX |
484 |
IF (DEXY(1,J,I).LT.EMIN) GO TO 710 |
IF (DEXY(1,J,I).GE.EMIN) THEN |
485 |
NCYL = NCYL + 1 |
NCYL = NCYL + 1 |
486 |
QCYL = QCYL + DEXY(1,J,I) |
QCYL = QCYL + DEXY(1,J,I) |
487 |
710 ENDDO |
ENDIF |
488 |
|
ENDDO |
489 |
ENDIF |
ENDIF |
490 |
NNY = IBAR(2,J) |
NNY = IBAR(2,J) |
491 |
IF (NNY.NE.-1) THEN |
IF (NNY.NE.-1) THEN |
494 |
INFY = NNY - 8 |
INFY = NNY - 8 |
495 |
ISUPY = NNY + 8 |
ISUPY = NNY + 8 |
496 |
DO I=INFY,ISUPY |
DO I=INFY,ISUPY |
497 |
IF (DEXY(2,J,I).LT.EMIN) GO TO 810 |
IF (DEXY(2,J,I).GE.EMIN) THEN |
498 |
NCYL = NCYL + 1 |
NCYL = NCYL + 1 |
499 |
QCYL = QCYL + DEXY(2,J,I) |
QCYL = QCYL + DEXY(2,J,I) |
500 |
810 ENDDO |
ENDIF |
501 |
|
ENDDO |
502 |
ENDIF |
ENDIF |
503 |
C |
C |
504 |
C QTR = DETECTED ENERGY AND NTR = NUMBER OF HIT STRIPS IN A CYLINDER oF |
C QTR = DETECTED ENERGY AND NTR = NUMBER OF HIT STRIPS IN A CYLINDER oF |
610 |
C |
C |
611 |
C CALCULATE NLAST AND QLAST |
C CALCULATE NLAST AND QLAST |
612 |
C |
C |
613 |
DO J = NPLA-4,NPLA |
MNPLA = NPLA -4 |
614 |
|
IF ( MNPLA .LT. 1 ) MNPLA = 1 |
615 |
|
DO J = MNPLA,NPLA |
616 |
NNX = IBAR(1,J) |
NNX = IBAR(1,J) |
617 |
IF (NNX.NE.-1) THEN |
IF (NNX.NE.-1) THEN |
618 |
IF (NNX.LT.5) NNX = 5 |
IF (NNX.LT.5) NNX = 5 |
650 |
ENDIF |
ENDIF |
651 |
ENDDO |
ENDDO |
652 |
C |
C |
|
EINF = EMIN |
|
|
ESUP = 15000. |
|
|
C |
|
|
IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN |
|
|
EINF = 50. |
|
|
ESUP = 15000. |
|
|
ENDIF |
|
653 |
C |
C |
654 |
C CALCULATE PLANETOT AND QMEAN |
C CALCULATE PLANETOT AND QMEAN |
655 |
C |
C |
660 |
NPIANI = 5 |
NPIANI = 5 |
661 |
QMEAN = 0. |
QMEAN = 0. |
662 |
INDEX = 0 |
INDEX = 0 |
663 |
CALL ELIO(RPIANO,NPIANI,QMEAN,NTOT,INDEX) |
C |
664 |
PLANETOT = RPIANO(1) + RPIANO(2) |
IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN |
665 |
|
EINF = 50. |
666 |
|
ESUP = 15000. |
667 |
|
CALL NUCLEI(RPIANO,NPIANI,QMEAN,NTOT,INDEX) |
668 |
|
PLANETOT = RPIANO(1) + RPIANO(2) |
669 |
|
ELSE |
670 |
|
EINF = EMIN |
671 |
|
ESUP = 15000. |
672 |
|
CALL ELIO(RPIANO,NPIANI,QMEAN,NTOT,INDEX) |
673 |
|
PLANETOT = RPIANO(1) + RPIANO(2) |
674 |
|
ENDIF |
675 |
C |
C |
676 |
50 CONTINUE |
50 CONTINUE |
677 |
C |
C |
678 |
|
c print *,' esco' |
679 |
RETURN |
RETURN |
680 |
END |
END |
681 |
|
|