Phase2L1ParticleFlow/python/l1TkEgAlgoEmulator_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003
0004 CompositeParametersTkEleVec = cms.VPSet(
0005     #Algorithm 0:  Elliptic ID
0006     cms.PSet( # this is not used by any CompositeID model!
0007         model=cms.string(""),
0008         loose_wp=cms.PSet(
0009             bins=cms.vdouble(0),
0010             values=cms.vdouble(-1)
0011             ),
0012         tight_wp=cms.PSet(
0013             bins=cms.vdouble(0),
0014             values=cms.vdouble(-1)
0015             ),
0016         dEta_max = cms.double(-1),
0017         dPhi_max = cms.double(-1),
0018     ),
0019     #Algorithm 1: compositeEE_v0
0020     cms.PSet(
0021         # NOTE: conifer BDT score is log(p/1-p)
0022         # the working points are cuts on BDT output logits [log(p/1-p)]/4 (range -1 to 1 to match the FW dataformat)
0023         model=cms.string("L1Trigger/Phase2L1ParticleFlow/data/egamma/compositeID_EE_v0.json"),
0024         loose_wp=cms.PSet(
0025             bins=cms.vdouble(0., 18., 28., 36.),
0026             values=cms.vdouble(-0.181641, -0.15, 0.075, -0.181641)
0027             ),
0028         tight_wp=cms.PSet(
0029             bins=cms.vdouble(0., 28., 40.),
0030             values=cms.vdouble(0.0527344, 0.3, 0.0527344)
0031             ),
0032         dPhi_max = cms.double(0.2),
0033         dEta_max = cms.double(0.2),
0034     ),
0035     #Algorithm 2: compositeEB_v0
0036     cms.PSet(
0037         model=cms.string("L1Trigger/Phase2L1ParticleFlow/data/egamma/compositeID_EB_v0.json"),
0038         loose_wp=cms.PSet(
0039             bins=cms.vdouble(0),
0040             values=cms.vdouble(-1)
0041             ),
0042         tight_wp=cms.PSet(
0043             bins=cms.vdouble(0, 5, 10, 20, 30, 50),
0044             values=cms.vdouble(0.19, 0.05, -0.35, -0.45, -0.5, -0.38),
0045         ),
0046         dPhi_max = cms.double(0.3),
0047         dEta_max = cms.double(0.03),
0048     ),
0049     #Algorithm 3: compositeEE_v1
0050     cms.PSet(
0051         model=cms.string("L1Trigger/Phase2L1ParticleFlow/data/egamma/compositeID_EE_v1.json"),
0052         loose_wp=cms.PSet(
0053             bins=cms.vdouble(0),
0054             values=cms.vdouble(-1)
0055             ),
0056         tight_wp=cms.PSet(
0057             bins=cms.vdouble(0, 5, 10, 20, 30, 50),
0058             values=cms.vdouble(0.309191, 0.0909913, -0.211824, -0.223846, -0.188443, -0.154760),
0059         ),
0060         dPhi_max = cms.double(0.3),
0061         dEta_max = cms.double(0.03),
0062     ),
0063     #Algorithm 4: compositeEB_v1
0064     cms.PSet(
0065         model=cms.string("L1Trigger/Phase2L1ParticleFlow/data/egamma/compositeID_EB_v1.json"),
0066         loose_wp=cms.PSet(
0067             bins=cms.vdouble(0),
0068             values=cms.vdouble(-1)
0069             ),
0070         tight_wp=cms.PSet(
0071             bins=cms.vdouble(0, 5, 10, 30),
0072             values=cms.vdouble(0.17, 0.018, -0.08, -0.11),
0073         ),
0074         dPhi_max = cms.double(0.3),
0075         dEta_max = cms.double(0.03),
0076     )
0077 )
0078
0079
0080
0081 tkEgAlgoParameters = cms.PSet(
0082     # debug=cms.untracked.uint32(4),
0083     nTRACK=cms.uint32(50),  # very large numbers for first test
0084     nTRACK_EGIN=cms.uint32(50),  # very large numbers for first test
0085     nEMCALO_EGIN=cms.uint32(50),  # very large numbers for first test
0086     nEM_EGOUT=cms.uint32(50),  # very large numbers for first test
0087     doBremRecovery=cms.bool(False),
0088     writeBeforeBremRecovery=cms.bool(True),
0089     filterHwQuality=cms.bool(False),
0090     caloHwQual=cms.int32(4),
0091     doEndcapHwQual=cms.bool(False),
0092     dEtaMaxBrem=cms.double(0.02),
0093     dPhiMaxBrem=cms.double(0.1),
0094     absEtaBoundaries=cms.vdouble(0.0, 0.9, 1.5),
0095     dEtaValues=cms.vdouble(0.025, 0.015, 0.01),  # last was  0.0075  in TDR
0096     dPhiValues=cms.vdouble(0.07, 0.07, 0.07),
0097     caloEtMin=cms.double(0.0),
0098     trkQualityPtMin=cms.double(10.0),
0099     writeEGSta=cms.bool(False),
0100     tkIsoParametersTkEm=cms.PSet(
0101         tkQualityPtMin=cms.double(2.),
0102         dZ=cms.double(0.6),
0103         dRMin=cms.double(0.07),
0104         dRMax=cms.double(0.30),
0105     ),
0106     tkIsoParametersTkEle=cms.PSet(
0107         tkQualityPtMin=cms.double(2.),
0108         dZ=cms.double(0.6),
0109         dRMin=cms.double(0.03),
0110         dRMax=cms.double(0.20),
0111     ),
0112     pfIsoParametersTkEm=cms.PSet(
0113         tkQualityPtMin=cms.double(1.),
0114         dZ=cms.double(0.6),
0115         dRMin=cms.double(0.07),
0116         dRMax=cms.double(0.30),
0117     ),
0118     pfIsoParametersTkEle=cms.PSet(
0119         tkQualityPtMin=cms.double(1.),
0120         dZ=cms.double(0.6),
0121         dRMin=cms.double(0.03),
0122         dRMax=cms.double(0.20),
0123     ),
0124     doTkIso=cms.bool(True),
0125     doPfIso=cms.bool(True),
0126     hwIsoTypeTkEle=cms.uint32(0),
0127     hwIsoTypeTkEm=cms.uint32(0),
0128     algorithm=cms.uint32(0), # 0 = elliptic , 1 = composite EE, 2 = composite EB
0129     nCompCandPerCluster=cms.uint32(3),
0130     compositeParametersTkEle=CompositeParametersTkEleVec,
0131 )
0132
0133 tkEgSorterParameters = cms.PSet(
0134     nObjToSort=cms.uint32(6),
0135     nObjSorted=cms.uint32(16),
0136 )