| da0e9ce3 |
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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame