BOSS 7.0.8
BESIII Offline Software System
Loading...
Searching...
No Matches
KarLud.h
Go to the documentation of this file.
1*//////////////////////////////////////////////////////////////////////////////
2*// //
3*// Pseudo-CLASS KarLud //
4*// //
5*// Purpose: //
6*// Top level Monte-Carlo event generator for ISR radiadion. //
7*// Administrates directly generation of v-variable //
8*// and indirectly of all other variables. //
9*//////////////////////////////////////////////////////////////////////////////
10* Class members:
11 INTEGER m_npmx ! maximum photon multiplicity, see also m_nmax
12 PARAMETER ( m_npmx = 100)
14 DOUBLE PRECISION m_vv, m_x1, m_x2
24*
25 COMMON /c_KarLud/
26 $ m_CMSene, ! Input CMS energy [GeV]
27 $ m_XXXene, ! CMS energy after beam-spread/beam-strahlung [GeV]
28 $ m_DelEne, ! Beam energy spread [GeV]
29 $ m_exe, ! z-boost due to beam spread
31* ff-pair spectrum
32 $ m_vv, ! v = 1-s'/s
33 $ m_x1, ! x=1-z for 1-st beam
34 $ m_x2, ! x=1-z for 2-nd beam
35 $ m_vvmin, ! minimum v, infrared cut
36 $ m_vvmax, ! maximum v
37 $ m_HadMin, ! minimum hadronization mass [GeV]
38* Basic QED
39 $ m_alfinv, ! Alpha_QED at Thomson limit
40* Normalization, MC technicalities
41 $ m_XCrude, ! Crude xsect. from Vesk1
42 $ m_xcgaus, ! from gauss integration
43 $ m_ErGaus, ! gauss precision
44 $ m_Xenph, ! crude photon multiplicity enhancement factor
45* ----------------- EVENT -------------------------
46 $ m_p1(4), ! beam momentum
47 $ m_p2(4), ! beam momentum
48 $ m_q1(4), ! final momentum UNUSED for FSR switched ON
49 $ m_q2(4), ! final momentum UNUSED for FSR switched ON
50 $ m_sphot(m_npmx,4), ! photon momenta
51 $ m_sphum(4), ! sum of photon momenta
52 $ m_PX(4), ! total 4-momentum minus ISR photons
53 $ m_yini(m_npmx), ! Sudakov y-variables
54 $ m_zini(m_npmx), ! Sudakov z-variables
55 $ m_WtBasic, ! Weight associated with V generation
56 $ m_Wt_KF, ! Weight associated with KF generation
57 $ m_WtIni, ! Weight associated with photon momenta generation
58 $ m_WtMass, ! WtIni component
59 $ m_WtDil, ! WtIni component m_WtIni=m_WtMass*m_WtDil*m_WtCut
60 $ m_WtCut, ! WtIni component
61 $ m_nphot, ! Photon multiplicity
62 $ m_KFfin, ! final fermion code
63* -------------------------------------------------
64* Test Switches
65 $ m_KeyZet, ! Z-boson on/off
66 $ m_KeyISR, ! ISR on/off
67 $ m_MltISR, ! fixed photon mult. VERY SPECIAL tests
68 $ m_KeyFix, ! fixed v=vvmav. VERY SPECIAL tests
69 $ m_KeyWtm, ! omitting mass terms in S-factors
70* Miscelaneous
71 $ m_nmax, ! maximum photon multiplicity
72 $ m_out, ! output unit number
73 $ m_NevGen, ! event counter
74 $ m_idyfs ! pointer for histograms
75*
76 SAVE /c_KarLud/
77*//////////////////////////////////////////////////////////////////////////////
78*// //
79*// End of CLASS KarLud //
80*//////////////////////////////////////////////////////////////////////////////
**********INTEGER nmxhep !maximum number of particles DOUBLE PRECISION vhep INTEGER jdahep COMMON hepevt $ !serial number $ !number of particles $ !status code $ !particle ident KF $ !parent particles $ !childreen particles $ !four momentum
double mass
************Class m_alfQCDMZ INTEGER m_KFfin INTEGER m_IVfin INTEGER m_ibox *COMMON c_DZface $ alphaQED at(Q^2=MZ^2) DIZET $ m_alfQCDMZ
@ electron
Definition: DstMdcDedx.h:9
#define PRECISION
double momenta[4][15]
double maximum
Definition: EvtPycont.cc:38
*******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
Definition: FoamA.h:89
*******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
Definition: FoamA.h:90
*******INTEGER m_nBinMax INTEGER m_NdiMax PARAMETER(m_NdiMax=5) ! maximum dimension for siplices INTEGER m_KdiMax PARAMETER(m_KdiMax
XmlRpcServer s
Definition: HelloServer.cpp:11
************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
Definition: KK2f.h:50
************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
Definition: KK2f.h:69
*********Class see also m_nmax DOUBLE PRECISION m_MasPhot DOUBLE PRECISION m_phsu DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_q2
Definition: KarFin.h:15
*********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_idyfs
Definition: KarFin.h:18
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyISR
Definition: KarLud.h:17
**********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_p2
Definition: KarLud.h:20
**********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 $ m_amel
Definition: KarLud.h:30
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_ErGaus
Definition: KarLud.h:16
**********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_q1
Definition: KarLud.h:20
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_vvmin
Definition: KarLud.h:15
**********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_NevGen
Definition: KarLud.h:18
**********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_sphum
Definition: KarLud.h:19
**********Class see also m_nmax DOUBLE PRECISION m_CMSene
Definition: KarLud.h:13
**********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 MC technicalities $ !Crude xsect from Vesk1 $ !from gauss integration $ !gauss precision $ !crude photon multiplicity enhancement factor *EVENT $ !beam momentum $ !beam momentum $ !final momentum UNUSED for FSR switched ON $ !final momentum UNUSED for FSR switched ON $ !photon momenta $ !sum of photon momenta $ !total momentum minus ISR photons $ !Sudakov y variables $ m_zini(m_npmx)
**********Class see also m_nmax DOUBLE PRECISION m_exe
Definition: KarLud.h:13
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_MltISR
Definition: KarLud.h:17
**********Class see also m_nmax DOUBLE PRECISION m_XXXene
Definition: KarLud.h:13
**********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 MC technicalities $ !Crude xsect from Vesk1 $ !from gauss integration $ !gauss precision $ !crude photon multiplicity enhancement factor *EVENT $ !beam momentum $ !beam momentum $ !final momentum UNUSED for FSR switched ON $ !final momentum UNUSED for FSR switched ON $ !photon momenta $ !sum of photon momenta $ !total momentum minus ISR photons $ !Sudakov y variables $ !Sudakov z variables $ !Weight associated with V generation $ m_Wt_KF
Definition: KarLud.h:56
**********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 MC technicalities $ !Crude xsect from Vesk1 $ !from gauss integration $ !gauss precision $ m_Xenph
Definition: KarLud.h:44
**********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_p1
Definition: KarLud.h:20
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_vv
Definition: KarLud.h:14
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_HadMin
Definition: KarLud.h:15
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyFix
Definition: KarLud.h:17
**********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_out
Definition: KarLud.h:18
**********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 MC technicalities $ !Crude xsect from Vesk1 $ !from gauss integration $ !gauss precision $ !crude photon multiplicity enhancement factor *EVENT $ !beam momentum $ !beam momentum $ !final momentum UNUSED for FSR switched ON $ !final momentum UNUSED for FSR switched ON $ !photon momenta $ !sum of photon momenta $ !total momentum minus ISR photons $ !Sudakov y variables $ !Sudakov z variables $ !Weight associated with V generation $ !Weight associated with KF generation $ !Weight associated with photon momenta generation $ !WtIni component $ !WtIni component !WtIni component $ !Photon multiplicity $ m_KFfin
Definition: KarLud.h:62
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_xcgaus
Definition: KarLud.h:16
**********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_yini
Definition: KarLud.h:19
**********Class members
Definition: KarLud.h:11
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyZet
Definition: KarLud.h:17
**********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
Definition: KarLud.h:35
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x1
Definition: KarLud.h:14
**********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_WtIni
Definition: KarLud.h:22
**********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
Definition: KarLud.h:40
**********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_PX
Definition: KarLud.h:19
**********Class see also m_nmax DOUBLE PRECISION m_DelEne
Definition: KarLud.h:13
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_XCrude
Definition: KarLud.h:16
**********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_WtMass
Definition: KarLud.h:22
**********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_WtBasic
Definition: KarLud.h:21
**********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_nmax
Definition: KarLud.h:18
**********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 $ m_alfinv
Definition: KarLud.h:39
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_vvmax
Definition: KarLud.h:15
**********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_nphot
Definition: KarLud.h:23
**********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_WtDil
Definition: KarLud.h:22
**********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_sphot
Definition: KarLud.h:19
****************DOUBLE PRECISION m_vvmin INTEGER m_KeyWtm
Definition: YFSini2.h:19
double y[1000]
IMPLICIT REAL *A O Z IMPLICIT INTEGER(I-N) C.... COMMONBALOCKS FOR MODE COMMON/MODEXS/ICH common/bornXS/xborn REAL *8 USER_XX(1000)