File indexing completed on 2022-11-09 11:29:57
0001 import FWCore.ParameterSet.Config as cms
0002 from PhysicsTools.NanoAOD.common_cff import *
0003 from PhysicsTools.NanoAOD.globalVariablesTableProducer_cfi import globalVariablesTableProducer
0004 from PhysicsTools.NanoAOD.simpleBeamspotFlatTableProducer_cfi import simpleBeamspotFlatTableProducer
0005 from PhysicsTools.NanoAOD.simpleGenEventFlatTableProducer_cfi import simpleGenEventFlatTableProducer
0006 from PhysicsTools.NanoAOD.simpleGenFilterFlatTableProducerLumi_cfi import simpleGenFilterFlatTableProducerLumi
0007
0008 beamSpotTable = simpleBeamspotFlatTableProducer.clone(
0009 src = cms.InputTag("offlineBeamSpot"),
0010 name = cms.string("BeamSpot"),
0011 doc = cms.string("offlineBeamSpot, the offline reconstructed beamspot"),
0012 variables = cms.PSet(
0013 type = Var("type()","int8",doc="BeamSpot type (Unknown = -1, Fake = 0, LHC = 1, Tracker = 2)"),
0014 z = Var("position().z()",float,doc="BeamSpot center, z coordinate (cm)",precision=-1),
0015 zError = Var("z0Error()",float,doc="Error on BeamSpot center, z coordinate (cm)",precision=-1),
0016 sigmaZ = Var("sigmaZ()",float,doc="Width of BeamSpot in z (cm)",precision=-1),
0017 sigmaZError = Var("sigmaZ0Error()",float,doc="Error on width of BeamSpot in z (cm)",precision=-1),
0018 ),
0019 )
0020
0021 rhoTable = globalVariablesTableProducer.clone(
0022 name = cms.string("Rho"),
0023 variables = cms.PSet(
0024 fixedGridRhoAll = ExtVar( cms.InputTag("fixedGridRhoAll"), "double", doc = "rho from all PF Candidates, no foreground removal (for isolation of prompt photons)" ),
0025 fixedGridRhoFastjetAll = ExtVar( cms.InputTag("fixedGridRhoFastjetAll"), "double", doc = "rho from all PF Candidates, used e.g. for JECs" ),
0026 fixedGridRhoFastjetCentralNeutral = ExtVar( cms.InputTag("fixedGridRhoFastjetCentralNeutral"), "double", doc = "rho from neutral PF Candidates with |eta| < 2.5, used e.g. for rho corrections of some lepton isolations" ),
0027 fixedGridRhoFastjetCentralCalo = ExtVar( cms.InputTag("fixedGridRhoFastjetCentralCalo"), "double", doc = "rho from calo towers with |eta| < 2.5, used e.g. egamma PFCluster isolation" ),
0028 fixedGridRhoFastjetCentral = ExtVar( cms.InputTag("fixedGridRhoFastjetCentral"), "double", doc = "rho from all PF Candidates for central region, used e.g. for JECs" ),
0029 fixedGridRhoFastjetCentralChargedPileUp = ExtVar( cms.InputTag("fixedGridRhoFastjetCentralChargedPileUp"), "double", doc = "rho from charged PF Candidates for central region, used e.g. for JECs" ),
0030 )
0031 )
0032
0033 puTable = cms.EDProducer("NPUTablesProducer",
0034 src = cms.InputTag("slimmedAddPileupInfo"),
0035 pvsrc = cms.InputTag("offlineSlimmedPrimaryVertices"),
0036 zbins = cms.vdouble( [0.0,1.7,2.6,3.0,3.5,4.2,5.2,6.0,7.5,9.0,12.0] ),
0037 savePtHatMax = cms.bool(False),
0038 )
0039
0040 genTable = simpleGenEventFlatTableProducer.clone(
0041 src = cms.InputTag("generator"),
0042 name= cms.string("Generator"),
0043 doc = cms.string("Generator information"),
0044 variables = cms.PSet(
0045 x1 = Var( "?hasPDF?pdf().x.first:-1", float, doc="x1 fraction of proton momentum carried by the first parton",precision=14 ),
0046 x2 = Var( "?hasPDF?pdf().x.second:-1", float, doc="x2 fraction of proton momentum carried by the second parton",precision=14 ),
0047 xpdf1 = Var( "?hasPDF?pdf().xPDF.first:-1", float, doc="x*pdf(x) for the first parton", precision=14 ),
0048 xpdf2 = Var( "?hasPDF?pdf().xPDF.second:-1", float, doc="x*pdf(x) for the second parton", precision=14 ),
0049 id1 = Var( "?hasPDF?pdf().id.first:-1", int, doc="id of first parton", precision=6 ),
0050 id2 = Var( "?hasPDF?pdf().id.second:-1", int, doc="id of second parton", precision=6 ),
0051 scalePDF = Var( "?hasPDF?pdf().scalePDF:-1", float, doc="Q2 scale for PDF", precision=14 ),
0052 binvar = Var("?hasBinningValues()?binningValues()[0]:-1", float, doc="MC generation binning value", precision=14),
0053 weight = Var("weight()", float,doc="MC generator weight", precision=14),
0054 ),
0055 )
0056
0057 genFilterTable = simpleGenFilterFlatTableProducerLumi.clone(
0058 src = cms.InputTag("genFilterEfficiencyProducer"),
0059 name= cms.string("GenFilter"),
0060 doc = cms.string("Generator filter information"),
0061 variables = cms.PSet(
0062 numEventsTotal = Var("numEventsTotal()", int, doc="generator filter: total number of events", precision=6),
0063 numEventsPassed = Var("numEventsPassed()", int, doc="generator filter: passed number of events", precision=6),
0064 filterEfficiency = Var("filterEfficiency()", float, doc="generator filter: efficiency", precision=14),
0065 filterEfficiencyError = Var("filterEfficiencyError()", float, doc="generator filter: efficiency error", precision=14),
0066 ),
0067 )
0068
0069 globalTablesTask = cms.Task(beamSpotTable, rhoTable)
0070 globalTablesMCTask = cms.Task(puTable,genTable,genFilterTable)