HERWIG/src/hwuant.f

   1
   2 CDECK  ID>, HWUANT.
   3
   4 *CMZ :-        -27/07/99  13.33.03  by  Mike Seymour
   5
   6 *-- Author :    Ian Knowles
   7
   8 C-----------------------------------------------------------------------
   9
  10       FUNCTION HWUANT(IPART)
  11
  12 C-----------------------------------------------------------------------
  13
  14 C     Returns the antiparticle of IPART; uses HERWIG numbering
  15
  16 C-----------------------------------------------------------------------
  17
  18       INCLUDE 'HERWIG61.INC'
  19
  20       INTEGER HWUANT,IPART,IPDG,IANTI,OLDERR
  21
  22       CHARACTER*8 CDUM
  23
  24       OLDERR=IERROR
  25
  26       IPDG=IDPDG(IPART)
  27
  28       IF (IPDG.EQ. 9.OR.IPDG.EQ.21.OR.IPDG.EQ.22.OR.IPDG.EQ.23.OR.
  29
  30      &    IPDG.EQ.25.OR.IPDG.EQ.26.OR.IPDG.EQ.32.OR.IPDG.EQ.35.OR.
  31
  32      &    IPDG.EQ.36.OR.IPDG.EQ.91.OR.IPDG.EQ.98.OR.IPDG.EQ.99.OR.
  33
  34      &    IPDG.EQ.130.OR.IPDG.EQ.310.OR.
  35
  36      &    IPDG.EQ.1000021.OR.IPDG.EQ.1000022.OR.IPDG.EQ.1000023.OR.
  37
  38      &    IPDG.EQ.1000025.OR.IPDG.EQ.1000035.OR.IPDG.EQ.1000039.OR.
  39
  40      &    (FLOAT(INT(RSPIN(IPART))).EQ.RSPIN(IPART).AND.
  41
  42      &     MOD(IPDG/100,10).EQ.MOD(IPDG/10,10).AND.
  43
  44      &     MOD(IPDG/10,10).NE.0)) THEN
  45
  46 C Self-conjugate boson
  47
  48         IANTI=IPART
  49
  50       ELSEIF(IPART.EQ.211.OR.IPART.EQ.212) THEN
  51
  52 C Fourth generation (anti-)quarks
  53
  54         IANTI=IPART+6
  55
  56       ELSEIF(IPART.EQ.217.OR.IPART.EQ.218) THEN
  57
  58         IANTI=IPART-6
  59
  60       ELSE
  61
  62 C Non-zero charge particle
  63
  64         CALL HWUIDT(1,-IPDG,IANTI,CDUM)
  65
  66       ENDIF
  67
  68       IF (IANTI.EQ.20) WRITE(6,10) RNAME(IPART)
  69
  70   10  FORMAT(1X,A8,' has no antiparticle'/)
  71
  72       HWUANT=IANTI
  73
  74       IERROR=OLDERR
  75
  76       END