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     multiClassPID = cms.PSet(
0038         isPUFilter = cms.bool(False),
0039         preselection = cms.string(""),
0040         model = cms.string("L1Trigger/Phase2L1ParticleFlow/data/multiclassID/multiclass_EE.json"),
0041         variables = cms.VPSet(
0042             cms.PSet(name = cms.string("showerlength"), value = cms.string("showerLength()")),
0043             cms.PSet(name = cms.string("coreshowerlength"), value = cms.string("coreShowerLength()")),
0044             cms.PSet(name = cms.string("eot"), value = cms.string("eot()")),
0045             cms.PSet(name = cms.string("eta"), value = cms.string("abs(eta())")),
0046             cms.PSet(name = cms.string("meanz"), value = cms.string("zBarycenter()")),
0047             cms.PSet(name = cms.string("seetot"), value = cms.string("sigmaEtaEtaTot()")),
0048             cms.PSet(name = cms.string("spptot"), value = cms.string("sigmaPhiPhiTot()")),
0049             cms.PSet(name = cms.string("szz"), value = cms.string("sigmaZZ()")),
0050         ) # variables are not used from here, but putting them here for future use when the code is more neat
0051     ),
0052     useEMInterpretation = cms.string("allKeepHad"), # for all clusters, use EM intepretation to redefine the EM part of the energy
0053     emOnly = cms.bool(False),
0054     etMin = cms.double(1.0),
0055     resol = cms.PSet(
0056         etaBins = cms.vdouble( 1.900,  2.200,  2.500,  2.800,  2.950),
0057         offset  = cms.vdouble( 2.593,  3.089,  2.879,  2.664,  2.947),
0058         scale   = cms.vdouble( 0.120,  0.098,  0.099,  0.098,  0.124),
0059         kind    = cms.string('calo')
0060     ),
0061 )
0062
0063
0064 from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
0065 from Configuration.Eras.Modifier_phase2_hgcalV11_cff import phase2_hgcalV11
0066 phase2_hgcalV10.toModify(l1tPFClustersFromHGC3DClusters,
0067     corrector = "L1Trigger/Phase2L1ParticleFlow/data/hadcorr_HGCal3D_TC_106X.root",
0068     resol = cms.PSet(
0069         etaBins = cms.vdouble( 1.700,  1.900,  2.200,  2.500,  2.800,  2.900),
0070         offset  = cms.vdouble(-0.819,  0.900,  2.032,  2.841,  2.865,  1.237),
0071         scale   = cms.vdouble( 0.320,  0.225,  0.156,  0.108,  0.119,  0.338),
0072         kind    = cms.string('calo')
0073     ),
0074 )
0075 phase2_hgcalV11.toModify(l1tPFClustersFromHGC3DClusters,
0076     corrector = "L1Trigger/Phase2L1ParticleFlow/data/hadcorr_HGCal3D_TC_110X.root",
0077     resol = cms.PSet(
0078         etaBins = cms.vdouble( 1.700,  1.900,  2.200,  2.500,  2.800,  2.900),
0079         offset  = cms.vdouble( 1.793,  1.827,  2.363,  2.538,  2.812,  2.642),
0080         scale   = cms.vdouble( 0.138,  0.137,  0.124,  0.115,  0.106,  0.121),
0081         kind    = cms.string('calo'),
0082     ),
0083 )