source: trunk/MagicSoft/Simulation/Corsika/Mmcs/box66.f@ 8266

Last change on this file since 8266 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 BOX66
2
3C-----------------------------------------------------------------------
4C
5C PION INTERACTIONS
6C LIGHT 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(' BOX66 : CURPAR=',1P,9E10.3)
57
58C SELECT LIGHT ISOBAR AND PION MASSES
59 CA = PAMA(8)
60 CB = C(3)
61
62C LIGHT ISOBAR BACKWARD
63 CALL ISOBAR( ECM,1,CB,CA,1 )
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.