TEvtGen/PHOTOS/photwo.F

   1       SUBROUTINE PHOTWO(MODE)
   2 C.----------------------------------------------------------------------
   3 C.
   4 C.    PHOTWO:   PHOtos but TWO mothers allowed
   5 C.
   6 C.    Purpose:  Combines two mothers into one in /PHOEVT/
   7 C.              necessary eg in case of g g (q qbar) --> t tbar
   8 C.
   9 C.    Input Parameters: Common /PHOEVT/ (/PHOCMS/)
  10 C.
  11 C.    Output Parameters:  Common /PHOEVT/, (stored mothers)
  12 C.
  13 C.    Author(s):  Z. Was                          Created at:  5/08/93
  14 C.                                                Last Update:10/08/93
  15 C.
  16 C.----------------------------------------------------------------------
  17       IMPLICIT NONE
  18       INTEGER NMXPHO
  19       PARAMETER (NMXPHO=10000)
  20       INTEGER IDPHO,ISTPHO,JDAPHO,JMOPHO,NEVPHO,NPHO
  21       REAL*8 PPHO,VPHO
  22       COMMON/PHOEVT/NEVPHO,NPHO,ISTPHO(NMXPHO),IDPHO(NMXPHO),
  23      &JMOPHO(2,NMXPHO),JDAPHO(2,NMXPHO),PPHO(5,NMXPHO),VPHO(4,NMXPHO)
  24       DOUBLE PRECISION BET(3),GAM
  25       COMMON /PHOCMS/ BET,GAM
  26       INTEGER I,MODE
  27       REAL*8 MPASQR
  28       LOGICAL IFRAD
  29 C logical IFRAD is used to tag cases when two mothers may be
  30 C merged to the sole one.
  31 C So far used in case:
  32 C                      1) of t tbar production
  33 C
  34 C t tbar case
  35       IF(MODE.EQ.0) THEN
  36        IFRAD=(IDPHO(1).EQ.21).AND.(IDPHO(2).EQ.21)
  37        IFRAD=IFRAD.OR.(IDPHO(1).EQ.-IDPHO(2).AND.ABS(IDPHO(1)).LE.6)
  38        IFRAD=IFRAD
  39      &       .AND.(ABS(IDPHO(3)).EQ.6).AND.(ABS(IDPHO(4)).EQ.6)
  40         MPASQR= (PPHO(4,1)+PPHO(4,2))**2-(PPHO(3,1)+PPHO(3,2))**2
  41      &          -(PPHO(2,1)+PPHO(2,2))**2-(PPHO(1,1)+PPHO(1,2))**2
  42        IFRAD=IFRAD.AND.(MPASQR.GT.0.0D0)
  43        IF(IFRAD) THEN
  44 c.....combining first and second mother
  45             DO I=1,4
  46             PPHO(I,1)=PPHO(I,1)+PPHO(I,2)
  47             ENDDO
  48             PPHO(5,1)=SQRT(MPASQR)
  49 c.....removing second mother,
  50             DO I=1,5
  51               PPHO(I,2)=0.0D0
  52             ENDDO
  53        ENDIF
  54       ELSE
  55 C boosting of the mothers to the reaction frame not implemented yet.
  56 C to do it in mode 0 original mothers have to be stored in new comon (?)
  57 C and in mode 1 boosted to cms.
  58       ENDIF
  59       END