18****
PARAMETER (m_fleps = 1d-35) ! original
19****
PARAMETER (m_fleps = 1d-45) ! enough???
21****
PARAMETER (m_fleps = 1d-200) ! enough!!!
25* 340-point grid, only 80pt
for NLC to be improved/tested in future
29 PARAMETER( m_poinQ = 4 ) ! No of QCD corrections
30*----------- Low energies and LEP1
33 DOUBLE
PRECISION m_WminLEP1, m_WmaxLEP1 ! LEP1 basic range (m_WminLEP1,m_WmaxLEP1)
34 PARAMETER( m_WminLEP1=0.010d0, m_WmaxLEP1= 95.000d0 ) ! LEP1 basic range (m_WminLEP1,m_WmaxLEP1)
35***
PARAMETER( m_WminLEP1=0.010d0, m_WmaxLEP1=100.001d0 ) ! LEP1 basic range (m_WminLEP1,m_WmaxLEP1)
36*----------- Z resonance
40 DOUBLE
PRECISION m_WminZ, m_WmaxZ, m_WdelZ ! Z range (amz + m_WdelZ)
41 PARAMETER( m_WdelZ = 5.000d0) ! Z range (amz + m_WdelZ)
42***
PARAMETER( m_WdelZ = 2d0 *2.495342704946d0) ! Old range 2*Gamma
48 DOUBLE
PRECISION m_WmaxLEP2 ! LEP2 interval (m_WmaxLEP1,m_WmaxLEP2)
49 PARAMETER( m_WmaxLEP2 =240.001d0 ) ! LEP2 interval (m_WmaxLEP1,m_WmaxLEP2)
50*----------- Linear Colliders
52 PARAMETER( m_poTh4 = 14 ) ! Cost(heta) spacing
53 DOUBLE
PRECISION m_WmaxNLC ! NLC range (m_WmaxLEP2,m_WmaxNLC)
54 PARAMETER( m_WmaxNLC =1040.001d0 ) ! NLC range (m_WmaxLEP2,m_WmaxNLC)
56* EW formfactors, all flavours!!
57 DOUBLE COMPLEX m_cyy, m_czz, m_ctt, m_clc ! Electroweak FFactors
58 DOUBLE
PRECISION m_syy, m_szz, m_stt, m_slc ! QCD corr.
66 DOUBLE
PRECISION m_swsq, m_gammz, m_MW, m_GammW
71 DOUBLE
PRECISION m_Qf, m_T3f, m_helic, m_amferm, m_auxpar !
74 INTEGER m_IsGenerated, m_KFferm, m_NCf
79 INTEGER m_KeyCNT, m_KeyPsi, m_KeyRho
82* Tables of EW formfactors
86 $ m_clc(
m_poin4+1,m_poTh4+1,
m_poinG,16), ! formfactor, table, costheta dependent
87 $ m_syy(
m_poin1+1, m_poinQ,16), ! QCD correction,
88 $ m_szz(
m_poin2+1, m_poinQ,16), ! QCD correction,
89 $ m_stt(
m_poin3+1, m_poinQ,16), ! QCD correction,
90 $ m_slc(
m_poin3+1, m_poinQ,16), ! QCD correction,
92 $ m_QCDcorR(m_poinQ), ! QCD correction,
at the actual
energy/angle
93 $ m_QCDcor, ! obsolete!!!!
97* --------------------
EVENT --------------------------
98 $
m_x1, ! 1-z1 = x1
for first beam(strahlung)
99 $ m_x2, ! 1-z2 = x2
for second beam(strahlung)
101 $ m_AvMult, ! Average photon multiplicity CRude
at given
v
102 $
m_YFSkon, ! YFS formfactor finite part
103 $
m_YFS_IR, ! YFS formfactor IR part
104* -----------------------------------------------------
109 $
m_alfinv, ! 1/alphaQED, Thomson limit (Q^2=0)
110 $ m_alfpi, ! alphaQED/
pi
111 $
m_Xenph, ! Enhancement factor
for Crude photon multiplicity
114 $ m_amh, ! Higgs
mass
115 $ m_amtop, ! Top
mass
116 $ m_swsq, !
sin(thetaW)**2
120 $ m_Gmu, ! Fermi constant (from
muon decay)
122* Table of fermion paramerets, quarks (1->6) and leptons (11->16)
123 $ m_KFferm(20), ! fermion KFcode (1->6) and (11->16)
124 $ m_NCf(20), ! number of colours
125 $ m_Qf(20), ! electric charge
126 $ m_T3f(20), ! isospin, L-hand component
127 $ m_helic(20), !
helicity or polarization
128 $ m_amferm(20), ! fermion
mass
130 $ m_IsGenerated(20), ! Generation
flag, only
for SAN !!!
133* Initial/final fermion types
137 $
m_KeyINT, ! ISR/FSR INTereference switch
138 $
m_KeyElw, ! Type of Electrowak Library
141 $ m_KeyRes, ! experim. R
for gamma* decays switch
143 $ m_KeyCNT, ! continuum generation
144 $ m_KeyPsi, ! psi, psi', psi" on / off
145 $ m_KeyRho, ! rho, omega, phi on/off
149* Formats
for writing EW tables onto disk
file.
150 CHARACTER*80 m_fmt0, m_fmt1, m_fmt2
152 $ m_fmt0 ='(4g20.13)', ! Mz,Mt,Mh etc.
153 $ m_fmt1 ='( a, i4, f10.5, i4, f10.5 )', ! header
154 $ m_fmt2 ='(6g13.7)' ) !
complex formfactors
*******DOUBLE PRECISION m_EGridB INTEGER m_out
*********DOUBLE PRECISION m_pi
************Class only for NLC to be improved tested in future INTEGER m_poinQ !INTEGER m_poTh2
************Class members
************Class only for NLC to be improved tested in future INTEGER m_poinQ !INTEGER m_poTh3
************Class only for NLC to be improved tested in future INTEGER m_poin1
************Class only for NLC to be improved tested in future INTEGER m_poinQ !INTEGER m_poTh1
************Class only for NLC to be improved tested in future INTEGER m_poin4
************Class only for NLC to be improved tested in future INTEGER m_poinG
************Class only for NLC to be improved tested in future INTEGER m_poin2
************Class only for NLC to be improved tested in future INTEGER m_poin3
************Class m_alfQCDMZ INTEGER m_KFfin INTEGER m_IVfin INTEGER m_ibox *COMMON c_DZface $ alphaQED at(Q^2=MZ^2) DIZET $ m_alfQCDMZ
*******INTEGER m_nBinMax INTEGER m_NdiMax !No of bins in histogram for cell exploration division $ !Last vertex $ !Last active cell $ !Last cell in buffer $ !No of sampling when dividing cell $ !No of function total $ !Flag for random ceel for $ !Flag for type of for WtMax $ !Flag which decides whether vertices are included in the sampling $ entire domain is hyp !Maximum effective eevents per saves r n generator level $ !Flag for chat level in output
*******INTEGER m_nBinMax INTEGER m_NdiMax !No of bins in histogram for cell exploration division $ !Last vertex $ !Last active cell $ !Last cell in buffer $ !No of sampling when dividing cell $ !No of function total $ !Flag for random ceel for $ !Flag for type of for WtMax $ !Flag which decides whether vertices are included in the sampling $ entire domain is hyp !Maximum effective eevents per saves r n generator level $ !Flag for chat level in !Latex Output unit
*******INTEGER m_nBinMax INTEGER m_NdiMax PARAMETER(m_NdiMax=5) ! maximum dimension for siplices INTEGER m_KdiMax PARAMETER(m_KdiMax
*********DOUBLE PRECISION m_pi INTEGER m_lenwt !max no of aux weights INTEGER m_phmax !maximum photon multiplicity ISR FSR *DOUBLE COMPLEX m_Pauli4 DOUBLE COMPLEX m_AmpBorn DOUBLE COMPLEX m_AmpBoxy DOUBLE COMPLEX m_AmpBorn1 DOUBLE COMPLEX m_AmpBorn2 DOUBLE COMPLEX m_AmpExpo2p DOUBLE COMPLEX m_Rmat DOUBLE COMPLEX m_BoxGZut !DOUBLE COMPLEX m_F1finPair2 !DOUBLE PRECISION m_Vcut DOUBLE PRECISION m_Alfinv DOUBLE PRECISION m_Lorin1 DOUBLE PRECISION m_Lorin2 DOUBLE PRECISION m_b3 DOUBLE PRECISION m_Gmu INTEGER m_out INTEGER m_Phel INTEGER m_KeyZet
*********DOUBLE PRECISION m_pi INTEGER m_lenwt !max no of aux weights INTEGER m_phmax !maximum photon multiplicity ISR FSR *DOUBLE COMPLEX m_Pauli4 DOUBLE COMPLEX m_AmpBorn DOUBLE COMPLEX m_AmpBoxy DOUBLE COMPLEX m_AmpBorn1 DOUBLE COMPLEX m_AmpBorn2 DOUBLE COMPLEX m_AmpExpo2p DOUBLE COMPLEX m_Rmat DOUBLE COMPLEX m_BoxGZut !DOUBLE COMPLEX m_F1finPair2 !DOUBLE PRECISION m_Vcut DOUBLE PRECISION m_Alfinv DOUBLE PRECISION m_Lorin1 DOUBLE PRECISION m_Lorin2 DOUBLE PRECISION m_b3 DOUBLE PRECISION m_MZ
*********DOUBLE PRECISION m_pi INTEGER m_lenwt !max no of aux weights INTEGER m_phmax !maximum photon multiplicity ISR FSR *DOUBLE COMPLEX m_Pauli4 DOUBLE COMPLEX m_AmpBorn DOUBLE COMPLEX m_AmpBoxy DOUBLE COMPLEX m_AmpBorn1 DOUBLE COMPLEX m_AmpBorn2 DOUBLE COMPLEX m_AmpExpo2p DOUBLE COMPLEX m_Rmat DOUBLE COMPLEX m_BoxGZut !DOUBLE COMPLEX m_F1finPair2 !DOUBLE PRECISION m_Vcut DOUBLE PRECISION m_Alfinv DOUBLE PRECISION m_Lorin1 DOUBLE PRECISION m_Lorin2 DOUBLE PRECISION m_b3 DOUBLE PRECISION m_Gmu INTEGER m_out INTEGER m_Phel INTEGER m_KeyElw
*********DOUBLE PRECISION m_pi INTEGER m_lenwt !max no of aux weights INTEGER m_phmax !maximum photon multiplicity ISR FSR *DOUBLE COMPLEX m_Pauli4 DOUBLE COMPLEX m_AmpBorn DOUBLE COMPLEX m_AmpBoxy DOUBLE COMPLEX m_AmpBorn1 DOUBLE COMPLEX m_AmpBorn2 DOUBLE COMPLEX m_AmpExpo2p DOUBLE COMPLEX m_Rmat DOUBLE COMPLEX m_BoxGZut !DOUBLE COMPLEX m_F1finPair2 !DOUBLE PRECISION m_Vcut DOUBLE PRECISION m_Alfinv DOUBLE PRECISION m_Lorin1 DOUBLE PRECISION m_Lorin2 DOUBLE PRECISION m_b3 DOUBLE PRECISION m_Gmu INTEGER m_out INTEGER m_Phel INTEGER m_KeyINT
double sin(const BesAngle a)
************Class m_ypar INTEGER m_KeyWgt INTEGER m_KeyIHVP INTEGER m_KeyGPS INTEGER m_IsBeamPolarized INTEGER m_EvtGenInterface DOUBLE PRECISION m_Emin DOUBLE PRECISION m_sphot DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_q2 DOUBLE PRECISION m_PolBeam2 DOUBLE PRECISION m_xErrPb *COMMON c_KK2f $ !CMS energy average $ !Spin Polarization vector first beam $ !Spin Polarization vector second beam $ !Beam energy spread[GeV] $ !minimum hadronization energy[GeV] $ !input parameters
************Class m_ypar INTEGER m_KeyWgt INTEGER m_KeyIHVP INTEGER m_KeyGPS INTEGER m_IsBeamPolarized INTEGER m_EvtGenInterface DOUBLE PRECISION m_Emin DOUBLE PRECISION m_sphot DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_q2 DOUBLE PRECISION m_PolBeam2 DOUBLE PRECISION m_xErrPb *COMMON c_KK2f $ !CMS energy average $ !Spin Polarization vector first beam $ !Spin Polarization vector second beam $ !Beam energy spread[GeV] $ !minimum hadronization energy[GeV] $ !input READ never touch them !$ !debug facility $ !maximum weight $ !inverse alfaQED $ !minimum real photon energy
************Class m_ypar INTEGER m_KeyWgt INTEGER m_KeyIHVP INTEGER m_KeyGPS INTEGER m_IsBeamPolarized INTEGER m_EvtGenInterface DOUBLE PRECISION m_Emin DOUBLE PRECISION m_sphot DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_q2 DOUBLE PRECISION m_PolBeam2 DOUBLE PRECISION m_xErrPb *COMMON c_KK2f $ !CMS energy average $ !Spin Polarization vector first beam $ !Spin Polarization vector second beam $ !Beam energy spread[GeV] $ !minimum hadronization energy[GeV] $ !input READ never touch them !$ !debug facility $ !maximum weight $ !inverse alfaQED $ !minimum real photon IR regulator $ !ficticious photon IR regulator $ !Enhancement factor for Crude photon multiplicity $ !technical cut on E Ebeam for non IR real contributions $ !output cross section available through getter $ !output crossxsection available through getter *EVENT $ !e beam $ !e beam $ !final fermion $ !final anti fermion $ !photon momenta $ !MAIN weight of KK2f $ !crude distr from ISR and FSR $ !complete list of weights $ !complete list of weights $ !crude in nanobarns $ !Crude Born $ for fsr $ !photon helicity
************Class m_ypar INTEGER m_KeyWgt INTEGER m_KeyIHVP INTEGER m_KeyGPS INTEGER m_KFini
************Class m_ypar INTEGER m_KeyWgt INTEGER m_KeyIHVP INTEGER m_KeyGPS INTEGER m_IsBeamPolarized INTEGER m_EvtGenInterface DOUBLE PRECISION m_Emin DOUBLE PRECISION m_sphot DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_q2 DOUBLE PRECISION m_PolBeam2 DOUBLE PRECISION m_xErrPb *COMMON c_KK2f $ !CMS energy average $ !Spin Polarization vector first beam $ !Spin Polarization vector second beam $ !Beam energy spread[GeV] $ !minimum hadronization energy[GeV] $ !input READ never touch them !$ !debug facility $ !maximum weight $ !inverse alfaQED $ !minimum real photon IR regulator $ !ficticious photon IR regulator $ !Enhancement factor for Crude photon multiplicity $ !technical cut on E Ebeam for non IR real contributions $ !output cross section available through getter $ !output crossxsection available through getter *EVENT $ !e beam $ !e beam $ !final fermion $ !final anti fermion $ !photon momenta $ !MAIN weight of KK2f $ !crude distr from ISR and FSR $ !complete list of weights $ !complete list of weights $ !crude in nanobarns $ !Crude Born $ for fsr $ !photon for
************Class m_ypar INTEGER m_KeyWgt INTEGER m_KeyIHVP INTEGER m_KeyGPS INTEGER m_IsBeamPolarized INTEGER m_EvtGenInterface DOUBLE PRECISION m_Emin DOUBLE PRECISION m_sphot DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_q2 DOUBLE PRECISION m_PolBeam2 DOUBLE PRECISION m_xErrPb *COMMON c_KK2f $ !CMS energy average $ !Spin Polarization vector first beam $ !Spin Polarization vector second beam $ !Beam energy spread[GeV] $ !minimum hadronization energy[GeV] $ !input READ never touch them !$ !debug facility $ !maximum weight $ !inverse alfaQED $ !minimum real photon IR regulator $ !ficticious photon IR regulator $ !Enhancement factor for Crude photon multiplicity $ !technical cut on E Ebeam for non IR real contributions $ !output cross section available through getter $ !output crossxsection available through getter *EVENT $ !e beam $ !e beam $ !final fermion $ !final anti fermion $ !photon momenta $ !MAIN weight of KK2f $ !crude distr from ISR and FSR $ !complete list of weights $ !complete list of weights $ !crude in nanobarns $ !Crude Born $ for fsr $ !photon $ !Total Photon multiplicity *$ !status of beam polarization $ !KF of beam fermion $ !hadronization switch $ !serial number of the event $ !maximum photon multiplicity $ !ISR switch $ !FSR switch $ !ISR FSR INTereference switch $ m_KeyIHVP
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at alphaQED at MZ $ !GeV nanobarns $ !minimum v $ !maximum v $ !propagator sprim $ !longit polarization st beam $ !longit polarization nd beam $ !longit polarization st final ferm $ !longit polarization nd final ferm $ !Enhancement factor for Crude photon multiplicity $ !type of Born $ !type of bremsstrahlung $ !KF flavour code st beam $ !KF flavour code st final fermion $ !togle swith for Forward or Backward $ !FSR on off $ m_KeyQCD
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_gnanob
*********Class see also m_nmax DOUBLE PRECISION m_MasPhot DOUBLE PRECISION m_phsu DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_r2 DOUBLE PRECISION m_YFSkon
*********Class see also m_nmax DOUBLE PRECISION m_alfinv
*********Class see also m_nmax DOUBLE PRECISION m_MasPhot DOUBLE PRECISION m_phsu DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_r2 DOUBLE PRECISION m_WtMass INTEGER m_nmax INTEGER m_Nevgen INTEGER m_IsFSR INTEGER m_MarTot *COMMON c_KarFin $ !Output file $ !Event serial number $ !alpha QED at Thomson limit $ !minimum energy at CMS for remooval $ !infrared dimensionless $ !dummy photon IR regulator $ m_Xenph
*********Class see also m_nmax DOUBLE PRECISION m_MasPhot DOUBLE PRECISION m_phsu DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_r2 DOUBLE PRECISION m_YFS_IR
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_vvmin
**********Class see also m_nmax DOUBLE PRECISION m_CMSene
**********Class see also m_nmax DOUBLE PRECISION m_XXXene
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_vv
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_HadMin
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x1
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum infrared cut $ !maximum v $ !minimum hadronization mass *[GeV] Basic QED $ !Alpha_QED at Thomson limit * Normalization
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_vvmax
****************DOUBLE PRECISION m_vvmin INTEGER m_KeyWtm
IMPLICIT REAL *A O Z IMPLICIT INTEGER(I-N) C.... COMMONBALOCKS FOR MODE COMMON/MODEXS/ICH common/bornXS/xborn REAL *8 USER_XX(1000)
double precision pisqo6 parameter(pi=3.14159265358979d0, pisq=pi *pi, pisqo6=pisq/6d0) double precision zip
1.9.6