Page MenuHomeHEPForge
Files F19250566

No OneTemporary
Actions

Size
80 KB
Referenced Files
None
Subscribers
None

View Options

Index: branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_M2Tree.F
===================================================================
--- branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_M2Tree.F (revision 1091)
+++ branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_M2Tree.F (revision 1092)
@@ -1,912 +1,784 @@
************************************************
function stst2qq_M2Tree()
#include "stst2QQ_Kinematics.h"
#include "stst2QQ_Model.h"
#include "stst2QQ_GenCouplings.h"
double precision findValW
integer i,j,k,l,m,n
double precision muupu(2,4,3), muuqu(2,4,3)
double complex stst2qq_M2Tree
complex*16 Den(4)
- do i=1,4 !hier gluino und neutralino unterscheiden
- Den(i)=(tman - MNeu(i)**2) !+ cI*MNeu(i)*wNeu(i)
- enddo
stst2qq_M2Tree = 0d0
-! write(*,*)'TEST'
-! write(*,*)'M2Tree', tman, MNeu(1),MNeu(2),MNeu(3),MNeu(4)
+
!Propagator
DenSV(1) = sqrtS**2
DenSV(2) = sqrtS**2
- DenSV(3) = (sqrtS**2-MZ**2)!+ cI*MZ*WZ
- DenSV(4) = (sqrtS**2-MW**2)
+ DenSV(3) = (sqrtS**2-MZ**2)+ cI*MZ*WZ
+ DenSV(4) = (sqrtS**2-MW**2)+ cI*MW*WW
do i=1,6
- DensH(i)=(sqrtS**2-MHk(i)**2)!+ cI*MHk(i)*wHk(i)
+ DensH(i)=(sqrtS**2-MHk(i)**2)+ cI*MHk(i)*wHk(i)
enddo
do i=1,4 !hier gluino und neutralino unterscheiden
- DenGauT(i)=(tman - MNeu(i)**2) !+ cI*MNeu(i)*wNeu(i)
+ DenGauT(i)=(tman - MNeu(i)**2) + cI*MNeu(i)*wNeu(i)
enddo
- DenGauT(5)=(tman - MCha(1)**2) !+ cI*MCha(1)*wCha(1)
- DenGauT(6)=(tman - MCha(2)**2) !+ cI*MCha(2)*wCha(2)
+ DenGauT(5)=(tman - MCha(1)**2) + cI*MCha(1)*wCha(1)
+ DenGauT(6)=(tman - MCha(2)**2) + cI*MCha(2)*wCha(2)
DenGlT = (tman-MGl**2)
-!Farbe:
- !N als NC, CF, CV ist CA
-! MSf(1,3,3) = 342.4197
+! ---------------------Amplitudensquare for s-channel with gluon-----------
-
-
-! write(*,*)'MGL', MGl
-
-! write(*,*) 'wNe',wNeu(1)
-! write(*,*) 'WSF',WSF(1,3,3)
-! write(*,*) 'MSF',MSf(1,3,3)
-! write(*,*) 'CF', CF, (Nc**2-1)/(2d0*Nc)
-! write(*,*) 'Zc11',USf(1,1,3,2)
-! write(*,*) 'Zc12', USf(1,2,3,2)
-! write(*,*) 'Zc21', USf(2,1,3,2)
-! write(*,*) 'Zc22', USf(2,2,3,2)
-
-
-! ! ---------------------Amplitudensquare for s-channel with gluon-----------
-!
n=1
k=1
stst2qq_M2Tree = stst2qq_M2Tree - 1/4d0*(Nc**2-1)/Nc**2*((sqrtS**4*(BL(k,ftt1,ftt2,fgen1,fgen2)*
- ((muu3**4 - muu4**2 + muu4**4 +
- muu3**2*(-1 + 2*muu4**2 - 4*tred) +
- 4*muu1**2*(muu2**2 - tred) + 4*tred -
- 4*muu2**2*tred - 4*muu4**2*tred + 4*tred**2)*
- BL(n,ftt1,ftt2,fgen1,fgen2) +
- 2*(-1 + 2*muu1**2 + 2*muu2**2)*muu3*muu4*
- BR(n,ftt1,ftt2,fgen1,fgen2)) +
- BR(k,ftt1,ftt2,fgen1,fgen2)*
- (2*(-1 + 2*muu1**2 + 2*muu2**2)*muu3*muu4*
- BL(n,ftt1,ftt2,fgen1,fgen2) +
- (muu3**4 - muu4**2 + muu4**4 +
- muu3**2*(-1 + 2*muu4**2 - 4*tred) +
- 4*muu1**2*(muu2**2 - tred) + 4*tred -
- 4*muu2**2*tred - 4*muu4**2*tred + 4*tred**2)*
- BR(n,ftt1,ftt2,fgen1,fgen2)))*
- dconjg(VSfSf(n,isf1,isf2,itt1,itt2,igen1,igen2))*
- VSfSf(k,isf1,isf2,itt1,itt2,igen1,igen2))/
- (dconjg(DenSV(n))*DenSV(k)))
-! ! !
-! ! ! ! ! !---------------------Amplitudensquare for s-channel with Vector-Boson-----------
-!
+
+!---------------------Amplitudensquare for s-channel with Vector-Boson-----------
+
do k=bmin,bmax
do n=cmin,cmax
-! write(*,*)k,n
-! k=3
-! n=3
-! k=4
-! n=4
-! !
+
stst2qq_M2Tree = stst2qq_M2Tree -((sqrtS**4*(BL(k,ftt1,ftt2,fgen1,fgen2)*
- ((muu3**4 - muu4**2 + muu4**4 +
- muu3**2*(-1 + 2*muu4**2 - 4*tred) +
- 4*muu1**2*(muu2**2 - tred) + 4*tred -
- 4*muu2**2*tred - 4*muu4**2*tred + 4*tred**2)*
- BL(n,ftt1,ftt2,fgen1,fgen2) +
- 2*(-1 + 2*muu1**2 + 2*muu2**2)*muu3*muu4*
- BR(n,ftt1,ftt2,fgen1,fgen2)) +
- BR(k,ftt1,ftt2,fgen1,fgen2)*
- (2*(-1 + 2*muu1**2 + 2*muu2**2)*muu3*muu4*
- BL(n,ftt1,ftt2,fgen1,fgen2) +
- (muu3**4 - muu4**2 + muu4**4 +
- muu3**2*(-1 + 2*muu4**2 - 4*tred) +
- 4*muu1**2*(muu2**2 - tred) + 4*tred -
- 4*muu2**2*tred - 4*muu4**2*tred + 4*tred**2)*
- BR(n,ftt1,ftt2,fgen1,fgen2)))*
- dconjg(VSfSf(n,isf1,isf2,itt1,itt2,igen1,igen2))*
- VSfSf(k,isf1,isf2,itt1,itt2,igen1,igen2))/
- (dconjg(DenSV(n))*DenSV(k)))
-!
-! write(*,*)'Den',DenSV(n)
-! write(*,*)'s-MW**2',((sqrtS**2-MW**2))
+
enddo
enddo
-!
-! !!---------------------------- S-channel Higgs-----------------------------
-!
+
+!---------------------------- S-channel Higgs-----------------------------
+
do k=imin,imax
do n=jmin,jmax
-! k=5
-! n=5
+
stst2qq_M2Tree = stst2qq_M2Tree -((sqrtS**2*dconjg(HSfSf(n,isf1,isf2,itt1,itt2,igen1,igen2))*
- (dconjg(DL(n,ftt1,ftt2,fgen1,fgen2))*
- ((-1 + muu3**2 + muu4**2)*
- DL(k,ftt1,ftt2,fgen1,fgen2) +
- 2*muu3*muu4*DR(k,ftt1,ftt2,fgen1,fgen2)) +
- dconjg(DR(n,ftt1,ftt2,fgen1,fgen2))*
- (2*muu3*muu4*DL(k,ftt1,ftt2,fgen1,fgen2) +
- (-1 + muu3**2 + muu4**2)*
- DR(k,ftt1,ftt2,fgen1,fgen2)))*
- HSfSf(k,isf1,isf2,itt1,itt2,igen1,igen2))/
- (dconjg(DenSH(n))*DenSH(k)))
-!
-! write(*,*)'HSfSf',HSfSf(5,1,1,3,4,3,3)
-! write(*,*)'HSfSf_hier', HSfSf(n,isf1,isf2,itt1,itt2,igen1,igen2)
-! write(*,*)n,isf1,isf2,itt1,itt2,igen1,igen2
-!
-! ! write(*,*)'Mb',Mf(4,3)
+
+
enddo
enddo
-!
-! !
-! !
-! ! !
-! ! ! ! !-----------------------T-channel Neutralino-------------------------------------------------------
-! ! !
-! write(*,*)pmin,pmax,qmin,qmax
+
+!-----------------------T-channel Neutralino-------------------------------------------------------
do n=pmin,pmax
do k=qmin,qmax
-! n=1
-! k=1
-! do n=2,4
-! do k=2,4
-! write(*,*)'muu4',Mb/sqrtS,muu4
-! write(*,*)'muu3',Mt/sqrtS,muu3
+
stst2qq_M2Tree = stst2qq_M2Tree -((sqrtS**4*(ALt2(k,isf2,itt2,ftt2,igen2,fgen2)*
- (ARt1(k,isf1,itt1,ftt1,igen1,fgen1)*
- (dconjg(ALt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- ((muu3**2*muu4**2 +
- muu1**2*
- (muu2**2 - muu4**2 - tred) +
- tred + tred**2 -
- muu2**2*(muu3**2 + tred))*
- dconjg(ARt1(n,isf1,itt1,ftt1,
- igen1,fgen1)) +
- muu3*(-muu2**2 + muu4**2 + tred)*
- dconjg(ALt1(n,isf1,itt1,ftt1,
- igen1,fgen1))*muupt(n)) +
- muu4*
- dconjg(ARt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (2*muu3*tred*
- dconjg(ALt1(n,isf1,itt1,ftt1,
- igen1,fgen1)) +
- (-muu1**2 + muu3**2 + tred)*
- dconjg(ARt1(n,isf1,itt1,ftt1,
- igen1,fgen1))*muupt(n))) +
- ALt1(k,isf1,itt1,ftt1,igen1,fgen1)*
- (dconjg(ALt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (muu3*(-muu2**2 + muu4**2 + tred)*
- dconjg(ARt1(n,isf1,itt1,ftt1,
- igen1,fgen1)) +
- (-1 + muu3**2 + muu4**2)*
- dconjg(ALt1(n,isf1,itt1,ftt1,
- igen1,fgen1))*muupt(n)) +
- muu4*
- dconjg(ARt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- ((-muu1**2 + muu3**2 + tred)*
- dconjg(ALt1(n,isf1,itt1,ftt1,
- igen1,fgen1)) +
- 2*muu3*
- dconjg(ARt1(n,isf1,itt1,ftt1,
- igen1,fgen1))*muupt(n)))*
- muuqt(k)) +
- ARt2(k,isf2,itt2,ftt2,igen2,fgen2)*
- (ALt1(k,isf1,itt1,ftt1,igen1,fgen1)*
- (muu4*
- dconjg(ALt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (2*muu3*tred*
- dconjg(ARt1(n,isf1,itt1,ftt1,
- igen1,fgen1)) +
- (-muu1**2 + muu3**2 + tred)*
- dconjg(ALt1(n,isf1,itt1,ftt1,
- igen1,fgen1))*muupt(n)) +
- dconjg(ARt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- ((muu3**2*muu4**2 +
- muu1**2*
- (muu2**2 - muu4**2 - tred) +
- tred + tred**2 -
- muu2**2*(muu3**2 + tred))*
- dconjg(ALt1(n,isf1,itt1,ftt1,
- igen1,fgen1)) +
- muu3*(-muu2**2 + muu4**2 + tred)*
- dconjg(ARt1(n,isf1,itt1,ftt1,
- igen1,fgen1))*muupt(n))) +
- ARt1(k,isf1,itt1,ftt1,igen1,fgen1)*
- (muu4*
- dconjg(ALt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- ((-muu1**2 + muu3**2 + tred)*
- dconjg(ARt1(n,isf1,itt1,ftt1,
- igen1,fgen1)) +
- 2*muu3*
- dconjg(ALt1(n,isf1,itt1,ftt1,
- igen1,fgen1))*muupt(n)) +
- dconjg(ARt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (muu3*(-muu2**2 + muu4**2 + tred)*
- dconjg(ALt1(n,isf1,itt1,ftt1,
- igen1,fgen1)) +
- (-1 + muu3**2 + muu4**2)*
- dconjg(ARt1(n,isf1,itt1,ftt1,
- igen1,fgen1))*muupt(n)))*
- muuqt(k))))/
- (dconjg(DenGauT(n))*DenGauT(k)))
- write(*,*)'stst2qq_M2Tree',stst2qq_M2Tree
-
-
-!Karols Version:
-!
-! stst2qq_M2Tree = stst2qq_M2Tree -((sqrtS**4*
-! - (ALt1(k,isf1,itt1,ftt1,igen1,fgen1)*
-! - (dconjg(ALt1(n,isf1,itt1,ftt1,igen1,fgen1))*
-! - (ARt2(k,isf2,itt2,ftt2,igen2,fgen2)*
-! - ((muu3**2*muu4**2 + muu1**2*(muu2**2 - muu4**2 - tred) + tred +
-! - tred**2 - muu2**2*(muu3**2 + tred))*
-! - dconjg(ARt2(n,isf2,itt2,ftt2,igen2,fgen2)) +
-! - muu4*(-muu1**2 + muu3**2 + tred)*
-! - dconjg(ALt2(n,isf2,itt2,ftt2,igen2,fgen2))*muupt(n)) +
-! - ALt2(k,isf2,itt2,ftt2,igen2,fgen2)*
-! - (muu4*(-muu1**2 + muu3**2 + tred)*
-! - dconjg(ARt2(n,isf2,itt2,ftt2,igen2,fgen2)) +
-! - (-1 + muu3**2 + muu4**2)*
-! - dconjg(ALt2(n,isf2,itt2,ftt2,igen2,fgen2))*muupt(n))*muuqt(k)
-! - ) + muu3*dconjg(ARt1(n,isf1,itt1,ftt1,igen1,fgen1))*
-! - (ARt2(k,isf2,itt2,ftt2,igen2,fgen2)*
-! - (2*muu4*tred*dconjg(ALt2(n,isf2,itt2,ftt2,igen2,fgen2)) +
-! - (-muu2**2 + muu4**2 + tred)*
-! - dconjg(ARt2(n,isf2,itt2,ftt2,igen2,fgen2))*muupt(n)) +
-! - ALt2(k,isf2,itt2,ftt2,igen2,fgen2)*
-! - ((-muu2**2 + muu4**2 + tred)*
-! - dconjg(ALt2(n,isf2,itt2,ftt2,igen2,fgen2)) +
-! - 2*muu4*dconjg(ARt2(n,isf2,itt2,ftt2,igen2,fgen2))*muupt(n))*
-! - muuqt(k))) + ARt1(k,isf1,itt1,ftt1,igen1,fgen1)*
-! - (muu3*dconjg(ALt1(n,isf1,itt1,ftt1,igen1,fgen1))*
-! - (ALt2(k,isf2,itt2,ftt2,igen2,fgen2)*
-! - (2*muu4*tred*dconjg(ARt2(n,isf2,itt2,ftt2,igen2,fgen2)) +
-! - (-muu2**2 + muu4**2 + tred)*
-! - dconjg(ALt2(n,isf2,itt2,ftt2,igen2,fgen2))*muupt(n)) +
-! - ARt2(k,isf2,itt2,ftt2,igen2,fgen2)*
-! - ((-muu2**2 + muu4**2 + tred)*
-! - dconjg(ARt2(n,isf2,itt2,ftt2,igen2,fgen2)) +
-! - 2*muu4*dconjg(ALt2(n,isf2,itt2,ftt2,igen2,fgen2))*muupt(n))*
-! - muuqt(k)) + dconjg(ARt1(n,isf1,itt1,ftt1,igen1,fgen1))*
-! - (ALt2(k,isf2,itt2,ftt2,igen2,fgen2)*
-! - ((muu3**2*muu4**2 + muu1**2*(muu2**2 - muu4**2 - tred) + tred +
-! - tred**2 - muu2**2*(muu3**2 + tred))*
-! - dconjg(ALt2(n,isf2,itt2,ftt2,igen2,fgen2)) +
-! - muu4*(-muu1**2 + muu3**2 + tred)*
-! - dconjg(ARt2(n,isf2,itt2,ftt2,igen2,fgen2))*muupt(n)) +
-! - ARt2(k,isf2,itt2,ftt2,igen2,fgen2)*
-! - (muu4*(-muu1**2 + muu3**2 + tred)*
-! - dconjg(ALt2(n,isf2,itt2,ftt2,igen2,fgen2)) +
-! - (-1 + muu3**2 + muu4**2)*
-! - dconjg(ARt2(n,isf2,itt2,ftt2,igen2,fgen2))*muupt(n))*muuqt(k)
-! - ))))/(dconjg(DenGauT(n))*DenGauT(k)))
-!
-! ! write(*,*)'DenGauT(n)',DenGauT(n),n
-! ! write(*,*)'Den',Den(n)
-! ! write(*,*)'(tman - MNeu(i)**2)',(tman - MNeu(n)**2)
-! write(*,*)'USF', USf(1,1,3,1)
-! write(*,*)'MC',Mf(3,2)
-!
+
enddo
enddo
-!
-! ! write(*,*)'AL',ALt2(k,isf2,itt2,ftt2,igen2,fgen2), ALt1(n,isf1,itt1,ftt1,igen1,fgen1)
-! ! write(*,*)'AR',ARt2(n,isf2,itt2,ftt2,igen2,fgen2), ARt1(k,isf1,itt1,ftt1,igen1,fgen1)
-!
-! ! -----------------------T-channel Gluino-------------------------------------------------------
+
+! -----------------------T-channel Gluino-------------------------------------------------------
stst2qq_M2Tree = stst2qq_M2Tree - 1/4d0*(Nc**2-1)/Nc**2*
- ((sqrtS**4*(dconjg(GluRt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2))*
- (muu3*
- GluLt1(isf1,itt1,ftt1,igen1,fgen1)*
- (muugl*(-muu2**2 + muu4**2 + tred)*
- GluLt2(isf2,itt2,ftt2,igen2,
- fgen2) +
- 2*muu4*tred*
- GluRt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- GluRt1(isf1,itt1,ftt1,igen1,fgen1)*
- ((muu3**2*muu4**2 +
- muu1**2*
- (muu2**2 - muu4**2 - tred) +
- tred + tred**2 -
- muu2**2*(muu3**2 + tred))*
- GluLt2(isf2,itt2,ftt2,igen2,
- fgen2) +
- muu4*muugl*
- (-muu1**2 + muu3**2 + tred)*
- GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))) +
- muugl*
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))*
- (GluRt1(isf1,itt1,ftt1,igen1,fgen1)*
- (muu4*(-muu1**2 + muu3**2 + tred)*
- GluLt2(isf2,itt2,ftt2,igen2,
- fgen2) +
- (-1 + muu3**2 + muu4**2)*muugl*
- GluRt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- muu3*
- GluLt1(isf1,itt1,ftt1,igen1,fgen1)*
- (2*muu4*muugl*
- GluLt2(isf2,itt2,ftt2,igen2,
- fgen2) +
- (-muu2**2 + muu4**2 + tred)*
- GluRt2(isf2,itt2,ftt2,igen2,
- fgen2)))) +
- dconjg(GluLt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (muugl*dconjg(GluLt2(isf2,itt2,ftt2,
- igen2,fgen2))*
- (muu3*
- GluRt1(isf1,itt1,ftt1,igen1,fgen1)*
- ((-muu2**2 + muu4**2 + tred)*
- GluLt2(isf2,itt2,ftt2,igen2,
- fgen2) +
- 2*muu4*muugl*
- GluRt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- GluLt1(isf1,itt1,ftt1,igen1,fgen1)*
- ((-1 + muu3**2 + muu4**2)*muugl*
- GluLt2(isf2,itt2,ftt2,igen2,
- fgen2) +
- muu4*(-muu1**2 + muu3**2 + tred)*
- GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))) +
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))*
- (muu3*
- GluRt1(isf1,itt1,ftt1,igen1,fgen1)*
- (2*muu4*tred*
- GluLt2(isf2,itt2,ftt2,igen2,
- fgen2) +
- muugl*
- (-muu2**2 + muu4**2 + tred)*
- GluRt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- GluLt1(isf1,itt1,ftt1,igen1,fgen1)*
- (muu4*muugl*
- (-muu1**2 + muu3**2 + tred)*
- GluLt2(isf2,itt2,ftt2,igen2,
- fgen2) +
- (muu3**2*muu4**2 +
- muu1**2*
- (muu2**2 - muu4**2 - tred) +
- tred + tred**2 -
- muu2**2*(muu3**2 + tred))*
- GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))))))/
- (DenGlT*dconjg(DenGlT)))
- write(*,*)'GluR',GluR(isf2,itt2,ftt2,igen2,fgen2)
- write(*,*)'GluL',GluL(isf2,itt2,ftt2,igen2,fgen2)
-
-
-
-
-
-
-
! ------------------------------INTERFERENZEN-----------------------------------------------
-!
+
! ---------------------Interferenz S-Kanal Gluon- Vectorboson------------------------------
-!
do k=bmin,bmax
n=1
stst2qq_M2Tree = stst2qq_M2Tree + 0d0*(-2*sqrtS**4*
- (muu3**4 - muu2**2*muu4**2 +
- muu3**2*(-1 + muu2**2 + muu4**2 - 3*tred) +
- 2*muu1**2*(muu2**2 - tred) + 2*tred -
- 2*muu2**2*tred - muu4**2*tred + 2*tred**2)*
- (BL(k,ftt1,ftt2,fgen1,fgen2)*
- BL(n,ftt1,ftt2,fgen1,fgen2) +
- BR(k,ftt1,ftt2,fgen1,fgen2)*
- BR(n,ftt1,ftt2,fgen1,fgen2))*
- dconjg(VSfSf(n,isf1,isf2,itt1,itt2,igen1,igen2))*
- VSfSf(k,isf1,isf2,itt1,itt2,igen1,igen2))/
- (dconjg(DenSV(n))*DenSV(k))
enddo
-! ! !
-! !
-! !
-! !
-! !
-! !
-! ! ! ! !---------------------Interferenz for s-channel with Vector-Boson and Higgs-----------
-! ! ! !
+
+!---------------------Interferenz for s-channel with Vector-Boson and Higgs-----------
+
do k=bmin,bmax
do n=jmin,jmax
stst2qq_M2Tree = stst2qq_M2Tree - (-2*sqrtS**3*(BL(k,ftt1,ftt2,fgen1,fgen2)*
- (muu3*(-1 + 2*muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(DL(n,ftt1,ftt2,fgen1,fgen2)) +
- muu4*(-1 + 2*muu1**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(DR(n,ftt1,ftt2,fgen1,fgen2))) +
- BR(k,ftt1,ftt2,fgen1,fgen2)*
- (muu4*(-1 + 2*muu1**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(DL(n,ftt1,ftt2,fgen1,fgen2)) +
- muu3*(-1 + 2*muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(DR(n,ftt1,ftt2,fgen1,fgen2))))*
- dconjg(HSfSf(n,isf1,isf2,itt1,itt2,igen1,
- igen2))*VSfSf(k,isf1,isf2,itt1,itt2,igen1,
- igen2))/(dconjg(DenSH(n))*DenSV(k))
enddo
enddo
-
-!
-!---------------------Interferenz for s-channel with Vector-Boson and t channel with Neutralino-----------
+!---------------------Interferenz for s-channel with Vector-Boson and t channel with Neutralino-----------
+!
do k=bmin,bmax
do n=qmin,qmax
stst2qq_M2Tree = stst2qq_M2Tree + 1/Nc*(2*sqrtS**4*(BL(k,ftt1,ftt2,fgen1,fgen2)*
- (dconjg(ALt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (((-2 + muu3**2 + muu4**2 - 2*tred)*
- tred +
- muu2**2*(muu3**2 + 2*tred) +
- muu1**2*
- (-2*muu2**2 + muu4**2 + 2*tred))*
- dconjg(ARt1(n,isf1,itt1,ftt1,igen1,
- fgen1)) +
- muu3*
- (-1 + 2*muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(ALt1(n,isf1,itt1,ftt1,igen1,
- fgen1))*muupt(n)) +
- muu4*dconjg(ARt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (muu3*(-muu1**2 - muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(ALt1(n,isf1,itt1,ftt1,igen1,
- fgen1)) +
- (-1 + 2*muu1**2 + muu3**2 + muu4**2 -
- 2*tred)*
- dconjg(ARt1(n,isf1,itt1,ftt1,igen1,
- fgen1))*muupt(n))) +
- BR(k,ftt1,ftt2,fgen1,fgen2)*
- (muu4*dconjg(ALt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (muu3*(-muu1**2 - muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(ARt1(n,isf1,itt1,ftt1,igen1,
- fgen1)) +
- (-1 + 2*muu1**2 + muu3**2 + muu4**2 -
- 2*tred)*
- dconjg(ALt1(n,isf1,itt1,ftt1,igen1,
- fgen1))*muupt(n)) +
- dconjg(ARt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (((-2 + muu3**2 + muu4**2 - 2*tred)*
- tred +
- muu2**2*(muu3**2 + 2*tred) +
- muu1**2*
- (-2*muu2**2 + muu4**2 + 2*tred))*
- dconjg(ALt1(n,isf1,itt1,ftt1,igen1,
- fgen1)) +
- muu3*
- (-1 + 2*muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(ARt1(n,isf1,itt1,ftt1,igen1,
- fgen1))*muupt(n))))*
- VSfSf(k,isf1,isf2,itt1,itt2,igen1,igen2))/
- (dconjg(DenGauT(n))*DenSV(k))
-
+!
enddo
enddo
-!
! !
-! ! !---------------------Interferenz for s-channel with Vector-Boson and t channel with Gluino-----------
+!---------------------Interferenz for s-channel with Vector-Boson and t channel with Gluino-----------
do k=bmin,bmax
-! k=3
+
stst2qq_M2Tree = stst2qq_M2Tree + (1/2d0*(Nc**2-1)/Nc**2)*(2*sqrtS**4*(BL(k,ftt1,ftt2,fgen1,fgen2)*
- (dconjg(GluRt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (((-2 + muu3**2 + muu4**2 - 2*tred)*
- tred +
- muu2**2*(muu3**2 + 2*tred) +
- muu1**2*
- (-2*muu2**2 + muu4**2 + 2*tred))*
- dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- muu4*muugl*
- (-1 + 2*muu1**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))) +
- muu3*dconjg(GluLt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (muugl*
- (-1 + 2*muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- muu4*
- (-muu1**2 - muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2)))) +
- BR(k,ftt1,ftt2,fgen1,fgen2)*
- (muu3*dconjg(GluRt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (muu4*(-muu1**2 - muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- muugl*
- (-1 + 2*muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))) +
- dconjg(GluLt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (muu4*muugl*
- (-1 + 2*muu1**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- ((-2 + muu3**2 + muu4**2 - 2*tred)*
- tred +
- muu2**2*(muu3**2 + 2*tred) +
- muu1**2*
- (-2*muu2**2 + muu4**2 + 2*tred))*
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2)))))*
- VSfSf(k,isf1,isf2,itt1,itt2,igen1,igen2))/
- (dconjg(DenGlT)*DenSV(k))
enddo
-! !
-!
-! ! !---------------------Interferenz for s-channel with Gluon and Higgs-----------
+
+!---------------------Interferenz for s-channel with Gluon and Higgs-----------
k=1
do n=jmin,jmax
-! n=5
stst2qq_M2Tree = stst2qq_M2Tree + 0d0*((sqrtS**3*
- (BL(k,ftt1,ftt2,fgen1,fgen2)*
- (muu3*
- (-1 + 2*muu2**2 +
- muu3**2 + muu4**2 -
- 2*tred)*
- dconjg(DL(n,ftt1,ftt2,
- fgen1,fgen2)) +
- muu4*
- (-1 + 2*muu1**2 +
- muu3**2 + muu4**2 -
- 2*tred)*
- dconjg(DR(n,ftt1,ftt2,
- fgen1,fgen2))) +
- BR(k,ftt1,ftt2,fgen1,
- fgen2)*
- (muu4*
- (-1 + 2*muu1**2 +
- muu3**2 + muu4**2 -
- 2*tred)*
- dconjg(DL(n,ftt1,ftt2,
- fgen1,fgen2)) +
- muu3*
- (-1 + 2*muu2**2 +
- muu3**2 + muu4**2 -
- 2*tred)*
- dconjg(DR(n,ftt1,ftt2,
- fgen1,fgen2))))*
- dconjg(HSfSf(n,isf1,isf2,
- itt1,itt2,igen1,igen2))*
- VSfSf(k,isf1,isf2,itt1,itt2,
- igen1,igen2))/
- (dconjg(DenSH(n))*DenSV(k)))
enddo
-! ! ! !
+
! !---------------------Interferenz for s-channel with Higgs and t-channel Neutralino-----------
!
do k=imin,imax
-! k=5
do n=qmin,qmax
stst2qq_M2Tree = stst2qq_M2Tree - 1/Nc*(2*sqrtS**3*HSfSf(k,isf1,isf2,itt1,itt2,igen1,igen2)*
- (dconjg(ALt2(n,isf2,itt2,ftt2,igen2,fgen2))*
- (dconjg(ARt1(n,isf1,itt1,ftt1,igen1,fgen1))*
- (muu3*(-muu2**2 + muu4**2 + tred)*
- DL(k,ftt1,ftt2,fgen1,fgen2) +
- muu4*(-muu1**2 + muu3**2 + tred)*
- DR(k,ftt1,ftt2,fgen1,fgen2)) +
- dconjg(ALt1(n,isf1,itt1,ftt1,igen1,fgen1))*
- ((-1 + muu3**2 + muu4**2)*
- DL(k,ftt1,ftt2,fgen1,fgen2) +
- 2*muu3*muu4*
- DR(k,ftt1,ftt2,fgen1,fgen2))*muupt(n)) +
- dconjg(ARt2(n,isf2,itt2,ftt2,igen2,fgen2))*
- (dconjg(ALt1(n,isf1,itt1,ftt1,igen1,fgen1))*
- (muu4*(-muu1**2 + muu3**2 + tred)*
- DL(k,ftt1,ftt2,fgen1,fgen2) +
- muu3*(-muu2**2 + muu4**2 + tred)*
- DR(k,ftt1,ftt2,fgen1,fgen2)) +
- dconjg(ARt1(n,isf1,itt1,ftt1,igen1,fgen1))*
- (2*muu3*muu4*
- DL(k,ftt1,ftt2,fgen1,fgen2) +
- (-1 + muu3**2 + muu4**2)*
- DR(k,ftt1,ftt2,fgen1,fgen2))*muupt(n))))/
- (dconjg(DenGauT(n))*DenSH(k))
enddo
enddo
-!
-!---------------------Interferenz for s-channel with Higgs and t-channel Gluino-----------
+
+! !---------------------Interferenz for s-channel with Higgs and t-channel Gluino-----------
do k=imin,imax
-! k=5
stst2qq_M2Tree = stst2qq_M2Tree - 1/(Nc)*CF*(2*sqrtS**3*(dconjg(GluRt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (muugl*dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))*
- (2*muu3*muu4*
- DL(k,ftt1,ftt2,fgen1,fgen2) +
- (-1 + muu3**2 + muu4**2)*
- DR(k,ftt1,ftt2,fgen1,fgen2)) +
- dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2))*
- (muu3*(-muu2**2 + muu4**2 + tred)*
- DL(k,ftt1,ftt2,fgen1,fgen2) +
- muu4*(-muu1**2 + muu3**2 + tred)*
- DR(k,ftt1,ftt2,fgen1,fgen2))) +
- dconjg(GluLt1(isf1,itt1,ftt1,igen1,fgen1))*
- (muugl*dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2))*
- ((-1 + muu3**2 + muu4**2)*
- DL(k,ftt1,ftt2,fgen1,fgen2) +
- 2*muu3*muu4*
- DR(k,ftt1,ftt2,fgen1,fgen2)) +
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))*
- (muu4*(-muu1**2 + muu3**2 + tred)*
- DL(k,ftt1,ftt2,fgen1,fgen2) +
- muu3*(-muu2**2 + muu4**2 + tred)*
- DR(k,ftt1,ftt2,fgen1,fgen2))))*
- HSfSf(k,isf1,isf2,itt1,itt2,igen1,igen2))/
- (dconjg(DenGlT)*DenSH(k))
enddo
-!
-
-!---------------------Interferenz for s-channel with Gluon and t channel with Neutralino-----------
-
+!
+! !---------------------Interferenz for s-channel with Gluon and t channel with Neutralino-----------
+!
k=1
do n=qmin,qmax
stst2qq_M2Tree = stst2qq_M2Tree + (1/2d0*(Nc**2-1)/Nc**2)*(2*sqrtS**4*(BL(k,ftt1,ftt2,fgen1,fgen2)*
- (dconjg(ALt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (((-2 + muu3**2 + muu4**2 - 2*tred)*
- tred +
- muu2**2*(muu3**2 + 2*tred) +
- muu1**2*
- (-2*muu2**2 + muu4**2 + 2*tred))*
- dconjg(ARt1(n,isf1,itt1,ftt1,igen1,
- fgen1)) +
- muu3*
- (-1 + 2*muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(ALt1(n,isf1,itt1,ftt1,igen1,
- fgen1))*muupt(n)) +
- muu4*dconjg(ARt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (muu3*(-muu1**2 - muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(ALt1(n,isf1,itt1,ftt1,igen1,
- fgen1)) +
- (-1 + 2*muu1**2 + muu3**2 + muu4**2 -
- 2*tred)*
- dconjg(ARt1(n,isf1,itt1,ftt1,igen1,
- fgen1))*muupt(n))) +
- BR(k,ftt1,ftt2,fgen1,fgen2)*
- (muu4*dconjg(ALt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (muu3*(-muu1**2 - muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(ARt1(n,isf1,itt1,ftt1,igen1,
- fgen1)) +
- (-1 + 2*muu1**2 + muu3**2 + muu4**2 -
- 2*tred)*
- dconjg(ALt1(n,isf1,itt1,ftt1,igen1,
- fgen1))*muupt(n)) +
- dconjg(ARt2(n,isf2,itt2,ftt2,igen2,
- fgen2))*
- (((-2 + muu3**2 + muu4**2 - 2*tred)*
- tred +
- muu2**2*(muu3**2 + 2*tred) +
- muu1**2*
- (-2*muu2**2 + muu4**2 + 2*tred))*
- dconjg(ALt1(n,isf1,itt1,ftt1,igen1,
- fgen1)) +
- muu3*
- (-1 + 2*muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(ARt1(n,isf1,itt1,ftt1,igen1,
- fgen1))*muupt(n))))*
- VSfSf(k,isf1,isf2,itt1,itt2,igen1,igen2))/
- (dconjg(DenGauT(n))*DenSV(k))
enddo
-!
!
!---------------------Interferenz for s-channel with Gluon and t channel with Gluino-----------
k=1
stst2qq_M2Tree = stst2qq_M2Tree - 1/(2d0*Nc**2)*CF*(2*sqrtS**4*(BL(k,ftt1,ftt2,fgen1,fgen2)*
- (dconjg(GluRt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (((-2 + muu3**2 + muu4**2 - 2*tred)*
- tred +
- muu2**2*(muu3**2 + 2*tred) +
- muu1**2*
- (-2*muu2**2 + muu4**2 + 2*tred))*
- dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- muu4*muugl*
- (-1 + 2*muu1**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))) +
- muu3*dconjg(GluLt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (muugl*
- (-1 + 2*muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- muu4*
- (-muu1**2 - muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2)))) +
- BR(k,ftt1,ftt2,fgen1,fgen2)*
- (muu3*dconjg(GluRt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (muu4*(-muu1**2 - muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- muugl*
- (-1 + 2*muu2**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2))) +
- dconjg(GluLt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (muu4*muugl*
- (-1 + 2*muu1**2 + muu3**2 +
- muu4**2 - 2*tred)*
- dconjg(GluLt2(isf2,itt2,ftt2,igen2,
- fgen2)) +
- ((-2 + muu3**2 + muu4**2 - 2*tred)*
- tred +
- muu2**2*(muu3**2 + 2*tred) +
- muu1**2*
- (-2*muu2**2 + muu4**2 + 2*tred))*
- dconjg(GluRt2(isf2,itt2,ftt2,igen2,
- fgen2)))))*
- VSfSf(k,isf1,isf2,itt1,itt2,igen1,igen2))/
- (dconjg(DenGlT)*DenSV(k))
-
!
!---------------------Interferenz for t-channel with Neutralino and t channel with Gluino-----------
do k=pmin,pmax
stst2qq_M2Tree = stst2qq_M2Tree + 0d0*(sqrtS**4*(ALt2(k,isf2,itt2,ftt2,igen2,fgen2)*
- (ARt1(k,isf1,itt1,ftt1,igen1,fgen1)*
- (muu3*dconjg(GluLt1(isf1,itt1,ftt1,
- igen1,fgen1))*
- (muugl*(muu2**2 - muu4**2 - tred)*
- dconjg(GluLt2(isf2,itt2,ftt2,
- igen2,fgen2)) -
- 2*muu4*tred*
- dconjg(GluRt2(isf2,itt2,ftt2,
- igen2,fgen2))) -
- dconjg(GluRt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- ((muu3**2*muu4**2 +
- muu1**2*
- (muu2**2 - muu4**2 - tred) +
- tred + tred**2 -
- muu2**2*(muu3**2 + tred))*
- dconjg(GluLt2(isf2,itt2,ftt2,
- igen2,fgen2)) +
- muu4*muugl*
- (-muu1**2 + muu3**2 + tred)*
- dconjg(GluRt2(isf2,itt2,ftt2,
- igen2,fgen2)))) -
- ALt1(k,isf1,itt1,ftt1,igen1,fgen1)*
- (muu3*dconjg(GluRt1(isf1,itt1,ftt1,
- igen1,fgen1))*
- ((-muu2**2 + muu4**2 + tred)*
- dconjg(GluLt2(isf2,itt2,ftt2,
- igen2,fgen2)) +
- 2*muu4*muugl*
- dconjg(GluRt2(isf2,itt2,ftt2,
- igen2,fgen2))) +
- dconjg(GluLt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- ((-1 + muu3**2 + muu4**2)*muugl*
- dconjg(GluLt2(isf2,itt2,ftt2,
- igen2,fgen2)) +
- muu4*(-muu1**2 + muu3**2 + tred)*
- dconjg(GluRt2(isf2,itt2,ftt2,
- igen2,fgen2))))*muuqt(k)) +
- ARt2(k,isf2,itt2,ftt2,igen2,fgen2)*
- (ALt1(k,isf1,itt1,ftt1,igen1,fgen1)*
- (-(muu3*
- dconjg(GluRt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (2*muu4*tred*
- dconjg(GluLt2(isf2,itt2,ftt2,
- igen2,fgen2)) +
- muugl*
- (-muu2**2 + muu4**2 + tred)*
- dconjg(GluRt2(isf2,itt2,ftt2,
- igen2,fgen2)))) +
- dconjg(GluLt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (muu4*muugl*
- (muu1**2 - muu3**2 - tred)*
- dconjg(GluLt2(isf2,itt2,ftt2,
- igen2,fgen2)) -
- (muu3**2*muu4**2 +
- muu1**2*
- (muu2**2 - muu4**2 - tred) +
- tred + tred**2 -
- muu2**2*(muu3**2 + tred))*
- dconjg(GluRt2(isf2,itt2,ftt2,
- igen2,fgen2)))) -
- ARt1(k,isf1,itt1,ftt1,igen1,fgen1)*
- (dconjg(GluRt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (muu4*(-muu1**2 + muu3**2 + tred)*
- dconjg(GluLt2(isf2,itt2,ftt2,
- igen2,fgen2)) +
- (-1 + muu3**2 + muu4**2)*muugl*
- dconjg(GluRt2(isf2,itt2,ftt2,
- igen2,fgen2))) +
- muu3*
- dconjg(GluLt1(isf1,itt1,ftt1,igen1,
- fgen1))*
- (2*muu4*muugl*
- dconjg(GluLt2(isf2,itt2,ftt2,
- igen2,fgen2)) +
- (-muu2**2 + muu4**2 + tred)*
- dconjg(GluRt2(isf2,itt2,ftt2,
- igen2,fgen2))))*muuqt(k))))/
- (dconjg(DenGlT)*DenGauT(k))
enddo
-!
+
end
Index: branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_Couplings.F
===================================================================
--- branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_Couplings.F (revision 1091)
+++ branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_Couplings.F (revision 1092)
@@ -1,745 +1,539 @@
subroutine stst2qq_SetCouplings
implicit none
#include "stst2QQ_Model.h"
#include "stst2QQ_Kinematics.h"
#include "stst2QQ_GenCouplings.h"
#include "stst2QQ_GenCounterterms.h"
integer i, j, k, l, m, n, xtype, ytype, a, b, o, p
double precision Krondelta
! ********************************
! Fermion-Vector boson couplings*
! ********************************
!e.g. BL(i,f3type,f4type,f3gen,f4gen), no generation-mixing, couplings multiplied by -i
!1=Gluon, 2= Photon, 3= Z-Boson, 4= W-Boson
do i = 1, 4
do j = 1, 4
do m = 1, 3
do n = 1, 3
if((i.gt.2).and.(j.gt.2)) then
BL(1,i,j,m,n) = -GS*Krondelta(i,j)*Krondelta(m,n)
BR(1,i,j,m,n) = -GS*Krondelta(i,j)*Krondelta(m,n)
else
BL(1,i,j,m,n) = (0d0,0d0)
BR(1,i,j,m,n) = (0d0,0d0)
endif
BL(2,i,j,m,n) = -ef_(i)*EL*Krondelta(i,j)*Krondelta(m,n)
BR(2,i,j,m,n) = -ef_(i)*EL*Krondelta(i,j)*Krondelta(m,n)
BV(i) = I3f_(i) - 2d0*SW**2*ef_(i)
BA(i) = I3f_(i)
BL(3,i,j,m,n) = - (EL/(2d0*CW*SW))*(BV(i) + BA(i))*Krondelta(i,j)*Krondelta(m,n)
BR(3,i,j,m,n) = - (EL/(2d0*CW*SW))*(BV(i) - BA(i))*Krondelta(i,j)*Krondelta(m,n)
if((i.eq.1).or.(i.eq.3)) then !incoming up type (->W+ -radiation)
BL(4,i,j,m,n) = - EL/(sqrt2*SW)*Krondelta(i,j-1)*Krondelta(m,n)
else !incoming down type (->W- -radiation)
BL(4,i,j,m,n) = - EL/(sqrt2*SW)*Krondelta(i,j+1)*Krondelta(m,n)
endif
BR(4,i,j,m,n) = (0d0,0d0)
-! do k=1,4
-!
-! if((k.gt.3).and.((xtype.eq.4).or.(xtype.eq.8)).and.((ytype.eq.4).or.(ytype.eq.8))) then
-!
-! if((i.eq.1).or.(i.eq.3)) then
-!
-! BLt(k,i,j,m,n) = BL(k,i,j,m,n)
-! BLu(k,i,j,m,n) = 0d0
-!
-! BRt(k,i,j,m,n) = BR(k,i,j,m,n)
-! BRu(k,i,j,m,n) = 0d0
-!
-! else
-!
-! BLu(k,i,j,m,n) = BL(k,i,j,m,n)
-! BLt(k,i,j,m,n) = 0d0
-!
-! BRu(k,i,j,m,n) = BR(k,i,j,m,n)
-! BRt(k,i,j,m,n) = 0d0
-!
-! endif
-!
-! else
-!
-! BLt(k,i,j,m,n) = BL(k,i,j,m,n)
-! BLu(k,i,j,m,n) = BL(k,i,j,m,n)
-!
-! BRt(k,i,j,m,n) = BR(k,i,j,m,n)
-! BRu(k,i,j,m,n) = BR(k,i,j,m,n)
-!
-! endif
-!
-! enddo
-
-
enddo
enddo
enddo
enddo
!**************************
! Fermion-Higgs couplings *
!**************************
!e.g. DL/R(i,f3type,f4type,f3gen,f4gen), couplings multiplied by -i
do i = 1, 4
do j = 1, 4
do m = 1, 3
do n = 1, 3
! DL(1,i,j,m,n) = 1d0
! DR(1,i,j,m,n) = 1d0
!
!--- incoming up-type fermion ---
! hf_ = EL*Mf(type,gen)/(sqrt2*MW*SW*SB)
if ((i.eq.1).or.(i.eq.3)) then
! h
DL(1,i,j,m,n) = -hf_(i,m)*CA/sqrt2*Krondelta(i,j)*Krondelta(m,n)
DR(1,i,j,m,n) = DL(1,i,j,m,n)
! H
DL(2,i,j,m,n) = -hf_(i,m)*SA/sqrt2*Krondelta(i,j)*Krondelta(m,n)
DR(2,i,j,m,n) = DL(2,i,j,m,n)
! H3
DL(3,i,j,m,n) = -cI*hf_(i,m)*CB/sqrt2*Krondelta(i,j)*Krondelta(m,n)
DR(3,i,j,m,n) = -DL(3,i,j,m,n)
! G0
DL(4,i,j,m,n) = -cI*hf_(i,m)*SB/sqrt2*Krondelta(i,j)*Krondelta(m,n)
DR(4,i,j,m,n) = -DL(4,i,j,m,n)
! H+(out)
DL(5,i,j,m,n) = hf_(i+1,m)*SB*Krondelta(i,j-1)*Krondelta(m,n)
DR(5,i,j,m,n) = hf_(i,m)*CB*Krondelta(i,j-1)*Krondelta(m,n)
! G+(out)
DL(6,i,j,m,n) = -hf_(i+1,m)*CB*Krondelta(i,j-1)*Krondelta(m,n)
DR(6,i,j,m,n) = hf_(i,m)*SB*Krondelta(i,j-1)*Krondelta(m,n)
-! do k =1,6
-!
-! if((k.gt.4).and.((xtype.eq.4).or.(xtype.eq.8)).and.((ytype.eq.4).or.(ytype.eq.8))) then
-!
-! DLt(k,i,j,m,n) = DL(k,i,j,m,n)
-! DLu(k,i,j,m,n) = 0d0
-!
-! DRt(k,i,j,m,n) = DR(k,i,j,m,n)
-! DRu(k,i,j,m,n) = 0d0
-!
-! else
-!
-! DLt(k,i,j,m,n) = DL(k,i,j,m,n)
-! DLu(k,i,j,m,n) = DL(k,i,j,m,n)
-!
-! DRt(k,i,j,m,n) = DR(k,i,j,m,n)
-! DRu(k,i,j,m,n) = DR(k,i,j,m,n)
-!
-! endif
-!
-!
-! enddo
-
+!
! --- incoming down-type fermion ---
else
! h
DL(1,i,j,m,n) = hf_(i,m)*SA/sqrt2*Krondelta(i,j)*Krondelta(m,n)
DR(1,i,j,m,n) = DL(1,i,j,m,n)
! H
DL(2,i,j,m,n) = -hf_(i,m)*CA/sqrt2*Krondelta(i,j)*Krondelta(m,n)
DR(2,i,j,m,n) = DL(2,i,j,m,n)
! H3
DL(3,i,j,m,n) = -cI*hf_(i,m)*SB/sqrt2*Krondelta(i,j)*Krondelta(m,n)
DR(3,i,j,m,n) = -DL(3,i,j,m,n)
! G0
DL(4,i,j,m,n) = cI*hf_(i,m)*CB/sqrt2*Krondelta(i,j)*Krondelta(m,n)
DR(4,i,j,m,n) = -DL(4,i,j,m,n)
! H-(out)
DL(5,i,j,m,n) = hf_(i-1,m)*CB*Krondelta(i,j+1)*Krondelta(m,n)
DR(5,i,j,m,n) = hf_(i,m)*SB*Krondelta(i,j+1)*Krondelta(m,n)
! G-(out)
DL(6,i,j,m,n) = hf_(i-1,m)*SB*Krondelta(i,j+1)*Krondelta(m,n)
DR(6,i,j,m,n) = -hf_(i,m)*CB*Krondelta(i,j+1)*Krondelta(m,n)
-! do k =1,6
-!
-!
-! if((k.gt.4).and.((xtype.eq.4).or.(xtype.eq.8)).and.((ytype.eq.4).or.(ytype.eq.8))) then
-!
-! DLu(k,i,j,m,n) = DL(k,i,j,m,n)
-! DLt(k,i,j,m,n) = 0d0
-!
-! DRu(k,i,j,m,n) = DR(k,i,j,m,n)
-! DRt(k,i,j,m,n) = 0d0
-!
-! else
-!
-! DLt(k,i,j,m,n) = DL(k,i,j,m,n)
-! DLu(k,i,j,m,n) = DL(k,i,j,m,n)
-!
-! DRt(k,i,j,m,n) = DR(k,i,j,m,n)
-! DRu(k,i,j,m,n) = DR(k,i,j,m,n)
-!
-! endif
-!
-! enddo
-
endif
enddo
enddo
enddo
enddo
!********************************************************
! Vertex Gluino - Squark (outgoing) - Quark (incoming) *
!********************************************************
! e.g. GluL/R(sfmeigenin,sftypein,ftypeout,sfgenin,fgenout), couplings multiplied by -i
do a=1, 2 !taken from NeuQ2qx + Sparticles p. 88
do b= 1, 4
do m= 1, 4
do i= 1, 3
do n= 1, 3
if((b.gt.2).and.(m.gt.2)) then
GluL(a,b,m,i,n) = - GS*sqrt2*USf(a,1,b,i)*Krondelta(b,m)*Krondelta(i,n)
GluR(a,b,m,i,n) = GS*sqrt2*USf(a,2,b,i)*Krondelta(b,m)*Krondelta(i,n)
else
GluL(a,b,m,i,n) = (0d0,0d0)
GluR(a,b,m,i,n) = (0d0,0d0)
endif
enddo
enddo
enddo
enddo
enddo
! ***********************************************
!Sfermion - Sfermion - Gluon / Photon-couplings *
!************************************************
! e.g. C(xtype,sftype1)
! e.g. VSfSf(xtype,sfmeigen1,sfmeigen2,sftype1,sftype2,sfgen1,sfgen2) / general structure (s-channel):VSfSf*(p1-p2), couplings multiplied by -i
!1=Gluon (=g), 2=Gamma (=G)
C(1,1) = (0d0,0d0)
C(1,2) = (0d0,0d0)
-c C(1,3) = -findValW("SC")
-c C(1,4) = -findValW("SC")
C(1,3) = -GS
C(1,4) = -GS
C(2,1) = -ef_(1)*EL
C(2,2) = -ef_(2)*EL
C(2,3) = -ef_(3)*EL
C(2,4) = -ef_(4)*EL
do i = 1, 2
do j = 1, 2
do k = 1, 4
do l = 1, 4
do m = 1, 3
do n = 1, 3
VSfSf(1,i,j,k,l,m,n) = C(1,k)*Krondelta(i,j)*Krondelta(k,l)*Krondelta(m,n)
VSfSf(2,i,j,k,l,m,n) = C(2,k)*Krondelta(i,j)*Krondelta(k,l)*Krondelta(m,n)
VSfSft(1,i,j,k,l,m,n) = VSfSf(1,i,j,k,l,m,n)
VSfSfu(1,i,j,k,l,m,n) = VSfSf(1,i,j,k,l,m,n)
VSfSft(2,i,j,k,l,m,n) = VSfSf(2,i,j,k,l,m,n)
VSfSfu(2,i,j,k,l,m,n) = VSfSf(2,i,j,k,l,m,n)
enddo
enddo
enddo
enddo
enddo
enddo
! *****************************************
! Sfermion - Sfermion - W / Z - couplings *
! *****************************************
! e.g. VSfSf(xtype,sfmeigenout,sfmeigenin,sftypeout,sftypein,sfgenout,sfgenin) / general structure (s-channel):VSfSf*(p1-p2), couplings multiplied by -i
!3=Z, 4=W+ (W-)
do i = 1, 2
do j = 1, 2
do k = 1, 4
do l = 1, 4
do m = 1, 3
do n = 1, 3
VSfSf(3,i,j,k,l,m,n)=-EL/(SW*CW)*(I3f_(k)*USf(j,1,l,n)*USf(i,1,k,m)-SW2*ef_(k)*Krondelta(j,i))*Krondelta(k,l)*Krondelta(m,n) !Z-Boson
VSfSft(3,i,j,k,l,m,n) = VSfSf(3,i,j,k,l,m,n)
VSfSfu(3,i,j,k,l,m,n) = VSfSf(3,i,j,k,l,m,n)
! --- incoming up-type sfermion ---
if ((l.eq.1).or.(l.eq.3)) then !W+/- -Boson
VSfSf(4,i,j,k,l,m,n) = - EL/(sqrt2*SW)*USf(j,1,l,n)*USf(i,1,k,m)*Krondelta(l+1,k)*Krondelta(m,n)
-! if(((xtype.eq.4).or.(xtype.eq.8)).and.((ytype.eq.4).or.(ytype.eq.8))) then !u-channel only needs to be switched of for '+-' final-states (e.g. W+H- etc.)
-! VSfSft(4,i,j,k,l,m,n) = VSfSf(4,i,j,k,l,m,n)
-! VSfSfu(4,i,j,k,l,m,n) = (0d0,0d0)
-! else
-! VSfSft(4,i,j,k,l,m,n) = VSfSf(4,i,j,k,l,m,n)
-! VSfSfu(4,i,j,k,l,m,n) = VSfSf(4,i,j,k,l,m,n)
-! endif
! --- incoming down-type sfermion ---
else
VSfSf(4,i,j,k,l,m,n) = - EL/(sqrt2*SW)*USf(j,1,l,n)*USf(i,1,k,m)*Krondelta(l-1,k)*Krondelta(m,n)
-! if(((xtype.eq.4).or.(xtype.eq.8)).and.((ytype.eq.4).or.(ytype.eq.8))) then !u-channel only needs to be switched of for '+-' final-states (e.g. W+H- etc.)
-! VSfSft(4,i,j,k,l,m,n) = (0d0,0d0)
-! VSfSfu(4,i,j,k,l,m,n) = VSfSf(4,i,j,k,l,m,n)
-! else
-! VSfSft(4,i,j,k,l,m,n) = VSfSf(4,i,j,k,l,m,n)
-! VSfSfu(4,i,j,k,l,m,n) = VSfSf(4,i,j,k,l,m,n)
-! endif
endif
enddo
enddo
enddo
enddo
enddo
enddo
!**************************
! Higgs-Sfermion-Sfermion *
!**************************
!e.g. HSfSf(xhtype,sfmeigenin,sfmeigenout,sftypein,sftypeout,sfgenin,sfgenout), couplings multiplied by -i
!From Higgs-Hunter (p.394 f)
do i= 1, 4
do j=1, 3
if ((i.eq.1).or.(i.eq.3)) then
GLR(1,1,1,i,j) = -sqrt2*hf_(i,j)*Mf(i,j)*CA + EL*MZ*(I3f_(i)-ef_(i)*SW2)*SAB/(SW*CW)
GLR(1,2,2,i,j) = -sqrt2*hf_(i,j)*Mf(i,j)*CA + EL*MZ*ef_(i)*SW2*SAB/(SW*CW)
GLR(1,2,1,i,j) = -hf_(i,j)/sqrt2*(Af(i,j)*CA + MUE*SA)
GLR(1,1,2,i,j) = -hf_(i,j)/sqrt2*(Af(i,j)*CA + MUE*SA)
GLR(2,1,1,i,j) = -sqrt2*hf_(i,j)*Mf(i,j)*SA - EL*MZ*(I3f_(i)-ef_(i)*SW2)*CAB/(SW*CW)
GLR(2,2,2,i,j) = -sqrt2*hf_(i,j)*Mf(i,j)*SA - EL*MZ*ef_(i)*SW2*CAB/(SW*CW)
GLR(2,2,1,i,j) = -hf_(i,j)/sqrt2*(Af(i,j)*SA - MUE*CA)
GLR(2,1,2,i,j) = -hf_(i,j)/sqrt2*(Af(i,j)*SA - MUE*CA)
GLR(3,1,1,i,j) = (0d0,0d0)
GLR(3,2,2,i,j) = (0d0,0d0)
GLR(3,2,1,i,j) = cI*hf_(i,j)/sqrt2*(Af(i,j)*CB + MUE*SB)
GLR(3,1,2,i,j) = -cI*hf_(i,j)/sqrt2*(Af(i,j)*CB + MUE*SB)
GLR(4,1,1,i,j) = (0d0,0d0)
GLR(4,2,2,i,j) = (0d0,0d0)
GLR(4,2,1,i,j) = cI*hf_(i,j)/sqrt2*(Af(i,j)*SB - MUE*CB)
GLR(4,1,2,i,j) = -cI*hf_(i,j)/sqrt2*(Af(i,j)*SB - MUE*CB)
! G_{LR}^{tb} where i = up-type
GLR(5,1,1,i,j) = hf_(i+1,j)*Mf(i+1,j)*SB + hf_(i,j)*Mf(i,j)*CB - EL*MW/(sqrt2*SW)*S2B
GLR(5,2,2,i,j) = hf_(i,j)*Mf(i+1,j)*CB + hf_(i+1,j)*Mf(i,j)*SB
GLR(5,2,1,i,j) = hf_(i,j)*(Af(i,j)*CB + MUE*SB)
GLR(5,1,2,i,j) = hf_(i+1,j)*(Af(i+1,j)*SB + MUE*CB)
GLR(6,1,1,i,j) = -hf_(i+1,j)*Mf(i+1,j)*CB + hf_(i,j)*Mf(i,j)*SB + EL*MW/(sqrt2*SW)*C2B
GLR(6,2,2,i,j) = hf_(i,j)*Mf(i+1,j)*SB - hf_(i+1,j)*Mf(i,j)*CB
GLR(6,2,1,i,j) = hf_(i,j)*(Af(i,j)*SB - MUE*CB)
GLR(6,1,2,i,j) = hf_(i+1,j)*(-Af(i+1,j)*CB + MUE*SB)
else
GLR(1,1,1,i,j) = sqrt2*hf_(i,j)*Mf(i,j)*SA + EL*MZ*(I3f_(i)-ef_(i)*SW2)*SAB/(SW*CW)
GLR(1,2,2,i,j) = sqrt2*hf_(i,j)*Mf(i,j)*SA + EL*MZ*ef_(i)*SW2*SAB/(SW*CW)
GLR(1,2,1,i,j) = hf_(i,j)/sqrt2*(Af(i,j)*SA + MUE*CA)
GLR(1,1,2,i,j) = hf_(i,j)/sqrt2*(Af(i,j)*SA + MUE*CA)
GLR(2,1,1,i,j) = -sqrt2*hf_(i,j)*Mf(i,j)*CA - EL*MZ*(I3f_(i)-ef_(i)*SW2)*CAB/(SW*CW)
GLR(2,2,2,i,j) = -sqrt2*hf_(i,j)*Mf(i,j)*CA - EL*MZ*ef_(i)*SW2*CAB/(SW*CW)
GLR(2,2,1,i,j) = -hf_(i,j)/sqrt2*(Af(i,j)*CA - MUE*SA)
GLR(2,1,2,i,j) = -hf_(i,j)/sqrt2*(Af(i,j)*CA - MUE*SA)
GLR(3,1,1,i,j) = (0d0,0d0)
GLR(3,2,2,i,j) = (0d0,0d0)
GLR(3,2,1,i,j) = cI*hf_(i,j)/sqrt2*(Af(i,j)*SB + MUE*CB)
GLR(3,1,2,i,j) = -cI*hf_(i,j)/sqrt2*(Af(i,j)*SB + MUE*CB)
GLR(4,1,1,i,j) = (0d0,0d0)
GLR(4,2,2,i,j) = (0d0,0d0)
GLR(4,2,1,i,j) = -cI*hf_(i,j)/sqrt2*(Af(i,j)*CB - MUE*SB)
GLR(4,1,2,i,j) = cI*hf_(i,j)/sqrt2*(Af(i,j)*CB - MUE*SB)
! G_{LR}^{bt} where i = down-type
GLR(5,1,1,i,j) = hf_(i,j)*Mf(i,j)*SB + hf_(i-1,j)*Mf(i-1,j)*CB - EL*MW/(sqrt2*SW)*S2B
GLR(5,2,2,i,j) = hf_(i-1,j)*Mf(i,j)*CB + hf_(i,j)*Mf(i-1,j)*SB
GLR(5,2,1,i,j) = hf_(i,j)*(Af(i,j)*SB + MUE*CB)
GLR(5,1,2,i,j) = hf_(i-1,j)*(Af(i-1,j)*CB + MUE*SB)
GLR(6,1,1,i,j) = -hf_(i,j)*Mf(i,j)*CB + hf_(i-1,j)*Mf(i-1,j)*SB + EL*MW/(sqrt2*SW)*C2B
GLR(6,2,2,i,j) = hf_(i-1,j)*Mf(i,j)*SB - hf_(i,j)*Mf(i-1,j)*CB
GLR(6,2,1,i,j) = hf_(i,j)*(-Af(i,j)*CB + MUE*SB)
GLR(6,1,2,i,j) = hf_(i-1,j)*(Af(i-1,j)*SB - MUE*CB)
endif
enddo
enddo
! neutral Higgs
do a=1,4
do i=1,2
do j=1,2
do m=1,4
do n=1,4
do o=1,3
do p=1,3
HSfSf(a,i,j,m,n,o,p) = (0d0,0d0)
do k=1,2
do l=1,2
HSfSf(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p) + GLR(a,k,l,m,o)*USf(i,k,m,o)*USf(j,l,n,p)*Krondelta(m,n)*Krondelta(o,p)
enddo
enddo
-!
-
-! HSfSft(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p)
-! HSfSfu(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p)
+
enddo
enddo
enddo
enddo
enddo
enddo
enddo
! charged Higgs
do a=5,6
do i=1,2
do j=1,2
do m=1,4
do n=1,4
do o=1,3
do p=1,3
HSfSf(a,i,j,m,n,o,p) = (0d0,0d0)
do k=1,2
do l=1,2
if((m.eq.1).or.(m.eq.3)) then !up-case
HSfSf(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p) + GLR(a,k,l,m,o)*USf(i,k,m,o)*USf(j,l,n,p)*Krondelta(m+1,n)*Krondelta(o,p)
else !down-case
HSfSf(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p) + GLR(a,k,l,m,o)*USf(i,k,m,o)*USf(j,l,n,p)*Krondelta(m-1,n)*Krondelta(o,p)
endif
enddo
enddo
-!
-! ! if((m.eq.1).or.(m.eq.3)) then
-! !
-! ! if(((xtype.eq.4).or.(xtype.eq.8)).and.((ytype.eq.4).or.(ytype.eq.8))) then !u-channel only needs to be switched of for +- final-states (e.g. W+H- etc.)
-! ! HSfSft(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p)
-! ! HSfSfu(a,i,j,m,n,o,p) = (0d0,0d0)
-! ! else
-! ! HSfSft(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p)
-! ! HSfSfu(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p)
-! ! endif
-! !
-! ! else
-! !
-! ! if(((xtype.eq.4).or.(xtype.eq.8)).and.((ytype.eq.4).or.(ytype.eq.8))) then !u-channel only needs to be switched of for +- final-states (e.g. W+H- etc.)
-! ! HSfSft(a,i,j,m,n,o,p) = (0d0,0d0)
-! ! HSfSfu(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p)
-! ! else
-! ! HSfSft(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p)
-! ! HSfSfu(a,i,j,m,n,o,p) = HSfSf(a,i,j,m,n,o,p)
-! ! endif
-! !
-! ! endif
+
enddo
enddo
enddo
enddo
enddo
enddo
enddo
!
!*************************************
! Gaugino-Fermion-Sfermion couplings *
!*************************************
! e.g. AL/R(chimeigen,sfmeigenin,sftypein,ftypeout,sfgenin,fgenout), couplings multiplied by -i
do a = 1, 2
do b = 1, 4
do m= 1, 4
do i = 1, 3
do n =1, 3
! ------------------------------ Vertex Neutralino - Sfermion (in) - Fermion (out) ------------------------------
do j = 1, 4 !taken from Sparticles p. 517 with GL->AR,GR->AL etc. (following the sfermion(in)-diagramm + ZneuC can be imaginary=> keep track of hermitian conjugation (see Sparticles p. 517))
! --- incoming up-type sfermion ---
if ((b.eq.1).or.(b.eq.3)) then
AR(j,a,b,m,i,n) = (- EL*Mf(b,i)/(SW*sqrt2*MW*SB)*ZNeu(j,4)*USf(a,2,b,i) -
& sqrt2*EL/SW*((ef_(b) - I3f_(b))*SW/CW*ZNeu(j,1) + I3f_(b)*ZNeu(j,2))*USf(a,1,b,i))*Krondelta(b,m)*Krondelta(i,n)
AL(j,a,b,m,i,n) = (- EL*Mf(b,i)/(SW*sqrt2*MW*SB)*ZNeuC(j,4)*USf(a,1,b,i) +
& sqrt2*EL*ef_(b)/CW*ZNeuC(j,1)*USf(a,2,b,i))*Krondelta(b,m)*Krondelta(i,n)
+!
- AR(j,1,b,m,i,n) = 1.0
- AL(j,1,b,m,i,n) = 1.0
ARt1(j,a,b,m,i,n) = AR(j,a,b,m,i,n)
ALt1(j,a,b,m,i,n) = AL(j,a,b,m,i,n)
ARt2(j,a,b,m,i,n) = AL(j,a,b,m,i,n)
ALt2(j,a,b,m,i,n) = AR(j,a,b,m,i,n)
-
+
! --- incoming down-type sfermion---
else
AR(j,a,b,m,i,n) = (- EL*Mf(b,i)/(SW*sqrt2*MW*CB)*ZNeu(j,3)*USf(a,2,b,i) -
& sqrt2*EL/SW*(I3f_(b)*ZNeu(j,2) + (ef_(b) - I3f_(b))*SW/CW*ZNeu(j,1))*USf(a,1,b,i))*Krondelta(b,m)*Krondelta(i,n)
AL(j,a,b,m,i,n) = (- EL*Mf(b,i)/(SW*sqrt2*MW*CB)*ZNeuC(j,3)*USf(a,1,b,i) +
& sqrt2*EL*ef_(b)/CW*ZNeuC(j,1)*USf(a,2,b,i))*Krondelta(b,m)*Krondelta(i,n)
-! AR(j,a,b,m,i,n) = 1.0
-!
-! AL(j,a,b,m,i,n) = 1.0
-
ARt1(j,a,b,m,i,n) = AR(j,a,b,m,i,n)
ALt1(j,a,b,m,i,n) = AL(j,a,b,m,i,n)
ARt2(j,a,b,m,i,n) = AL(j,a,b,m,i,n)
ALt2(j,a,b,m,i,n) = AR(j,a,b,m,i,n)
endif
enddo
! ------------------------------ Vertex Chargino - Sfermion (in) - Fermion (out) ------------------------------ !differs from couplings implemented in CalcHEP->not comparable
do j= 5, 6
if ((b.eq.1).or.(b.eq.3)) then
-
-! KL(j,a,b,i) = -EL/SW*dconjg(UCha(j-4,1))*dconjg(USf(a,1,b+1,i)) + EL*Mf(b+1,i)/(SW*sqrt2*MW*CB)*dconjg(UCha(j-4,2))*dconjg(USf(a,2,b+1,i)) != C_L in Sparticles p. 515
-! KR(j,a,b,i) = EL*Mf(b,i)/(SW*sqrt2*MW*SB)*VCha(j-4,2)*dconjg(USf(a,1,b+1,i)) != F_R in Sparticles p. 515
+
! --- incoming up-type sfermion ---
AL(j,a,b,m,i,n) = EL*Mf(b+1,i)/(SW*sqrt2*MW*CB)*dconjg(UCha(j-4,2))*USf(a,1,b,i)*Krondelta(b,m-1)*Krondelta(i,n) != E_R in Sparticles p. 515 , no t_i-e-coupl. etc. and no generation-mixing
AR(j,a,b,m,i,n) = (-EL/SW*VCha(j-4,1)*USf(a,1,b,i) + EL*Mf(b,i)/(SW*sqrt2*MW*SB)*VCha(j-4,2)*USf(a,2,b,i))*Krondelta(b,m-1)*Krondelta(i,n) != D_L in Sparticles p. 515
! --- incoming up-type sfermion ---
ARt1(j,a,b,m,i,n) = AR(j,a,b,m,i,n)
ALt1(j,a,b,m,i,n) = AL(j,a,b,m,i,n)
! --- outgoing up-type sfermion ---
ARt2(j,a,b,m,i,n) = AL(j,a,b,m,i,n)
ALt2(j,a,b,m,i,n) = AR(j,a,b,m,i,n)
else
-! KL(j,a,b,i) = -EL/SW*dconjg(VCha(j-4,1))*dconjg(USf(a,1,b-1,i)) + EL*Mf(b-1,i)/(SW*sqrt2*MW*SB)*dconjg(VCha(j-4,2))*dconjg(USf(a,2,b-1,i)) != D_L in Sparticles p. 515
-! KR(j,a,b,i) = EL*Mf(b,i)/(SW*sqrt2*MW*CB)*UCha(j-4,2)*dconjg(USf(a,1,b-1,i)) != E_R in Sparticles p. 515
! --- incoming down-type sfermion ---
AL(j,a,b,m,i,n) = EL*Mf(b-1,i)/(SW*sqrt2*MW*SB)*dconjg(VCha(j-4,2))*USf(a,1,b,i)*Krondelta(b,m+1)*Krondelta(i,n) != F_R in Sparticles p. 515
AR(j,a,b,m,i,n) = (-EL/SW*UCha(j-4,1)*USf(a,1,b,i) + EL*Mf(b,i)/(SW*sqrt2*MW*CB)*UCha(j-4,2)*USf(a,2,b,i))*Krondelta(b,m+1)*Krondelta(i,n) != C_L in Sparticles p. 515
! --- incoming down-type sfermion ---
ARt1(j,a,b,m,i,n) = AR(j,a,b,m,i,n)
ALt1(j,a,b,m,i,n) = AL(j,a,b,m,i,n)
! --- outgoing down-type sfermion ---
ARt2(j,a,b,m,i,n) = AL(j,a,b,m,i,n)
ALt2(j,a,b,m,i,n) = AR(j,a,b,m,i,n)
endif
enddo
enddo
enddo
enddo
enddo
enddo
!********************************************************
! Vertex Gluino - Squark (outgoing) - Quark (incoming) *
!********************************************************
! e.g. GluL/R(sfmeigenin,sftypein,ftypeout,sfgenin,fgenout), couplings multiplied by -i
do a=1, 2 !taken from NeuQ2qx + Sparticles p. 88
do b= 1, 4
do m= 1, 4
do i= 1, 3
do n= 1, 3
if((b.gt.2).and.(m.gt.2)) then
GluL(a,b,m,i,n) = - GS*sqrt2*USf(a,1,b,i)*Krondelta(b,m)*Krondelta(i,n)
GluR(a,b,m,i,n) = GS*sqrt2*USf(a,2,b,i)*Krondelta(b,m)*Krondelta(i,n)
GluLt2(a,b,m,i,n) = GluL(a,b,m,i,n)
GluRt2(a,b,m,i,n) = GluR(a,b,m,i,n)
GluLt1(a,b,m,i,n) = GluR(a,b,m,i,n)
GluRt1(a,b,m,i,n) = GluL(a,b,m,i,n)
else
GluL(a,b,m,i,n) = (0d0,0d0)
GluR(a,b,m,i,n) = (0d0,0d0)
endif
enddo
enddo
enddo
enddo
enddo
- end
-
-
-!
-!
-! !functions used in couplings
-!
-! ! ===============================================================================================================
-! ! ============================================== Fermion charges ================================================
-! ! ===============================================================================================================
-!
-!
-! double precision function ef_(type)
-! integer type
-!
-! if (type.eq.1) then
-!
-! ef_ = 0d0
-!
-! elseif (type.eq.2) then
-!
-! ef_ = -1d0
-!
-! elseif (type.eq.3) then
-!
-! ef_ = 2/3d0
-!
-! elseif (type.eq.4) then
-!
-! ef_ = -1/3d0
-!
-! else
-! endif
-! end
-!
-!
-!
-!
-! ! ===============================================================================================================
-! ! ============================================== Fermion isospin ================================================
-! ! ===============================================================================================================
-!
-! double precision function I3f_(type)
-! integer type
-!
-! if ((type.eq.1).or.(type.eq.3)) then
-!
-! I3f_ = 1/2d0
-!
-! else
-!
-! I3f_ = -1/2d0
-!
-! endif
-! end
-!
-!
-!
-!
-! ! ===============================================================================================================
-! ! ============================================== Yukawa coupling ================================================
-! ! ===============================================================================================================
-!
-! function hf_(type, gen)
-! #include "stst2QQ_GenCouplings.h"
-! #include "stst2QQ_Model.h"
-! integer type, gen
-!
-! if ((type.eq.1).or.(type.eq.3)) then
-!
-! hf_ = EL*Mf(type,gen)/(sqrt2*MW*SW*SB)
-!
-! else
-!
-! hf_ = EL*Mf(type,gen)/(sqrt2*MW*SW*CB)
-!
-! endif
-! end
\ No newline at end of file
+ end
\ No newline at end of file
Index: branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_Main.F
===================================================================
--- branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_Main.F (revision 1091)
+++ branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_Main.F (revision 1092)
@@ -1,515 +1,487 @@
* Template for a new process @ NLO
* DM@NLO Project - B.Herrmann, K.Kovarik, M.Klasen, Q.Le Boulc`h, J.Harz, P.Steppeler, M. Meinecke, S. Schmiemann
#include "stst2QQ_Couplings.F"
#include "stst2QQ_Kinematics.F"
********************************
subroutine stst2qq_cs(renscale,scalevarpara,Pcm,iflag,myresult)
implicit none
********************************
#include "../util/DMNLO_LoopIntegrals.h"
#include "stst2QQ_Model.h"
#include "stst2QQ_Kinematics.h"
#include "stst2QQ_GenCouplings.h"
#include "stst2QQ_GenCounterterms.h"
- integer ndim, ncomp, iflag(15), inptype, inpgen, iflux, i, j, k, platzhalter
+ integer ndim, ncomp, iflag(15), inptype, inpgen, iflux, i, j, k, platzhalter, u
c parameter (ndim = 2)
parameter (ncomp = 1)
double precision Pcm,myresult(9), sinth, costh, avgfac, massth
double precision renscale, scalevarpara(5)
double precision result(9), error(9)
double precision result23(ncomp),error23(ncomp)
double precision cut
double precision theta,thetastep
double precision AS,zetaDC(2),ASDR
+
+ character buff(4)
+ double precision p(0:3,4)
* subroutines
external DMNLO_ModelPara, DMNLO_ModelIni, DMNLO_ModelDigest, Init_RenScheme
external DMNLO_CalcPifth, DMNLO_CalcDeltaMb, DMNLO_CalcPiftest
external ASzetaDecCoe
QScale = renscale
Qscalealphas = Qscale
ischeme = 1
** flags assignment
c number of Sfermion1
isf1 = iflag(1)
c type of initial state partical
itt1 = iflag(2)
! number of family of initial state partical
igen1 = iflag(3)
! number of Sfermion2
isf2 = iflag(4)
! type of initial state anti-partical
itt2 = iflag(5)
! number of family of initial state anti-partical
igen2 = iflag(6)
! type of final state partical
ftt1 = iflag(7)
! number of family of final state partical
fgen1 = iflag(8)
! type of final state anti-partical
ftt2 = iflag(9)
! number of family of final state anti-partical
fgen2 = iflag(10)
! imass = 0- take masses from MicrOmegas, 1-diagonalize here
imass = iflag(11)
! iflux = 0 - calculate v.sigma, 1- standard cross-section sigma
iflux = iflag(12)
! itree = 0- full one loop, 1- only tree level (quicker)
itree = iflag(13)
-! write(*,*)'MT', Mf(3,3), 'MB', Mf(4,3), 'MSt',MSf(1,3,3), 'MSb',MSf(1,4,3)
-! write(*,*)'MNEU1',MNeu(1)
+
if ((itt1.eq.itt2).AND.(ftt1.eq.ftt2))then
!Higgs
imin=1
imax=4
jmin=1
jmax=4
!Vektorboson
bmin=2
bmax=3
cmin=2
cmax=3
-! write(*,*)'This do work by now'
-!
+
if(itt1.eq.ftt1) then
! write(*,*)'S-channel and t-channel with neutral particles are activated'
!Gaugino
pmin=1
pmax=4
qmin=1
qmax=4
else
!write(*,*)'S-channel with neutral particles is activated and and t-channel with Chargino is activated'
!Gaugino
pmin=5
pmax=6
qmin=5
qmax=6
endif
elseif ((itt1.eq.ftt1).AND.(itt2.eq.ftt2)) then
if(itt1.eq.itt2-1) then
!write(*,*)'S-channel with charged particles is activated and and t-channel with neutral particles is activated'
elseif (itt2.eq.itt1-1) then
!write(*,*)'Changing up and downtype quark as particle and antiparticle in both state'
platzhalter = itt1
itt1 = itt2
itt2 = platzhalter
platzhalter = ftt1
ftt1 = ftt2
ftt2=platzhalter
endif
!setzten der Propagator-Indizes
!Higgs
imin=5
imax=6
jmin=5
jmax=6
!Vekotrboson
bmin=4
bmax=4
cmin=4
cmax=4
!Gaugino
pmin=1
pmax=4
qmin=1
qmax=4
elseif ((itt2.eq.ftt1).AND.(itt2.eq.ftt1)) then
if(itt2.eq.itt1-1) then
!write(*,*)'Changing up and downtype quark as particle and antiparticle in initial state'
platzhalter = itt1
itt1 = itt2
itt2 = platzhalter
elseif (ftt2.eq.ftt1-1) then
!write(*,*)'Changing up and downtype quark as particle and antiparticle in final state'
platzhalter = ftt1
ftt1 = ftt2
ftt2=platzhalter
endif
!setzten der Propagator-Indizes
!Higgs
imin=5
imax=6
jmin=5
jmax=6
!Vekotrboson
bmin=4
bmax=4
cmin=4
cmax=4
!Gaugino
pmin=1
pmax=4
qmin=1
qmax=4
else
write(*,*)'This do not work - please enter first the uptype-quark and then the downtype quark'
endif
-! write(*,*)'TEST'
-! write(*,*)'MTos',MTos
-!Check ob Prozess richtig angegeben wurde:
+
! ===== Set model parameters =====
! reading the MSSM parameters from micrOmegas and setting the masses & mixings & all shortforms e.g. TB,SA
call DMNLO_ModelPara()
call DMNLO_ModelIni()
! call DMNLO_ModelDigest
if (scalevarpara(1).eq.1) then
! print*,"scalevariation"
QScale = scalevarpara(2)
! switch for alpha_s scale
Qscalealphas = scalevarpara(2)
! Qscalealphas = renscale
Af(3,3) = scalevarpara(3)
Af(4,3) = scalevarpara(4)
endif
-! write(*,*)'TEST2'
! ===== set sfermion-widths to zero/10 ===========
do i = 1,2
do j = 3,4
do k = 1,3
if(WSf(i,j,k).lt.10.0) then
!see Warning below for purpose of this
!print*,'WARNING: MicrOMEGAs sQuark-width small (=',WSf(i,j,k),'GeV )=> set width WSf(',i,',',j,',',k,') to 10 GeV'
WSf(i,j,k)=10.0
endif
enddo
enddo
enddo
-! write(*,*)'TEST3'
-! write(*,*)itt3, igen3, itt4,igen4
-
-! Define the CMS \sqrt{s} and the scale Q2
-! if (icha.eq.0) then
- sqrtS = dsqrt(Pcm**2 + MSf2(isf1,itt1,igen1)) + dsqrt(Pcm**2 + MSf2(isf2,itt2,igen2))
-! else if (icha.eq.1) then
-! sqrtS = dsqrt(Pcm**2 + MCha2(indi)) + dsqrt(Pcm**2 + MNeu2(indj))
-! else if (icha.eq.2) then
-! sqrtS = dsqrt(Pcm**2 + MCha2(indi)) + dsqrt(Pcm**2 + MCha2(indj))
-! endif
-! write(*,*)'TEST4'
-! write(*,*)itt3, igen3, itt4,igen4
- write(*,*)'sqrtS',sqrtS,Pcm,MSf(isf1,itt1,igen1),MSf(isf2,itt2,igen2)
+
+ sqrtS = dsqrt(Pcm**2 + MSf2(isf1,itt1,igen1)) + dsqrt(Pcm**2 + MSf2(isf2,itt2,igen2))
! ===== Definition of Renormalization Scheme, Scale & Input =====
call Init_RenScheme(99)
-! write(*,*) 'QScale',QScale
-! write(*,*)'ALPHAMZ', AlfasMZ
-! write(*,*)'GS',GS
-! write(*,*)'AS', AS(sqrtS**2, MZ**2, AlfasMZ/(4d0*pi), 5, 3, 5)
+
! call ASzetaDecCoe(sqrtS**2,zetaDC)
! write(*,*) 'zeta',zetaDC
! write(*,*)'ASDR',ASDR(sqrtS**2)
+
! ===== Renormalization Scales for Loop integrals =====
! Divergent UV & IR poles for Loop integrals
UVdiv = 0d0
IRdiv = 0d0
! Renormalization scale \mu & switch between dim.reg and mass regularization for IR divergence - xiIR
! xiIR enables a check with mass regularization if set to 0d0 -> DEFAULT is xiIR = 1d0 !!!!
xiIR = 1d0
muSc = QScale
! EpsPoles = -2,-1,0,1,... steers LoopFunctions
! e.g. EpsPoles = -2/-1 lets all LoopFunctions return the IR diveregent coefficient corresponding to double/single pole
! checking UV divergence requires setting EpsPole=1 and UVdiv=1 (loop amplitude should be zero or extremely small)
EpsPole = 0d0
! flux factor (in mykinematics.h)
fluxflag = iflux
! write(*,*)'TEST5'
! write(*,*)itt3, igen3, itt4,igen4
! Generic couplings (here the type and generation is set !!! )
call stst2qq_SetCouplings
! write(*,*)'TEST6'
! write(*,*)itt3, igen3, itt4,igen4
! FA kinematics & reduced masses
call stst2qq_SetKinematics
! self-energies & counterterms
! call Template_SetCounterTerms()
*******************************************
! Mass threshold
massth = sqrtS*(muu3+muu4)
-! print*,'sqrtS=',sqrtS
-! print*,'massth=',massth
-
if (sqrtS.le.massth) then
do i=1,9
result(i) = 0d0
enddo
result23(1) = 0d0
else
! Compute integrated cross-section
call stst2qq_IntCS(result,error,itree)
-! call stst2qq_DiffCS(result, dcos(1.16025d0))
-! write(*,*)'TEST'
+
+
! Compute differential cross-section
cc call stst2qq_DiffCS(result, dcos(pi/4d0))
! open(5, FILE='Diff_CShphm.txt',status='REPLACE')
! ! write(5,*)'sqrtS', sqrtS
! ! write(5,*)'result', 'dcos(...)'
! do i=0,20
! call stst2qq_DiffCS(result, dcos(pi/10d0*REAL(i)))
! ! write(5,*)result, dcos(REAL(i)/10d0*pi)
! ! write(5,*)result,i,REAL(i),dcos((pi/10d0))
! ! write(5,*) i
! ! write(5,*) REAL(i)
! write(5,*) dcos(pi/10d0*REAL(i)),result(1)
! ! write(5,*) result(1)
! ! write(5,*)'--------------------------'
! end do
! close(5)
result23(1) = 0d0
! Copy to output vectors
if (itree.eq.0) then
do i=1,6
myresult(i) = result(i)
enddo
myresult(7) = result23(1)
myresult(8) = error23(1)
myresult(9) = result(6) + result23(1)
else
do i=1,9
myresult(i) = result(i)
enddo
endif
endif
end
**************************************************************
subroutine stst2qq_IntCS(result,error,mytree)
implicit none
#include "stst2QQ_Kinematics.h"
#include "stst2QQ_Model.h"
double precision result(9), error(9)
double precision costhmin, costhmax, relaccuracy, absaccuracy
double precision kappa
integer c, neval, fail, ncomp, mytree
external stst2qq_TreeDiffCS, stst2qq_DiffCS, Patterson
ncomp = 7
! Integration limits for cos(theta)
costhmin = -1d0
costhmax = 1d0
! Set requested accuracy
absaccuracy = 1d-12
relaccuracy = 1d-5
! Full calculation
if (mytree.eq.0) then
call Patterson(ncomp, costhmin, costhmax, stst2qq_DiffCS,
& relaccuracy, absaccuracy, neval, fail, result, error)
! Only tree-level calculation
else
call Patterson(ncomp, costhmin, costhmax, stst2qq_TreeDiffCS,
& relaccuracy, absaccuracy, neval, fail, result, error)
endif
c if( fail .ne. 0 )
c & print *, "Failed to reach the desired accuracy."
end
*************************************
subroutine stst2qq_DiffCS(result, costh)
implicit none
-#include "stst2QQ_Kinematics.h"f
+#include "stst2QQ_Kinematics.h"
#include "stst2QQ_Model.h"
#include "stst2QQ_GenCouplings.h"
double precision costh, intfac, intfacMO, avgfac, testmin, testmax
double precision result(9)
double precision kappa
double complex stst2qq_M2Tree
double precision test
-! write(*,*) 'cos(theta)', costh
+
! Mandelstam variables t=(p1-k1)^2 , u=(p1-k2)^2
tman = sqrtS**2*(muu1**2 + muu3**2 - 0.5d0*(1d0+muu1**2-muu2**2)*(1d0+muu3**2-muu4**2) +
& 0.5d0*kappa(1d0,muu1**2,muu2**2)*kappa(1d0,muu3**2,muu4**2)*costh)
tred = (muu1**2 + muu3**2 - 0.5d0*(1d0+muu1**2-muu2**2)*(1d0+muu3**2-muu4**2) +
& 0.5d0*kappa(1d0,muu1**2,muu2**2)*kappa(1d0,muu3**2,muu4**2)*costh)
uman = sqrtS**2*(muu1**2 + muu2**2 + muu3**2 + muu4**2 - 1d0 - tred)
ured = muu1**2 + muu2**2 + muu3**2 + muu4**2 - 1d0 - tred
! Integration factor for 2->2 phase-space integration in CMS
! (the azimuthal integration constant 2pi is inlcuded in intfac)
intfac = 2*pi*(kappa(1d0,muu3**2,muu4**2)/2d0)/(4*(2*pi)**2)
intfacMO = 2*pi*(kappa(sqrtS**2,Mfm(4)**2,Mfm(4)**2)/(2d0*sqrtS))/(4*(2*pi)**2*sqrtS)
! Flux factors for different units of cross-section
if (fluxflag.eq.1) then
! in units of pb
cc flux = hbar_c2/(4d0*(kappa(sqrtS**2,muu1**2*sqrtS**2,muu2**2*sqrtS**2)/(2d0*sqrtS))*sqrtS)
! in units of GeV^-2
flux = 1d0/(4d0*(kappa(1d0,muu1**2,muu2**2)/2d0)*sqrtS**2)
else
! in units of cm^3/sec
c flux = 2.9979d-26*hbar_c2/sqrtS**2
! in units of GeV^-2
flux = 1d0/sqrtS**2
endif
!!distinguish between identical and non identical particels
! if(indk.ge.indl) then
! if(indk.lt.5)then
!! Average factor+Symmetryfactor from equal final particels
! avgfac = 1d0/4d0*1d0/2d0
! end if
! else
! Average factor
avgfac = 1d0
! end if
! ***************************************
! contributions to the matrix element
! ***************************************
! (1-tree, 2-full one loop virtual, 3-vertex, 4-propagator, 5-box, 6-dipole)
result(1) = intfac*flux*avgfac*stst2qq_M2Tree()
end
*************************************
subroutine stst2qq_TreeDiffCS(result, costh)
implicit none
#include "stst2QQ_Kinematics.h"
#include "stst2QQ_Model.h"
#include "stst2QQ_GenCouplings.h"
double precision costh, intfac, intfacMO, avgfac
double precision result(9)
double precision kappa, stst2qq_M2TreeM0, stst2qq_M2TreeHO, stst2qq_M2TreePerc
double complex stst2qq_M2Tree
-! write(*,*)'TEST'
-
-
! Mandelstam variables t=(p2-k1)**2 , u=(p1-k1)**2
! Scattering angle theta between p1 and k1
tman = sqrtS**2*(muu2**2 + muu3**2 - 0.5d0*(1d0-muu1**2+muu2**2)*(1d0-muu4**2+muu3**2) -
& 0.5d0*kappa(1d0,muu1**2,muu2**2)*kappa(1d0,muu3**2,muu4**2)*costh)
tred = (muu2**2 + muu3**2 - 0.5d0*(1d0-muu1**2+muu2**2)*(1d0-muu4**2+muu3**2) -
& 0.5d0*kappa(1d0,muu1**2,muu2**2)*kappa(1d0,muu3**2,muu4**2)*costh)
uman = sqrtS**2*(muu1**2 + muu2**2 + muu3**2 + muu4**2 - 1d0 - tred)
ured = muu1**2 + muu2**2 + muu3**2 + muu4**2 - 1d0 - tred
! Integration factor for 2->2 phase-space integration in CMS
! the azimuthal integration constant 2pi is inlcuded in intfac
intfac = 2*pi*(kappa(1d0,muu3**2,muu4**2)/2d0)/(4*(2*pi)**2)
! Flux factor
if (fluxflag.eq.1) then
! in units of pb
cc flux = hbar_c2/(4d0*(kappa(sqrtS**2,muu1**2*sqrtS**2,muu2**2*sqrtS**2)/(2d0*sqrtS))*sqrtS)
! in units of GeV^-2
flux = 1d0/(4d0*(kappa(1d0,muu1**2,muu2**2)/2d0)*sqrtS**2)
else
! in units of cm^3/sec
c flux = 2.9979d-26*hbar_c2/sqrtS**2
! in units of GeV^-2
flux = 1d0/sqrtS**2
endif
avgfac = 1d0
! no flux and averaging factor!
result(1) = intfac*flux*avgfac*stst2qq_M2Tree()
-! write(*,*)'t=', tman
-! write(*,*)'ZL=', ZL(1,1,5)
-! write(*,*)'ZR=', ZR(1,1,5)
-! write(*,*)'MNE1=', MNeu(1)
-! write(*,*)'MC1=',MCha(1)
-! write(*,*)'MHc=',Mhk(5)
-! write(*,*)'s=', SqrtS**2
+
end
#include "stst2QQ_M2Tree.F"
Index: branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_Kinematics.F
===================================================================
--- branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_Kinematics.F (revision 1091)
+++ branches/DMNLO/run_dmnlo/stst2QQ/stst2QQ_Kinematics.F (revision 1092)
@@ -1,56 +1,28 @@
!**********************************
!**** subroutine to set the kinematics of the ingoing 4-momenta and magnitude of
!**** in/out-going momenta
!***********************************
subroutine stst2qq_SetKinematics
#include "stst2QQ_Kinematics.h"
#include "stst2QQ_Model.h"
-! write(*,*) ftt1, fgen1, ftt2,fgen2
-! write(*,*) iSf1,iSf2
muu1 = MSf(iSf1,itt1,igen1)/sqrtS
muu2 = MSf(iSf2,itt2,igen2)/sqrtS
muu3 = Mf(ftt1,fgen1)/sqrtS
muu4 = Mf(ftt2,fgen2)/sqrtS
-! muu4 = 2.750102/sqrtS
-! write(*,*) MT/sqrtS,muu3,muu4
!Reduced masses of propagator
do i=1,4
muuqt(i) = MNeu(i)/sqrtS
muupt(i) = MNeu(i)/sqrtS
enddo
muuqt(5) = MCha(1)/sqrtS
muupt(5) = MCha(1)/sqrtS
muuqt(6) = MCha(2)/sqrtS
muupt(6) = MCha(2)/sqrtS
muugl = MGl/sqrtS
-! !Propagator
-! DenSV(1) = sqrtS**2
-! DenSV(2) = sqrtS**2
-! DenSV(3) = (sqrtS**2-MZ**2)!+ cI*MZ*WZ
-! DenSV(4) = (sqrtS**2-MW**2)
-!
-! do i=1,6
-! DensH(i)=(sqrtS**2-MHk(i)**2)!+ cI*MHk(i)*wHk(i)
-! enddo
-! write(*,*)'kinematics', tman, MNeu(1),MNeu(2),MNeu(3),MNeu(4)
-! wNeu(1)=findValW("wNE1")
-! wNeu(2)=findValW("wNE2")
-! wNeu(3)=findValW("wNE3")
-! wNeu(4)=findValW("wNE4")
-
-! do i=1,4 !hier gluino und neutralino unterscheiden
-! DenGauT(i)=(tman - MNeu(i)**2) !+ cI*MNeu(i)*wNeu(i)
-! enddo
-! DenGauT(5)=(tman - MCha(1)**2) !+ cI*MCha(1)*wCha(1)
-! DenGauT(6)=(tman - MCha(2)**2) !+ cI*MCha(2)*wCha(2)
-! DenGlT = (tman-MGl**2)
-
-
-
end

File Metadata

Mime Type
text/x-diff
Expires
Tue, Sep 30, 5:42 AM (1 h, 9 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
6566272
Default Alt Text
(80 KB)

Event Timeline

HEPForge.org · .