ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
profile.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file profile.f
1 C
2 C
3 C THIS FUNCTION IS TO CALCULATE THE NUCLEAR PROFILE FUNCTION
4 C OF THE COLLIDERING SYSTEM (IN UNITS OF 1/mb)
5 C
6  FUNCTION profile(XB)
7  common/pact/bb,b1,phi,z1
8  SAVE /pact/
9  common/hiparnt/hipr1(100),ihpr2(50),hint1(100),ihnt2(50)
10  SAVE /hiparnt/
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