8* Version 2-nd
August 1997, by S. Jadach.
10* ISR O(alf3) + FSR O(alf2).
12* 1. Kinematical reduction
procedure kicked out.
13* 2. Numerical stability
for collinear photons improved up to
BHLUMI level.
17* All that should be organized as separate Class
18* 1. wtset should be accessible through getter
19* 2. ISR+FSR and ISR parts beter separated,
if possible.
21* The auxiliary weights
for variable-
weight runs are:
22* wt=wtcru1*wtcru2*wtset(i) where
23* wtset(i) = ypar(500+i) are auxiliary weights
25* wtset(71) =
zero-th order initial+
final
26* wtset(72) = first order initial+
final
27* wtset(73) = second order initial+
final
28* wtset(74) = third order initial+
final
29* and the following provide the corresponding
components of
x-section.
30* wtset(80-84) = first order beta
's
31* wtset(90-97) = second order beta's
32* wtset(100-111) = third order beta
's
33* furthermore, for the initial state alone we provide:
34* wtset( 1) = zero-th order initial
35* wtset( 2) = first order initial
36* wtset( 3) = second order initial
37* wtset( 4) = third order initial
38* and the corresponding components
39* wtset(10-11) = first order beta's
40* wtset(20-22) = second order beta
's
41* wtset(30-33) = third order beta's
45* Below limits m_vlim1/2 the contribution from beta1/2 is not calculated,
55 $ m_xSfac,m_ySfac ! Soft factors
59 DOUBLE
PRECISION m_xBet10, m_yBet10 ! O(alf1)
60 DOUBLE
PRECISION m_Beta02, m_xBet11, m_yBet11 ! O(alf2)
61 DOUBLE
PRECISION m_xxBet20, m_xyBet20, m_yyBet20 ! O(alf2)
62 DOUBLE
PRECISION m_Beta03, m_xBet12, m_yBet12 ! O(alf3)
63 DOUBLE
PRECISION m_xxBet21, m_xyBet21, m_yyBet21 ! O(alf3)
64 DOUBLE
PRECISION m_xxxBet30,m_xxyBet30,m_xyyBet30 ! O(alf3)
67 DOUBLE
PRECISION m_beti01,m_sbti10 ! O(alf1)
68 DOUBLE
PRECISION m_beti02,m_sbti11,m_sbti20 ! O(alf2)
69 DOUBLE
PRECISION m_beti03,m_sbti12,m_sbti21,m_sbti30 ! O(alf2)
80 $ m_xSfac(npmx), ! ISR soft factors
81 $ m_ySfac(npmx), ! FSR soft factors
82 $ m_Beta00, ! beta0 O(alf0) ISR+FSR
83 $ m_Beta01, m_xBet10, m_yBet10, ! beta1 O(alf1) ISR+FSR
84 $ m_Beta02, ! beta0 O(alf2) ISR+FSR
85 $ m_xBet11, m_yBet11, ! beta1 O(alf2) ISR+FSR
86 $ m_xxBet20, m_xyBet20, m_yyBet20, ! beta2 O(alf2) ISR+FSR
87 $ m_Beta03, ! beta0 O(alf3) ISR+FSR
88 $ m_xBet12, m_yBet12, ! beta1 O(alf3) ISR+FSR
89 $ m_xxBet21, m_xyBet21, m_yyBet21, ! beta2 O(alf3) ISR+FSR
90 $ m_xxxBet30,m_xxyBet30,m_xyyBet30, ! beta3 O(alf3) ISR+FSR
91 $ m_beti00, ! O(alf0) pure ISR
92 $ m_beti01,m_sbti10, ! O(alf1) pure ISR
93 $ m_beti02,m_sbti11,m_sbti20, ! O(alf2) pure ISR
94 $ m_beti03,m_sbti12,m_sbti21,m_sbti30, ! O(alf2) pure ISR
95 $ m_dis0,m_dis1,m_dis2, ! UNEXP O(alf2) ISR
96 $ m_dig1, ! UNEXP auxiliary
#define BHLUMI(MODE, XPAR, NPAR)
*******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_KeyISR
*********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_WtBest
*********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_WtSet
*********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 $ !crude photon multiplicity enhancement factor *EVENT $ !MC crude volume of PhhSpace *Sfactors $ !YFS formfactor IR part only $ !YFS formfactor non IR finite part $ !mass weight
*********Class see also m_nmax DOUBLE PRECISION m_alfinv
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_vvmin
********Version nd August
********Version nd by S Jadach *Matrix element for fermion pair production
*************DOUBLE PRECISION m_pi *DOUBLE PRECISION m_HvecTau2 DOUBLE PRECISION m_HvClone2 DOUBLE PRECISION m_gamma1 DOUBLE PRECISION m_gamma2 DOUBLE PRECISION m_thet1 DOUBLE PRECISION m_thet2 INTEGER m_IFPHOT *COMMON c_Taupair $ !Spin Polarimeter vector first Tau $ !Spin Polarimeter vector second Tau $ !Clone Spin Polarimeter vector first Tau $ !Clone Spin Polarimeter vector second Tau $ !Random Euler angle for cloning st tau $ !Random Euler angle for cloning st tau $ !Random Euler angle for cloning st tau $ !Random Euler angle for cloning nd tau $ !Random Euler angle for cloning nd tau $ !Random Euler angle for cloning nd tau $ !phi of HvecTau1 $ !theta of HvecTau1 $ !phi of HvecTau2 $ !theta of HvecTau2 $ !super for inhibiting all tauola activity $ !key for PHOTOS $ !switch for cloning procedure
IMPLICIT REAL *A O Z IMPLICIT INTEGER(I-N) C.... COMMONBALOCKS FOR MODE COMMON/MODEXS/ICH common/bornXS/xborn REAL *8 USER_XX(1000)
uint32_t components(const eformat::FullEventFragment< TPointer > &f)
1.9.6