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 #define BLANKET_SAVE
   9 #include "pact_hijing.inc"
  10 #include "hiparnt.inc"
  11         SAVE
  12         EXTERNAL FLAP, FLAP1, FLAP2
  13 C
  14         BB=XB
  15         PROFILE=1.0
  16         IF(IHNT2(1).GT.1 .AND. IHNT2(3).GT.1) THEN
  17            PROFILE=float(IHNT2(1))*float(IHNT2(3))*0.1*
  18      &          GAUSS1(FLAP,0.0,HIPR1(34),0.01)
  19         ELSE IF(IHNT2(1).EQ.1 .AND. IHNT2(3).GT.1) THEN
  20            PROFILE=0.2*float(IHNT2(3))*
  21      &          GAUSS1(FLAP2,0.0,HIPR1(35),0.001)
  22         ELSE IF(IHNT2(1).GT.1 .AND. IHNT2(3).EQ.1) THEN
  23            PROFILE=0.2*float(IHNT2(1))*
  24      &          GAUSS1(FLAP1,0.0,HIPR1(34),0.001)
  25         ENDIF
  26         RETURN
  27         END