pscut_min = 'PhaseSpace:mHatMin = ' + str(1800*0.7) # pscut_max = 'PhaseSpace:mHatMax = ' + str(1800*1.3) SLHA_TABLE = """## Important note! BLOCK SPINFO 1 FeynHiggs 2 2.10.3 2 built on Mar 12, 2015 BLOCK MODSEL 1 0 # Model 2 1 # GridPts 3 0 # Content 4 0 # RPV 5 0 # CPV 6 0 # FV BLOCK SMINPUTS 1 1.28952896E+02 # invAlfaMZ 2 1.16637000E-05 # GF 3 1.19000000E-01 # AlfasMZ 4 9.11876000E+01 # MZ 5 4.16000000E+00 # Mb 6 1.72500000E+02 # Mt 7 1.77703000E+00 # Mtau 11 5.10998902E-04 # Me 13 1.05658357E-01 # Mmu 21 6.00000000E-03 # Md 22 3.00000000E-03 # Mu 23 9.50000000E-02 # Ms 24 1.28600000E+00 # Mc BLOCK MINPAR 3 3.00000000E+01 # TB BLOCK EXTPAR 0 0.00000000E+00 # Q 1 9.54716519E+01 # M1 2 2.00000000E+02 # M2 3 1.50000000E+03 # M3 11 1.50666667E+03 # At 12 1.50666667E+03 # Ab 13 1.50666667E+03 # Atau 23 2.00000000E+02 # MUE 25 3.00000000E+01 # TB 26 1.80000000E+03 # MA0 27 1.80179462E+03 # MHp 31 5.00000000E+02 # MSL(1) 32 5.00000000E+02 # MSL(2) 33 1.00000000E+03 # MSL(3) 34 5.00000000E+02 # MSE(1) 35 5.00000000E+02 # MSE(2) 36 1.00000000E+03 # MSE(3) 41 1.50000000E+03 # MSQ(1) 42 1.50000000E+03 # MSQ(2) 43 1.00000000E+03 # MSQ(3) 44 1.50000000E+03 # MSU(1) 45 1.50000000E+03 # MSU(2) 46 1.00000000E+03 # MSU(3) 47 1.50000000E+03 # MSD(1) 48 1.50000000E+03 # MSD(2) 49 1.00000000E+03 # MSD(3) BLOCK MASS 1000012 4.95834286E+02 # MSf(1,1,1) 1000011 5.01843810E+02 # MSf(1,2,1) 2000011 5.02295880E+02 # MSf(2,2,1) 1000002 1.49902739E+03 # MSf(1,3,1) 2000002 1.49958945E+03 # MSf(2,3,1) 1000001 1.50020509E+03 # MSf(1,4,1) 2000001 1.50117722E+03 # MSf(2,4,1) 1000014 4.95834286E+02 # MSf(1,1,2) 1000013 5.01398885E+02 # MSf(1,2,2) 2000013 5.02740033E+02 # MSf(2,2,2) 1000004 1.49902793E+03 # MSf(1,3,2) 2000004 1.49959001E+03 # MSf(2,3,2) 1000003 1.50017147E+03 # MSf(1,4,2) 2000003 1.50121083E+03 # MSf(2,4,2) 1000016 9.97923664E+02 # MSf(1,1,3) 1000015 9.97040283E+02 # MSf(1,2,3) 2000015 1.00502007E+03 # MSf(2,2,3) 1000006 8.76886935E+02 # MSf(1,3,3) 2000006 1.13421049E+03 # MSf(2,3,3) 1000005 9.92457871E+02 # MSf(1,4,3) 2000005 1.00955661E+03 # MSf(2,4,3) 25 1.26034276E+02 # Mh0 35 1.79937794E+03 # MHH 36 1.80000000E+03 # MA0 37 1.80140236E+03 # MHp 1000022 8.83213889E+01 # MNeu(1) 1000023 1.52049632E+02 # MNeu(2) 1000025 -2.10470056E+02 # MNeu(3) 1000035 2.65570688E+02 # MNeu(4) 1000024 1.48692309E+02 # MCha(1) 1000037 2.66117031E+02 # MCha(2) 1000021 1.50000000E+03 # MGl BLOCK DMASS 0 1.72500000E+02 # Q 25 1.37334112E+00 # Delta Mh0 35 2.54539060E-02 # Delta MHH 36 0.00000000E+00 # Delta MA0 37 1.33583318E-02 # Delta MHp BLOCK NMIX 1 1 9.33161766E-01 # ZNeu(1,1) 1 2 -1.11637757E-01 # ZNeu(1,2) 1 3 3.09461555E-01 # ZNeu(1,3) 1 4 -1.44843625E-01 # ZNeu(1,4) 2 1 -3.14536695E-01 # ZNeu(2,1) 2 2 -6.93472054E-01 # ZNeu(2,2) 2 3 5.12603053E-01 # ZNeu(2,3) 2 4 -3.96738310E-01 # ZNeu(2,4) 3 1 9.74524166E-02 # ZNeu(3,1) 3 2 -1.35722538E-01 # ZNeu(3,2) 3 3 -6.77905232E-01 # ZNeu(3,3) 3 4 -7.15909852E-01 # ZNeu(3,4) 4 1 -1.44148578E-01 # ZNeu(4,1) 4 2 6.98722344E-01 # ZNeu(4,2) 4 3 4.26516298E-01 # ZNeu(4,3) 4 4 -5.55960539E-01 # ZNeu(4,4) BLOCK UMIX 1 1 -6.06292879E-01 # UCha(1,1) 1 2 7.95241438E-01 # UCha(1,2) 2 1 7.95241438E-01 # UCha(2,1) 2 2 6.06292879E-01 # UCha(2,2) BLOCK VMIX 1 1 -7.95241438E-01 # VCha(1,1) 1 2 6.06292879E-01 # VCha(1,2) 2 1 6.06292879E-01 # VCha(2,1) 2 2 7.95241438E-01 # VCha(2,2) BLOCK STAUMIX 1 1 6.96989496E-01 # USf(1,1) 1 2 7.17081336E-01 # USf(1,2) 2 1 7.17081336E-01 # USf(2,1) 2 2 -6.96989496E-01 # USf(2,2) BLOCK STOPMIX 1 1 7.08257287E-01 # USf(1,1) 1 2 -7.05954401E-01 # USf(1,2) 2 1 7.05954401E-01 # USf(2,1) 2 2 7.08257287E-01 # USf(2,2) BLOCK SBOTMIX 1 1 6.76310146E-01 # USf(1,1) 1 2 7.36616988E-01 # USf(1,2) 2 1 7.36616988E-01 # USf(2,1) 2 2 -6.76310146E-01 # USf(2,2) BLOCK ALPHA -3.36184710E-02 # Alpha BLOCK DALPHA 1.07290538E-05 # Delta Alpha BLOCK HMIX Q= -0.99900000E+03 1 2.00000000E+02 # MUE 2 3.00000000E+01 # TB BLOCK MSOFT Q= 0.00000000E+00 1 9.54716519E+01 # M1 2 2.00000000E+02 # M2 3 1.50000000E+03 # M3 31 5.00000000E+02 # MSL(1) 32 5.00000000E+02 # MSL(2) 33 1.00000000E+03 # MSL(3) 34 5.00000000E+02 # MSE(1) 35 5.00000000E+02 # MSE(2) 36 1.00000000E+03 # MSE(3) 41 1.50000000E+03 # MSQ(1) 42 1.50000000E+03 # MSQ(2) 43 1.00000000E+03 # MSQ(3) 44 1.50000000E+03 # MSU(1) 45 1.50000000E+03 # MSU(2) 46 1.00000000E+03 # MSU(3) 47 1.50000000E+03 # MSD(1) 48 1.50000000E+03 # MSD(2) 49 1.00000000E+03 # MSD(3) BLOCK AE Q= 0.00000000E+00 1 1 0.00000000E+00 # Af(1,1) 2 2 0.00000000E+00 # Af(2,2) 3 3 1.50666667E+03 # Af(3,3) BLOCK AU Q= 0.00000000E+00 1 1 0.00000000E+00 # Af(1,1) 2 2 0.00000000E+00 # Af(2,2) 3 3 1.50666667E+03 # Af(3,3) BLOCK AD Q= 0.00000000E+00 1 1 0.00000000E+00 # Af(1,1) 2 2 0.00000000E+00 # Af(2,2) 3 3 1.50666667E+03 # Af(3,3) BLOCK YE Q= 0.00000000E+00 1 1 8.80994160E-05 # Yf(1,1) 2 2 1.82161635E-02 # Yf(2,2) 3 3 3.06371119E-01 # Yf(3,3) BLOCK YU Q= 0.00000000E+00 1 1 1.72406273E-05 # Yf(1,1) 2 2 7.39048222E-03 # Yf(2,2) 3 3 9.91336068E-01 # Yf(3,3) BLOCK YD Q= 0.00000000E+00 1 1 1.00363958E-03 # Yf(1,1) 2 2 1.58899417E-02 # Yf(2,2) 3 3 6.54339723E-01 # Yf(3,3) BLOCK VCKMIN 1 2.25300000E-01 # lambda 2 8.08000000E-01 # A 3 1.32000000E-01 # rhobar 4 3.41000000E-01 # etabar BLOCK MSL2 Q= 0.00000000E+00 1 1 2.50000000E+05 # MSL2(1,1) 2 2 2.50000000E+05 # MSL2(2,2) 3 3 1.00000000E+06 # MSL2(3,3) BLOCK MSE2 Q= 0.00000000E+00 1 1 2.50000000E+05 # MSE2(1,1) 2 2 2.50000000E+05 # MSE2(2,2) 3 3 1.00000000E+06 # MSE2(3,3) BLOCK MSQ2 Q= 0.00000000E+00 1 1 2.25000000E+06 # MSQ2(1,1) 2 2 2.25000000E+06 # MSQ2(2,2) 3 3 1.00000000E+06 # MSQ2(3,3) BLOCK MSU2 Q= 0.00000000E+00 1 1 2.25000000E+06 # MSU2(1,1) 2 2 2.25000000E+06 # MSU2(2,2) 3 3 1.00000000E+06 # MSU2(3,3) BLOCK MSD2 Q= 0.00000000E+00 1 1 2.25000000E+06 # MSD2(1,1) 2 2 2.25000000E+06 # MSD2(2,2) 3 3 1.00000000E+06 # MSD2(3,3) BLOCK TE Q= 0.00000000E+00 1 1 0.00000000E+00 # Tf(1,1) 2 2 0.00000000E+00 # Tf(2,2) 3 3 4.61599152E+02 # Tf(3,3) BLOCK TU Q= 0.00000000E+00 1 1 0.00000000E+00 # Tf(1,1) 2 2 0.00000000E+00 # Tf(2,2) 3 3 1.49361301E+03 # Tf(3,3) BLOCK TD Q= 0.00000000E+00 1 1 0.00000000E+00 # Tf(1,1) 2 2 0.00000000E+00 # Tf(2,2) 3 3 9.85871849E+02 # Tf(3,3) BLOCK CVHMIX 1 1 9.99999992E-01 # UH(1,1) 1 2 1.26515337E-04 # UH(1,2) 1 3 0.00000000E+00 # UH(1,3) 2 1 -1.26515337E-04 # UH(2,1) 2 2 9.99999992E-01 # UH(2,2) 2 3 0.00000000E+00 # UH(2,3) 3 1 0.00000000E+00 # UH(3,1) 3 2 0.00000000E+00 # UH(3,2) 3 3 1.00000000E+00 # UH(3,3) BLOCK PRECOBS 1 4.17666251E-05 # DeltaRho 2 8.03586103E+01 # MWMSSM 3 8.03562586E+01 # MWSM 4 2.31530215E-01 # SW2effMSSM 5 2.31543247E-01 # SW2effSM 11 3.36059481E-09 # gminus2mu 21 0.00000000E+00 # EDMeTh 22 0.00000000E+00 # EDMn 23 0.00000000E+00 # EDMHg 31 4.95603571E-04 # bsgammaMSSM 32 3.98954070E-04 # bsgammaSM 33 2.15529640E+01 # DeltaMsMSSM 34 2.12958950E+01 # DeltaMsSM 35 4.01186251E-08 # BsmumuMSSM 36 3.50847021E-09 # BsmumuSM DECAY 25 4.00576356E-03 # Gamma(h0) 2.37905735E-03 2 22 22 # BR(h0 -> photon photon) 1.61122214E-03 2 22 23 # BR(h0 -> photon Z) 3.04464978E-02 2 23 23 # BR(h0 -> Z Z) 2.46470530E-01 2 -24 24 # BR(h0 -> W W) 7.24243403E-02 2 21 21 # BR(h0 -> gluon gluon) 5.10701398E-09 2 -11 11 # BR(h0 -> Electron electron) 2.27169573E-04 2 -13 13 # BR(h0 -> Muon muon) 6.52922542E-02 2 -15 15 # BR(h0 -> Tau tau) 2.02575352E-07 2 -2 2 # BR(h0 -> Up up) 2.80593958E-02 2 -4 4 # BR(h0 -> Charm charm) 8.36383821E-07 2 -1 1 # BR(h0 -> Down down) 2.10043405E-04 2 -3 3 # BR(h0 -> Strange strange) 5.52878445E-01 2 -5 5 # BR(h0 -> Bottom bottom) DECAY 35 4.11780753E+01 # Gamma(HH) 1.59871837E-07 2 22 22 # BR(HH -> photon photon) 3.98501533E-07 2 22 23 # BR(HH -> photon Z) 3.02693637E-06 2 23 23 # BR(HH -> Z Z) 1.86779843E-06 2 -24 24 # BR(HH -> W W) 1.34397398E-06 2 21 21 # BR(HH -> gluon gluon) 6.16660505E-09 2 -11 11 # BR(HH -> Electron electron) 2.74522962E-04 2 -13 13 # BR(HH -> Muon muon) 7.74685751E-02 2 -15 15 # BR(HH -> Tau tau) 1.71035551E-13 2 -2 2 # BR(HH -> Up up) 2.36777321E-08 2 -4 4 # BR(HH -> Charm charm) 1.35480570E-03 2 -6 6 # BR(HH -> Top top) 6.16083771E-07 2 -1 1 # BR(HH -> Down down) 1.54714491E-04 2 -3 3 # BR(HH -> Strange strange) 3.55868040E-01 2 -5 5 # BR(HH -> Bottom bottom) 1.31905624E-01 2 -1000024 1000024 # BR(HH -> Chargino1 chargino1) 8.58839688E-02 2 -1000037 1000024 # BR(HH -> Chargino2 chargino1) 8.58839688E-02 2 -1000024 1000037 # BR(HH -> Chargino1 chargino2) 3.73915555E-02 2 -1000037 1000037 # BR(HH -> Chargino2 chargino2) 1.22102946E-02 2 1000022 1000022 # BR(HH -> neutralino1 neutralino1) 3.76159462E-02 2 1000022 1000023 # BR(HH -> neutralino1 neutralino2) 2.41975173E-02 2 1000022 1000025 # BR(HH -> neutralino1 neutralino3) 2.77989462E-05 2 1000022 1000035 # BR(HH -> neutralino1 neutralino4) 2.35386566E-02 2 1000023 1000023 # BR(HH -> neutralino2 neutralino2) 1.32098313E-02 2 1000023 1000025 # BR(HH -> neutralino2 neutralino3) 4.92692320E-03 2 1000023 1000035 # BR(HH -> neutralino2 neutralino4) 5.71730426E-03 2 1000025 1000025 # BR(HH -> neutralino3 neutralino3) 6.64679026E-02 2 1000025 1000035 # BR(HH -> neutralino3 neutralino4) 3.12304174E-02 2 1000035 1000035 # BR(HH -> neutralino4 neutralino4) 2.14935120E-05 2 25 25 # BR(HH -> h0 h0) 2.17756432E-07 2 -1000011 1000011 # BR(HH -> Selectron1 selectron1) 3.36555206E-13 2 -1000011 2000011 # BR(HH -> Selectron1 selectron2) 3.36555206E-13 2 -2000011 1000011 # BR(HH -> Selectron2 selectron1) 3.37773448E-07 2 -2000011 2000011 # BR(HH -> Selectron2 selectron2) 2.16839982E-07 2 -1000013 1000013 # BR(HH -> Smuon1 smuon1) 1.63486596E-09 2 -1000013 2000013 # BR(HH -> Smuon1 smuon2) 1.63486596E-09 2 -2000013 1000013 # BR(HH -> Smuon2 smuon1) 3.35958314E-07 2 -2000013 2000013 # BR(HH -> Smuon2 smuon2) 4.63937548E-03 2 -1000006 1000006 # BR(HH -> Stop1 stop1) DECAY 36 4.04818085E+01 # Gamma(A0) -2.24381563E-07 2 22 22 # BR(A0 -> photon photon) -4.51000258E-07 2 22 23 # BR(A0 -> photon Z) -2.04162948E-05 2 21 21 # BR(A0 -> gluon gluon) -6.08894146E-09 2 -11 11 # BR(A0 -> Electron electron) -2.71065595E-04 2 -13 13 # BR(A0 -> Muon muon) 7.64931726E-02 2 -15 15 # BR(A0 -> Tau tau) -1.73861937E-13 2 -2 2 # BR(A0 -> Up up) -2.40693245E-08 2 -4 4 # BR(A0 -> Charm charm) -1.55798844E-03 2 -6 6 # BR(A0 -> Top top) -6.08295221E-07 2 -1 1 # BR(A0 -> Down down) -1.52758692E-04 2 -3 3 # BR(A0 -> Strange strange) -3.51445787E-01 2 -5 5 # BR(A0 -> Bottom bottom) -1.49167991E-01 2 -1000024 1000024 # BR(A0 -> Chargino1 chargino1) -7.26880042E-02 2 -1000037 1000024 # BR(A0 -> Chargino2 chargino1) -7.26880042E-02 2 -1000024 1000037 # BR(A0 -> Chargino1 chargino2) -5.25193379E-02 2 -1000037 1000037 # BR(A0 -> Chargino2 chargino2) -1.33575327E-02 2 1000022 1000022 # BR(A0 -> neutralino1 neutralino1) -4.25759923E-02 2 1000022 1000023 # BR(A0 -> neutralino1 neutralino2) -2.01966186E-02 2 1000022 1000025 # BR(A0 -> neutralino1 neutralino3) -1.38808325E-04 2 1000022 1000035 # BR(A0 -> neutralino1 neutralino4) -2.73507116E-02 2 1000023 1000023 # BR(A0 -> neutralino2 neutralino2) -1.02318929E-02 2 1000023 1000025 # BR(A0 -> neutralino2 neutralino3) -5.34289221E-03 2 1000023 1000035 # BR(A0 -> neutralino2 neutralino4) -5.34276333E-03 2 1000025 1000025 # BR(A0 -> neutralino3 neutralino3) -5.69928978E-02 2 1000025 1000035 # BR(A0 -> neutralino3 neutralino4) -4.14599572E-02 2 1000035 1000035 # BR(A0 -> neutralino4 neutralino4) -4.09015046E-06 2 23 25 # BR(A0 -> Z h0) -9.85665589E-38 2 25 25 # BR(A0 -> h0 h0) -3.77010562E-14 2 -1000011 2000011 # BR(A0 -> Selectron1 selectron2) -3.77010562E-14 2 -2000011 1000011 # BR(A0 -> Selectron2 selectron1) -1.61183779E-09 2 -1000013 2000013 # BR(A0 -> Smuon1 smuon2) -1.61183779E-09 2 -2000013 1000013 # BR(A0 -> Smuon2 smuon1) DECAY 37 4.14366100E+01 # Gamma(Hp) 6.47558820E-09 2 -11 12 # BR(Hp -> Electron nu_e) 2.76851571E-04 2 -13 14 # BR(Hp -> Muon nu_mu) 7.83119365E-02 2 -15 16 # BR(Hp -> Tau nu_tau) 5.89157662E-07 2 -1 2 # BR(Hp -> Down up) 6.86020640E-06 2 -3 2 # BR(Hp -> Strange up) 3.86760818E-06 2 -5 2 # BR(Hp -> Bottom up) 2.85549295E-08 2 -1 4 # BR(Hp -> Down charm) 1.47455532E-04 2 -3 4 # BR(Hp -> Strange charm) 5.41603203E-04 2 -5 4 # BR(Hp -> Bottom charm) 1.33727573E-07 2 -1 6 # BR(Hp -> Down top) 3.12976011E-06 2 -3 6 # BR(Hp -> Strange top) 3.65297970E-01 2 -5 6 # BR(Hp -> Bottom top) 5.75347532E-02 2 1000022 1000024 # BR(Hp -> neutralino1 chargino1) 2.16304621E-03 2 1000022 1000037 # BR(Hp -> neutralino1 chargino2) 1.04565650E-02 2 1000023 1000024 # BR(Hp -> neutralino2 chargino1) 2.03074496E-01 2 1000023 1000037 # BR(Hp -> neutralino2 chargino2) 7.18193727E-02 2 1000024 1000025 # BR(Hp -> chargino1 neutralino3) 8.40655159E-02 2 1000025 1000037 # BR(Hp -> neutralino3 chargino2) 1.24571218E-01 2 1000024 1000035 # BR(Hp -> chargino1 neutralino4) 1.71796755E-03 2 1000035 1000037 # BR(Hp -> neutralino4 chargino2) 4.01237112E-06 2 24 25 # BR(Hp -> W h0) 1.11559246E-12 2 24 35 # BR(Hp -> W HH) 1.78010125E-13 2 24 36 # BR(Hp -> W A0) 5.56338939E-11 2 -1000011 1000012 # BR(Hp -> Selectron1 snu_e1) 1.30963838E-06 2 -2000011 1000012 # BR(Hp -> Selectron2 snu_e1) 3.75230861E-07 2 -1000013 1000014 # BR(Hp -> Smuon1 snu_mu1) 9.35546835E-07 2 -2000013 1000014 # BR(Hp -> Smuon2 snu_mu1) DECAY 6 1.35195755E+00 # Gamma(top) 1.00000000E+00 2 5 24 # BR(top -> bottom W) """ import FWCore.ParameterSet.Config as cms from Configuration.Generator.Pythia8CommonSettings_cfi import * from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import * generator = cms.EDFilter("Pythia8GeneratorFilter", pythiaPylistVerbosity = cms.untracked.int32(1), filterEfficiency = cms.untracked.double(1.0), pythiaHepMCVerbosity = cms.untracked.bool(False), comEnergy = cms.double(13000.0), crossSection = cms.untracked.double(518.3), maxEventsToPrint = cms.untracked.int32(1), SLHATableForPythia8 = cms.string('%s' % SLHA_TABLE), PythiaParameters = cms.PSet( pythia8CommonSettingsBlock, pythia8CUEP8M1SettingsBlock, processParameters = cms.vstring('Higgs:useBSM = on', 'HiggsBSM:gg2A3 = on', 'SLHA:allowUserOverride = off', 'SLHA:minMassSM = 100.', pscut_min, # pscut_max # '36:m0 = 90', # '36:onMode = off', # '36:onIfAny = 15' ), parameterSets = cms.vstring('pythia8CommonSettings', 'pythia8CUEP8M1Settings', 'processParameters') ) )