source: trunk/MagicSoft/Simulation/Corsika/Mmcs/box68.f@ 10146

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