Prompt skims for PbPb and pp at 5 TeV

Configuration/Skimming/python/HI_BJetSkim_cff.py

@@ -0,0 +1,34 @@
+import FWCore.ParameterSet.Config as cms
+
+# HLT jet trigger
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltJetHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltJetHI.HLTPaths = ["HLT_PuAK4CaloBJetCSV80_Eta2p1_v*"]
+hltJetHI.throw = False
+hltJetHI.andOr = True
+
+# selection of valid vertex
+primaryVertexFilterForBJets = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("hiSelectedVertex"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"), 
+    filter = cms.bool(True),   # otherwise it won't filter the events
+    )
+
+# jet energy correction (L2+L3)
+from JetMETCorrections.Configuration.JetCorrectionServicesAllAlgos_cff import *
+akVs4PFJetsL2L3 = cms.EDProducer('PFJetCorrectionProducer',
+    src = cms.InputTag('akVs4PFJets'),
+    correctors = cms.vstring('ak4PFL2L3')
+    )
+
+hiPtBJet = cms.EDFilter("PFJetSelector",
+    src = cms.InputTag("akVs4PFJetsL2L3"),
+    cut = cms.string("pt > 110")
+    )
+
+# dijet skim sequence
+bJetSkimSequence = cms.Sequence(hltJetHI
+                                 * primaryVertexFilterForBJets
+                                 * akVs4PFJetsL2L3
+                                 * hiPtBJet
+                                 )

Configuration/Skimming/python/HI_D0MesonSkim_cff.py

@@ -0,0 +1,20 @@
+import FWCore.ParameterSet.Config as cms
+
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltDmeson60 = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltDmeson60.HLTPaths = ["HLT_DmesonHITrackingGlobal_Dpt60_v*"]
+hltDmeson60.throw = False
+hltDmeson60.andOr = True
+
+# selection of valid vertex                                                                                                                                
+primaryVertexFilterForD0Meson = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("hiSelectedVertex"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"),
+    filter = cms.bool(True),   # otherwise it won't filter the events                                                                                       
+
+    )
+
+d0MesonSkimSequence = cms.Sequence(
+    primaryVertexFilterForD0Meson 
+    *hltDmeson60
+)

Configuration/Skimming/python/HI_HighPtJetSkim_cff.py

@@ -0,0 +1,22 @@
+import FWCore.ParameterSet.Config as cms
+
+# HLT jet trigger
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltHIJet150 = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltHIJet150.HLTPaths = ["HLT_PuAK4CaloJet150_Eta2p1_v*"]
+hltHIJet150.throw = False
+hltHIJet150.andOr = True
+
+# selection of valid vertex
+primaryVertexFilterForHighPtJets = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("hiSelectedVertex"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"), 
+    filter = cms.bool(True),   # otherwise it won't filter the events
+    )
+
+highPtJetSkimSequence = cms.Sequence(
+        hltHIJet150*
+        primaryVertexFilterForHighPtJets 
+)
+
+

Configuration/Skimming/python/HI_MinBiasSkim_cff.py

@@ -0,0 +1,21 @@
+import FWCore.ParameterSet.Config as cms
+
+# HLT dimuon trigger
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltMinBiasHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltMinBiasHI.HLTPaths = ["HLT_L1MinimumBiasHF2AND_v*"]
+hltMinBiasHI.throw = False
+hltMinBiasHI.andOr = True
+
+# selection of valid vertex
+primaryVertexFilterForMinBias = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("hiSelectedVertex"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"), 
+    filter = cms.bool(True),   # otherwise it won't filter the events
+    )
+
+# MinBias skim sequence
+minBiasSkimSequence = cms.Sequence(
+    hltMinBiasHI *
+    primaryVertexFilterForMinBias
+    )

Configuration/Skimming/python/HI_OniaCentralSkim_cff.py

