hijing/src/profile.f

0001 C
0002 C
0003 C       THIS FUNCTION IS TO CALCULATE THE NUCLEAR PROFILE FUNCTION
0004 C       OF THE  COLLIDERING SYSTEM (IN UNITS OF 1/mb)
0005 C
0006         FUNCTION  PROFILE(XB)
0007         COMMON/PACT/BB,B1,PHI,Z1
0008         SAVE  /PACT/
0009         COMMON/HIPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)
0010         SAVE  /HIPARNT/
0011         EXTERNAL FLAP, FLAP1, FLAP2
0012 C
0013         BB=XB
0014         PROFILE=1.0
0015         IF(IHNT2(1).GT.1 .AND. IHNT2(3).GT.1) THEN
0016            PROFILE=float(IHNT2(1))*float(IHNT2(3))*0.1*
0017      &          GAUSS1(FLAP,0.0,HIPR1(34),0.01)
0018         ELSE IF(IHNT2(1).EQ.1 .AND. IHNT2(3).GT.1) THEN
0019            PROFILE=0.2*float(IHNT2(3))*
0020      &          GAUSS1(FLAP2,0.0,HIPR1(35),0.001)
0021         ELSE IF(IHNT2(1).GT.1 .AND. IHNT2(3).EQ.1) THEN
0022            PROFILE=0.2*float(IHNT2(1))*
0023      &          GAUSS1(FLAP1,0.0,HIPR1(34),0.001)
0024         ENDIF
0025         RETURN
0026         END