import FWCore.ParameterSet.Config as cms from Configuration.Generator.Pythia8CommonSettings_cfi import * from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import * from Configuration.Generator.Pythia8aMCatNLOSettings_cfi import * baseSLHATable=""" BLOCK SPINFO 1 FeynHiggs 2 2.10.0 2 built on Mar 10, 2014 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 4.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.50500000E+03 # At 12 1.50500000E+03 # Ab 13 1.50500000E+03 # Atau 23 2.00000000E+02 # MUE 25 4.00000000E+01 # TB 26 3.00000000E+01 # MA0 27 8.58128103E+01 # 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.95830217E+02 # MSf(1,1,1) 1000011 5.01845584E+02 # MSf(1,2,1) 2000011 5.02298124E+02 # MSf(2,2,1) 1000002 1.49902645E+03 # MSf(1,3,1) 2000002 1.49958905E+03 # MSf(2,3,1) 1000001 1.50020519E+03 # MSf(1,4,1) 2000001 1.50117846E+03 # MSf(2,4,1) 1000014 4.95830217E+02 # MSf(1,1,2) 1000013 5.01199509E+02 # MSf(1,2,2) 2000013 5.02942810E+02 # MSf(2,2,2) 1000004 1.49902699E+03 # MSf(1,3,2) 2000004 1.49958961E+03 # MSf(2,3,2) 1000003 1.50014799E+03 # MSf(1,4,2) 2000003 1.50123564E+03 # MSf(2,4,2) 1000016 9.97921642E+02 # MSf(1,1,3) 1000015 9.95256404E+02 # MSf(1,2,3) 2000015 1.00678866E+03 # MSf(2,2,3) 1000006 8.76885784E+02 # MSf(1,3,3) 2000006 1.13420960E+03 # MSf(2,3,3) 1000005 9.88986250E+02 # MSf(1,4,3) 2000005 1.01295893E+03 # MSf(2,4,3) 25 1.25000000E+02 # Mh0 35 1.26348865E+02 # MHH 36 %MA0% # MA0 37 9.47952678E+01 # MHp 1000022 8.85696665E+01 # MNeu(1) 1000023 1.52403239E+02 # MNeu(2) 1000025 -2.10646406E+02 # MNeu(3) 1000035 2.65145152E+02 # MNeu(4) 1000024 1.49281133E+02 # MCha(1) 1000037 2.65787171E+02 # MCha(2) 1000021 1.50000000E+03 # MGl BLOCK DMASS 0 1.72500000E+02 # Q 25 0.00000000E+00 # Delta Mh0 35 1.35783674E+00 # Delta MHH 36 0.00000000E+00 # Delta MA0 37 6.04629261E-01 # Delta MHp BLOCK NMIX 1 1 9.34616473E-01 # ZNeu(1,1) 1 2 -1.08169841E-01 # ZNeu(1,2) 1 3 3.07485524E-01 # ZNeu(1,3) 1 4 -1.42281363E-01 # ZNeu(1,4) 2 1 -3.10165758E-01 # ZNeu(2,1) 2 2 -6.93221305E-01 # ZNeu(2,2) 2 3 5.14891896E-01 # ZNeu(2,3) 2 4 -3.97652814E-01 # ZNeu(2,4) 3 1 9.81909017E-02 # ZNeu(3,1) 3 2 -1.36771096E-01 # ZNeu(3,2) 3 3 -6.77763154E-01 # ZNeu(3,3) 3 4 -7.15743893E-01 # ZNeu(3,4) 4 1 -1.43693415E-01 # ZNeu(4,1) 4 2 6.99311930E-01 # ZNeu(4,2) 4 3 4.25412853E-01 # ZNeu(4,3) 4 4 -5.56182463E-01 # ZNeu(4,4) BLOCK UMIX 1 1 -6.05392853E-01 # UCha(1,1) 1 2 7.95926814E-01 # UCha(1,2) 2 1 7.95926814E-01 # UCha(2,1) 2 2 6.05392853E-01 # UCha(2,2) BLOCK VMIX 1 1 -7.95926814E-01 # VCha(1,1) 1 2 6.05392853E-01 # VCha(1,2) 2 1 6.05392853E-01 # VCha(2,1) 2 2 7.95926814E-01 # VCha(2,2) BLOCK STAUMIX 1 1 7.00114782E-01 # USf(1,1) 1 2 7.14030316E-01 # USf(1,2) 2 1 7.14030316E-01 # USf(2,1) 2 2 -7.00114782E-01 # USf(2,2) BLOCK STOPMIX 1 1 7.08258405E-01 # USf(1,1) 1 2 -7.05953279E-01 # USf(1,2) 2 1 7.05953279E-01 # USf(2,1) 2 2 7.08258405E-01 # USf(2,2) BLOCK SBOTMIX 1 1 6.85260003E-01 # USf(1,1) 1 2 7.28298516E-01 # USf(1,2) 2 1 7.28298516E-01 # USf(2,1) 2 2 -6.85260003E-01 # USf(2,2) BLOCK ALPHA -1.54441184E+00 # Alpha BLOCK DALPHA 9.83428934E-04 # Delta Alpha BLOCK HMIX Q= -0.99900000E+03 1 2.00000000E+02 # MUE 2 4.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.50500000E+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.50500000E+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.50500000E+03 # Af(3,3) BLOCK YE Q= 0.00000000E+00 1 1 1.17437365E-04 # Yf(1,1) 2 2 2.42823205E-02 # Yf(2,2) 3 3 4.08395636E-01 # Yf(3,3) BLOCK YU Q= 0.00000000E+00 1 1 1.72364410E-05 # Yf(1,1) 2 2 7.38868770E-03 # Yf(2,2) 3 3 9.91095356E-01 # Yf(3,3) BLOCK YD Q= 0.00000000E+00 1 1 1.32471393E-03 # Yf(1,1) 2 2 2.09728535E-02 # Yf(2,2) 3 3 8.47502026E-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 6.14635432E+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.49159851E+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 1.27549055E+03 # Tf(3,3) BLOCK CVHMIX 1 1 9.99989125E-01 # UH(1,1) 1 2 4.66359301E-03 # UH(1,2) 1 3 0.00000000E+00 # UH(1,3) 2 1 -4.66359301E-03 # UH(2,1) 2 2 9.99989125E-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.17304250E-05 # DeltaRho 2 8.03584639E+01 # MWMSSM 3 8.03561143E+01 # MWSM 4 2.31531417E-01 # SW2effMSSM 5 2.31544437E-01 # SW2effSM 11 5.58086942E-09 # gminus2mu 21 0.00000000E+00 # EDMeTh 22 0.00000000E+00 # EDMn 23 0.00000000E+00 # EDMHg 31 8.04820993E-04 # bsgammaMSSM 32 3.98981472E-04 # bsgammaSM 33 2.16046139E+01 # DeltaMsMSSM 34 2.13349597E+01 # DeltaMsSM 35 4.00124270E-08 # BsmumuMSSM 36 3.51490606E-09 # BsmumuSM DECAY 25 4.00000000E-03 # Gamma(h0) 1.00000000E+00 2 36 36 # BR(h0 -> A0 A0) DECAY 35 1.39295461E-01 # Gamma(HH) 6.41895580E-05 2 22 22 # BR(HH -> photon photon) 1.97232862E-04 2 22 23 # BR(HH -> photon Z) 9.17994878E-04 2 23 23 # BR(HH -> Z Z) 7.38657150E-03 2 -24 24 # BR(HH -> W W) 2.55885394E-03 2 21 21 # BR(HH -> gluon gluon) 1.47272074E-10 2 -11 11 # BR(HH -> Electron electron) 6.55087625E-06 2 -13 13 # BR(HH -> Muon muon) 1.74659059E-03 2 -15 15 # BR(HH -> Tau tau) 5.79395931E-09 2 -2 2 # BR(HH -> Up up) 8.01266697E-04 2 -4 4 # BR(HH -> Charm charm) 1.78029768E-08 2 -1 1 # BR(HH -> Down down) 4.46929265E-06 2 -3 3 # BR(HH -> Strange strange) 8.05657718E-03 2 -5 5 # BR(HH -> Bottom bottom) 1.73933172E-04 2 23 36 # BR(HH -> Z A0) 4.82185990E-01 2 25 25 # BR(HH -> h0 h0) 2.12714798E-37 2 25 36 # BR(HH -> h0 A0) 6.60852658E-36 2 25 36 # BR(HH -> h0 A0) 4.95899756E-01 2 36 36 # BR(HH -> A0 A0) DECAY 36 6.00000000E-03 # Gamma(A0) 0.50000000E+00 2 -13 13 # BR(A0 -> muon muon) 0.50000000E+00 2 -5 5 # BR(A0 -> b b_bar) DECAY 37 3.09600270E-01 # Gamma(Hp) 8.10803314E-08 2 -11 12 # BR(Hp -> Electron nu_e) 3.46642819E-03 2 -13 14 # BR(Hp -> Muon nu_mu) 9.79850491E-01 2 -15 16 # BR(Hp -> Tau nu_tau) 1.14219676E-05 2 -1 2 # BR(Hp -> Down up) 1.26170695E-04 2 -3 2 # BR(Hp -> Strange up) 6.82337775E-05 2 -5 2 # BR(Hp -> Bottom up) 5.09808557E-07 2 -1 4 # BR(Hp -> Down charm) 2.85719255E-03 2 -3 4 # BR(Hp -> Strange charm) 9.55163990E-03 2 -5 4 # BR(Hp -> Bottom charm) 9.78756246E-10 2 -5 6 # BR(Hp -> Bottom top) 2.00794797E-03 2 24 25 # BR(Hp -> W h0) 2.05988251E-03 2 24 36 # BR(Hp -> W A0) DECAY 6 1.59404260E+00 # Gamma(top) 8.48131384E-01 2 5 24 # BR(top -> bottom W) 1.51868616E-01 2 5 37 # BR(top -> bottom Hp) """ ma0 = 17 slhatable = baseSLHATable.replace(' %MA0%','%e' % ma0) generator = cms.EDFilter("Pythia8ConcurrentHadronizerFilter", eventsToPrint = cms.untracked.uint32(5), UseExternalGenerators = cms.untracked.bool(True), maxEventsToPrint = cms.untracked.int32(5), pythiaPylistVerbosity = cms.untracked.int32(1), filterEfficiency = cms.untracked.double(1.0), pythiaHepMCVerbosity = cms.untracked.bool(False), comEnergy = cms.double(13000.), ## set up mA mass and force h125->A0A0 and A->mumu, bb decays SLHATableForPythia8 = cms.string('%s' % slhatable), PythiaParameters = cms.PSet( pythia8CommonSettingsBlock, pythia8CP5SettingsBlock, pythia8aMCatNLOSettingsBlock, processParameters = cms.vstring( 'Main:timesAllowErrors = 10000', 'ParticleDecays:limitTau0 = on', 'ParticleDecays:tau0Max = 10', 'Check:epTolErr = 0.01', 'SLHA:keepSM = on', 'SLHA:minMassSM = 10.', 'TimeShower:nPartonsInBorn = 0', #number of coloured particles (before resonance decays) in LO ), parameterSets = cms.vstring('pythia8CommonSettings', 'pythia8CP5Settings', 'pythia8aMCatNLOSettings', 'processParameters',) ) ) ProductionFilterSequence = cms.Sequence(generator)