@@ -0,0 +1,49 @@
+import FWCore.ParameterSet.Config as cms
+
+# HLT dimuon trigger
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltCentralOniaMMHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltCentralOniaMMHI.HLTPaths = ["HLT_HIL3DoubleMu0_Cent30_OS_m2p5to4p5_v*","HLT_HIL3DoubleMu0_Cent30_OS_m7to14_v*"] 
+hltCentralOniaMMHI.throw = False
+hltCentralOniaMMHI.andOr = True
+
+# selection of valid vertex
+primaryVertexFilterForOniaMMCentral = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("hiSelectedVertex"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"), 
+    filter = cms.bool(True),   # otherwise it won't filter the events
+    )
+
+# selection of dimuons with mass in Jpsi or 
+muonSelectorForOniaMMCentral = cms.EDFilter("MuonSelector",
+    src = cms.InputTag("muons"),
+    cut = cms.string("(isTrackerMuon && isGlobalMuon) && pt > 1.5"),
+    filter = cms.bool(True)
+    )
+
+muonFilterForOniaMMCentral = cms.EDFilter("MuonCountFilter",
+    src = cms.InputTag("muonSelectorForOniaMMCentral"),
+    minNumber = cms.uint32(2)
+    )
+
+# opposite charge only 
+dimuonMassCutForOniaMMCentral = cms.EDProducer("CandViewShallowCloneCombiner",
+    checkCharge = cms.bool(True),
+    cut = cms.string(' (2.6 < mass < 3.5) || (7.0 < mass < 14.0)'),
+    decay = cms.string("muonSelectorForOniaMMCentral@+ muonSelectorForOniaMMCentral@-")
+    )
+
+dimuonMassCutFilterForOniaMMCentral = cms.EDFilter("CandViewCountFilter",
+    src = cms.InputTag("dimuonMassCutForOniaMMCentral"),
+    minNumber = cms.uint32(1)
+    )
+
+# onia skim sequence
+oniaCentralSkimSequence = cms.Sequence(
+    hltCentralOniaMMHI *
+    primaryVertexFilterForOniaMMCentral *
+    muonSelectorForOniaMMCentral *
+    muonFilterForOniaMMCentral *
+    dimuonMassCutForOniaMMCentral *
+    dimuonMassCutFilterForOniaMMCentral
+    )

Configuration/Skimming/python/HI_OniaPeripheralSkim_cff.py

@@ -0,0 +1,50 @@
+import FWCore.ParameterSet.Config as cms
+
+# HLT dimuon trigger
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltPeripheralOniaMMHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltPeripheralOniaMMHI.HLTPaths = ["HLT_HIL1DoubleMu0_2HF_Cent30100_v*","HLT_HIL1DoubleMu0_2HF0_Cent30100_v*"] 
+hltPeripheralOniaMMHI.throw = False
+hltPeripheralOniaMMHI.andOr = True
+
+# selection of valid vertex
+primaryVertexFilterForOniaMMPeripheral = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("hiSelectedVertex"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"), 
+    filter = cms.bool(True),   # otherwise it won't filter the events
+    )
+
+# selection of dimuons with mass in Jpsi or 
+muonSelectorForOniaMMPeripheral = cms.EDFilter("MuonSelector",
+    src = cms.InputTag("muons"),
+    cut = cms.string("(isTrackerMuon && isGlobalMuon) && pt > 1.5"),
+    filter = cms.bool(True)
+    )
+
+muonFilterForOniaMMPeripheral = cms.EDFilter("MuonCountFilter",
+    src = cms.InputTag("muonSelectorForOniaMMPeripheral"),
+    minNumber = cms.uint32(2)
+    )
+
+# opposite charge only 
+dimuonMassCutForOniaMMPeripheral = cms.EDProducer("CandViewShallowCloneCombiner",
+    checkCharge = cms.bool(True),
+    cut = cms.string(' (2.6 < mass < 3.5) || (7.0 < mass < 14.0)'),
+    decay = cms.string("muonSelectorForOniaMMPeripheral@+ muonSelectorForOniaMMPeripheral@-")
+    )
+
+dimuonMassCutFilterForOniaMMPeripheral = cms.EDFilter("CandViewCountFilter",
+    src = cms.InputTag("dimuonMassCutForOniaMMPeripheral"),
+    minNumber = cms.uint32(1)
+    )
+
+# onia skim sequence
+oniaPeripheralSkimSequence = cms.Sequence(
+    hltPeripheralOniaMMHI *
+    primaryVertexFilterForOniaMMPeripheral *
+    muonSelectorForOniaMMPeripheral *
+    muonFilterForOniaMMPeripheral *
+    dimuonMassCutForOniaMMPeripheral *
+    dimuonMassCutFilterForOniaMMPeripheral
+    )
+

