RecoTrack/python/GenParticleSelectionsForEfficiency_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 GenParticleSelectionForEfficiency = cms.PSet(
0004     lipGP = cms.double(30.0),
0005     chargedOnlyGP = cms.bool(True),
0006     pdgIdGP = cms.vint32(),
0007     minRapidityGP = cms.double(-2.5),
0008     ptMinGP = cms.double(0.005),
0009     maxRapidityGP = cms.double(2.5),
0010     tipGP = cms.double(60),
0011     statusGP = cms.int32(1)
0012 )
0013 from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel
0014 from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
0015 phase1Pixel.toModify(GenParticleSelectionForEfficiency,minRapidityGP = -3.0, maxRapidityGP = 3.0)
0016 phase2_tracker.toModify(GenParticleSelectionForEfficiency,minRapidityGP = -4.5, maxRapidityGP = 4.5)
0017
0018 generalGpSelectorBlock = cms.PSet(
0019     status = cms.int32(1),
0020     lip = cms.double(30.0),
0021     chargedOnly = cms.bool(True),
0022     pdgId = cms.vint32(),
0023     minRapidity = cms.double(-2.5),
0024     ptMin = cms.double(0.9),
0025     maxRapidity = cms.double(2.5),
0026     tip = cms.double(3.5),
0027     invertRapidityCut = cms.bool(False),
0028     maxPhi = cms.double(3.2),
0029     minPhi = cms.double(-3.2)
0030 )
0031
0032
0033 GpSelectorForEfficiencyVsEtaBlock = cms.PSet(
0034     status = cms.int32(1),
0035     lip = cms.double(30.0),
0036     chargedOnly = cms.bool(True),
0037     pdgId = cms.vint32(),
0038     minRapidity = cms.double(-2.5),
0039     ptMin = cms.double(0.9),
0040     maxRapidity = cms.double(2.5),
0041     tip = cms.double(3.5),
0042     invertRapidityCut = cms.bool(False),
0043     maxPhi = cms.double(3.2),
0044     minPhi = cms.double(-3.2)
0045 )
0046
0047 GpSelectorForEfficiencyVsPhiBlock = cms.PSet(
0048     status = cms.int32(1),
0049     lip = cms.double(30.0),
0050     chargedOnly = cms.bool(True),
0051     pdgId = cms.vint32(),
0052     minRapidity = cms.double(-2.5),
0053     ptMin = cms.double(0.9),
0054     maxRapidity = cms.double(2.5),
0055     tip = cms.double(3.5),
0056     invertRapidityCut = cms.bool(False),
0057     maxPhi = cms.double(3.2),
0058     minPhi = cms.double(-3.2)
0059 )
0060
0061 GpSelectorForEfficiencyVsPtBlock = cms.PSet(
0062     status = cms.int32(1),
0063     chargedOnly = cms.bool(True),
0064     pdgId = cms.vint32(),
0065     minRapidity = cms.double(-2.5),
0066     maxRapidity = cms.double(2.5),
0067     ptMin = cms.double(0.050),
0068     tip = cms.double(3.5),
0069     lip = cms.double(30.0),
0070     invertRapidityCut = cms.bool(False),
0071     maxPhi = cms.double(3.2),
0072     minPhi = cms.double(-3.2)
0073 )
0074
0075 GpSelectorForEfficiencyVsVTXRBlock = cms.PSet(
0076     status = cms.int32(1),
0077     chargedOnly = cms.bool(True),
0078     pdgId = cms.vint32(),
0079     minRapidity = cms.double(-2.5),
0080     ptMin = cms.double(0.9),
0081     maxRapidity = cms.double(2.5),
0082     lip = cms.double(30.0),
0083     tip = cms.double(30.0),
0084     invertRapidityCut = cms.bool(False),
0085     maxPhi = cms.double(3.2),
0086     minPhi = cms.double(-3.2)
0087 )
0088
0089 GpSelectorForEfficiencyVsVTXZBlock = cms.PSet(
0090     status = cms.int32(1),
0091     chargedOnly = cms.bool(True),
0092     pdgId = cms.vint32(),
0093     minRapidity = cms.double(-2.5),
0094     ptMin = cms.double(0.9),
0095     maxRapidity = cms.double(2.5),
0096     lip = cms.double(35.0),
0097     tip = cms.double(3.5),
0098     invertRapidityCut = cms.bool(False),
0099     maxPhi = cms.double(3.2),
0100     minPhi = cms.double(-3.2)
0101 )
0102
0103 def _modifyForPhase1(pset):
0104     pset.minRapidity = -3
0105     pset.maxRapidity = 3
0106     pset.tip = 2.5 # beampipe is around 2.0, BPIX1 is at 2.9
0107 from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel
0108 phase1Pixel.toModify(generalGpSelectorBlock,             _modifyForPhase1)
0109 phase1Pixel.toModify(GpSelectorForEfficiencyVsEtaBlock,  _modifyForPhase1)
0110 phase1Pixel.toModify(GpSelectorForEfficiencyVsPhiBlock,  _modifyForPhase1)
0111 phase1Pixel.toModify(GpSelectorForEfficiencyVsPtBlock,   _modifyForPhase1)
0112 phase1Pixel.toModify(GpSelectorForEfficiencyVsVTXRBlock, _modifyForPhase1)
0113 phase1Pixel.toModify(GpSelectorForEfficiencyVsVTXZBlock, _modifyForPhase1)
0114
0115 def _modifyForPhase2(pset):
0116     pset.minRapidity = -4.5
0117     pset.maxRapidity = 4.5
0118     pset.tip = 2.5 # IT1 will be around 3.0 (as in Phase1)
0119 from Configuration.Eras.Modifier_phase2_tracker_cff import phase2_tracker
0120 phase2_tracker.toModify(generalGpSelectorBlock,             _modifyForPhase2)
0121 phase2_tracker.toModify(GpSelectorForEfficiencyVsEtaBlock,  _modifyForPhase2)
0122 phase2_tracker.toModify(GpSelectorForEfficiencyVsPhiBlock,  _modifyForPhase2)
0123 phase2_tracker.toModify(GpSelectorForEfficiencyVsPtBlock,   _modifyForPhase2)
0124 phase2_tracker.toModify(GpSelectorForEfficiencyVsVTXRBlock, _modifyForPhase2)
0125 phase2_tracker.toModify(GpSelectorForEfficiencyVsVTXZBlock, _modifyForPhase2)