profile.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file profile.f

C
C
C       THIS FUNCTION IS TO CALCULATE THE NUCLEAR PROFILE FUNCTION
C       OF THE  COLLIDERING SYSTEM (IN UNITS OF 1/mb)
C
        FUNCTION  profile(XB)
        common/pact/bb,b1,phi,z1
        SAVE  /pact/
        common/hiparnt/hipr1(100),ihpr2(50),hint1(100),ihnt2(50)
        SAVE  /hiparnt/
        EXTERNAL flap, flap1, flap2
C
        bb=xb
        profile=1.0
        IF(ihnt2(1).GT.1 .AND. ihnt2(3).GT.1) THEN
           profile=float(ihnt2(1))*float(ihnt2(3))*0.1*
     &          gauss1(flap,0.0,hipr1(34),0.01)
        ELSE IF(ihnt2(1).EQ.1 .AND. ihnt2(3).GT.1) THEN
           profile=0.2*float(ihnt2(3))*
     &          gauss1(flap2,0.0,hipr1(35),0.001)
        ELSE IF(ihnt2(1).GT.1 .AND. ihnt2(3).EQ.1) THEN
           profile=0.2*float(ihnt2(1))*
     &          gauss1(flap1,0.0,hipr1(34),0.001)
        ENDIF
        RETURN
        END