/[PAMELA software]/DarthVader/CalorimeterLevel2/src/calol2tr.for
ViewVC logotype

Diff of /DarthVader/CalorimeterLevel2/src/calol2tr.for

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.10 by mocchiut, Mon Jan 22 10:02:27 2007 UTC revision 1.15 by mocchiut, Mon Apr 2 14:42:14 2007 UTC
# Line 142  c               print *,' CALORIMETER - Line 142  c               print *,' CALORIMETER -
142              TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4))              TY = TAN(ASIN(AL_P(3))) * SIN(AL_P(4))
143              DO I = 1, NPLA              DO I = 1, NPLA
144                 NN = 0                 NN = 0
145  C?????               IF (M.EQ.2) NN = 1                 IF (M.EQ.2) NN = 1
146                 IF (MOD(I,2).EQ.NN) THEN                 IF (MOD(I,2).EQ.NN) THEN
147                    SHIFT = +0.5                    SHIFT = +0.5
148                 ELSE                 ELSE
# Line 253  C Line 253  C
253  C      C    
254           DO M = 1,2           DO M = 1,2
255  C  C
256    c            print *,' ax ',ax,' ay ',ay
257    c            print *,' bx ',bx,' by ',by
258              IF (NPCFIT(M).GE.2) THEN              IF (NPCFIT(M).GE.2) THEN
259                 IF (M.EQ.1) THEN                 IF (M.EQ.1) THEN
260                    VARCFIT(1) = CHI2X                    VARCFIT(1) = CHI2X
261                    IMPX = 10. * ( AX + 12.1 )                    IMPX = AX ! PAMELA REF
262                    TANX = BX                    TANX = BX
263                 ELSE                 ELSE
264                    VARCFIT(2) = CHI2Y                    VARCFIT(2) = CHI2Y
265                    IMPY = 10. * ( AY + 12.1 )                    IMPY = AY ! PAMELA REF
266                    TANY = BY                    TANY = BY
267                 ENDIF                 ENDIF
268  C  C
269                 DO I = 1,NPLA                     DO I = 1,NPLA    
270                    NN = 0                    NN = 0
271  C?????                  IF (M.EQ.2) NN = 1                    IF (M.EQ.2) NN = 1
272                    IF (MOD(I,2).EQ.NN) THEN                    IF (MOD(I,2).EQ.NN) THEN
273                       SHIFT = +0.5                       SHIFT = +0.5
274                    ELSE                    ELSE
# Line 275  C?????                  IF (M.EQ.2) NN = Line 277  C?????                  IF (M.EQ.2) NN =
277  C      C    
278                    IF (M.EQ.1) THEN                    IF (M.EQ.1) THEN
279                       DISTX = PIANO(I) - 5.81                       DISTX = PIANO(I) - 5.81
280                       Y(I) = DISTX * TANX + IMPX                       Y(I) = (DISTX * TANX) +  AX - XALIG
281  c                     CBAR(M,I) = Y(I)  c                     CBAR(M,I) = Y(I)
282                       BAR(M,I) = Y(I)                       BAR(M,I) = Y(I)
283                       CBAR(M,I) = (Y(I) - XALIG)/10.                               CBAR(M,I) = (Y(I) + XALIG)/10.
284                       IF (I.EQ.22) MX=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22))                       IF (I.EQ.22) MX=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22))
285  C      C    
286                    ELSE                    ELSE
287                       DISTY = PIANO(I)                                       DISTY = PIANO(I)                
288                       YY(I) = DISTY * TANY + IMPY                       YY(I) = (DISTY * TANY) + AY - YALIG
289  c                     CBAR(M,I) = YY(I)  c                     CBAR(M,I) = YY(I)
290                       BAR(M,I) = YY(I)                       BAR(M,I) = YY(I)
291                       CBAR(M,I) = (-YALIG + YY(I))/10.                       CBAR(M,I) = (YY(I) + YALIG)/10.
292                       IF (I.EQ.22) MY=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22))                       IF (I.EQ.22) MY=ABS(Y(1)-Y(22))/ABS(ZIN(1)-ZIN(22))
293  C      C    
294                    ENDIF                    ENDIF
# Line 370  C Line 372  C
372        RNSS = 0.        RNSS = 0.
373        QTOTT = 0.        QTOTT = 0.
374  C  C
375        PPLANEMAX = 1.01*(LOG(ABS(RIG)/0.0081)-1.)        PPLANEMAX = 1.01*(LOG(ABS(RIG)/0.0081)-1.) / 0.74
376  C  C
377        IPLANE = INT(ANINT(PPLANEMAX)) + 5        IPLANE = INT(ANINT(PPLANEMAX)) + 5
378  C  C
# Line 399  C     8 STRIPS ARE 2.88 cm , A MOLIERE R Line 401  C     8 STRIPS ARE 2.88 cm , A MOLIERE R
401  C      C    
402        DO J = 1,IPLANE        DO J = 1,IPLANE
403           NNX = IBAR(1,J)           NNX = IBAR(1,J)
404             RNSS = 0.               ! BACO!!
405             QTOTT = 0.               ! BACO!!
406           IF (NNX.NE.-1) THEN           IF (NNX.NE.-1) THEN
407              IF (NNX.LT.9) NNX = 9              IF (NNX.LT.9) NNX = 9
408              IF (NNX.GT.88) NNX = 88              IF (NNX.GT.88) NNX = 88
# Line 614  c            ISUPY = NNY + 8 Line 618  c            ISUPY = NNY + 8
618           ENDIF           ENDIF
619        ENDDO        ENDDO
620  C  C
       EINF = EMIN  
       ESUP = 150.  
 C  
       IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN  
          EINF = 50.  
          ESUP = 1500.  
       ENDIF  
621  C  C
622  C     CALCULATE PLANETOT AND QMEAN  C     CALCULATE PLANETOT AND QMEAN
623  C  C
# Line 631  C Line 628  C
628        NPIANI = 5        NPIANI = 5
629        QMEAN = 0.        QMEAN = 0.
630        INDEX = 0        INDEX = 0
631        CALL ELIO(RPIANO,NPIANI,QMEAN,NTOT,INDEX)  C
632        PLANETOT = RPIANO(1) + RPIANO(2)          IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN
633             EINF = 50.
634             ESUP = 15000.
635             CALL NUCLEI(RPIANO,NPIANI,QMEAN,NTOT,INDEX)
636             PLANETOT = RPIANO(1) + RPIANO(2)  
637          ELSE
638             EINF = EMIN
639             ESUP = 15000.
640             CALL ELIO(RPIANO,NPIANI,QMEAN,NTOT,INDEX)
641             PLANETOT = RPIANO(1) + RPIANO(2)  
642          ENDIF
643  C  C
644   50   CONTINUE   50   CONTINUE
645  C  C
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.23