source: trunk/MagicSoft/Simulation/Corsika/Mmcs/box67.f@ 502

Last change on this file since 502 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 BOX67
2
3C-----------------------------------------------------------------------
4C
5C PION INTERACTIONS
6C HEAVY MESON FORWARD, NUCLEON BACKWARD
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(' BOX67 : CURPAR=',1P,9E10.3)
57
58C SELECT HEAVY MESON AND NUCLEON MASSES
59 CA = C(5)
60 CB = PAMA(14)
61
62C HEAVY MESON FORWARD
63 CALL HMESON( ECM,CA,CB )
64
65C NUCLEON BACKWARD
66 CALL SINGLE( ECM,1,CB,CA )
67
68 RETURN
69 END
Note: See TracBrowser for help on using the repository browser.