Added MG/MW patches and plot macros

BrieucF · Jan 16, 2015 · 6a94240 · 6a94240
1 parent b581100
commit 6a94240
Show file tree

Hide file tree

Showing 19 changed files with 1,093 additions and 0 deletions.
diff --git a/MG_interference_patches/MG_prodonly/patch_MG.sh b/MG_interference_patches/MG_prodonly/patch_MG.sh
@@ -0,0 +1,27 @@
+#!/bin/bash
+
+MGDIR=$HOME/scratch/Madgraph/MG5_aMC_v2_2_1/
+
+MGWORKDIR=$MGDIR/$1/SubProcesses/
+
+SUBPROCS=("gg" "ccx" "ddx" "dsx" "sdx" "ssx" "uux")
+DECAYS=("t_bwp_wp_vlep_tx_bxwm_wm_vlem" "t_bwp_wp_vlep_tx_bxwm_wm_vlmum" "t_bwp_wp_vlmup_tx_bxwm_wm_vlem" "t_bwp_wp_vlmup_tx_bxwm_wm_vlmum") 
+
+for ((i=0;i<${#SUBPROCS[*]};i++)); 
+do
+	proc=${SUBPROCS[$i]}
+	for ((j=0;j<${#DECAYS[*]};j++)); 
+	do
+		decay=${DECAYS[$j]}
+		dir=${MGWORKDIR}/P0_${proc}_ttx_${decay}
+		if [ -d "$dir" ]
+		then
+			echo "Back-uping and patching ${dir}/matrix1.f"
+			cp ${dir}/matrix1.f ${dir}/matrix1.f.old
+			patch ${dir}/matrix1.f patch_${proc}.f
+		else
+			echo "Couldn't find directory $dir".
+		fi
+	done
+done
+
diff --git a/MG_interference_patches/MG_prodonly/patch_ccx.f b/MG_interference_patches/MG_prodonly/patch_ccx.f
@@ -0,0 +1,40 @@
+490c490
+<       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,NAMPSO)
+---
+>       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,2)
+562,565c562,567
+<       JAMP(1,1)=+AMP(1)-1D0/6D0*AMP(2)-1D0/6D0*AMP(3)-1D0/6D0*AMP(4)
+<      $ +AMP(5)+AMP(6)-1D0/6D0*AMP(7)-1D0/6D0*AMP(8)+AMP(9)+AMP(10)
+<      $ +AMP(11)+AMP(12)
+<       JAMP(2,1)=+1D0/2D0*(+AMP(2)+AMP(3)+AMP(4)+AMP(7)+AMP(8))
+---
+>       JAMP(1,1)=+AMP(6)-1D0/6D0*AMP(8)+AMP(10)+AMP(12)
+>       JAMP(1,2)=+AMP(1)-1D0/6D0*AMP(2)-1D0/6D0*AMP(3)-1D0/6D0*AMP(4)
+>      $ +AMP(5)-1D0/6D0*AMP(7)+AMP(9)+AMP(11)
+>       
+>       JAMP(2,1)=+1D0/2D0*(AMP(8))
+>       JAMP(2,2)=+1D0/2D0*(+AMP(2)+AMP(3)+AMP(4)+AMP(7))
+568c570
+<       DO M = 1, NAMPSO
+---
+>       DO M = 1, 2 
+574,576c576,578
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               MATRIX1 = MATRIX1 + ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+---
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               MATRIX1=MATRIX1+ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+591,594c593,596
+<         DO M = 1, NAMPSO
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               JAMP2(I)=JAMP2(I)+JAMP(I,M)*DCONJG(JAMP(I,N))
+---
+>         DO M = 1, 2
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               JAMP2(I)=JAMP2(I)+ABS(JAMP(I,M)*DCONJG(JAMP(I,N)))
+599a602
+> 
diff --git a/MG_interference_patches/MG_prodonly/patch_ddx.f b/MG_interference_patches/MG_prodonly/patch_ddx.f
@@ -0,0 +1,40 @@
+490c490
+<       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,NAMPSO)
+---
+>       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,2)
+565,569c565,571
+<       JAMP(1,1)=+AMP(1)-1D0/6D0*AMP(2)-1D0/6D0*AMP(3)-1D0/6D0*AMP(4)
+<      $ +AMP(5)+AMP(6)-1D0/6D0*AMP(7)-1D0/6D0*AMP(8)+AMP(9)+AMP(10)
+<      $ +AMP(11)+AMP(12)
+<       JAMP(2,1)=+1D0/2D0*AMP(2)+1D0/2D0*AMP(3)+1D0/2D0*AMP(4)
+<      $ +1D0/2D0*AMP(7)+1D0/2D0*AMP(8)-AMP(13)
+---
+>       JAMP(1,1)=+AMP(6)-1D0/6D0*AMP(8)+AMP(10)+AMP(12)
+>       JAMP(1,2)=+AMP(1)-1D0/6D0*AMP(2)-1D0/6D0*AMP(3)-1D0/6D0*AMP(4)
+>      $ +AMP(5)-1D0/6D0*AMP(7)+AMP(9)+AMP(11)
+>       
+>       JAMP(2,1)=+1D0/2D0*AMP(8)-AMP(13)
+>       JAMP(2,2)=+1D0/2D0*AMP(2)+1D0/2D0*AMP(3)+1D0/2D0*AMP(4)
+>      $ +1D0/2D0*AMP(7)
+572c574
+<       DO M = 1, NAMPSO
+---
+>       DO M = 1, 2 
+578,580c580,582
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               MATRIX1 = MATRIX1 + ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+---
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               MATRIX1=MATRIX1+ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+596,599c598,601
+<         DO M = 1, NAMPSO
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               JAMP2(I)=JAMP2(I)+JAMP(I,M)*DCONJG(JAMP(I,N))
+---
+>         DO M = 1, 2
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               JAMP2(I)=JAMP2(I)+ABS(JAMP(I,M)*DCONJG(JAMP(I,N)))
diff --git a/MG_interference_patches/MG_prodonly/patch_dsx.f b/MG_interference_patches/MG_prodonly/patch_dsx.f
@@ -0,0 +1,28 @@
+533,545c533,545
+<       DO M = 1, NAMPSO
+<         DO I = 1, NCOLOR
+<           ZTEMP = (0.D0,0.D0)
+<           DO J = 1, NCOLOR
+<             ZTEMP = ZTEMP + CF(J,I)*JAMP(J,M)
+<           ENDDO
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               MATRIX1 = MATRIX1 + ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+<             ENDIF
+<           ENDDO
+<         ENDDO
+<       ENDDO
+---
+> C     DO M = 1, NAMPSO
+> C       DO I = 1, NCOLOR
+> C         ZTEMP = (0.D0,0.D0)
+> C         DO J = 1, NCOLOR
+> C           ZTEMP = ZTEMP + CF(J,I)*JAMP(J,M)
+> C         ENDDO
+> C         DO N = 1, NAMPSO
+> C           IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+> C             MATRIX1 = MATRIX1 + ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+> C           ENDIF
+> C         ENDDO
+> C       ENDDO
+> C     ENDDO
diff --git a/MG_interference_patches/MG_prodonly/patch_gg.f b/MG_interference_patches/MG_prodonly/patch_gg.f
@@ -0,0 +1,41 @@
+490c490
+<       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,NAMPSO)
+---
+>       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,2)
+562,566c562,569
+<       JAMP(1,1)=-IMAG1*AMP(1)-IMAG1*AMP(3)-IMAG1*AMP(4)-IMAG1*AMP(5)
+<      $ +AMP(6)+AMP(7)+AMP(8)
+<       JAMP(2,1)=+IMAG1*AMP(1)+IMAG1*AMP(3)+IMAG1*AMP(4)+IMAG1*AMP(5)
+<      $ +AMP(9)+AMP(10)+AMP(11)
+<       JAMP(3,1)=+2D0*(+AMP(2))
+---
+>       JAMP(1,1)=-IMAG1*AMP(5)+AMP(8)
+>       JAMP(1,2)=-IMAG1*AMP(1)-IMAG1*AMP(3)-IMAG1*AMP(4)+AMP(6)+AMP(7)
+> 
+>       JAMP(2,1)=+IMAG1*AMP(5)+AMP(11)
+>       JAMP(2,2)=+IMAG1*AMP(1)+IMAG1*AMP(3)+IMAG1*AMP(4)+AMP(9)+AMP(10)
+> 
+>       JAMP(3,1)=0.D0
+>       JAMP(3,2)=+2D0*(+AMP(2))
+569c572
+<       DO M = 1, NAMPSO
+---
+>       DO M = 1, 2 
+575,577c578,580
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               MATRIX1 = MATRIX1 + ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+---
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               MATRIX1=MATRIX1+ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+591,594c594,597
+<         DO M = 1, NAMPSO
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               JAMP2(I)=JAMP2(I)+JAMP(I,M)*DCONJG(JAMP(I,N))
+---
+>         DO M = 1, 2
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               JAMP2(I)=JAMP2(I)+ABS(JAMP(I,M)*DCONJG(JAMP(I,N)))
diff --git a/MG_interference_patches/MG_prodonly/patch_sdx.f b/MG_interference_patches/MG_prodonly/patch_sdx.f
@@ -0,0 +1,28 @@
+533,545c533,545
+<       DO M = 1, NAMPSO
+<         DO I = 1, NCOLOR
+<           ZTEMP = (0.D0,0.D0)
+<           DO J = 1, NCOLOR
+<             ZTEMP = ZTEMP + CF(J,I)*JAMP(J,M)
+<           ENDDO
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               MATRIX1 = MATRIX1 + ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+<             ENDIF
+<           ENDDO
+<         ENDDO
+<       ENDDO
+---
+> C     DO M = 1, NAMPSO
+> C       DO I = 1, NCOLOR
+> C         ZTEMP = (0.D0,0.D0)
+> C         DO J = 1, NCOLOR
+> C           ZTEMP = ZTEMP + CF(J,I)*JAMP(J,M)
+> C         ENDDO
+> C         DO N = 1, NAMPSO
+> C           IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+> C             MATRIX1 = MATRIX1 + ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+> C           ENDIF
+> C         ENDDO
+> C       ENDDO
+> C     ENDDO
diff --git a/MG_interference_patches/MG_prodonly/patch_ssx.f b/MG_interference_patches/MG_prodonly/patch_ssx.f
@@ -0,0 +1,40 @@
+490c490
+<       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,NAMPSO)
+---
+>       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,2)
+565,569c565,571
+<       JAMP(1,1)=+AMP(1)-1D0/6D0*AMP(2)-1D0/6D0*AMP(3)-1D0/6D0*AMP(4)
+<      $ +AMP(5)+AMP(6)-1D0/6D0*AMP(7)-1D0/6D0*AMP(8)+AMP(9)+AMP(10)
+<      $ +AMP(11)+AMP(12)
+<       JAMP(2,1)=+1D0/2D0*AMP(2)+1D0/2D0*AMP(3)+1D0/2D0*AMP(4)
+<      $ +1D0/2D0*AMP(7)+1D0/2D0*AMP(8)-AMP(13)
+---
+>       JAMP(1,1)=+AMP(6)-1D0/6D0*AMP(8)+AMP(10)+AMP(12)
+>       JAMP(1,2)=+AMP(1)-1D0/6D0*AMP(2)-1D0/6D0*AMP(3)-1D0/6D0*AMP(4)
+>      $ +AMP(5)-1D0/6D0*AMP(7)+AMP(9)+AMP(11)
+>       
+>       JAMP(2,1)=+1D0/2D0*AMP(8)-AMP(13)
+>       JAMP(2,2)=+1D0/2D0*AMP(2)+1D0/2D0*AMP(3)+1D0/2D0*AMP(4)
+>      $ +1D0/2D0*AMP(7)
+572c574
+<       DO M = 1, NAMPSO
+---
+>       DO M = 1, 2 
+578,580c580,582
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               MATRIX1 = MATRIX1 + ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+---
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               MATRIX1=MATRIX1+ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+596,599c598,601
+<         DO M = 1, NAMPSO
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               JAMP2(I)=JAMP2(I)+JAMP(I,M)*DCONJG(JAMP(I,N))
+---
+>         DO M = 1, 2
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               JAMP2(I)=JAMP2(I)+ABS(JAMP(I,M)*DCONJG(JAMP(I,N)))
diff --git a/MG_interference_patches/MG_prodonly/patch_uux.f b/MG_interference_patches/MG_prodonly/patch_uux.f
@@ -0,0 +1,38 @@
+490c490
+<       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,NAMPSO)
+---
+>       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,2)
+562,565c562,567
+<       JAMP(1,1)=+AMP(1)-1D0/6D0*AMP(2)-1D0/6D0*AMP(3)-1D0/6D0*AMP(4)
+<      $ +AMP(5)+AMP(6)-1D0/6D0*AMP(7)-1D0/6D0*AMP(8)+AMP(9)+AMP(10)
+<      $ +AMP(11)+AMP(12)
+<       JAMP(2,1)=+1D0/2D0*(+AMP(2)+AMP(3)+AMP(4)+AMP(7)+AMP(8))
+---
+>       JAMP(1,1)=+AMP(6)-1D0/6D0*AMP(8)+AMP(10)+AMP(12)
+>       JAMP(1,2)=+AMP(1)-1D0/6D0*AMP(2)-1D0/6D0*AMP(3)-1D0/6D0*AMP(4)
+>      $ +AMP(5)-1D0/6D0*AMP(7)+AMP(9)+AMP(11)
+>       
+>       JAMP(2,1)=+1D0/2D0*(AMP(8))
+>       JAMP(2,2)=+1D0/2D0*(+AMP(2)+AMP(3)+AMP(4)+AMP(7))
+568c570
+<       DO M = 1, NAMPSO
+---
+>       DO M = 1, 2 
+574,576c576,578
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               MATRIX1 = MATRIX1 + ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+---
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               MATRIX1=MATRIX1+ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+591,594c593,596
+<         DO M = 1, NAMPSO
+<           DO N = 1, NAMPSO
+<             IF (CHOSEN_SO_CONFIGS(SQSOINDEX1(M,N))) THEN
+<               JAMP2(I)=JAMP2(I)+JAMP(I,M)*DCONJG(JAMP(I,N))
+---
+>         DO M = 1, 2
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               JAMP2(I)=JAMP2(I)+ABS(JAMP(I,M)*DCONJG(JAMP(I,N)))
diff --git a/MG_interference_patches/MW_prodonly/patch_1_gg.f b/MG_interference_patches/MW_prodonly/patch_1_gg.f
@@ -0,0 +1,42 @@
+357c357
+<       INTEGER I,J
+---
+>       INTEGER I,J,M,N
+360c360
+<       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
+---
+>       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,2)
+424,428c424,431
+<       JAMP(1)=-IMAG1*AMP(1)-IMAG1*AMP(3)-IMAG1*AMP(4)-IMAG1*AMP(5)
+<      $ +AMP(6)+AMP(7)+AMP(8)
+<       JAMP(2)=+IMAG1*AMP(1)+IMAG1*AMP(3)+IMAG1*AMP(4)+IMAG1*AMP(5)
+<      $ +AMP(9)+AMP(10)+AMP(11)
+<       JAMP(3)=+2D0*(+AMP(2))
+---
+>       JAMP(1,1)=-IMAG1*AMP(5)+AMP(8)
+>       JAMP(1,2)=-IMAG1*AMP(1)-IMAG1*AMP(3)-IMAG1*AMP(4)+AMP(6)+AMP(7)
+> 
+>       JAMP(2,1)=+IMAG1*AMP(5)+AMP(11)
+>       JAMP(2,2)=+IMAG1*AMP(1)+IMAG1*AMP(3)+IMAG1*AMP(4)+AMP(9)+AMP(10)
+> 
+>       JAMP(3,1)=0.D0
+>       JAMP(3,2)=+2D0*(+AMP(2))
+431,434c434,444
+<       DO I = 1, NCOLOR
+<         ZTEMP = (0.D0,0.D0)
+<         DO J = 1, NCOLOR
+<           ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
+---
+>       DO M = 1, 2 
+>         DO I = 1, NCOLOR
+>           ZTEMP = (0.D0,0.D0)
+>           DO J = 1, NCOLOR
+>             ZTEMP = ZTEMP + CF(J,I)*JAMP(J,M)
+>           ENDDO
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               MATRIX1=MATRIX1+ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+>             ENDIF
+>           ENDDO
+436d445
+<         MATRIX1 = MATRIX1+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
diff --git a/MG_interference_patches/MW_prodonly/patch_2_uux.f b/MG_interference_patches/MW_prodonly/patch_2_uux.f
@@ -0,0 +1,39 @@
+357c357
+<       INTEGER I,J
+---
+>       INTEGER I,J,M,N
+360c360
+<       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
+---
+>       COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR,2)
+424,427c424,429
+<       JAMP(1)=+AMP(1)-1D0/6D0*AMP(2)-1D0/6D0*AMP(3)-1D0/6D0*AMP(4)
+<      $ +AMP(5)+AMP(6)-1D0/6D0*AMP(7)-1D0/6D0*AMP(8)+AMP(9)+AMP(10)
+<      $ +AMP(11)+AMP(12)
+<       JAMP(2)=+1D0/2D0*(+AMP(2)+AMP(3)+AMP(4)+AMP(7)+AMP(8))
+---
+>       JAMP(1,1)=+AMP(6)-1D0/6D0*AMP(8)+AMP(10)+AMP(12)
+>       JAMP(1,2)=+AMP(1)-1D0/6D0*AMP(2)-1D0/6D0*AMP(3)-1D0/6D0*AMP(4)
+>      $ +AMP(5)-1D0/6D0*AMP(7)+AMP(9)+AMP(11)
+>       
+>       JAMP(2,1)=+1D0/2D0*(AMP(8))
+>       JAMP(2,2)=+1D0/2D0*(+AMP(2)+AMP(3)+AMP(4)+AMP(7))
+430,433c432,442
+<       DO I = 1, NCOLOR
+<         ZTEMP = (0.D0,0.D0)
+<         DO J = 1, NCOLOR
+<           ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
+---
+>       DO M = 1, 2 
+>         DO I = 1, NCOLOR
+>           ZTEMP = (0.D0,0.D0)
+>           DO J = 1, NCOLOR
+>             ZTEMP = ZTEMP + CF(J,I)*JAMP(J,M)
+>           ENDDO
+>           DO N = 1, 2
+>             IF ((M+N).EQ.3) THEN
+>               MATRIX2=MATRIX2+ZTEMP*DCONJG(JAMP(I,N))/DENOM(I)
+>             ENDIF
+>           ENDDO
+435d443
+<         MATRIX2 = MATRIX2+ZTEMP*DCONJG(JAMP(I))/DENOM(I)