source: trunk/MagicSoft/Simulation/Corsika/Mmcs/box69.f@ 6783

Last change on this file since 6783 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.4 KB
Line 
1 SUBROUTINE BOX69
2
3C-----------------------------------------------------------------------
4C
5C PION INTERACTIONS
6C HEAVY ISOBAR BACKWARD, HEAVY MESON FORWARD
7C THIS SUBROUTINE IS CALLED FROM NUCINT
8C-----------------------------------------------------------------------
9
10 IMPLICIT DOUBLE PRECISION (A-H,O-Z)
11*KEEP,PARPAR.
12 COMMON /PARPAR/ CURPAR,SECPAR,PRMPAR,OUTPAR,C,
13 * E00,E00PN,PTOT0,PTOT0N,THICKH,ITYPE,LEVL
14 DOUBLE PRECISION CURPAR(14),SECPAR(14),PRMPAR(14),OUTPAR(14),
15 * C(50),E00,E00PN,PTOT0,PTOT0N,THICKH
16 INTEGER ITYPE,LEVL
17*KEEP,PARPAE.
18 DOUBLE PRECISION GAMMA,COSTHE,PHI,H,T,X,Y,CHI,BETA,GCM,ECM
19 EQUIVALENCE (CURPAR(2),GAMMA), (CURPAR(3),COSTHE),
20 * (CURPAR(4), PHI ), (CURPAR(5), H ),
21 * (CURPAR(6), T ), (CURPAR(7), X ),
22 * (CURPAR(8), Y ), (CURPAR(9), CHI ),
23 * (CURPAR(10),BETA), (CURPAR(11),GCM ),
24 * (CURPAR(12),ECM )
25*KEEP,RUNPAR.
26 COMMON /RUNPAR/ FIXHEI,THICK0,HILOECM,HILOELB,
27 * STEPFC,NRRUN,NSHOW,PATAPE,MONIIN,
28 * MONIOU,MDEBUG,NUCNUC,
29 * CETAPE,
30 * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
31 * N1STTR,MDBASE,
32 * DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
33 * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
34 * ,GHEISH,GHESIG
35 COMMON /RUNPAC/ DSN,HOST,USER
36 DOUBLE PRECISION FIXHEI,THICK0,HILOECM,HILOELB
37 REAL STEPFC
38 INTEGER NRRUN,NSHOW,PATAPE,MONIIN,MONIOU,MDEBUG,NUCNUC,
39 * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
40 * N1STTR,MDBASE
41 INTEGER CETAPE
42 CHARACTER*79 DSN
43 CHARACTER*20 HOST,USER
44
45 LOGICAL DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
46 * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
47 * ,GHEISH,GHESIG
48*KEND.
49
50C-----------------------------------------------------------------------
51
52 IF ( DEBUG ) WRITE(MDEBUG,444) (CURPAR(I),I=1,9)
53 444 FORMAT(' BOX69 : CURPAR=',1P,9E10.3)
54
55C SELECT HEAVY MESON AND HEAVY ISOBAR MASSES
56 CA = C(5)
57 CB = C(4)
58
59C HEAVY ISOBAR BACKWARD
60 CALL ISOBAR( ECM,1,CB,CA,3 )
61
62C HEAVY MESON FORWARD
63 CALL HMESON( ECM,CA,CB )
64
65 RETURN
66 END
Note: See TracBrowser for help on using the repository browser.