stopp.f

计算能量损失,重离子经过介质后的能量沉积LET,射程

*$declare
*$debug


	SUBROUTINE STOPP(NT)
	integer	NT

	integer	i,j

	external	RANGE,DEDX
	real		RANGE,DEDX

****************	real	sum
**********************	real	z1,a1

	real	DIADJ(10)
	COMMON	/GEN/DIADJ

	include	'stpwtbl.i'

****** SUBROUTINE FOR CALCULATING ENERGY-RANGE AND ENERGY-DEDX **************
****** TABLES FOR THE PRIMARY PARTICLES *************************************

*		Original code	L. Sihver
*		1995-10-31	H. Tomura

      AIADJ=0.0
      MAT=TGNAME(NT)
      M1=TG1(1,NT)
      RHO=TG1(2,NT)
      AIADJ=TG1(3,NT)
*??????      NELECT=0.0
**************      SUM=0.0
      DO 10 I=1,M1
          SZ(I)=TG2(1,I,NT)
          SA(I)=TG2(2,I,NT)
          FM(I)=TG2(3,I,NT)
          DIADJ(I)=TG2(4,I,NT)
***************         SUM=SUM+FM(I)*SA(I)
  10  CONTINUE
*************************      Z1=PZ
**************************      A1=PA
      CALL MATPAR(M1,SZ,SA,FM)
      IF(DIADJ(1).GT.0.0) CALL MAT2

        ER1(22)=10.
        ER1(23)=10.4
        I=23
      DO 408 J=12,50,2
        I=I+1
  408   ER1(I)=J
      DO 409 J=55,200,5
        I=I+1
  409   ER1(I)=J
      DO 410 J=210,1000,10
        I=I+1
  410   ER1(I)=J
        NR=I
        ER1(1)=.01
        ER1(2)=.015
      DO 519 I=2,9
  519   ER1(I+1)=.01*FLOAT(I)
        ER1(11)=.1
        ER1(12)=.15
        J=12
      DO 524 I=2,9
        J=J+1
  524   ER1(J)=0.1*FLOAT(I)
        J=20
      DO 525 I=2,20
        J=J+1
  525   ER1(J)=0.5*FLOAT(I)
        ER1(40)=10.4
        I=40
      DO 504 J=12,50,2
        I=I+1
  504   ER1(I)=FLOAT(J)
      DO 505 J=55,200,5
        I=I+1
  505   ER1(I)=FLOAT(J)
      DO 506 J=210,1000,10
        I=I+1
  506   ER1(I)=FLOAT(J)
        NR=I

      CALL STNEW

      DO 30 I=1,170
***********	write(8,*) i,er1(i)
         RR(I)=RANGE(pz,pa,ER1(I))
***********	write(8,*) rr(i)
	 DEDX1R(I)=DEDX(pZ,pA,ER1(I))
***********	write(8,*) dedx1r(i)
 30   CONTINUE

	RETURN
	END