Back to home page

Project CMSSW displayed by LXR

 
 

    


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)