25 PARAMETER ( m_NpairMax = m_NdiMax*(m_NdiMax+1)/2 ) !
no. of pairs of vertices
27 PARAMETER ( m_NdiviMax = m_NpairMax+m_KdiMax ) ! Total
no. of division samplings
29 INTEGER m_CeStat, m_CePare, m_CeDau1, m_CeDau2, m_CeSamp, m_CeVert, m_CeBest !
30 DOUBLE
PRECISION m_CeIntg, m_CeSum, m_CeDriv, m_CePrim, m_CePrCu, m_CeXave, m_CeVolu !
38 DOUBLE
PRECISION m_MCresult, m_MCerror, m_MCwt, m_MCvector !
43 $ m_CeStat(m_nBufMax), ! Cell member: status=0 inactive, =1 active
44 $ m_CePare(m_nBufMax), ! Cell member: parent cell pointer
45 $ m_CeDau1(m_nBufMax), ! Cell member: daughter1 cell pointer
46 $ m_CeDau2(m_nBufMax), ! Cell member: daughter2 cell pointer
47 $ m_CeSamp(m_nBufMax), ! Cell member: No of MC events in exploration
48 $ m_CeVert(m_nBufMax,m_NdiMax+1),! Cell member: simplex vertex pointers
49 $ m_CeVer1(m_nBufMax,m_KdiMax), ! Cell member: hypercubic vertex position
50 $ m_CeVer2(m_nBufMax,m_KdiMax), ! Cell member: hypercubic vertex size
51 $ m_CeIntg(m_nBufMax), ! Cell member: integral estimator
52 $ m_CeDriv(m_nBufMax), ! Cell member:
Drive integral estimate, from build-up
53 $ m_CePrim(m_nBufMax), ! Cell member: Primary integral estimate, MC generation
54 $ m_CePrCu(0:m_nBufMax), ! Cell member: Cumulative Primary
55 $ m_CeVolu(m_nBufMax), ! Cell member: Cartesian volume
56 $ m_CeXave(m_nBufMax), ! Cell member: Average best X
57 $ m_CeBest(m_nBufMax), ! Cell member: Best pair of vertices, pointer
58 $ m_CeSum( m_nBufMax,m_sMax), ! Cell member:
weight summaries
59 $ m_VerX( m_vMax, m_NdiMax), ! List of all VERTEX positions
60 $ m_ActC(
m_cMax), ! List of all pointers to ACTIVE cells
61 $ m_VolTot, ! Estimate of Volume total, without error
62 $ m_Drive, ! M.
C. generation
Drive value of integral
63 $ m_SumWt, ! M.
C. generation sum of Wt
64 $ m_SumWt2, ! M.
C. generation sum of Wt**2
66 $ m_WtMax, ! M.
C. generation
maximum wt
67 $ m_WtMin, ! M.
C. generation minimum wt
68 $ m_MCresult, ! M.
C. generation Final value of INTEGRAL
69 $ m_MCerror, ! M.
C. generation Final walue of
ERROR
70 $ m_MCwt, ! M.
C. generation current event
weight
71 $ m_MCvector(m_DimMax),! M.
C. generated vector, dimension=m_NdiMax+ m_KdiMax
72 $ m_Ndim, ! dimension of the siplices
73 $ m_Kdim, ! dimension of the hypercubics
74 $ m_Dimen, ! total dimension of the problem =m_Ndim+m_Kdim
75 $ m_nBuf, ! Actual dynamic lenth of the buffer m_nBuf<m_nBufMax
76 $ m_nBin, ! No of
bins in histogram
for cell exploration/division
80 $
m_nSampl, ! No. of sampling when dividing cell
85 $
m_OptOrd, ! Flag =0 default, entire domain is hyp-
cubic, =1 it is simplex
86 $ m_EvPerBin, ! Maximum effective eevents per
bin, saves
CPU, =0 inactive
87 $ m_OptRanIni, ! Flag =0 default, =1 r.n. generator not initialized in PianA
88 $ m_OptRanLux, ! Flag =-1,0,1,2,3,4 r.n. generator level
92 $ m_MagicInit ! Magic cookie of initialization (global variable)
*******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
*******INTEGER m_nBinMax INTEGER m_NdiMax !No of bins in histogram for cell exploration division $ !Last vertex $ !Last active cell $ m_LastCe
*******INTEGER m_nBinMax INTEGER m_NdiMax !No of bins in histogram for cell exploration division $ !Last vertex $ !Last active cell $ !Last cell in buffer $ m_nSampl
*******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 $ m_Ncalls
*******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 CPU
*******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 choice
*******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 $ m_OptPeek
*******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 TrueVol
*******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 for debug $ m_Out
*******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 $ m_OptOrd
*******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 $ m_OptEdge
*******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 bin
*******INTEGER m_nBinMax INTEGER m_NdiMax !No of bins in histogram for cell exploration division $ m_LastVe
*******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 cubic
*******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 calls
*******INTEGER m_nBinMax INTEGER m_NdiMax !No of bins in histogram for cell exploration division $ !Last vertex $ m_LastAc
*******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 Sigma
*******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 Drive
*******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 $ m_Chat
*******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 Chat
*******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
*******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 $ m_OptDrive
************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
****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_cMax
*********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_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 $ m_NevGen
***************************************************************************************Pseudo Class RRes *****************************************************************************************Parameters and physical constants **Maarten sept ************************************************************************DOUBLE PRECISION xsmu **************************************************************************PARTICLE DATA all others are from PDG *Only resonances with known widths into electron pairs are sept ************************************************************************C Declarations C
***************************************************************************************Pseudo Class RRes *****************************************************************************************Parameters and physical constants **Maarten sept ************************************************************************DOUBLE PRECISION xsmu **************************************************************************PARTICLE DATA all others are from PDG *Only resonances with known widths into electron pairs are * included(others are taken out by setting zero width into electrons). **Maarten Boonekamp
IMPLICIT REAL *A O Z IMPLICIT INTEGER(I-N) C.... COMMONBALOCKS FOR MODE COMMON/MODEXS/ICH common/bornXS/xborn REAL *8 USER_XX(1000)
1.9.6