RWGCF converted to native cmake
[u/mrichter/AliRoot.git] / THbtp / common_coulomb.inc
1 CCC   common/coulmb/
2 C
3 C     This common block includes the following data:
4 C
5 C          o input q mesh (GeV/c) 
6 C          o interpolated Coulomb correction for a specified
7 C            Coulomb source size (input as parameter Q0 when
8 C            switch_coulomb = 3) for both like and unlike pairs.
9 C        
10 C          The correction corresponds to John Cramer's tabulation
11 C          of Scott Pratt's Coulomb wave function integration code
12 C          results for finite size spherical sources.  The present
13 C          code reads these tables and interpolates to the requested
14 C          source radius.
15 C
16 C          The input source radii, the number of input files of data
17 C          and the q-space grid are assumed static and are therefore
18 C          hard-coded as follows:
19
20       integer*4 max_c2_coul         ! # q-mesh pts for input Coul. correc.
21       integer*4 ncoulradsteps       ! # input radii values (= # files)
22       parameter (max_c2_coul = 288)
23       parameter (ncoulradsteps = 7)
24
25       real*4 coulradmin     ! min source radius (fm) for Coul Correc. tables
26       real*4 coulradmax     ! max source radius (fm) for Coul Correc. tables
27       real*4 coulradstep    ! Source radius step size for Coul Correc. tables
28       parameter (coulradmin  =  6.0)
29       parameter (coulradmax  = 18.0)
30       parameter (coulradstep =  2.0)
31
32       common/coulmb/
33      1   c2_coul_like(max_c2_coul),c2_coul_unlike(max_c2_coul),
34      2   q_coul(max_c2_coul)
35
36       real*4 c2_coul_like   ! Finite source size Coul. correc for like pairs
37       real*4 c2_coul_unlike ! Finite source size Coul. correc for unlike pairs
38       real*4 q_coul         ! Momentum transfer grid (GeV/c).
39
40 CCC   END common/coulmb/