HIJING/hijing1_36/profile.F

   1 * $Id$
   2 C
   3 C
   4 C       THIS FUNCTION IS TO CALCULATE THE NUCLEAR PROFILE FUNCTION
   5 C       OF THE  COLLIDERING SYSTEM (IN UNITS OF 1/mb)
   6 C
   7         FUNCTION  PROFILE(XB)
   8 #include "pact_hijing.inc"
   9 #include "hiparnt.inc"
  10         SAVE
  11         EXTERNAL FLAP, FLAP1, FLAP2
  12 C
  13         BB=XB
  14         PROFILE=1.0
  15         IF(IHNT2(1).GT.1 .AND. IHNT2(3).GT.1) THEN
  16            PROFILE=float(IHNT2(1))*float(IHNT2(3))*0.1*
  17      &          GAUSS1(FLAP,0.0,HIPR1(34),0.01)
  18         ELSE IF(IHNT2(1).EQ.1 .AND. IHNT2(3).GT.1) THEN
  19            PROFILE=0.2*float(IHNT2(3))*
  20      &          GAUSS1(FLAP2,0.0,HIPR1(35),0.001)
  21         ELSE IF(IHNT2(1).GT.1 .AND. IHNT2(3).EQ.1) THEN
  22            PROFILE=0.2*float(IHNT2(1))*
  23      &          GAUSS1(FLAP1,0.0,HIPR1(34),0.001)
  24         ENDIF
  25         RETURN
  26         END