Configuration/Skimming/python/HI_OniaUPCSkim_cff.py

@@ -0,0 +1,13 @@
+import FWCore.ParameterSet.Config as cms
+
+# HLT dimuon trigger
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltOniaUPCHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltOniaUPCHI.HLTPaths = ["HLT_HIUPCDoubleMuNotHF2Pixel_SingleTrack_v*"]
+hltOniaUPCHI.throw = False
+hltOniaUPCHI.andOr = True
+
+# UPC double mu skim sequence
+oniaUPCSkimSequence = cms.Sequence(
+    hltOniaUPCHI 
+    )

Configuration/Skimming/python/HI_PhotonSkim_cff.py

@@ -0,0 +1,29 @@
+import FWCore.ParameterSet.Config as cms
+
+# HLT photon trigger
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltPhotonHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltPhotonHI.HLTPaths = ["HLT_HISinglePhoton50_Eta1p5_v*"]
+hltPhotonHI.throw = False
+hltPhotonHI.andOr = True
+
+# photon selection
+goodPhotons = cms.EDFilter("PhotonSelector",
+    src = cms.InputTag("photons"),
+    cut = cms.string('et > 60 && hadronicOverEm < 0.1 && r9 > 0.8 && sigmaIetaIeta > 0.002')
+)
+
+
+# leading photon E_T filter
+photonFilter = cms.EDFilter("EtMinPhotonCountFilter",
+    src = cms.InputTag("goodPhotons"),
+    etMin = cms.double(60.0),
+    minNumber = cms.uint32(1)
+)
+
+# photon skim sequence
+photonSkimSequence = cms.Sequence(hltPhotonHI
+                                  * goodPhotons
+                                  * photonFilter
+                                  )
+

Configuration/Skimming/python/HI_SingleTrackSkim_cff.py

@@ -0,0 +1,20 @@
+import FWCore.ParameterSet.Config as cms
+
+# HLT jet trigger
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltTrackHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltTrackHI.HLTPaths = ["HLT_HIFullTrack45_v*"]
+hltTrackHI.throw = False
+hltTrackHI.andOr = True
+
+# selection of valid vertex                                                                                                                                             
+primaryVertexFilterForTrack = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("hiSelectedVertex"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"),
+    filter = cms.bool(True),   # otherwise it won't filter the events                                                                                                    
+    )
+
+singleTrackSkimSequence = cms.Sequence(
+    primaryVertexFilterForTrack 
+    *hltTrackHI
+)

Configuration/Skimming/python/HI_ZEESkim_cff.py

@@ -0,0 +1,21 @@
+import FWCore.ParameterSet.Config as cms
+
+# HLT dimuon trigger
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltZEEHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltZEEHI.HLTPaths = ["HLT_HIDoublePhoton15_Eta1p5_Mass50_1000_R9HECut_v*"]
+hltZEEHI.throw = False
+hltZEEHI.andOr = True
+
+# selection of valid vertex
+primaryVertexFilterForZEE = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("hiSelectedVertex"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"), 
+    filter = cms.bool(True),   # otherwise it won't filter the events
+    )
+
+# Z->ee skim sequence
+zEESkimSequence = cms.Sequence(
+    hltZEEHI *
+    primaryVertexFilterForZEE
+    )

HeavyIonsAnalysis/Configuration/python/HI_ZMMSkim_cff.py → Configuration/Skimming/python/HI_ZMMSkim_cff.py

@@ -3,7 +3,7 @@
 # HLT dimuon trigger
 import HLTrigger.HLTfilters.hltHighLevel_cfi
 hltZMMHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
