Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / GEANT321 / giface / gmunui.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.1 1995/10/24 10:21:15 cernlib
6* Geant
7*
8*
9#include "geant321/pilot.h"
10*CMZ : 3.21/02 29/03/94 15.41.38 by S.Giani
11*-- Author :
12 SUBROUTINE GMUNUI
13C
14C *** CALCULATION OF MUON-NUCLEUS CROSS SECTIONS IN CURRENT MATERIAL ***
15C *** NVE 16-MAR-1988 CERN GENEVA ***
16C
17C CALLED BY : GPHYSI
18C ORIGIN : F.CARMINATI, H.FESEFELDT
19C
20#include "geant321/gcbank.inc"
21#include "geant321/gcjloc.inc"
22#include "geant321/gcmulo.inc"
23#include "geant321/gcmate.inc"
24#include "geant321/gconsp.inc"
25#include "geant321/gctrak.inc"
26C
27C
28C --- GET MUON MASS FROM GEANT PARTICLE DATA BANK ---
29 LNVE=LQ(JPART-5)
30 HMAMU=Q(LNVE+7)
31C
32C --- Infinite mean free patch for EKIN < 1 GeV (HJT/NVE 27-sep-1990) ---
33 SIG=0.
34 IF(JMUNU.EQ.0)GO TO 9999
35 IF (ELOW(IEKBIN) .LT. 1.) GO TO 90
36C
37 IF(Z.LT.1.)GO TO 90
38 EN = ELOW(IEKBIN)+HMAMU
39 SIG = 0.0003
40 IF(EN.GT.30.) SIG = SIG*(EN/30.)**0.25
41 SIG = SIG*DENS*AVO*1E-3
42C
43 90 IF(SIG.GT.0.)THEN
44 Q(JMUNU+IEKBIN)=1./SIG
45 ELSE
46 Q(JMUNU+IEKBIN)=BIG
47 ENDIF
48C
49 9999 CONTINUE
50 END