RecoTrack/python/TrackingParticleSelectionsForEfficiency_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel
0003 from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
0004
0005 generalTpSelectorBlock = cms.PSet(
0006     lip = cms.double(30.0),
0007     chargedOnly = cms.bool(True),
0008     pdgId = cms.vint32(),
0009     signalOnly = cms.bool(True),
0010     intimeOnly = cms.bool(False),
0011     stableOnly = cms.bool(False),
0012     minRapidity = cms.double(-2.5),
0013     minHit = cms.int32(0),
0014     ptMin = cms.double(0.9),
0015     ptMax = cms.double(1e100),
0016     maxRapidity = cms.double(2.5),
0017     tip = cms.double(3.5),
0018     minPhi = cms.double(-3.2),
0019     maxPhi = cms.double(3.2),
0020     invertRapidityCut = cms.bool(False)
0021 )
0022
0023 from Configuration.Eras.Modifier_fastSim_cff import fastSim
0024 fastSim.toModify(generalTpSelectorBlock, stableOnly = True)
0025
0026 generalTpSelectorForEfficiencyBlock = generalTpSelectorBlock.clone() # TP selector block for efficiency (with additional selections)
0027
0028 def _modifyForPhase1Efficiency(pset):
0029     pset.tip = 2.5 # beampipe is around 2.0, BPIX1 is at 2.9
0030
0031 phase1Pixel.toModify(generalTpSelectorForEfficiencyBlock, _modifyForPhase1Efficiency)
0032
0033 def _modifyForPhase2Efficiency(pset):
0034     pset.minRapidity = -3.5 # within efficient eta range in phase-2
0035     pset.maxRapidity = 3.5 # within efficient eta range in phase-2
0036     pset.tip = 2.5 # IT1 will be around 3.0 (as in Phase1)
0037
0038 phase2_tracker.toModify(generalTpSelectorForEfficiencyBlock, _modifyForPhase2Efficiency)
0039
0040 TpSelectorForEfficiencyVsEtaBlock = generalTpSelectorForEfficiencyBlock.clone()
0041 TpSelectorForEfficiencyVsPhiBlock = generalTpSelectorForEfficiencyBlock.clone()
0042 TpSelectorForEfficiencyVsPtBlock = generalTpSelectorForEfficiencyBlock.clone(ptMin = 0.050 )
0043 TpSelectorForEfficiencyVsVTXRBlock = generalTpSelectorForEfficiencyBlock.clone(tip = 60.0)
0044 TpSelectorForEfficiencyVsVTXZBlock = generalTpSelectorForEfficiencyBlock.clone()
0045
0046 def _modifyForPhase1(pset):
0047     pset.minRapidity = -3
0048     pset.maxRapidity = 3
0049     pset.tip = 2.5 # beampipe is around 2.0, BPIX1 is at 2.9
0050
0051 phase1Pixel.toModify(generalTpSelectorBlock, _modifyForPhase1) # for general TP selector, extend eta to full acceptance
0052 phase1Pixel.toModify(TpSelectorForEfficiencyVsEtaBlock, _modifyForPhase1) # for efficiency vs eta, also extend eta to full acceptance
0053
0054 def _modifyForPhase2(pset):
0055     pset.minRapidity = -4.5
0056     pset.maxRapidity = 4.5
0057     pset.tip = 2.5 # IT1 will be around 3.0 (as in Phase1)
0058
0059 phase2_tracker.toModify(generalTpSelectorBlock, _modifyForPhase2) # for general TP selector, extend eta to full acceptance
0060 phase2_tracker.toModify(TpSelectorForEfficiencyVsEtaBlock, _modifyForPhase2) # for efficiency vs eta, also extend eta to full acceptance