-hltZMMHI.HLTPaths = ["HLT_HIL2DoubleMu3"]
+hltZMMHI.HLTPaths = ["HLT_HIL1DoubleMu10_v*","HLT_HIL3Mu20_2HF_v*","HLT_HIL3Mu20_2HF0_v*"]
 hltZMMHI.throw = False
 hltZMMHI.andOr = True
 
@@ -15,34 +15,35 @@
     )
 
 # selection of dimuons (at least STA+STA) with mass in Z range
-muonSelector = cms.EDFilter("MuonSelector",
+muonSelectorForZMM = cms.EDFilter("MuonSelector",
     src = cms.InputTag("muons"),
-    cut = cms.string("(isStandAloneMuon || isGlobalMuon) && pt > 1."),
+    cut = cms.string("(isTrackerMuon && isGlobalMuon) && pt > 10."),
     filter = cms.bool(True)
     )
 
-muonFilter = cms.EDFilter("MuonCountFilter",
-    src = cms.InputTag("muonSelector"),
-    minNumber = cms.uint32(1)
+muonFilterForZMM = cms.EDFilter("MuonCountFilter",
+    src = cms.InputTag("muonSelectorForZMM"),
+    minNumber = cms.uint32(2)
     )
 
-dimuonMassCut = cms.EDProducer("CandViewShallowCloneCombiner",
+dimuonMassCutForZMM = cms.EDProducer("CandViewShallowCloneCombiner",
     checkCharge = cms.bool(True),
-    cut = cms.string(' 60 < mass < 120'),
-    decay = cms.string("muonSelector@+ muonSelector@-")
+    cut = cms.string(' 80 < mass < 110'),
+    decay = cms.string("muonSelectorForZMM@+ muonSelectorForZMM@-")
     )
 
-dimuonMassCutFilter = cms.EDFilter("CandViewCountFilter",
-    src = cms.InputTag("dimuonMassCut"),
+dimuonMassCutFilterForZMM = cms.EDFilter("CandViewCountFilter",
+    src = cms.InputTag("dimuonMassCutForZMM"),
     minNumber = cms.uint32(1)
     )
 
 # Z->mumu skim sequence
 zMMSkimSequence = cms.Sequence(
     hltZMMHI *
     primaryVertexFilterForZMM *
-    muonSelector *
-    muonFilter *
-    dimuonMassCut *
-    dimuonMassCutFilter
+    muonSelectorForZMM *
+    muonFilterForZMM *
+    dimuonMassCutForZMM *
+    dimuonMassCutFilterForZMM
     )
+

Configuration/Skimming/python/PP_D0MesonSkim_cff.py

@@ -0,0 +1,19 @@
+import FWCore.ParameterSet.Config as cms
+
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltDmeson60 = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltDmeson60.HLTPaths = ["HLT_DmesonPPTrackingGlobal_Dpt60_v*"]
+hltDmeson60.throw = False
+hltDmeson60.andOr = True
+
+# selection of valid vertex                                                                                                                                             
+primaryVertexFilterForD0Meson = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("offlinePrimaryVerticesWithBS"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"),
+    filter = cms.bool(True),   # otherwise it won't filter the events                                                                                         
+    )
+
+d0MesonSkimSequence = cms.Sequence(
+    primaryVertexFilterForD0Meson 
+    *hltDmeson60
+)

Configuration/Skimming/python/PP_HighPtJetSkim_cff.py

@@ -0,0 +1,20 @@
+import FWCore.ParameterSet.Config as cms
+
+# HLT jet trigger
+import HLTrigger.HLTfilters.hltHighLevel_cfi
+hltJet150 = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
+hltJet150.HLTPaths = ["HLT_AK4CaloJet150_Eta2p1_v*"]
+hltJet150.throw = False
+hltJet150.andOr = True
+
+# selection of valid vertex
+primaryVertexFilterForHighPtJets = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("offlinePrimaryVerticesWithBS"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"), 
+    filter = cms.bool(True),   # otherwise it won't filter the events
+    )
+
+highPtJetSkimSequence = cms.Sequence(
+        hltJet150*
+        primaryVertexFilterForHighPtJets 
+)