Phase2L1ParticleFlow/python/l1tPFClustersFromHGC3DClusters_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 l1tPFClustersFromHGC3DClusters = cms.EDProducer("PFClusterProducerFromHGC3DClusters",
0004     src = cms.InputTag("l1tHGCalBackEndLayer2Producer","HGCalBackendLayer2Processor3DClustering"),
0005     corrector = cms.string("L1Trigger/Phase2L1ParticleFlow/data/hadcorr_HGCal3D_TC.root"),
0006     correctorEmfMax = cms.double(1.125),
0007     preEmId  = cms.string("hOverE < 0.3 && hOverE >= 0"),
0008     emVsPionID = cms.PSet(
0009         isPUFilter = cms.bool(False),
0010         preselection = cms.string(""),
0011         method = cms.string("BDT"), # "" to be disabled, "BDT" to be enabled
0012         variables = cms.VPSet(
0013             cms.PSet(name = cms.string("fabs(eta)"), value = cms.string("abs(eta())")),
0014             cms.PSet(name = cms.string("eMax"), value = cms.string("eMax()")),
0015             cms.PSet(name = cms.string("sigmaPhiPhiTot"), value = cms.string("sigmaPhiPhiTot()")),
0016             cms.PSet(name = cms.string("sigmaZZ"), value = cms.string("sigmaZZ()")),
0017             cms.PSet(name = cms.string("layer50percent"), value = cms.string("layer50percent()")),
0018             cms.PSet(name = cms.string("triggerCells67percent"), value = cms.string("triggerCells67percent()")),
0019         ),
0020         weightsFile = cms.string("L1Trigger/Phase2L1ParticleFlow/data/hgcal_egID/Photon_vs_Pion_BDTweights_1116.xml.gz"),
0021         wp = cms.string("0.05")
0022     ),
0023     emVsPUID = cms.PSet(
0024         isPUFilter = cms.bool(True),
0025         preselection = cms.string(""),
0026         method = cms.string("BDT"), # "" to be disabled, "BDT" to be enabled
0027         variables = cms.VPSet(
0028             cms.PSet(name = cms.string("eMax"), value = cms.string("eMax()")),
0029             cms.PSet(name = cms.string("eMaxOverE"), value = cms.string("eMax()/energy()")),
0030             cms.PSet(name = cms.string("sigmaPhiPhiTot"), value = cms.string("sigmaPhiPhiTot()")),
0031             cms.PSet(name = cms.string("sigmaRRTot"), value = cms.string("sigmaRRTot()")),
0032             cms.PSet(name = cms.string("triggerCells90percent"), value = cms.string("triggerCells90percent()")),
0033         ),
0034         weightsFile = cms.string("L1Trigger/Phase2L1ParticleFlow/data/hgcal_egID/Photon_Pion_vs_Neutrino_BDTweights_1116.xml.gz"),
0035         wp = cms.string("0.15")
0036     ),
0037     useEMInterpretation = cms.string("allKeepHad"), # for all clusters, use EM intepretation to redefine the EM part of the energy
0038     emOnly = cms.bool(False),
0039     etMin = cms.double(1.0),
0040     resol = cms.PSet(
0041         etaBins = cms.vdouble( 1.900,  2.200,  2.500,  2.800,  2.950),
0042         offset  = cms.vdouble( 2.593,  3.089,  2.879,  2.664,  2.947),
0043         scale   = cms.vdouble( 0.120,  0.098,  0.099,  0.098,  0.124),
0044         kind    = cms.string('calo')
0045     ),
0046 )
0047
0048
0049 from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
0050 from Configuration.Eras.Modifier_phase2_hgcalV11_cff import phase2_hgcalV11
0051 phase2_hgcalV10.toModify(l1tPFClustersFromHGC3DClusters,
0052     corrector = "L1Trigger/Phase2L1ParticleFlow/data/hadcorr_HGCal3D_TC_106X.root",
0053     resol = cms.PSet(
0054         etaBins = cms.vdouble( 1.700,  1.900,  2.200,  2.500,  2.800,  2.900),
0055         offset  = cms.vdouble(-0.819,  0.900,  2.032,  2.841,  2.865,  1.237),
0056         scale   = cms.vdouble( 0.320,  0.225,  0.156,  0.108,  0.119,  0.338),
0057         kind    = cms.string('calo')
0058     ),
0059 )
0060 phase2_hgcalV11.toModify(l1tPFClustersFromHGC3DClusters,
0061     corrector = "L1Trigger/Phase2L1ParticleFlow/data/hadcorr_HGCal3D_TC_110X.root",
0062     resol = cms.PSet(
0063         etaBins = cms.vdouble( 1.700,  1.900,  2.200,  2.500,  2.800,  2.900),
0064         offset  = cms.vdouble( 1.793,  1.827,  2.363,  2.538,  2.812,  2.642),
0065         scale   = cms.vdouble( 0.138,  0.137,  0.124,  0.115,  0.106,  0.121),
0066         kind    = cms.string('calo'),
0067     ),
0068 )