Merge pull request #29 from peiffer/master

tprime and bprime fragments with Pythia8
CmsHI · Nov 13, 2013 · 16d3aea · 16d3aea
2 parents 1fcc1a6 + 27ba5bb
commit 16d3aea
Show file tree

Hide file tree

Showing 6 changed files with 309 additions and 0 deletions.
diff --git a/...eenTeV/BprimeBprimeToBHBZTWinc_M_1000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py b/...eenTeV/BprimeBprimeToBHBZTWinc_M_1000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py
@@ -0,0 +1,52 @@
+import FWCore.ParameterSet.Config as cms
+from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
+
+generator = cms.EDFilter("Pythia8HadronizerFilter",
+    maxEventsToPrint = cms.untracked.int32(1),
+    pythiaPylistVerbosity = cms.untracked.int32(1),
+    filterEfficiency = cms.untracked.double(1.0),
+    pythiaHepMCVerbosity = cms.untracked.bool(False),
+    comEnergy = cms.double(13000.),
+    ExternalDecays = cms.PSet(
+        Tauola = cms.untracked.PSet(
+            TauolaPolar,
+            TauolaDefaultInputCards
+        ),
+        parameterSets = cms.vstring('Tauola')
+    ),
+    UseExternalGenerators = cms.untracked.bool(True),
+    jetMatching = cms.untracked.PSet(
+        scheme = cms.string("Madgraph"),
+        mode = cms.string("auto"),# soup, or "inclusive" / "exclusive"
+        MEMAIN_etaclmax = cms.double(-1),
+        MEMAIN_qcut = cms.double(-1),
+        MEMAIN_minjets = cms.int32(-1),
+        MEMAIN_maxjets = cms.int32(-1),
+        MEMAIN_showerkt = cms.double(0), # use 1=yes only for pt-ordered showers !
+        MEMAIN_nqmatch = cms.int32(5), #PID of the flavor until which the QCD radiation are kept in the matching procedure;
+                                         # if nqmatch=4, then all showered partons from b's are NOT taken into account
+         # Note (JY): I think the default should be 5 (b); anyway, don't try -1 as it'll result in a throw...
+        MEMAIN_excres = cms.string(""),
+        outTree_flag = cms.int32(0) # 1=yes, write out the tree for future sanity check
+    ),
+    PythiaParameters = cms.PSet(
+        processParameters = cms.vstring(
+         'Main:timesAllowErrors = 10000',
+        'ParticleDecays:limitTau0 = on',
+         'ParticleDecays:tauMax = 10',
+        'Tune:ee 3',
+        'Tune:pp 5',
+        '7:doForceWidth = 1',
+        '7:m0 = 1000 ! bprime mass',
+        '7:mWidth = 10 ! bprime width',
+        '7:mMin = 900 ! minimal bprime mass',
+        '7:mMax = 1100 ! maximal bprime mass',
+        '6:m0 = 172.5 ! top mass',
+        '7:oneChannel = 1 0.5 100 24 6 ! add channel bprime to W t with BR=0.5',
+        '7:addChannel = 1 0.25 100 23 5 ! add channel bprime to Z b with BR=0.25',
+        '7:addChannel = 1 0.25 100 25 5 ! add channel bprime to H b with BR=0.25',
+
+         ),
+        parameterSets = cms.vstring('processParameters')
+    )
+)
diff --git a/...eenTeV/BprimeBprimeToBHBZTWinc_M_2000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py b/...eenTeV/BprimeBprimeToBHBZTWinc_M_2000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py
@@ -0,0 +1,52 @@
+import FWCore.ParameterSet.Config as cms
+from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
+
+generator = cms.EDFilter("Pythia8HadronizerFilter",
+    maxEventsToPrint = cms.untracked.int32(1),
+    pythiaPylistVerbosity = cms.untracked.int32(1),
+    filterEfficiency = cms.untracked.double(1.0),
+    pythiaHepMCVerbosity = cms.untracked.bool(False),
+    comEnergy = cms.double(13000.),
+    ExternalDecays = cms.PSet(
+        Tauola = cms.untracked.PSet(
+            TauolaPolar,
+            TauolaDefaultInputCards
+        ),
+        parameterSets = cms.vstring('Tauola')
+    ),
+    UseExternalGenerators = cms.untracked.bool(True),
+    jetMatching = cms.untracked.PSet(
+        scheme = cms.string("Madgraph"),
+        mode = cms.string("auto"),# soup, or "inclusive" / "exclusive"
+        MEMAIN_etaclmax = cms.double(-1),
+        MEMAIN_qcut = cms.double(-1),
+        MEMAIN_minjets = cms.int32(-1),
+        MEMAIN_maxjets = cms.int32(-1),
+        MEMAIN_showerkt = cms.double(0), # use 1=yes only for pt-ordered showers !
+        MEMAIN_nqmatch = cms.int32(5), #PID of the flavor until which the QCD radiation are kept in the matching procedure;
+                                         # if nqmatch=4, then all showered partons from b's are NOT taken into account
+         # Note (JY): I think the default should be 5 (b); anyway, don't try -1 as it'll result in a throw...
+        MEMAIN_excres = cms.string(""),
+        outTree_flag = cms.int32(0) # 1=yes, write out the tree for future sanity check
+    ),
+    PythiaParameters = cms.PSet(
+        processParameters = cms.vstring(
+         'Main:timesAllowErrors = 10000',
+        'ParticleDecays:limitTau0 = on',
+         'ParticleDecays:tauMax = 10',
+        'Tune:ee 3',
+        'Tune:pp 5',
+        '7:doForceWidth = 1',
+        '7:m0 = 2000 ! bprime mass',
+        '7:mWidth = 20 ! bprime width',
+        '7:mMin = 1900 ! minimal bprime mass',
+        '7:mMax = 2100 ! maximal bprime mass',
+        '6:m0 = 172.5 ! top mass',
+        '7:oneChannel = 1 0.5 100 24 6 ! add channel bprime to W t with BR=0.5',
+        '7:addChannel = 1 0.25 100 23 5 ! add channel bprime to Z b with BR=0.25',
+        '7:addChannel = 1 0.25 100 25 5 ! add channel bprime to H b with BR=0.25',
+
+         ),
+        parameterSets = cms.vstring('processParameters')
+    )
+)
diff --git a/...eenTeV/BprimeBprimeToBHBZTWinc_M_3000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py b/...eenTeV/BprimeBprimeToBHBZTWinc_M_3000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py
@@ -0,0 +1,52 @@
+import FWCore.ParameterSet.Config as cms
+from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
+
+generator = cms.EDFilter("Pythia8HadronizerFilter",
+    maxEventsToPrint = cms.untracked.int32(1),
+    pythiaPylistVerbosity = cms.untracked.int32(1),
+    filterEfficiency = cms.untracked.double(1.0),
+    pythiaHepMCVerbosity = cms.untracked.bool(False),
+    comEnergy = cms.double(13000.),
+    ExternalDecays = cms.PSet(
+        Tauola = cms.untracked.PSet(
+            TauolaPolar,
+            TauolaDefaultInputCards
+        ),
+        parameterSets = cms.vstring('Tauola')
+    ),
+    UseExternalGenerators = cms.untracked.bool(True),
+    jetMatching = cms.untracked.PSet(
+        scheme = cms.string("Madgraph"),
+        mode = cms.string("auto"),# soup, or "inclusive" / "exclusive"
+        MEMAIN_etaclmax = cms.double(-1),
+        MEMAIN_qcut = cms.double(-1),
+        MEMAIN_minjets = cms.int32(-1),
+        MEMAIN_maxjets = cms.int32(-1),
+        MEMAIN_showerkt = cms.double(0), # use 1=yes only for pt-ordered showers !
+        MEMAIN_nqmatch = cms.int32(5), #PID of the flavor until which the QCD radiation are kept in the matching procedure;
+                                         # if nqmatch=4, then all showered partons from b's are NOT taken into account
+         # Note (JY): I think the default should be 5 (b); anyway, don't try -1 as it'll result in a throw...
+        MEMAIN_excres = cms.string(""),
+        outTree_flag = cms.int32(0) # 1=yes, write out the tree for future sanity check
+    ),
+    PythiaParameters = cms.PSet(
+        processParameters = cms.vstring(
+         'Main:timesAllowErrors = 10000',
+        'ParticleDecays:limitTau0 = on',
+         'ParticleDecays:tauMax = 10',
+        'Tune:ee 3',
+        'Tune:pp 5',
+        '7:doForceWidth = 1',
+        '7:m0 = 3000 ! bprime mass',
+        '7:mWidth = 30 ! bprime width',
+        '7:mMin = 2900 ! minimal bprime mass',
+        '7:mMax = 3100 ! maximal bprime mass',
+        '6:m0 = 172.5 ! top mass',
+        '7:oneChannel = 1 0.5 100 24 6 ! add channel bprime to W t with BR=0.5',
+        '7:addChannel = 1 0.25 100 23 5 ! add channel bprime to Z b with BR=0.25',
+        '7:addChannel = 1 0.25 100 25 5 ! add channel bprime to H b with BR=0.25',
+
+         ),
+        parameterSets = cms.vstring('processParameters')
+    )
+)
diff --git a/...eenTeV/TprimeTprimeToTHTZBWinc_M_1000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py b/...eenTeV/TprimeTprimeToTHTZBWinc_M_1000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py
@@ -0,0 +1,51 @@
+import FWCore.ParameterSet.Config as cms
+from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
+
+generator = cms.EDFilter("Pythia8HadronizerFilter",
+    maxEventsToPrint = cms.untracked.int32(1),
+    pythiaPylistVerbosity = cms.untracked.int32(1),
+    filterEfficiency = cms.untracked.double(1.0),
+    pythiaHepMCVerbosity = cms.untracked.bool(False),
+    comEnergy = cms.double(13000.),
+    ExternalDecays = cms.PSet(
+        Tauola = cms.untracked.PSet(
+            TauolaPolar,
+            TauolaDefaultInputCards
+        ),
+        parameterSets = cms.vstring('Tauola')
+    ),
+    UseExternalGenerators = cms.untracked.bool(True),
+    jetMatching = cms.untracked.PSet(
+        scheme = cms.string("Madgraph"),
+        mode = cms.string("auto"),# soup, or "inclusive" / "exclusive"
+        MEMAIN_etaclmax = cms.double(-1),
+        MEMAIN_qcut = cms.double(-1),
+        MEMAIN_minjets = cms.int32(-1),
+        MEMAIN_maxjets = cms.int32(-1),
+        MEMAIN_showerkt = cms.double(0), # use 1=yes only for pt-ordered showers !
+        MEMAIN_nqmatch = cms.int32(5), #PID of the flavor until which the QCD radiation are kept in the matching procedure;
+                                         # if nqmatch=4, then all showered partons from b's are NOT taken into account
+         # Note (JY): I think the default should be 5 (b); anyway, don't try -1 as it'll result in a throw...
+        MEMAIN_excres = cms.string(""),
+        outTree_flag = cms.int32(0) # 1=yes, write out the tree for future sanity check
+    ),
+    PythiaParameters = cms.PSet(
+        processParameters = cms.vstring(
+         'Main:timesAllowErrors = 10000',
+        'ParticleDecays:limitTau0 = on',
+         'ParticleDecays:tauMax = 10',
+        'Tune:ee 3',
+        'Tune:pp 5',
+        '8:doForceWidth = 1',
+        '8:m0 = 1000 ! tprime mass',
+        '8:mWidth = 10 ! tprime width',
+        '8:mMin = 900 ! minimal tprime mass',
+        '8:mMax = 1100 ! maximal tprime mass',
+        '6:m0 = 172.5 ! top mass',
+        '8:oneChannel = 1 0.5 100 24 5 ! add channel tprime to W b with BR=0.5',
+        '8:addChannel = 1 0.25 100 23 6 ! add channel tprime to Z t with BR=0.25',
+        '8:addChannel = 1 0.25 100 25 6 ! add channel tprime to H t with BR=0.25',
+         ),
+        parameterSets = cms.vstring('processParameters')
+    )
+)
diff --git a/...eenTeV/TprimeTprimeToTHTZBWinc_M_2000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py b/...eenTeV/TprimeTprimeToTHTZBWinc_M_2000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py
@@ -0,0 +1,51 @@
+import FWCore.ParameterSet.Config as cms
+from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
+
+generator = cms.EDFilter("Pythia8HadronizerFilter",
+    maxEventsToPrint = cms.untracked.int32(1),
+    pythiaPylistVerbosity = cms.untracked.int32(1),
+    filterEfficiency = cms.untracked.double(1.0),
+    pythiaHepMCVerbosity = cms.untracked.bool(False),
+    comEnergy = cms.double(13000.),
+    ExternalDecays = cms.PSet(
+        Tauola = cms.untracked.PSet(
+            TauolaPolar,
+            TauolaDefaultInputCards
+        ),
+        parameterSets = cms.vstring('Tauola')
+    ),
+    UseExternalGenerators = cms.untracked.bool(True),
+    jetMatching = cms.untracked.PSet(
+        scheme = cms.string("Madgraph"),
+        mode = cms.string("auto"),# soup, or "inclusive" / "exclusive"
+        MEMAIN_etaclmax = cms.double(-1),
+        MEMAIN_qcut = cms.double(-1),
+        MEMAIN_minjets = cms.int32(-1),
+        MEMAIN_maxjets = cms.int32(-1),
+        MEMAIN_showerkt = cms.double(0), # use 1=yes only for pt-ordered showers !
+        MEMAIN_nqmatch = cms.int32(5), #PID of the flavor until which the QCD radiation are kept in the matching procedure;
+                                         # if nqmatch=4, then all showered partons from b's are NOT taken into account
+         # Note (JY): I think the default should be 5 (b); anyway, don't try -1 as it'll result in a throw...
+        MEMAIN_excres = cms.string(""),
+        outTree_flag = cms.int32(0) # 1=yes, write out the tree for future sanity check
+    ),
+    PythiaParameters = cms.PSet(
+        processParameters = cms.vstring(
+         'Main:timesAllowErrors = 10000',
+        'ParticleDecays:limitTau0 = on',
+         'ParticleDecays:tauMax = 10',
+        'Tune:ee 3',
+        'Tune:pp 5',
+        '8:doForceWidth = 1',
+        '8:m0 = 2000 ! tprime mass',
+        '8:mWidth = 20 ! tprime width',
+        '8:mMin = 1900 ! minimal tprime mass',
+        '8:mMax = 2100 ! maximal tprime mass',
+        '6:m0 = 172.5 ! top mass',
+        '8:oneChannel = 1 0.5 100 24 5 ! add channel tprime to W b with BR=0.5',
+        '8:addChannel = 1 0.25 100 23 6 ! add channel tprime to Z t with BR=0.25',
+        '8:addChannel = 1 0.25 100 25 6 ! add channel tprime to H t with BR=0.25',
+         ),
+        parameterSets = cms.vstring('processParameters')
+    )
+)
diff --git a/...eenTeV/TprimeTprimeToTHTZBWinc_M_3000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py b/...eenTeV/TprimeTprimeToTHTZBWinc_M_3000_MgmMatchTune4C_13TeV_madgraph_pythia8_tauola_cff.py
@@ -0,0 +1,51 @@
+import FWCore.ParameterSet.Config as cms
+from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
+
+generator = cms.EDFilter("Pythia8HadronizerFilter",
+    maxEventsToPrint = cms.untracked.int32(1),
+    pythiaPylistVerbosity = cms.untracked.int32(1),
+    filterEfficiency = cms.untracked.double(1.0),
+    pythiaHepMCVerbosity = cms.untracked.bool(False),
+    comEnergy = cms.double(13000.),
+    ExternalDecays = cms.PSet(
+        Tauola = cms.untracked.PSet(
+            TauolaPolar,
+            TauolaDefaultInputCards
+        ),
+        parameterSets = cms.vstring('Tauola')
+    ),
+    UseExternalGenerators = cms.untracked.bool(True),
+    jetMatching = cms.untracked.PSet(
+        scheme = cms.string("Madgraph"),
+        mode = cms.string("auto"),# soup, or "inclusive" / "exclusive"
+        MEMAIN_etaclmax = cms.double(-1),
+        MEMAIN_qcut = cms.double(-1),
+        MEMAIN_minjets = cms.int32(-1),
+        MEMAIN_maxjets = cms.int32(-1),
+        MEMAIN_showerkt = cms.double(0), # use 1=yes only for pt-ordered showers !
+        MEMAIN_nqmatch = cms.int32(5), #PID of the flavor until which the QCD radiation are kept in the matching procedure;
+                                         # if nqmatch=4, then all showered partons from b's are NOT taken into account
+         # Note (JY): I think the default should be 5 (b); anyway, don't try -1 as it'll result in a throw...
+        MEMAIN_excres = cms.string(""),
+        outTree_flag = cms.int32(0) # 1=yes, write out the tree for future sanity check
+    ),
+    PythiaParameters = cms.PSet(
+        processParameters = cms.vstring(
+         'Main:timesAllowErrors = 10000',
+        'ParticleDecays:limitTau0 = on',
+         'ParticleDecays:tauMax = 10',
+        'Tune:ee 3',
+        'Tune:pp 5',
+        '8:doForceWidth = 1',
+        '8:m0 = 3000 ! tprime mass',
+        '8:mWidth = 30 ! tprime width',
+        '8:mMin = 2900 ! minimal tprime mass',
+        '8:mMax = 3100 ! maximal tprime mass',
+        '6:m0 = 172.5 ! top mass',
+        '8:oneChannel = 1 0.5 100 24 5 ! add channel tprime to W b with BR=0.5',
+        '8:addChannel = 1 0.25 100 23 6 ! add channel tprime to Z t with BR=0.25',
+        '8:addChannel = 1 0.25 100 25 6 ! add channel tprime to H t with BR=0.25',
+         ),
+        parameterSets = cms.vstring('processParameters')
+    )
+)