File: EX2FPE.FT of Tape: Various/ETH/eth11-3
(Source file text)
DIMENSION INPUT(1024),OUTPUT(10),EMU(4),SIGMA(4),SIZE(4)
DIMENSION ITABLE(68),VTYPE(2,2)
DATA VTYPE/' VA','LLEY','BASE','LINE'/
DATA EMU/20.,70.,600.,1000./
DATA SIGMA/20.,10.,200.,100./
DATA SIZE/950.,400.,300.,200./
DATA ITABLE/1,2,3,1,1,63*0/
DATA INLAST,INPTR,IDIMO,NPEAKS/1024,0,1,0/
L=0
TYPE 1000
1000 FORMAT(1H1,'PEAK NO.',8X,'AREA',4X,'P HEIGHT',6X,'P TIME',4X,
A 'L HEIGHT',6X,'L TIME',/,11X,'HALF WIDTH',4X,'T HEIGHT',6X,
B 'T TIME',8X,'TYPE',8X,'RATE'//)
DO 1 I=1,1024
A=0.
X=I
DO 2 J=1,4
2 A=A+SIZE(J)*EXP(-.5*((X-EMU(J))/SIGMA(J))**2)
1 INPUT(I)=A
3 CALL PEAK(ITABLE,INPUT,INLAST,INPTR,OUTPUT,IDIMO,NPEAKS)
IF(INPTR.LT.0.AND.NPEAKS.EQ.0) STOP
L=L+1
KK=OUTPUT(9)+1
4 TYPE 2000,(L,(OUTPUT(I),I=1,8),(VTYPE(K,KK),K=1,2),OUTPUT(10))
2000 FORMAT(I9,5F12.0,/,9X,3F12.0,4X,2A4,F12.0)
IF(INPTR.GE.0) GO TO 3
STOP
END