00001 *
00002 *
00003 *
00004 *
00005 *
00006 *
00007 *
00008 *
00009 *
00010 DOUBLE PRECISION m_pi
00011 PARAMETER( m_pi=3.1415926535897932d0)
00012 INTEGER m_lenwt ! max. no. of aux weights
00013 PARAMETER (m_lenwt=200)
00014 INTEGER m_phmax ! maximum photon multiplicity ISR+FSR
00015 PARAMETER ( m_phmax = 100)
00016 *
00017 DOUBLE COMPLEX m_Pauli, m_Pauli4
00018 DOUBLE COMPLEX m_AmpBorn
00019 DOUBLE COMPLEX m_AmpBoxy
00020 DOUBLE COMPLEX m_AmpBorn1
00021 DOUBLE COMPLEX m_AmpBorn2
00022 DOUBLE COMPLEX m_AmpExpo0, m_AmpExpo1, m_AmpExpo2, m_AmpExpo2p
00023 DOUBLE COMPLEX m_Rmat
00024 DOUBLE COMPLEX m_IntReson, m_BoxGGtu, m_BoxGZtu, m_BoxGGut, m_BoxGZut !
00025 DOUBLE COMPLEX m_F1ini1, m_F1fin1, m_F1ini2, m_F1fin2, m_F1iniPair2, m_F1finPair2 !
00026 DOUBLE PRECISION m_WtSet, m_WtBest, m_Debg, m_Vcut
00027 DOUBLE PRECISION m_e_QED, m_Alfpi , m_Alfinv
00028 DOUBLE PRECISION m_Loren1, m_Lorin1
00029 DOUBLE PRECISION m_Loren2, m_Lorin2
00030 DOUBLE PRECISION m_xi, m_eta, m_b,m_b1, m_b2, m_b3
00031 DOUBLE PRECISION m_Sw2, m_MZ, m_GammZ, m_Gmu
00032 INTEGER m_KeyArb, m_isr, m_out
00033 INTEGER m_Phel
00034 INTEGER m_KeyISR, m_KeyFSR, m_KeyINT, m_KeyElw, m_KeyZet, m_KeyGPS
00035 INTEGER m_HasFSR
00036 DOUBLE PRECISION m_PolBeam1, m_PolBeam2
00037 DOUBLE PRECISION m_HvecFer1, m_HvecFer2
00038 DOUBLE PRECISION m_RhoCrud, m_ExpoNorm
00039 DOUBLE PRECISION m_RhoExp0, m_RhoExp1, m_RhoExp2, m_RhoExp2p !
00040 DOUBLE COMPLEX m_SDMat1, m_SDMat2, m_SDMat3, m_SDMat4
00041 *
00042 DOUBLE COMPLEX sAph1(2),sBph1(2),sAph2(2),sBph2(2),sAph1o(2),sBph1o(2),sAph2o(2),sBph2o(2)
00043 COMMON /c_GPS/
00044 $ m_Pauli( 0:3, 1:2, 1:2),! Pauli matrices
00045 $ m_Pauli4(1:4, 1:2, 1:2),! Pauli matrices, other vector index numbering
00046 $ m_IntReson, ! Virtual: resonance part of correction
00047 $ m_BoxGGtu, ! Virtual: part of box
00048 $ m_BoxGZtu, ! Virtual: part of box
00049 $ m_BoxGGut, ! Virtual: part of box
00050 $ m_BoxGZut, ! Virtual: part of box
00051 $ m_F1ini1, ! Virtual: vertex ISR, O(alf1) corr.
00052 $ m_F1fin1, ! Virtual: vertex ISR, O(alf2) corr.
00053 $ m_F1ini2, ! Virtual: vertex FSR, O(alf1) corr.
00054 $ m_F1fin2, ! Virtual: vertex FSR, O(alf2) corr.
00055 $ m_F1iniPair2, ! Virtual: IS Pairs in vertex, O(alf2) corr.
00056 $ m_F1finPair2, ! Virtual: FS Pairs in vertex, O(alf2) corr.
00057 $ m_Rmat( 0:3, 0:3), ! Spin correlation matrix for final particles (test)
00058 $ m_AmpBorn( 2,2,2,2), ! Born spin amplitudes, final mass kept exactly
00059 $ m_AmpBoxy( 2,2,2,2), ! Box spin amplitudes, gamma-gamma and gamma-Z
00060 $ m_AmpBorn1(2,2,2,2), ! Born spin amplitudes, massles
00061 $ m_AmpBorn2(2,2,2,2), ! Born spin amplitudes, massles
00062 $ m_AmpExpo0(2,2,2,2), ! CEEX spin amplitudes O(alf0)ceex
00063 $ m_AmpExpo1(2,2,2,2), ! CEEX spin amplitudes O(alf1)ceex
00064 $ m_AmpExpo2(2,2,2,2), ! CEEX spin amplitudes O(alf2)ceex
00065 $ m_AmpExpo2p(2,2,2,2), ! CEEX spin amplitudes O(alf2)ceex with PAIRS
00066 $ m_PolBeam1(4), ! POLARIZATION vector 1-st beam PRIVATE no getter!!!
00067 $ m_PolBeam2(4), ! POLARIZATION vector 2-nd beam PRIVATE no getter!!!
00068 $ m_HvecFer1(4), ! POLARIMETER vector 1-st fin.ferm. PRIVATE no getter!!!
00069 $ m_HvecFer2(4), ! POLARIMETER vector 2-nd fin.ferm. PRIVATE no getter!!!
00070 $ m_SDMat1(2,2), ! Spin Density Matrix 1-st beam
00071 $ m_SDMat2(2,2), ! Spin Density Matrix 2-nd beam
00072 $ m_SDMat3(2,2), ! Polarimeter Density Matrix 1-st fin.ferm
00073 $ m_SDMat4(2,2), ! Polarimeter Density Matrix 2-nd fin.ferm
00074 $ m_RhoCrud, ! Diff.distr. Crude MC, normalized to LIPS
00075 $ m_RhoExp0, ! Diff.distr. O(alf0)ceex, normalized to LIPS
00076 $ m_RhoExp1, ! Diff.distr. O(alf1)ceex, normalized to LIPS
00077 $ m_RhoExp2, ! Diff.distr. O(alf2)ceex, normalized to LIPS
00078 $ m_RhoExp2p, ! Diff.distr. O(alf2)ceex+PAIRS, normalized to LIPS
00079 $ m_ExpoNorm, ! Normalized factor to LIPS for CEEX amplitudes
00080 $ m_Debg(0:200), ! Debug facility
00081 $ m_WtSet(m_lenwt), ! List of weights
00082 $ m_WtBest, ! The best of the weights
00083 $ m_Vcut(3), ! technical cut on E/Ebeam for non-IR real contributions
00084 $ m_xi(4), ! GPS axial vector, lightlike xi**2=0
00085 $ m_eta(4), ! GPS auxiliar vector, spacelike eta**2=-1
00086 $ m_b(4), ! Beta vector for photon polarization (Axial gauge)
00087 $ m_b1(4), ! Beta vector as above for tests
00088 $ m_b2(4), ! Beta vector as above for tests
00089 $ m_b3(4), ! Beta vector as above for tests
00090 $ m_Loren1(4,4), ! Lorenz tranfrom ONE for TRALOR (down to CMS)
00091 $ m_Lorin1(4,4), ! and its inverse (for tests)
00092 $ m_Loren2(4,4), ! Lorenz tranfrom TWO for TRALOR (down to CMS)
00093 $ m_Lorin2(4,4), ! and its inverse (for tests)
00094 $ m_Sw2, ! Electroweak: mixing angle sin(theta_W)^2
00095 $ m_MZ, ! Electroweak: Z mass
00096 $ m_GammZ, ! Electroweak: Z width
00097 $ m_Gmu, ! Electroweak: Fermi constant
00098 $ m_e_QED, ! Electroweak: Electric coupling constant
00099 $ m_Alfpi, ! Electroweak: alfQED/pi
00100 $ m_Alfinv, ! Electroweak: 1/alfQED
00101 $ m_KeyArb, ! Key for switching on/off the use of m_b
00102 $ m_isr( m_phmax), ! partition =1,0 for isr/fsr
00103 $ m_Phel(m_phmax), ! photon helicity =1,0 for +,-
00104 $ m_KeyElw, ! Electroweak switch
00105 $ m_KeyZet, ! Z switch
00106 $ m_KeyISR, ! ISR static switch
00107 $ m_KeyFSR, ! FSR static switch
00108 $ m_KeyINT, ! INT static interference switch
00109 $ m_KeyGPS, ! CEEX perturb. level switch
00110 $ m_HasFSR, ! FSR dynamic switch
00111 $ m_out ! output unit number
00112
00113 SAVE /c_GPS/
00114 *
00115 *
00116 *
00117 *
00118 *
00119
00120
00121
