This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / gphys / gbflos.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:21:21  cernlib
6 * Geant
7 *
8 *
9 #include "geant321/pilot.h"
10 *CMZ :  3.21/02 29/03/94  15.41.21  by  S.Giani
11 *-- Author :
12       FUNCTION GBFLOS(T,C)
13 C.    ******************************************************************
14 C.    *                                                                *
15 C.    *  It calculates a weight factor for the positron continuous     *
16 C.    *       bremsstrahlung energy loss.                              *
17 C.    *                                                                *
18 C.    *  Input:T kinetic energy (GeV)                                  *
19 C.    *        C cut for bremsstrahlung (BCUTE)                        *
20 C.    *        ETA ratio of positron to electron brems loss            *
21 C.    *  Output: weight factor                                         *
22 C.    *                                                                *
23 C.    *  ===>   positron loss = FLOS * electron loss   <===            *
24 C.    *                                                                *
25 C.    *  Author   L. Urban                                             *
26 C.    *                                                                *
27 C.    ******************************************************************
28 C.
29       ETA = GCBREM(T)
30       IF(ETA.LT.0.0001) THEN
31          F=1.E-10
32       ELSEIF(ETA.GT.0.9999) THEN
33          F=1.
34       ELSE
35 C
36          E0=C/T
37          IF(E0.GT.1.) E0=1.
38 C
39          IF(E0.LT.1.E-8) THEN
40             F=1.
41          ELSE
42             ALFA=1./ETA
43             E1=1.-E0
44             F=ETA*(1.-E1**ALFA)/E0
45          ENDIF
46 C
47       ENDIF
48 *
49       GBFLOS=F
50       END