source: trunk/MagicSoft/Simulation/Corsika/Mmcs/box73.f@ 4438

Last change on this file since 4438 was 286, checked in by harald, 25 years ago
This is the start point for further developments of the Magic Monte Carlo Simulation written by Jose Carlos Gonzales. Now it is under control of one CVS repository for the whole collaboration. Everyone should use this CVS repository for further developments.
File size: 2.5 KB
Line 
1 SUBROUTINE BOX73
2
3C-----------------------------------------------------------------------
4C
5C KAON INTERACTIONS
6C HEAVY ISOBAR BACKWARD, KAON FORWARD
7C THIS SUBROUTINE IS CALLED FROM NUCINT
8C-----------------------------------------------------------------------
9
10 IMPLICIT DOUBLE PRECISION (A-H,O-Z)
11*KEEP,PAM.
12 COMMON /PAM/ PAMA,SIGNUM
13 DOUBLE PRECISION PAMA(6000),SIGNUM(6000)
14*KEEP,PARPAR.
15 COMMON /PARPAR/ CURPAR,SECPAR,PRMPAR,OUTPAR,C,
16 * E00,E00PN,PTOT0,PTOT0N,THICKH,ITYPE,LEVL
17 DOUBLE PRECISION CURPAR(14),SECPAR(14),PRMPAR(14),OUTPAR(14),
18 * C(50),E00,E00PN,PTOT0,PTOT0N,THICKH
19 INTEGER ITYPE,LEVL
20*KEEP,PARPAE.
21 DOUBLE PRECISION GAMMA,COSTHE,PHI,H,T,X,Y,CHI,BETA,GCM,ECM
22 EQUIVALENCE (CURPAR(2),GAMMA), (CURPAR(3),COSTHE),
23 * (CURPAR(4), PHI ), (CURPAR(5), H ),
24 * (CURPAR(6), T ), (CURPAR(7), X ),
25 * (CURPAR(8), Y ), (CURPAR(9), CHI ),
26 * (CURPAR(10),BETA), (CURPAR(11),GCM ),
27 * (CURPAR(12),ECM )
28*KEEP,RUNPAR.
29 COMMON /RUNPAR/ FIXHEI,THICK0,HILOECM,HILOELB,
30 * STEPFC,NRRUN,NSHOW,PATAPE,MONIIN,
31 * MONIOU,MDEBUG,NUCNUC,
32 * CETAPE,
33 * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
34 * N1STTR,MDBASE,
35 * DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
36 * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
37 * ,GHEISH,GHESIG
38 COMMON /RUNPAC/ DSN,HOST,USER
39 DOUBLE PRECISION FIXHEI,THICK0,HILOECM,HILOELB
40 REAL STEPFC
41 INTEGER NRRUN,NSHOW,PATAPE,MONIIN,MONIOU,MDEBUG,NUCNUC,
42 * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
43 * N1STTR,MDBASE
44 INTEGER CETAPE
45 CHARACTER*79 DSN
46 CHARACTER*20 HOST,USER
47
48 LOGICAL DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
49 * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
50 * ,GHEISH,GHESIG
51*KEND.
52
53C-----------------------------------------------------------------------
54
55 IF ( DEBUG ) WRITE(MDEBUG,444) (CURPAR(I),I=1,9)
56 444 FORMAT(' BOX73 : CURPAR=',1P,9E10.3)
57
58C SELECT HEAVY ISOBAR AND KAON MASSES
59 CA = PAMA(ITYPE)
60 CB = C(4)
61
62C HEAVY ISOBAR BACKWARD
63 CALL ISOBAR( ECM,1,CA,CB,3 )
64
65C KAON FORWARD
66 CA = PAMA(ITYPE)
67 CALL SINGLE( ECM,0,CA,CB )
68
69 RETURN
70 END
Note: See TracBrowser for help on using the repository browser.