Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-08-05 08:15:44

0001     
0002 C*********************************************************************  
0003     
0004       BLOCK DATA PYHIDATA
0005     
0006 C...Give sensible default values to all status codes and parameters.    
0007       COMMON/PYHISUBS/MSEL,MSUB(200),KFIN(2,-40:40),CKIN(200) 
0008       SAVE /PYHISUBS/ 
0009       COMMON/PYHIPARS/MSTP(200),PARP(200),MSTI(200),PARI(200) 
0010       SAVE /PYHIPARS/ 
0011       COMMON/PYHIINT1/MINT(400),VINT(400) 
0012       SAVE /PYHIINT1/ 
0013       COMMON/PYHIINT2/ISET(200),KFPR(200,2),COEF(200,20),ICOL(40,4,2) 
0014       SAVE /PYHIINT2/ 
0015       COMMON/PYHIINT3/XSFX(2,-40:40),ISIG(1000,3),SIGH(1000)  
0016       SAVE /PYHIINT3/ 
0017       COMMON/PYHIINT4/WIDP(21:40,0:40),WIDE(21:40,0:40),WIDS(21:40,3) 
0018       SAVE /PYHIINT4/ 
0019       COMMON/PYHIINT5/NGEN(0:200,3),XSEC(0:200,3) 
0020       SAVE /PYHIINT5/ 
0021       COMMON/PYHIINT6/PROC(0:200) 
0022       CHARACTER PROC*28 
0023       SAVE /PYHIINT6/ 
0024     
0025 C...Default values for allowed processes and kinematics constraints.    
0026       DATA MSEL/1/  
0027       DATA MSUB/200*0/  
0028       DATA ((KFIN(I,J),J=-40,40),I=1,2)/40*1,0,80*1,0,40*1/ 
0029       DATA CKIN/    
0030      &   2.0, -1.0,  0.0, -1.0,  1.0,  1.0, -10.,  10., -10.,  10., 
0031      1  -10.,  10., -10.,  10., -10.,  10., -1.0,  1.0, -1.0,  1.0, 
0032      2   0.0,  1.0,  0.0,  1.0, -1.0,  1.0, -1.0,  1.0,   0.,   0., 
0033      3   2.0, -1.0,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0034      4   160*0./    
0035     
0036 C...Default values for main switches and parameters. Reset information. 
0037       DATA (MSTP(I),I=1,100)/   
0038      &     3,    1,    2,    0,    0,    0,    0,    0,    0,    0, 
0039      1     0,    0,    0,    0,    0,    0,    0,    0,    0,    0, 
0040      2     0,    0,    0,    0,    0,    0,    0,    0,    0,    0, 
0041      3     1,    2,    0,    0,    0,    2,    0,    0,    0,    0, 
0042      4     1,    0,    3,    7,    1,    0,    0,    0,    0,    0, 
0043      5     1,    1,   20,    6,    0,    0,    0,    0,    0,    0, 
0044      6     1,    2,    2,    2,    1,    0,    0,    0,    0,    0, 
0045      7     1,    0,    0,    0,    0,    0,    0,    0,    0,    0, 
0046      8     1,    1,  100,    0,    0,    0,    0,    0,    0,    0, 
0047      9     1,    4,    0,    0,    0,    0,    0,    0,    0,    0/ 
0048       DATA (MSTP(I),I=101,200)/ 
0049      &     1,    0,    0,    0,    0,    0,    0,    0,    0,    0, 
0050      1     1,    1,    1,    0,    0,    0,    0,    0,    0,    0, 
0051      2     0,    1,    2,    1,    1,   20,    0,    0,    0,    0, 
0052      3     0,    4,    0,    1,    0,    0,    0,    0,    0,    0, 
0053      4     0,    0,    0,    0,    0,    0,    0,    0,    0,    0, 
0054      5     0,    0,    0,    0,    0,    0,    0,    0,    0,    0, 
0055      6     0,    0,    0,    0,    0,    0,    0,    0,    0,    0, 
0056      7     0,    0,    0,    0,    0,    0,    0,    0,    0,    0, 
0057      8     5,    3, 1989,   11,   24,    0,    0,    0,    0,    0, 
0058      9     0,    0,    0,    0,    0,    0,    0,    0,    0,    0/ 
0059       DATA (PARP(I),I=1,100)/   
0060      &  0.25,  10.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0061      1    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0062      2    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0063      3   1.5,  2.0, 0.075,  0.,  0.2,   0.,   0.,   0.,   0.,   0., 
0064      4    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0065      5   1.0, 2.26, 1.E4, 1.E-4,  0.,   0.,   0.,   0.,   0.,   0., 
0066      6  0.25,  1.0, 0.25,  1.0,  2.0, 1.E-3, 4.0,   0.,   0.,   0., 
0067      7   4.0,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0068      8   1.6, 1.85,  0.5,  0.2, 0.33, 0.66,  0.7,  0.5,   0.,   0., 
0069      9  0.44, 0.44,  2.0,  1.0,   0.,  3.0,  1.0, 0.75,   0.,   0./ 
0070       DATA (PARP(I),I=101,200)/ 
0071      & -0.02,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0072      1   2.0,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0073      2   0.4,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0074      3  0.01,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0075      4    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0076      5    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0077      6    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0078      7    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0079      8    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0., 
0080      9    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0./ 
0081       DATA MSTI/200*0/  
0082       DATA PARI/200*0./ 
0083       DATA MINT/400*0/  
0084       DATA VINT/400*0./ 
0085     
0086 C...Constants for the generation of the various processes.  
0087       DATA (ISET(I),I=1,100)/   
0088      &    1,    1,    1,   -1,    3,   -1,   -1,    3,   -2,   -2,  
0089      1    2,    2,    2,    2,    2,    2,   -1,    2,    2,    2,  
0090      2   -1,    2,    2,    2,    2,    2,   -1,    2,    2,    2,  
0091      3    2,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,  
0092      4   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,  
0093      5   -1,   -1,    2,   -1,   -1,   -1,   -1,   -1,   -1,   -1,  
0094      6   -1,   -1,   -1,   -1,   -1,   -1,   -1,    2,   -1,   -1,  
0095      7    4,    4,    4,   -1,   -1,    4,    4,   -1,   -1,   -2,  
0096      8    2,    2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,  
0097      9    0,    0,    0,   -1,    0,    5,   -2,   -2,   -2,   -2/  
0098       DATA (ISET(I),I=101,200)/ 
0099      &   -1,    1,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,  
0100      1    2,    2,    2,    2,   -1,   -1,   -1,   -2,   -2,   -2,  
0101      2   -1,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,  
0102      3   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,  
0103      4    1,    1,    1,   -2,   -2,   -2,   -2,   -2,   -2,   -2,  
0104      5   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,  
0105      6    2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,  
0106      7   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,  
0107      8   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,  
0108      9   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2/  
0109       DATA ((KFPR(I,J),J=1,2),I=1,50)/  
0110      &   23,    0,   24,    0,   25,    0,   24,    0,   25,    0,  
0111      &   24,    0,   23,    0,   25,    0,    0,    0,    0,    0,  
0112      1    0,    0,    0,    0,   21,   21,   21,   22,   21,   23,  
0113      1   21,   24,   21,   25,   22,   22,   22,   23,   22,   24,  
0114      2   22,   25,   23,   23,   23,   24,   23,   25,   24,   24,  
0115      2   24,   25,   25,   25,    0,   21,    0,   22,    0,   23,  
0116      3    0,   24,    0,   25,    0,   21,    0,   22,    0,   23,  
0117      3    0,   24,    0,   25,    0,   21,    0,   22,    0,   23,  
0118      4    0,   24,    0,   25,    0,   21,    0,   22,    0,   23,  
0119      4    0,   24,    0,   25,    0,   21,    0,   22,    0,   23/  
0120       DATA ((KFPR(I,J),J=1,2),I=51,100)/    
0121      5    0,   24,    0,   25,    0,    0,    0,    0,    0,    0,  
0122      5    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0123      6    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0124      6    0,    0,    0,    0,   21,   21,   24,   24,   22,   24,  
0125      7   23,   23,   24,   24,   23,   24,   23,   25,   22,   22,  
0126      7   23,   23,   24,   24,   24,   25,   25,   25,    0,    0,  
0127      8    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0128      8    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0129      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0130      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0/  
0131       DATA ((KFPR(I,J),J=1,2),I=101,150)/   
0132      &   23,    0,   25,    0,    0,    0,    0,    0,    0,    0,  
0133      &    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0134      1   21,   25,    0,   25,   21,   25,   22,   22,   22,   23,  
0135      1   23,   23,   24,   24,    0,    0,    0,    0,    0,    0,  
0136      2    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0137      2    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0138      3    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0139      3    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0140      4   32,    0,   37,    0,   40,    0,    0,    0,    0,    0,  
0141      4    0,    0,    0,    0,    0,    0,    0,    0,    0,    0/  
0142       DATA ((KFPR(I,J),J=1,2),I=151,200)/   
0143      5    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0144      5    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0145      6    0,   37,    0,    0,    0,    0,    0,    0,    0,    0,  
0146      6    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0147      7    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0148      7    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0149      8    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0150      8    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0151      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,  
0152      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0/  
0153       DATA COEF/4000*0./    
0154       DATA (((ICOL(I,J,K),K=1,2),J=1,4),I=1,40)/    
0155      1 4,0,3,0,2,0,1,0,3,0,4,0,1,0,2,0,2,0,0,1,4,0,0,3,3,0,0,4,1,0,0,2, 
0156      2 3,0,0,4,1,4,3,2,4,0,0,3,4,2,1,3,2,0,4,1,4,0,2,3,4,0,3,4,2,0,1,2, 
0157      3 3,2,1,0,1,4,3,0,4,3,3,0,2,1,1,0,3,2,1,4,1,0,0,2,2,4,3,1,2,0,0,1, 
0158      4 3,2,1,4,1,4,3,2,4,2,1,3,4,2,1,3,3,4,4,3,1,2,2,1,2,0,3,1,2,0,0,0, 
0159      5 4,2,1,0,0,0,1,0,3,0,0,3,1,2,0,0,4,0,0,4,0,0,1,2,2,0,0,1,4,4,3,3, 
0160      6 2,2,1,1,4,4,3,3,3,3,4,4,1,1,2,2,3,2,1,3,1,2,0,0,4,2,1,4,0,0,1,2, 
0161      7 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0162      8 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0163      9 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0164      & 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0/ 
0165     
0166 C...Character constants: name of processes. 
0167       DATA PROC(0)/                    'All included subprocesses   '/  
0168       DATA (PROC(I),I=1,20)/    
0169      1'f + fb -> gamma*/Z0         ',  'f + fb'' -> W+/-             ', 
0170      2'f + fb -> H0                ',  'gamma + W+/- -> W+/-        ',  
0171      3'Z0 + Z0 -> H0               ',  'Z0 + W+/- -> W+/-           ',  
0172      4'                            ',  'W+ + W- -> H0               ',  
0173      5'                            ',  '                            ',  
0174      6'f + f'' -> f + f''            ','f + fb -> f'' + fb''          ',    
0175      7'f + fb -> g + g             ',  'f + fb -> g + gamma         ',  
0176      8'f + fb -> g + Z0            ',  'f + fb'' -> g + W+/-         ', 
0177      9'f + fb -> g + H0            ',  'f + fb -> gamma + gamma     ',  
0178      &'f + fb -> gamma + Z0        ',  'f + fb'' -> gamma + W+/-     '/ 
0179       DATA (PROC(I),I=21,40)/   
0180      1'f + fb -> gamma + H0        ',  'f + fb -> Z0 + Z0           ',  
0181      2'f + fb'' -> Z0 + W+/-        ', 'f + fb -> Z0 + H0           ',  
0182      3'f + fb -> W+ + W-           ',  'f + fb'' -> W+/- + H0        ', 
0183      4'f + fb -> H0 + H0           ',  'f + g -> f + g              ',  
0184      5'f + g -> f + gamma          ',  'f + g -> f + Z0             ',  
0185      6'f + g -> f'' + W+/-          ', 'f + g -> f + H0             ',  
0186      7'f + gamma -> f + g          ',  'f + gamma -> f + gamma      ',  
0187      8'f + gamma -> f + Z0         ',  'f + gamma -> f'' + W+/-      ', 
0188      9'f + gamma -> f + H0         ',  'f + Z0 -> f + g             ',  
0189      &'f + Z0 -> f + gamma         ',  'f + Z0 -> f + Z0            '/  
0190       DATA (PROC(I),I=41,60)/   
0191      1'f + Z0 -> f'' + W+/-         ', 'f + Z0 -> f + H0            ',  
0192      2'f + W+/- -> f'' + g          ', 'f + W+/- -> f'' + gamma      ', 
0193      3'f + W+/- -> f'' + Z0         ', 'f + W+/- -> f'' + W+/-       ', 
0194      4'f + W+/- -> f'' + H0         ', 'f + H0 -> f + g             ',  
0195      5'f + H0 -> f + gamma         ',  'f + H0 -> f + Z0            ',  
0196      6'f + H0 -> f'' + W+/-         ', 'f + H0 -> f + H0            ',  
0197      7'g + g -> f + fb             ',  'g + gamma -> f + fb         ',  
0198      8'g + Z0 -> f + fb            ',  'g + W+/- -> f + fb''         ', 
0199      9'g + H0 -> f + fb            ',  'gamma + gamma -> f + fb     ',  
0200      &'gamma + Z0 -> f + fb        ',  'gamma + W+/- -> f + fb''     '/ 
0201       DATA (PROC(I),I=61,80)/   
0202      1'gamma + H0 -> f + fb        ',  'Z0 + Z0 -> f + fb           ',  
0203      2'Z0 + W+/- -> f + fb''        ', 'Z0 + H0 -> f + fb           ',  
0204      3'W+ + W- -> f + fb           ',  'W+/- + H0 -> f + fb''        ', 
0205      4'H0 + H0 -> f + fb           ',  'g + g -> g + g              ',  
0206      5'gamma + gamma -> W+ + W-    ',  'gamma + W+/- -> gamma + W+/-',  
0207      6'Z0 + Z0 -> Z0 + Z0          ',  'Z0 + Z0 -> W+ + W-          ',  
0208      7'Z0 + W+/- -> Z0 + W+/-      ',  'Z0 + Z0 -> Z0 + H0          ',  
0209      8'W+ + W- -> gamma + gamma    ',  'W+ + W- -> Z0 + Z0          ',  
0210      9'W+/- + W+/- -> W+/- + W+/-  ',  'W+/- + H0 -> W+/- + H0      ',  
0211      &'H0 + H0 -> H0 + H0          ',  '                            '/  
0212       DATA (PROC(I),I=81,100)/  
0213      1'q + qb -> Q + QB, massive   ',  'g + g -> Q + QB, massive    ',  
0214      2'                            ',  '                            ',  
0215      3'                            ',  '                            ',  
0216      4'                            ',  '                            ',  
0217      5'                            ',  '                            ',  
0218      6'Elastic scattering          ',  'Single diffractive          ',  
0219      7'Double diffractive          ',  'Central diffractive         ',  
0220      8'Low-pT scattering           ',  'Semihard QCD 2 -> 2         ',  
0221      9'                            ',  '                            ',  
0222      &'                            ',  '                            '/  
0223       DATA (PROC(I),I=101,120)/ 
0224      1'g + g -> gamma*/Z0          ',  'g + g -> H0                 ',  
0225      2'                            ',  '                            ',  
0226      3'                            ',  '                            ',  
0227      4'                            ',  '                            ',  
0228      5'                            ',  '                            ',  
0229      6'f + fb -> g + H0            ',  'q + g -> q + H0             ',  
0230      7'g + g -> g + H0             ',  'g + g -> gamma + gamma      ',  
0231      8'g + g -> gamma + Z0         ',  'g + g -> Z0 + Z0            ',  
0232      9'g + g -> W+ + W-            ',  '                            ',  
0233      &'                            ',  '                            '/  
0234       DATA (PROC(I),I=121,140)/ 
0235      1'g + g -> f + fb + H0        ',  '                            ',  
0236      2'                            ',  '                            ',  
0237      3'                            ',  '                            ',  
0238      4'                            ',  '                            ',  
0239      5'                            ',  '                            ',  
0240      6'                            ',  '                            ',  
0241      7'                            ',  '                            ',  
0242      8'                            ',  '                            ',  
0243      9'                            ',  '                            ',  
0244      &'                            ',  '                            '/  
0245       DATA (PROC(I),I=141,160)/ 
0246      1'f + fb -> gamma*/Z0/Z''0     ', 'f + fb'' -> H+/-             ', 
0247      2'f + fb -> R                 ',  '                            ',  
0248      3'                            ',  '                            ',  
0249      4'                            ',  '                            ',  
0250      5'                            ',  '                            ',  
0251      6'                            ',  '                            ',  
0252      7'                            ',  '                            ',  
0253      8'                            ',  '                            ',  
0254      9'                            ',  '                            ',  
0255      &'                            ',  '                            '/  
0256       DATA (PROC(I),I=161,180)/ 
0257      1'f + g -> f'' + H+/-          ', '                            ',  
0258      2'                            ',  '                            ',  
0259      3'                            ',  '                            ',  
0260      4'                            ',  '                            ',  
0261      5'                            ',  '                            ',  
0262      6'                            ',  '                            ',  
0263      7'                            ',  '                            ',  
0264      8'                            ',  '                            ',  
0265      9'                            ',  '                            ',  
0266      &'                            ',  '                            '/  
0267       DATA (PROC(I),I=181,200)/     20*'                            '/  
0268     
0269       END