/[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.8 by mocchiut, Mon Jan 22 09:17:01 2007 UTC revision 1.12 by mocchiut, Tue Jan 23 11:52:26 2007 UTC
# Line 260  C Line 260  C
260                    TANX = BX                    TANX = BX
261                 ELSE                 ELSE
262                    VARCFIT(2) = CHI2Y                    VARCFIT(2) = CHI2Y
263                    IMPY = 10. * ( AY + 12.2 )                    IMPY = 10. * ( AY + 12.1 )
264                    TANY = BY                    TANY = BY
265                 ENDIF                 ENDIF
266  C  C
# Line 275  C?????                  IF (M.EQ.2) NN = Line 275  C?????                  IF (M.EQ.2) NN =
275  C      C    
276                    IF (M.EQ.1) THEN                    IF (M.EQ.1) THEN
277                       DISTX = PIANO(I) - 5.81                       DISTX = PIANO(I) - 5.81
278                       Y(I) = DISTX * TANX + CX                       Y(I) = DISTX * TANX + IMPX
279    c                     CBAR(M,I) = Y(I)
280                       BAR(M,I) = Y(I)                       BAR(M,I) = Y(I)
281                       CBAR(M,I) = Y(I)                               CBAR(M,I) = (Y(I)-XALIG)/10.        
282                       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))
283  C      C    
284                    ELSE                    ELSE
285                       DISTY = PIANO(I)                                       DISTY = PIANO(I)                
286                       YY(I) = DISTY * TANY + CY                       YY(I) = DISTY * TANY + IMPY
287    c                     CBAR(M,I) = YY(I)
288                       BAR(M,I) = YY(I)                       BAR(M,I) = YY(I)
289                       CBAR(M,I) = YY(I)                       CBAR(M,I) = (YY(I)-YALIG)/10.
290                       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))
291  C      C    
292                    ENDIF                    ENDIF
# Line 612  c            ISUPY = NNY + 8 Line 614  c            ISUPY = NNY + 8
614           ENDIF           ENDIF
615        ENDDO        ENDDO
616  C  C
       EINF = EMIN  
       ESUP = 150.  
 C  
       IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN  
          EINF = 50.  
          ESUP = 1500.  
       ENDIF  
617  C  C
618  C     CALCULATE PLANETOT AND QMEAN  C     CALCULATE PLANETOT AND QMEAN
619  C  C
# Line 629  C Line 624  C
624        NPIANI = 5        NPIANI = 5
625        QMEAN = 0.        QMEAN = 0.
626        INDEX = 0        INDEX = 0
627        CALL ELIO(RPIANO,NPIANI,QMEAN,NTOT,INDEX)  C
628        PLANETOT = RPIANO(1) + RPIANO(2)          IF (TRIGTY.GE.2.AND.HZN.NE.0) THEN
629             EINF = 50.
630             ESUP = 15000.
631             CALL NUCLEI(RPIANO,NPIANI,QMEAN,NTOT,INDEX)
632             PLANETOT = RPIANO(1) + RPIANO(2)  
633          ELSE
634             EINF = EMIN
635             ESUP = 15000.
636             CALL ELIO(RPIANO,NPIANI,QMEAN,NTOT,INDEX)
637             PLANETOT = RPIANO(1) + RPIANO(2)  
638          ENDIF
639  C  C
640   50   CONTINUE   50   CONTINUE
641  C  C
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.23