File indexing completed on 2021-07-29 01:21:40
0001 import FWCore.ParameterSet.Config as cms
0002 from PhysicsTools.NanoAOD.common_cff import *
0003
0004
0005
0006
0007
0008 mergedGenParticles = cms.EDProducer("MergedGenParticleProducer",
0009 inputPruned = cms.InputTag("prunedGenParticles"),
0010 inputPacked = cms.InputTag("packedGenParticles"),
0011 )
0012
0013 genParticles2HepMC = cms.EDProducer("GenParticles2HepMCConverter",
0014 genParticles = cms.InputTag("mergedGenParticles"),
0015 genEventInfo = cms.InputTag("generator"),
0016 signalParticlePdgIds = cms.vint32(),
0017 )
0018
0019 genParticles2HepMCHiggsVtx = cms.EDProducer("GenParticles2HepMCConverter",
0020 genParticles = cms.InputTag("mergedGenParticles"),
0021 genEventInfo = cms.InputTag("generator"),
0022 signalParticlePdgIds = cms.vint32(25),
0023 )
0024
0025
0026 particleLevel = cms.EDProducer("ParticleLevelProducer",
0027 src = cms.InputTag("genParticles2HepMC:unsmeared"),
0028
0029 doJetClustering = cms.bool(False),
0030 usePromptFinalStates = cms.bool(True),
0031 excludePromptLeptonsFromJetClustering = cms.bool(False),
0032 excludeNeutrinosFromJetClustering = cms.bool(True),
0033
0034 particleMinPt = cms.double(0.),
0035 particleMaxEta = cms.double(5.),
0036
0037 lepConeSize = cms.double(0.1),
0038 lepMinPt = cms.double(1.),
0039 lepMaxEta = cms.double(2.5),
0040
0041 jetConeSize = cms.double(0.4),
0042 jetMinPt = cms.double(10.),
0043 jetMaxEta = cms.double(999.),
0044
0045 fatJetConeSize = cms.double(0.8),
0046 fatJetMinPt = cms.double(170.),
0047 fatJetMaxEta = cms.double(999.),
0048
0049 phoIsoConeSize = cms.double(0.4),
0050 phoMaxRelIso = cms.double(0.5),
0051 phoMinPt = cms.double(1.),
0052 phoMaxEta = cms.double(2.5),
0053 )
0054
0055 rivetProducerHTXS = cms.EDProducer('HTXSRivetProducer',
0056 HepMCCollection = cms.InputTag('genParticles2HepMCHiggsVtx','unsmeared'),
0057 LHERunInfo = cms.InputTag('externalLHEProducer'),
0058 ProductionMode = cms.string('AUTO'),
0059 )
0060
0061
0062
0063 rivetLeptonTable = cms.EDProducer("SimpleCandidateFlatTableProducer",
0064 src = cms.InputTag("particleLevel:leptons"),
0065 cut = cms.string("pt > 15"),
0066 name= cms.string("GenDressedLepton"),
0067 doc = cms.string("Dressed leptons from Rivet-based ParticleLevelProducer"),
0068 singleton = cms.bool(False),
0069 extension = cms.bool(False),
0070 externalVariables = cms.PSet(
0071 hasTauAnc = ExtVar(cms.InputTag("tautagger"),bool, doc="true if Dressed lepton has a tau as ancestor"),
0072 ),
0073 variables = cms.PSet(
0074 P4Vars,
0075 pdgId = Var("pdgId", int, doc="PDG id"),
0076 )
0077 )
0078
0079 rivetPhotonTable = cms.EDProducer("SimpleCandidateFlatTableProducer",
0080 src = cms.InputTag("particleLevel:photons"),
0081 cut = cms.string("pt > 10"),
0082 name= cms.string("GenIsolatedPhoton"),
0083 doc = cms.string("Isolated photons from Rivet-based ParticleLevelProducer"),
0084 singleton = cms.bool(False),
0085 extension = cms.bool(False),
0086 variables = cms.PSet(
0087 P4Vars
0088 )
0089 )
0090
0091 tautagger = cms.EDProducer("GenJetTauTaggerProducer",
0092 src = rivetLeptonTable.src,
0093 )
0094
0095
0096
0097
0098
0099
0100
0101
0102
0103
0104
0105
0106
0107
0108
0109
0110
0111
0112
0113
0114
0115
0116
0117
0118
0119
0120
0121
0122
0123
0124
0125
0126
0127
0128
0129
0130
0131
0132
0133
0134
0135 rivetMetTable = cms.EDProducer("SimpleCandidateFlatTableProducer",
0136 src = cms.InputTag("particleLevel:mets"),
0137 name = cms.string("MET"),
0138 doc = cms.string("MET from Rivet-based ParticleLevelProducer in fiducial volume abs(eta)<5"),
0139 singleton = cms.bool(True),
0140 extension = cms.bool(True),
0141 variables = cms.PSet(
0142 fiducialGenPt = Var("pt", float, precision=10),
0143 fiducialGenPhi = Var("phi", float, precision=10),
0144 ),
0145 )
0146
0147 HTXSCategoryTable = cms.EDProducer("SimpleHTXSFlatTableProducer",
0148 src = cms.InputTag("rivetProducerHTXS","HiggsClassification"),
0149 cut = cms.string(""),
0150 name = cms.string("HTXS"),
0151 doc = cms.string("HTXS classification"),
0152 singleton = cms.bool(True),
0153 extension = cms.bool(False),
0154 variables=cms.PSet(
0155 stage_0 = Var("stage0_cat",int, doc="HTXS stage-0 category"),
0156 stage_1_pTjet30 = Var("stage1_cat_pTjet30GeV",int, doc="HTXS stage-1 category (jet pt>30 GeV)"),
0157 stage_1_pTjet25 = Var("stage1_cat_pTjet25GeV",int, doc="HTXS stage-1 category (jet pt>25 GeV)"),
0158 stage1_1_cat_pTjet30GeV = Var("stage1_1_cat_pTjet30GeV",int,doc="HTXS stage-1.1 category(jet pt>30 GeV)"),
0159 stage1_1_cat_pTjet25GeV = Var("stage1_1_cat_pTjet25GeV",int,doc="HTXS stage-1.1 category(jet pt>25 GeV)"),
0160 stage1_1_fine_cat_pTjet30GeV = Var("stage1_1_fine_cat_pTjet30GeV",int,doc="HTXS stage-1.1-fine category(jet pt>30 GeV)"),
0161 stage1_1_fine_cat_pTjet25GeV = Var("stage1_1_fine_cat_pTjet25GeV",int,doc="HTXS stage-1.1-fine category(jet pt>25 GeV)"),
0162 stage1_2_cat_pTjet30GeV = Var("stage1_2_cat_pTjet30GeV",int,doc="HTXS stage-1.2 category(jet pt>30 GeV)"),
0163 stage1_2_cat_pTjet25GeV = Var("stage1_2_cat_pTjet25GeV",int,doc="HTXS stage-1.2 category(jet pt>25 GeV)"),
0164 stage1_2_fine_cat_pTjet30GeV = Var("stage1_2_fine_cat_pTjet30GeV",int,doc="HTXS stage-1.2-fine category(jet pt>30 GeV)"),
0165 stage1_2_fine_cat_pTjet25GeV = Var("stage1_2_fine_cat_pTjet25GeV",int,doc="HTXS stage-1.2-fine category(jet pt>25 GeV)"),
0166 Higgs_pt = Var("higgs.Pt()",float, doc="pt of the Higgs boson as identified in HTXS", precision=14),
0167 Higgs_y = Var("higgs.Rapidity()",float, doc="rapidity of the Higgs boson as identified in HTXS", precision=12),
0168 njets30 = Var("jets30.size()","uint8", doc="number of jets with pt>30 GeV as identified in HTXS"),
0169 njets25 = Var("jets25.size()","uint8", doc="number of jets with pt>25 GeV as identified in HTXS"),
0170 )
0171 )
0172
0173 lheInfoTable = cms.EDProducer("LHETablesProducer",
0174 lheInfo = cms.VInputTag(cms.InputTag("externalLHEProducer"), cms.InputTag("source")),
0175 precision = cms.int32(14),
0176 storeLHEParticles = cms.bool(True)
0177 )
0178
0179 particleLevelTask = cms.Task(mergedGenParticles,genParticles2HepMC,particleLevel,tautagger,genParticles2HepMCHiggsVtx,rivetProducerHTXS)
0180 particleLevelTablesTask = cms.Task(rivetLeptonTable,rivetPhotonTable,rivetMetTable,HTXSCategoryTable,lheInfoTable)