CGEM BOSS 6.6.5.g
BESIII Offline Software System
Loading...
Searching...
No Matches
QED3.h
Go to the documentation of this file.
1*/////////////////////////////////////////////////////////////////////////////////////
2*// //
3*// Pseudo-CLASS QED3 //
4*// //
5*// Purpose: QED bremsstrahlung matrix element up to O(alfa^3) //
6*// //
7*/////////////////////////////////////////////////////////////////////////////////////
8* Version 2-nd August 1997, by S. Jadach.
9* Matrix element for fermion pair production,
10* ISR O(alf3) + FSR O(alf2).
11* New features:
12* 1. Kinematical reduction procedure kicked out.
13* 2. Numerical stability for collinear photons improved up to BHLUMI level.
14* 3. O(alf3)LL for ISR
15*****
16* To be done:
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.
20*/////////////////////////////////////////////////////////////////////////////////////
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
24*
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
42*/////////////////////////////////////////////////////////////////////////////////////
43 DOUBLE PRECISION pi
44 PARAMETER( pi=3.1415926535897932d0 )
45* Below limits m_vlim1/2 the contribution from beta1/2 is not calculated,
46 DOUBLE PRECISION m_vlim1,m_vlim2
47 PARAMETER(m_vlim1= 1.d-9, m_vlim2= 1.d-9)
48 INTEGER npmx ! maximum number of photons
49 PARAMETER( npmx=100)
50 INTEGER m_lenwt ! max. no. of aux weights
51 PARAMETER (m_lenwt=200)
52*/////////////////////////////////////////////////////////////////////////////////////
53* amplitude components
54 DOUBLE PRECISION
55 $ m_xSfac,m_ySfac ! Soft factors
56* ISR+FSR
57 DOUBLE PRECISION m_Beta00 ! O(alf0)
58 DOUBLE PRECISION m_Beta01 ! O(alf1)
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)
65* pure ISR
66 DOUBLE PRECISION m_beti00 ! O(alf0)
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)
70* UNEXP
71 DOUBLE PRECISION m_dis0,m_dis1,m_dis2
72 DOUBLE PRECISION m_dig1
73*
74 DOUBLE PRECISION m_WtSet,m_WtBest,m_alfinv,m_vvmin
75 INTEGER m_IdeWgt,m_KeyISR
76*
77 COMMON /c_QED3/
78 $ m_WtSet(m_lenwt),
79 $ m_WtBest,
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
97 $ m_alfinv,
98 $ m_vvmin,
99 $ m_IdeWgt,
100 $ m_KeyISR
101*
102 SAVE /c_QED3/
103*
104*/////////////////////////////////////////////////////////////////////////////////////
105*// //
106*// End of Pseudo-CLASS QED3 //
107*// //
108*/////////////////////////////////////////////////////////////////////////////////////
BHLUMI
#define BHLUMI(MODE, XPAR, NPAR)
Definition: Bhlumi.cxx:61
x
Double_t x[10]
Definition: DataBase/tau_mode.c:57
zero
const DifNumber zero
PRECISION
#define PRECISION
Definition: EvtDalitzReso.cc:31
maximum
*******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 maximum
Definition: FoamA.h:82
PARAMETER
*******INTEGER m_nBinMax INTEGER m_NdiMax PARAMETER(m_NdiMax=5) ! maximum dimension for siplices INTEGER m_KdiMax PARAMETER(m_KdiMax
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_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
Definition: GPS.h:34
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_WtBest
Definition: GPS.h:26
m_WtSet
*********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
Definition: GPS.h:26
s
XmlRpcServer s
Definition: HelloServer.cpp:11
weight
*********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
Definition: KarFin.h:34
m_alfinv
*********Class see also m_nmax DOUBLE PRECISION m_alfinv
Definition: KarFin.h:12
m_vvmin
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_vvmin
Definition: KarLud.h:15
August
********Version nd August
Definition: QED3.h:8
production
********Version nd by S Jadach *Matrix element for fermion pair production
Definition: QED3.h:9
procedure
*************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
Definition: Taupair.h:44
INTEGER
IMPLICIT REAL *A O Z IMPLICIT INTEGER(I-N) C.... COMMONBALOCKS FOR MODE COMMON/MODEXS/ICH common/bornXS/xborn REAL *8 USER_XX(1000)
components
uint32_t components(const eformat::FullEventFragment< TPointer > &f)
Definition: validate.cxx:29
pi
const float pi
Definition: vector3.h:133