python/slimming/isolatedTracks_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002 from TrackingTools.TrackAssociator.default_cfi import TrackAssociatorParameterBlock
0003
0004 tkAssocParamBlock = TrackAssociatorParameterBlock.clone()
0005 tkAssocParamBlock.TrackAssociatorParameters.useMuon = cms.bool(False)
0006 tkAssocParamBlock.TrackAssociatorParameters.useCalo = cms.bool(False)
0007 tkAssocParamBlock.TrackAssociatorParameters.useHO = cms.bool(False)
0008 tkAssocParamBlock.TrackAssociatorParameters.usePreshower = cms.bool(False)
0009 tkAssocParamBlock.TrackAssociatorParameters.EERecHitCollectionLabel = cms.InputTag("reducedEcalRecHitsEE")
0010 tkAssocParamBlock.TrackAssociatorParameters.EBRecHitCollectionLabel = cms.InputTag("reducedEcalRecHitsEB")
0011 tkAssocParamBlock.TrackAssociatorParameters.HBHERecHitCollectionLabel = cms.InputTag("reducedHcalRecHits","hbhereco")
0012 tkAssocParamBlock.TrackAssociatorParameters.HORecHitCollectionLabel = cms.InputTag("reducedHcalRecHits","horeco")
0013
0014 _susySoftDisappearingTrackCut = (
0015     "pt > 10. && (pt > 15. || hitPattern().pixelLayersWithMeasurement() == hitPattern().trackerLayersWithMeasurement())  && "+
0016     " abs(dxy) < 0.02 && abs(dz) < 0.1 && "+
0017     " (miniPFIsolation().chargedHadronIso()/pt < 0.2) && "+
0018     " !(pfLepOverlap) && "+
0019     " pfNeutralSum/pt < 0.2"
0020 )
0021 _exoHighPtTrackCut = (
0022     "pt > 50 && "+
0023     " isHighPurityTrack && "+
0024     " abs(dxy) < 0.5 && abs(dz) < 0.5 && "+
0025     " (miniPFIsolation().chargedHadronIso()/pt < 1.0 || pt > 100)"
0026 )
0027 _exoDisappearingTrackCut = (
0028     "pt > 30 && "+
0029     "abs(dxy) < 0.05 && abs(dz) < 1.0 &&"+
0030     "lostInnerLayers == 0 &&"+
0031     "lostLayers == 0 &&"+
0032     "lostOuterLayers >= 1 &&"+
0033     "(pfIsolationDR03().chargedHadronIso)/pt < 0.1"
0034 )
0035 isolatedTracks = cms.EDProducer("PATIsolatedTrackProducer",
0036     tkAssocParamBlock,
0037     packedPFCandidates = cms.InputTag("packedPFCandidates"),
0038     lostTracks = cms.InputTag("lostTracks"),
0039     generalTracks = cms.InputTag("generalTracks"),
0040     primaryVertices = cms.InputTag("offlinePrimaryVertices"),
0041     caloJets = cms.InputTag("ak4CaloJets"),
0042     dEdxDataStrip = cms.InputTag("dedxHarmonic2"),
0043     dEdxDataPixel = cms.InputTag("dedxPixelHarmonic2"),
0044     dEdxHitInfo = cms.InputTag("dedxHitInfo"),
0045     dEdxHitInfoPrescale = cms.InputTag("dedxHitInfo","prescale"),
0046     addPrescaledDeDxTracks = cms.bool(False),
0047     usePrecomputedDeDxStrip = cms.bool(True),        # if these are set to True, will get estimated DeDx from DeDxData branches
0048     usePrecomputedDeDxPixel = cms.bool(True),        # if set to False, will manually compute using dEdxHitInfo
0049     pT_cut = cms.double(5.0),         # save tracks above this pt
0050     pT_cut_noIso = cms.double(20.0),  # for tracks with at least this pT, don't apply any iso cut
0051     pfIsolation_DR = cms.double(0.3),
0052     pfIsolation_DZ = cms.double(0.1),
0053     miniIsoParams = cms.vdouble(0.05, 0.2, 10.0), # (minDR, maxDR, kT)
0054                                                   # dR for miniiso is max(minDR, min(maxDR, kT/pT))
0055     absIso_cut = cms.double(5.0),
0056     relIso_cut = cms.double(0.2),
0057     miniRelIso_cut = cms.double(0.2),
0058
0059     caloJet_DR = cms.double(0.3),
0060
0061     pflepoverlap_DR = cms.double(0.001),
0062     pflepoverlap_pTmin = cms.double(5.0),
0063
0064     pcRefNearest_DR = cms.double(0.3),
0065     pcRefNearest_pTmin = cms.double(5.0),
0066
0067     pfneutralsum_DR = cms.double(0.05),
0068
0069     useHighPurity = cms.bool(False),
0070
0071     saveDeDxHitInfo = cms.bool(True),
0072     saveDeDxHitInfoCut = cms.string("(%s) || (%s) || (%s)"
0073       % (_susySoftDisappearingTrackCut,_exoHighPtTrackCut,_exoDisappearingTrackCut)),
0074 )
0075
0076 from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
0077 pp_on_AA.toModify(isolatedTracks, useHighPurity = True)
0078
0079 def miniAOD_customizeIsolatedTracksFastSim(process):
0080     """Switch off dE/dx hit info on fast sim, as it's not available"""
0081     process.isolatedTracks.saveDeDxHitInfo = False
0082     return process