! ! PROTON STRUCTURE FUNCTION SELECTION FROM PDFLIB ! ! SET MSTP(52) = 2 FOR PDFLIB CHOICE ! SET MSTP(51) = 1000 * NGROUP + NSET ! EXAMPLES: ! MRSD-=3031 ; MRSG=3041 ! CTEQ2M=4018 ; CTEQ2PM=4024 ; CTEQ3M=4030 ! GRVH0=5002 ; GRV94H0=5006 ! MSTP(51)=4030 MSTP(52)=2 ! ! CHANGING RANDOM NUMBER GENERATOR SEED ! MRPY(1) = RANSEED2 ! RANDOM SEED FIELD ! MSEL = 5 ! b bbar ! ! Kinematical cuts ! CKIN(3)=PTLOW ! Minimum Pt for hard 2->2 process. CKIN(4)=PTHIGH ! Maximum Pt for hard 2->2 process. CKIN(9)=-4.2 ! Rapidity range for product with CKIN(10)=4.2 ! largest rapidity in hard 2->2 process. ! Inhibit Jpsi decays MDME(938,1) = 1 ! 0 0.060200 e- e+ MDME(939,1) = 1 ! 0 0.060100 mu- mu+ MDME(940,1) = 0 ! 12 0.879700 rndmflav rnd ! Inhibit B0 and B+ decays - but we want b->j/psi with bbar unforced ! Neutral B decays MDME(943,1) = 2 ! 42 0.020000 nu_e e+ D- MDME(944,1) = 2 ! 42 0.055000 nu_e e+ D*- MDME(945,1) = 2 ! 42 0.005000 nu_e e+ D_1- MDME(946,1) = 2 ! 42 0.005000 nu_e e+ D*_0 MDME(947,1) = 2 ! 42 0.008000 nu_e e+ D*_1 MDME(948,1) = 2 ! 42 0.012000 nu_e e+ D*_2 MDME(949,1) = 2 ! 42 0.020000 nu_mu mu+ D- MDME(950,1) = 2 ! 42 0.055000 nu_mu mu+ D*- MDME(951,1) = 2 ! 42 0.005000 nu_mu mu+ D_1- MDME(952,1) = 2 ! 42 0.005000 nu_mu mu+ D*_0 MDME(953,1) = 2 ! 42 0.008000 nu_mu mu+ D*_1 MDME(954,1) = 2 ! 42 0.012000 nu_mu mu+ D*_2 MDME(955,1) = 2 ! 42 0.010000 nu_tau tau+ D- MDME(956,1) = 2 ! 42 0.030000 nu_tau tau+ D*- MDME(957,1) = 2 ! 0 0.003500 D- pi+ MDME(958,1) = 2 ! 0 0.011000 D- rho+ MDME(959,1) = 2 ! 0 0.005500 D- a_1+ MDME(960,1) = 2 ! 0 0.004200 D*- pi+ MDME(961,1) = 2 ! 0 0.009000 D*- rho+ MDME(962,1) = 2 ! 0 0.018000 D*- a_1+ MDME(963,1) = 2 ! 0 0.015000 D- D_s+ MDME(964,1) = 2 ! 0 0.018500 D- D*_s+ MDME(965,1) = 2 ! 0 0.013500 D*- D_s+ MDME(966,1) = 2 ! 0 0.025000 D*- D*_s+ MDME(967,1) = 2 ! 0 0.000400 eta_c K0 MDME(968,1) = 2 ! 0 0.000700 eta_c K*0 MDME(969,1) = 1 ! 0 0.000800 J/psi K0 MDME(970,1) = 1 ! 0 0.001400 J/psi K*0 MDME(971,1) = 2 ! 0 0.001900 chi_1c K0 MDME(972,1) = 2 ! 0 0.002500 chi_1c K*0 MDME(973,1) = 2 ! 48 0.429100 u dbar cbar MDME(974,1) = 2 ! 13 0.080000 u cbar dbar MDME(975,1) = 2 ! 13 0.070000 c sbar cbar MDME(976,1) = 2 ! 13 0.020000 c cbar sbar MDME(977,1) = 2 ! 42 0.015000 u dbar ubar MDME(978,1) = 2 ! 42 0.005000 c sbar ubar ! Charged B decays MDME(988,1) = 2 ! 42 0.020000 nu_e e+ Dbar MDME(989,1) = 2 ! 42 0.055000 nu_e e+ D*ba MDME(990,1) = 2 ! 42 0.005000 nu_e e+ D_1b MDME(991,1) = 2 ! 42 0.005000 nu_e e+ D*_0 MDME(992,1) = 2 ! 42 0.008000 nu_e e+ D*_1 MDME(993,1) = 2 ! 42 0.012000 nu_e e+ D*_2 MDME(994,1) = 2 ! 42 0.020000 nu_mu mu+ Dbar MDME(995,1) = 2 ! 42 0.055000 nu_mu mu+ D*ba MDME(996,1) = 2 ! 42 0.005000 nu_mu mu+ D_1b MDME(997,1) = 2 ! 42 0.005000 nu_mu mu+ D*_0 MDME(998,1) = 2 ! 42 0.008000 nu_mu mu+ D*_1 MDME(999,1) = 2 ! 42 0.012000 nu_mu mu+ D*_2 MDME(1000,1) = 2 ! 42 0.010000 nu_tau tau+ Dbar MDME(1001,1) = 2 ! 42 0.030000 nu_tau tau+ D*ba MDME(1002,1) = 2 ! 0 0.003500 Dbar0 pi+ MDME(1003,1) = 2 ! 0 0.011000 Dbar0 rho+ MDME(1004,1) = 2 ! 0 0.005500 Dbar0 a_1+ MDME(1005,1) = 2 ! 0 0.004200 D*bar0 pi+ MDME(1006,1) = 2 ! 0 0.009000 D*bar0 rho+ MDME(1007,1) = 2 ! 0 0.018000 D*bar0 a_1+ MDME(1008,1) = 2 ! 0 0.015000 Dbar0 D_s+ MDME(1009,1) = 2 ! 0 0.018500 Dbar0 D*_s+ MDME(1010,1) = 2 ! 0 0.013500 D*bar0 D_s+ MDME(1011,1) = 2 ! 0 0.025000 D*bar0 D*_s+ MDME(1012,1) = 2 ! 0 0.000400 eta_c K+ MDME(1013,1) = 2 ! 0 0.000700 eta_c K*+ MDME(1014,1) = 1 ! 0 0.000800 J/psi K+ MDME(1015,1) = 1 ! 0 0.001400 J/psi K*+ MDME(1016,1) = 2 ! 0 0.001900 chi_1c K+ MDME(1017,1) = 2 ! 0 0.002500 chi_1c K*+ MDME(1018,1) = 2 ! 48 0.429100 u dbar cbar MDME(1019,1) = 2 ! 13 0.080000 u cbar dbar MDME(1020,1) = 2 ! 13 0.070000 c sbar cbar MDME(1021,1) = 2 ! 13 0.020000 c cbar sbar MDME(1022,1) = 2 ! 42 0.015000 u dbar ubar MDME(1023,1) = 2 ! 42 0.005000 c sbar ubar MSTP(81)=0 !TURN OFF MULTIPLE INTERACTIONS ! MDCY(C111,1)=0 !SETS PI ZERO STABLE (LIMITS SIZE OF OUTPUT) MDCY(15,1)=0 ! sets tau stable - this is important (for tauola)! ! The next line is _mandatory_ for RunII production MSTJ(22)=2 !A particle is decayed only if its proper !lifetime is smaller than PARJ(71) ! default should be ok. (strange particles stable, charm & bottom ! decay).