← Back to branch summary

~maddevelopers/mg5amcnlo/2.9.4

« back to all changes in this revision

Viewing changes to vendor/StdHEP/src/stdhep/stdquarks.F

pass to v2.0.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
vendor/StdHEP/src/stdhep/stdquarks.F
 
1
      subroutine stdquarks(it,kq1,kq2,kq3,kql,kqj,kqr,kqx)
 
2
 
 
3
C...Get constituent quarks of this particle
 
4
C
 
5
C      IT = index to HEPEVT common block
 
6
C      For particle ID, +/- HIJKLMN
 
7
C            KQX = H = flag tentative and special assignments
 
8
C            KQR = I = radial excitations
 
9
C            KQL = J = spin and orbital angular momentum
 
10
C            KQ3 = K = quark
 
11
C            KQ2 = L = quark
 
12
C            KQ1 = M = quark
 
13
C            KQJ = N = 2*Jspin + 1
 
14
 
 
15
      implicit none
 
16
 
 
17
      integer it,kq,kqa,kqx,kq3,kq2,kq1,kqj,kql,kqr,irt
 
18
 
 
19
#include "stdhep.inc"
 
20
 
 
21
      kq=idhep(it)
 
22
      kqa=iabs(kq)
 
23
      kqx=mod(kqa/1000000,10)
 
24
      kqr=mod(kqa/100000,10)
 
25
      kql=mod(kqa/10000,10)
 
26
      kq3=mod(kqa/1000,10)
 
27
      kq2=mod(kqa/100,10)
 
28
      kq1=mod(kqa/10,10)
 
29
      kqj=mod(kqa,10)
 
30
      irt = mod(kqa,10000)
 
31
C... set quarks to zero if not meson, baryon, or diquark
 
32
      if(kqa.le.100 .or. kqa.ge.10000000) then
 
33
        kq1=0
 
34
        kq2=0
 
35
        kq3=0
 
36
      elseif(kqx.gt.0 .and. irt.le.100) then
 
37
        kq1=0
 
38
        kq2=0
 
39
        kq3=0
 
40
      endif
 
41
      return
 
42
      end