← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.1

« back to all changes in this revision

Viewing changes to MG4_DECAY/decay.inc

move ./decay to ./mg5decay; resolve unit tests (n.b. __init__ does not check keys of input dictionaries, followed last revision)

Show diffs side-by-side

added added

removed removed

Lines of Context:
MG4_DECAY/decay.inc
1
 
c*****************************************************
2
 
c     PARAMETERS
3
 
c*****************************************************          
4
 
c
5
 
c     VEGAS MAX DIMENSIONS
6
 
c
7
 
      integer    mxdim          !to be the same as in vegas.f 
8
 
      parameter (mxdim=10)      !rebin.f
9
 
c                               
10
 
      integer    MaxParticles    ! max number of particles in 
11
 
      parameter (MaxParticles=25)! the event
12
 
      integer    MaxDecPart      ! max particles in the decay
13
 
      parameter (MaxDecPart=5)   ! 1->2,..,MaxDecPart  
14
 
      integer    MaxiEvent       ! max nr of processes
15
 
      parameter (MaxiEvent=100)  ! in the proc_card
16
 
c
17
 
      integer MaxPoints         ! maxpoints used in the integration
18
 
      parameter(MaxPoints=500000)
19
 
c
20
 
c     Numbers
21
 
c
22
 
      double precision  Zero, One, Two, Three, Four, Half,Rt2
23
 
      parameter( Zero = 0.0d0, One = 1.0d0, Two = 2.0d0 )
24
 
      parameter( Three = 3.0d0, Four = 4.0d0, Half = 0.5d0 )
25
 
      parameter( Rt2   = 1.414213562d0 )
26
 
      double precision  Pi, Fourpi
27
 
      parameter( Pi = 3.14159265358979323846d0 )
28
 
      parameter( Fourpi = Four * Pi )
29
 
c
30
 
c     Units     
31
 
c
32
 
      integer luni,lunr,lunw,lunp,lunrnd,lunout
33
 
      parameter (luni=15,lunr=16,lunw=17,lunp=18,lunrnd=19,lunout=20)   
34
 
c
35
 
c*****************************************************
36
 
c     GLOBAL
37
 
c*****************************************************          
38
 
c
39
 
c     IP    = PDG CODE FOR THE PARTICLE TO DECAY
40
 
c     IMODE = DECAY MODE
41
 
c     ND    = NUMBER OF DECAY PRODUCTS
42
 
c
43
 
      INTEGER       IP,IMODE,ND 
44
 
      COMMON/DEC_ID/IP,IMODE,ND
45
 
c
46
 
c     masses and couplings for the decay routines
47
 
c
48
 
      REAL*8          M1,M2,M3,M4,M5,MV,GV
49
 
      COMMON/m_decays/M1,M2,M3,M4,M5,MV,GV
50
 
      DOUBLE COMPLEX       GX,GXX(2),GXX1(2)  
51
 
      COMMON/WEAK_COUPL/   GX,GXX   ,GXX1
52
 
c
53
 
c    MC width, analytic width and best branching ratio
54
 
c       
55
 
      real*8          MC_width,calc_width,bratio,calc_br
56
 
      common/to_width/MC_width,calc_width,bratio,calc_br
57
 
c
58
 
c    Particle's ID
59
 
c       
60
 
      character*3 id(-100:100)
61
 
      real*8 pdgmass(-100:100)
62
 
      common/to_id/pdgmass,id
63
 
c
64
 
c    Color label
65
 
c       
66
 
      integer cindex 
67
 
      common/c_ind/cindex
68
 
c
69
 
c    Vegas      
70
 
c
71
 
      real*4          x(MXDIM)
72
 
      COMMON/TO_VEGAS/x
73
 
      integer                                 n_wei,n_ove
74
 
      real*8          mxwgt,s_wei,s_unw,s_ove
75
 
      common/to_event/mxwgt,s_wei,s_unw,s_ove,n_wei,n_ove
76
 
c
77
 
c     HDECAY branching ratios
78
 
c
79
 
      integer           ihiggs, nnlo, ipole
80
 
      common /flag/     ihiggs, nnlo, ipole  
81
 
      integer           ionsh, ionwz, iofsusy
82
 
      common /onshell/  ionsh, ionwz, iofsusy 
83
 
      double precision  gf, alph, almass, ammuon, amz, amw
84
 
      common /param/    gf, alph, almass, ammuon, amz, amw 
85
 
      double precision  ams, amc, amb, amt
86
 
      common /masses/   ams, amc, amb, amt 
87
 
      double precision  gamc0, gamt0, gamt1, gamw, gamz
88
 
      common /wzwdth/   gamc0, gamt0, gamt1, gamw, gamz 
89
 
      double precision  SMBRB,SMBRL,SMBRM,SMBRS,SMBRC,SMBRT,
90
 
     &                  SMBRG,SMBRGA,SMBRZGA,SMBRW,SMBRZ,SMWDTH
91
 
      common /widthsm/  SMBRB,SMBRL,SMBRM,SMBRS,SMBRC,SMBRT,
92
 
     &                  SMBRG,SMBRGA,SMBRZGA,SMBRW,SMBRZ,SMWDTH
93
 
c
94
 
      double precision  hdals, hdmhbeg, hdmhend, tgbet
95
 
      common /hdparms/  hdals, hdmhbeg, hdmhend, tgbet
96
 
c
97
 
c     PDG branching ratios
98
 
c
99
 
      double precision br_ta_lv,br_ta_pi,br_ta_ro
100
 
      double precision br_w_lv,br_w_jj
101
 
      double precision br_z_ll,br_z_vv,br_z_cc,br_z_bb,br_z_jj 
102
 
      common/brs/      br_ta_lv,br_ta_pi,br_ta_ro,
103
 
     &                 br_w_lv,br_w_jj,
104
 
     &                 br_z_ll,br_z_vv,br_z_cc,br_z_bb,br_z_jj 
105