Adding BPH skim module

Configuration/Skimming/python/PDWG_BPHSkim_cff.py

@@ -0,0 +1,88 @@
+import FWCore.ParameterSet.Config as cms
+
+from Configuration.StandardSequences.Reconstruction_cff import *
+
+# muons with trigger info
+import PhysicsTools.PatAlgos.producersLayer1.muonProducer_cfi
+oniaPATMuonsWithoutTrigger = PhysicsTools.PatAlgos.producersLayer1.muonProducer_cfi.patMuons.clone(
+    muonSource = 'muons',
+    embedTrack          = True,
+    embedCombinedMuon   = True,
+    embedStandAloneMuon = True,
+    embedPFCandidate    = False,
+    embedCaloMETMuonCorrs = cms.bool(False),
+    embedTcMETMuonCorrs   = cms.bool(False),
+    embedPfEcalEnergy     = cms.bool(False),
+    embedPickyMuon = False,
+    embedTpfmsMuon = False, 
+    userIsolation = cms.PSet(),   # no extra isolation beyond what's in reco::Muon itself
+    isoDeposits = cms.PSet(),     # no heavy isodeposits
+    addGenMatch = False,          # no mc
+)
+
+oniaSelectedMuons = cms.EDFilter('PATMuonSelector',
+   src = cms.InputTag('oniaPATMuonsWithoutTrigger'),
+   cut = cms.string('muonID(\"TMOneStationTight\")'
+                    ' && abs(innerTrack.dxy) < 0.3'
+                    ' && abs(innerTrack.dz)  < 20.'
+                    ' && innerTrack.hitPattern.trackerLayersWithMeasurement > 5'
+                    ' && innerTrack.hitPattern.pixelLayersWithMeasurement > 0'
+                    ' && innerTrack.quality(\"highPurity\")'
+                    ' && ((abs(eta) <= 0.9 && pt > 2.5) || (0.9 < abs(eta) <= 2.4 && pt > 1.5))'
+   ),
+   filter = cms.bool(True)
+)
+
+# tracks
+oniaSelectedTracks=cms.EDFilter("TrackSelector",
+     src = cms.InputTag("generalTracks"),
+     cut = cms.string('pt > 0.7 && abs(eta) <= 3.0'
+                      '&& charge !=0'
+                      '&& quality(\"highPurity\")')     
+)
+
+# dimuon = Onia2MUMU
+from HeavyFlavorAnalysis.Onia2MuMu.onia2MuMuPAT_cfi import *
+onia2MuMuPAT.muons=cms.InputTag('oniaSelectedMuons')
+onia2MuMuPAT.primaryVertexTag=cms.InputTag('offlinePrimaryVertices')
+onia2MuMuPAT.beamSpotTag=cms.InputTag('offlineBeamSpot')
+onia2MuMuPAT.dimuonSelection=cms.string("0.2 < mass && abs(daughter('muon1').innerTrack.dz - daughter('muon2').innerTrack.dz) < 25")
+onia2MuMuPAT.addMCTruth = cms.bool(False)
+
+onia2MuMuPATCounter = cms.EDFilter('CandViewCountFilter',
+      src = cms.InputTag('onia2MuMuPAT'),
+      minNumber = cms.uint32(1),
+      filter = cms.bool(True)
+   )
+
+# make photon candidate conversions for P-wave studies
+from HeavyFlavorAnalysis.Onia2MuMu.OniaPhotonConversionProducer_cfi import PhotonCandidates as oniaPhotonCandidates
+
+# add v0 with tracks embed
+from HeavyFlavorAnalysis.Onia2MuMu.OniaAddV0TracksProducer_cfi import *
+
+# Pick branches you want to keep
+BPHSkim_EventContent = cms.PSet(
+     outputCommands = cms.untracked.vstring(
+                     'drop *',
+                     'keep recoVertexs_offlinePrimaryVertices_*_*',
+                     'keep *_offlineBeamSpot_*_*',
+                     'keep *_TriggerResults_*_HLT',
+                     'keep *_gtDigis_*_RECO',
+                     'keep *_oniaSelectedTracks_*_*',
+                     'keep *_oniaPhotonCandidates_*_*',
+                     'keep *_onia2MuMuPAT_*_*',
+                     'keep *_oniaV0Tracks_*_*',
+                     'keep PileupSummaryInfos_*_*_*'
+     )
+)
+
+BPHSkimSequence = cms.Sequence(
+            oniaPATMuonsWithoutTrigger *
+	    oniaSelectedMuons *
+            onia2MuMuPAT *
+	    onia2MuMuPATCounter *
+	    oniaPhotonCandidates *
+	    oniaV0Tracks *
+	    oniaSelectedTracks
+)

Configuration/Skimming/python/Skims_PDWG_cff.py

@@ -9,6 +9,19 @@
 
 #####################
 
+from Configuration.Skimming.PDWG_BPHSkim_cff import *
+BPHSkimPath = cms.Path(BPHSkimSequence)
+SKIMStreamBPHSkim = cms.FilteredStream(
+    responsible = 'PDWG',
+    name = 'BPHSkim',
+    paths = (BPHSkimPath),
+    content = BPHSkim_EventContent.outputCommands,
+    selectEvents = cms.untracked.PSet(),
+    dataTier = cms.untracked.string('USER')
+    )
+
+#####################
+
 from Configuration.Skimming.PDWG_DiJetAODSkim_cff import *
 diJetAveSkimPath = cms.Path(DiJetAveSkim_Trigger)
 SKIMStreamDiJet = cms.FilteredStream(