00001 *
00002 *
00003 *
00004 *
00005 *
00006 *
00007 *
00008 *
00009 *
00010 *
00011 *
00012 * Class members:
00013 *
00014 *
00015 DOUBLE PRECISION m_pi
00016 PARAMETER (m_pi =3.1415926535897932d0)
00017 DOUBLE PRECISION m_fleps
00018 **** PARAMETER (m_fleps = 1d-35) ! original
00019 **** PARAMETER (m_fleps = 1d-45) ! enough???
00020 PARAMETER (m_fleps = 1d-100) ! enough!!!
00021 **** PARAMETER (m_fleps = 1d-200) ! enough!!!
00022 *
00023 *
00024 *
00025 * 340-point grid, only 80pt for NLC to be improved/tested in future
00026 INTEGER m_poin1, m_poin2, m_poin3, m_poin4, m_poinG , m_poinQ !
00027 INTEGER m_poTh1, m_poTh2, m_poTh3, m_poTh4
00028 PARAMETER( m_poinG = 7 ) ! No of EW formfactors
00029 PARAMETER( m_poinQ = 4 ) ! No of QCD corrections
00030 *----------- Low energies and LEP1
00031 PARAMETER( m_poin1 = 100 ) ! LEP1 LOG(SQRT(s)) spacing
00032 PARAMETER( m_poTh1 = 0 ) ! Cost(heta) spacing
00033 DOUBLE PRECISION m_WminLEP1, m_WmaxLEP1 ! LEP1 basic range (m_WminLEP1,m_WmaxLEP1)
00034 PARAMETER( m_WminLEP1=0.010d0, m_WmaxLEP1= 95.000d0 ) ! LEP1 basic range (m_WminLEP1,m_WmaxLEP1)
00035 *** PARAMETER( m_WminLEP1=0.010d0, m_WmaxLEP1=100.001d0 ) ! LEP1 basic range (m_WminLEP1,m_WmaxLEP1)
00036 *----------- Z resonance
00037 PARAMETER( m_poin2 = 20 ) ! Z range sqrt(s) spacing
00038 *** PARAMETER( m_poTh2 = 0 ) ! Cost(heta) spacing
00039 PARAMETER( m_poTh2 = 14 ) ! =14 is overkill?
00040 DOUBLE PRECISION m_WminZ, m_WmaxZ, m_WdelZ ! Z range (amz + m_WdelZ)
00041 PARAMETER( m_WdelZ = 5.000d0) ! Z range (amz + m_WdelZ)
00042 *** PARAMETER( m_WdelZ = 2d0 *2.495342704946d0) ! Old range 2*Gamma
00043 *----------- LEP2
00044 PARAMETER( m_poTh3 = 30 ) ! Overkill, bit lets kkep it
00045 *** PARAMETER( m_poTh3 = 14 ) ! Cost(heta) spacing
00046 PARAMETER( m_poin3 = 145 ) ! LEP2 interval sqrt(s) spacing
00047 *** PARAMETER( m_poin3 = 140 ) ! as in KORALZ, also enough
00048 DOUBLE PRECISION m_WmaxLEP2 ! LEP2 interval (m_WmaxLEP1,m_WmaxLEP2)
00049 PARAMETER( m_WmaxLEP2 =240.001d0 ) ! LEP2 interval (m_WmaxLEP1,m_WmaxLEP2)
00050 *----------- Linear Colliders
00051 PARAMETER( m_poin4 = 80 ) ! NLC range sqrt(s) spacing
00052 PARAMETER( m_poTh4 = 14 ) ! Cost(heta) spacing
00053 DOUBLE PRECISION m_WmaxNLC ! NLC range (m_WmaxLEP2,m_WmaxNLC)
00054 PARAMETER( m_WmaxNLC =1040.001d0 ) ! NLC range (m_WmaxLEP2,m_WmaxNLC)
00055 *
00056 * EW formfactors, all flavours!!
00057 DOUBLE COMPLEX m_cyy, m_czz, m_ctt, m_clc ! Electroweak FFactors
00058 DOUBLE PRECISION m_syy, m_szz, m_stt, m_slc ! QCD corr.
00059 DOUBLE COMPLEX m_GSW
00060 DOUBLE PRECISION m_QCDcor ! obsolete
00061 DOUBLE PRECISION m_QCDcorR
00062 *
00063 * EW parameters
00064 DOUBLE PRECISION m_Gmu
00065 DOUBLE PRECISION m_MZ, m_amh, m_amtop
00066 DOUBLE PRECISION m_swsq, m_gammz, m_MW, m_GammW
00067 *
00068 DOUBLE PRECISION m_CMSene, m_XXXene, m_HadMin, m_vvmin, m_vvmax !
00069 DOUBLE PRECISION m_AvMult, m_YFSkon, m_YFS_IR, m_alfinv, m_alfpi, m_Xenph !
00070 DOUBLE PRECISION m_vv, m_x1, m_x2
00071 DOUBLE PRECISION m_Qf, m_T3f, m_helic, m_amferm, m_auxpar !
00072 DOUBLE PRECISION m_gnanob
00073
00074 INTEGER m_IsGenerated, m_KFferm, m_NCf
00075 INTEGER m_KFini, m_KeyINT, m_KeyQCD, m_KeyRes
00076 INTEGER m_KeyElw, m_KeyZet, m_KeyWtm
00077 INTEGER m_out
00078 CC+++ hekl
00079 INTEGER m_KeyCNT, m_KeyPsi, m_KeyRho
00080 CC--- hekl
00081 COMMON /c_BornV/
00082 * Tables of EW formfactors
00083 $ m_cyy(m_poin1+1, m_poinG,16), ! formfactor, table
00084 $ m_czz(m_poin2+1,m_poTh2+1,m_poinG,16), ! formfactor, table
00085 $ m_ctt(m_poin3+1,m_poTh3+1,m_poinG,16), ! formfactor, table, costheta dependent
00086 $ m_clc(m_poin4+1,m_poTh4+1,m_poinG,16), ! formfactor, table, costheta dependent
00087 $ m_syy(m_poin1+1, m_poinQ,16), ! QCD correction,
00088 $ m_szz(m_poin2+1, m_poinQ,16), ! QCD correction,
00089 $ m_stt(m_poin3+1, m_poinQ,16), ! QCD correction,
00090 $ m_slc(m_poin3+1, m_poinQ,16), ! QCD correction,
00091 $ m_GSW( m_poinG), ! form-factors, at the actual energy/angle
00092 $ m_QCDcorR(m_poinQ), ! QCD correction, at the actual energy/angle
00093 $ m_QCDcor, ! obsolete!!!!
00094 *
00095 $ m_CMSene, ! Initial value of CMS energy
00096 $ m_XXXene, ! CMS energy after beamsstrahlung or beam spread
00097 * -------------------- EVENT --------------------------
00098 $ m_x1, ! 1-z1 = x1 for first beam(strahlung)
00099 $ m_x2, ! 1-z2 = x2 for second beam(strahlung)
00100 $ m_vv, ! v = 1-sprim/s
00101 $ m_AvMult, ! Average photon multiplicity CRude at given v
00102 $ m_YFSkon, ! YFS formfactor finite part
00103 $ m_YFS_IR, ! YFS formfactor IR part
00104 * -----------------------------------------------------
00105 $ m_vvmin, ! minimum v, infrared cut
00106 $ m_vvmax, ! maximum v
00107 $ m_HadMin, ! minimum hadronization mass [GeV]
00108 * Basic QED and QCD
00109 $ m_alfinv, ! 1/alphaQED, Thomson limit (Q^2=0)
00110 $ m_alfpi, ! alphaQED/pi
00111 $ m_Xenph, ! Enhancement factor for Crude photon multiplicity
00112 * EW parameters
00113 $ m_MZ, ! Z mass
00114 $ m_amh, ! Higgs mass
00115 $ m_amtop, ! Top mass
00116 $ m_swsq, ! sin(thetaW)**2
00117 $ m_gammz, ! Z width
00118 $ m_MW, ! W mass
00119 $ m_GammW, ! W width
00120 $ m_Gmu, ! Fermi constant (from muon decay)
00121
00122 * Table of fermion paramerets, quarks (1->6) and leptons (11->16)
00123 $ m_KFferm(20), ! fermion KFcode (1->6) and (11->16)
00124 $ m_NCf(20), ! number of colours
00125 $ m_Qf(20), ! electric charge
00126 $ m_T3f(20), ! isospin, L-hand component
00127 $ m_helic(20), ! helicity or polarization
00128 $ m_amferm(20), ! fermion mass
00129 $ m_auxpar(20), ! auxiliary parameter
00130 $ m_IsGenerated(20), ! Generation flag, only for SAN !!!
00131 * Normalization
00132 $ m_gnanob, ! GeV^(-2) to nanobarns
00133 * Initial/final fermion types
00134 $ m_KFini, ! KF code of beam
00135 * Test switches
00136 $ m_KeyQCD, ! QCD FSR corr. switch
00137 $ m_KeyINT, ! ISR/FSR INTereference switch
00138 $ m_KeyElw, ! Type of Electrowak Library
00139 $ m_KeyZet, ! Z-boson on/off
00140 $ m_KeyWtm, ! Photon emission without mass terms
00141 $ m_KeyRes, ! experim. R for gamma* decays switch
00142 $ m_out, ! output unit for printouts
00143 $ m_KeyCNT, ! continuum generation
00144 $ m_KeyPsi, ! psi, psi', psi" on / off
00145 $ m_KeyRho ! rho, omega, phi on/off
00146 SAVE /c_BornV/
00147 *
00148 * Formats for writing EW tables onto disk file.
00149 CHARACTER*80 m_fmt0, m_fmt1, m_fmt2
00150 PARAMETER (
00151 $ m_fmt0 ='(4g20.13)', ! Mz,Mt,Mh etc.
00152 $ m_fmt1 ='( a, i4, f10.5, i4, f10.5 )', ! header
00153 $ m_fmt2 ='(6g13.7)' ) ! complex formfactors
00154 *
00155 * Class procedures:
00156 *
00157 *
00158 *
00159 *
00160 *