]>
Commit | Line | Data |
---|---|---|
1 | * | |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1995/10/24 10:19:59 cernlib | |
6 | * Geant | |
7 | * | |
8 | * | |
9 | #include "geant321/pilot.h" | |
10 | *CMZ : 3.21/02 29/03/94 15.41.44 by S.Giani | |
11 | *-- Author : | |
12 | *$ CREATE CALUMO.FOR | |
13 | *COPY CALUMO | |
14 | * | |
15 | *=== calumo ===========================================================* | |
16 | * | |
17 | SUBROUTINE CALUMO (N,ITTA) | |
18 | ||
19 | #include "geant321/dblprc.inc" | |
20 | #include "geant321/dimpar.inc" | |
21 | #include "geant321/iounit.inc" | |
22 | C*** C.M.S.-ENERGY AND REACTION CHANNEL THRESHOLD CALCULATION | |
23 | * | |
24 | #include "geant321/reac.inc" | |
25 | #include "geant321/redver.inc" | |
26 | #include "geant321/split.inc" | |
27 | * | |
28 | **** *,IEII,IKII,NURE | |
29 | COMMON / FKABLT / AM (110), GA (110), TAU (110), ICH (110), | |
30 | & IBAR (110), K1 (110), K2 (110) | |
31 | IRE=NURE(N,ITTA/8+1) | |
32 | IEO=IEII (IRE)+1 | |
33 | IEE=IEII (IRE+1) | |
34 | AM1 = AM (N) | |
35 | AM12 = AM1**2 | |
36 | AM2 = AM (ITTA) | |
37 | AM22 = AM2**2 | |
38 | DO 1 IE=IEO,IEE | |
39 | PLAB2=PLABF(IE)**2 | |
40 | ELAB=SQRT(AM12+AM22+2.D+00*SQRT(PLAB2+AM12)*AM2) | |
41 | UMO(IE)=ELAB | |
42 | 1 CONTINUE | |
43 | IKO=IKII(IRE)+1 | |
44 | IKE=IKII(IRE +1) | |
45 | UMOO=UMO(IEO) | |
46 | DO 2 IK=IKO,IKE | |
47 | IF (NRK(2,IK).GT.0) GO TO 2 | |
48 | IKI=NRK(1,IK) | |
49 | AMSS=5.D+00 | |
50 | K11=K1(IKI) | |
51 | K22=K2(IKI) | |
52 | DO 3 IK1=K11,K22 | |
53 | IN=NZK(IK1,1) | |
54 | AMS=AM(IN) | |
55 | IN=NZK(IK1,2) | |
56 | IF (IN.GT.0) AMS=AMS+AM(IN) | |
57 | IN=NZK(IK1,3) | |
58 | IF (IN.GT.0) AMS=AMS+AM(IN) | |
59 | IF (AMS.LT.AMSS) AMSS=AMS | |
60 | 3 CONTINUE | |
61 | IF ( UMOO .LT. AMSS ) UMOO=AMSS | |
62 | THRESH (IK) = UMOO | |
63 | 2 CONTINUE | |
64 | RETURN | |
65 | END |