Merge pull request #46447 from flodamas/PbPbSkims2024_142X

Update the PbPb skims for 2024 data taking

Author: cmsbuild

Configuration/Skimming/python/PbPb_EMuSkim_cff.py

@@ -1,9 +1,9 @@
 import FWCore.ParameterSet.Config as cms
 
-# HLT dimuon trigger
+# Electron or muon HLT trigger
 import HLTrigger.HLTfilters.hltHighLevel_cfi
 hltEMuHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
-hltEMuHI.HLTPaths = ["HLT_HIEle*Gsf_v*","HLT_HIL3SingleMu*_v*"]
+hltEMuHI.HLTPaths = ["HLT_HIEle*Gsf_v*","HLT_HIL2SingleMu*_v*","HLT_HIL3SingleMu*_v*"]
 hltEMuHI.throw = False
 hltEMuHI.andOr = True
 
@@ -22,7 +22,7 @@
 
 muonSelectorForEMu = cms.EDFilter("MuonSelector",
                                   src = cms.InputTag("muons"),
-                                  cut = cms.string("isPFMuon && isGlobalMuon && pt > 20."),
+                                  cut = cms.string("isPFMuon && isGlobalMuon && pt > 20. && abs(eta)<2.4 && (globalTrack().normalizedChi2() < 10) && (globalTrack().hitPattern().numberOfValidMuonHits()>0) && (numberOfMatchedStations() > 1)"),
                                   filter = cms.bool(True)
     )
 

Configuration/Skimming/python/PbPb_ZEESkim_cff.py

@@ -4,7 +4,7 @@
 # HLT dimuon trigger
 import HLTrigger.HLTfilters.hltHighLevel_cfi
 hltZEEHI = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
-hltZEEHI.HLTPaths = ["HLT_HIEle*Gsf_v*"]
+hltZEEHI.HLTPaths = ["HLT_HIEle*Gsf_v*","HLT_HIDoubleEle*"]
 hltZEEHI.throw = False
 hltZEEHI.andOr = True
 

Configuration/Skimming/python/PbPb_ZMMSkim_cff.py

@@ -3,7 +3,7 @@
 # HLT dimuon trigger
 import HLTrigger.HLTfilters.hltHighLevel_cfi
 hltZMMPbPb = HLTrigger.HLTfilters.hltHighLevel_cfi.hltHighLevel.clone()
-hltZMMPbPb.HLTPaths = ["HLT_HIL3SingleMu*_v*"]
+hltZMMPbPb.HLTPaths = ["HLT_HIL2SingleMu*_v*","HLT_HIL3SingleMu*_v*"]
 hltZMMPbPb.throw = False
 hltZMMPbPb.andOr = True
 
@@ -17,7 +17,7 @@
 # selection of dimuons with mass in Z range
 muonSelectorForZMM = cms.EDFilter("MuonSelector",
     src = cms.InputTag("muons"),
-    cut = cms.string("(isPFMuon && isGlobalMuon) && pt > 20."),
+    cut = cms.string("(isPFMuon && isGlobalMuon) && pt > 15. && abs(eta) < 2.4"),
     filter = cms.bool(True)
     )
 

Configuration/Skimming/python/PbPb_ZMuSkimMuonDPG_cff.py

@@ -7,15 +7,22 @@
 
 PbPbZMuHLTFilter = copy.deepcopy(hltHighLevel)
 PbPbZMuHLTFilter.throw = cms.bool(False)
-PbPbZMuHLTFilter.HLTPaths = ["HLT_HIL3Mu*"]
+PbPbZMuHLTFilter.HLTPaths = ["HLT_HIL2SingleMu*","HLT_HIL3SingleMu*"]
+
+# selection of valid vertex
+primaryVertexFilterForPbPbZMuSkim = cms.EDFilter("VertexSelector",
+    src = cms.InputTag("offlinePrimaryVertices"),
+    cut = cms.string("!isFake && abs(z) <= 25 && position.Rho <= 2"), 
+    filter = cms.bool(True),
+    )
 
 ### Z -> MuMu candidates
 # Get muons of needed quality for Zs
 
 ###create a track collection with generic kinematic cuts
 looseMuonsForPbPbZMuSkim = cms.EDFilter("TrackSelector",
                              src = cms.InputTag("generalTracks"),
-                             cut = cms.string('pt > 10 &&  abs(eta)<2.4 &&  (charge!=0)'),
+                             cut = cms.string('pt > 10 && abs(eta)<2.4 && (charge!=0)'),
                              filter = cms.bool(True)                                
                              )
 
@@ -86,7 +93,7 @@
 
 tightMuonsForPbPbZMuSkim = cms.EDFilter("MuonSelector",
     src = cms.InputTag("muons"),
-    cut = cms.string("(isGlobalMuon) && pt > 25. &&  (abs(eta)<2.4) &&  (isPFMuon>0) && (globalTrack().normalizedChi2() < 10) && (globalTrack().hitPattern().numberOfValidMuonHits()>0)&& (numberOfMatchedStations() > 1)&& (innerTrack().hitPattern().numberOfValidPixelHits() > 0)&& (innerTrack().hitPattern().trackerLayersWithMeasurement() > 5) && ((isolationR03().sumPt/pt)<0.1)"),
+    cut = cms.string("(isGlobalMuon) && pt > 15. && (abs(eta)<2.4) && (isPFMuon>0) && (globalTrack().normalizedChi2() < 10) && (globalTrack().hitPattern().numberOfValidMuonHits()>0) && (numberOfMatchedStations() > 1) && (innerTrack().hitPattern().numberOfValidPixelHits() > 0) && (innerTrack().hitPattern().trackerLayersWithMeasurement() > 5)"),
     filter = cms.bool(True)
     )
 
@@ -96,8 +103,8 @@
 # build Z-> MuMu candidates
 dimuonsForPbPbZMuSkim = cms.EDProducer("CandViewShallowCloneCombiner",
                          checkCharge = cms.bool(False),
-                         cut = cms.string('(mass > 60) &&  (charge=0)'),       
-                         decay = cms.string("tightMuonsForPbPbZMuSkim looseIsoMuonsForPbPbZMuSkim")
+                         cut = cms.string('(mass > 60)'),       
+                         decay = cms.string("tightMuonsForPbPbZMuSkim looseMuonsForPbPbZMuSkim")
                          )                                    
 
 
@@ -111,11 +118,12 @@
 
 diMuonSelSeqForPbPbZMuSkim = cms.Sequence(
                             PbPbZMuHLTFilter *
+                            primaryVertexFilterForPbPbZMuSkim *
                             looseMuonsForPbPbZMuSkim *
-                            ConcretelooseMuonsForPbPbZMuSkim *
-                            tkIsoDepositTkForPbPbZMuSkim *
-                            allPatTracksForPbPbZMuSkim *
-                            looseIsoMuonsForPbPbZMuSkim * 
+                            #ConcretelooseMuonsForPbPbZMuSkim *
+                            #tkIsoDepositTkForPbPbZMuSkim *
+                            #allPatTracksForPbPbZMuSkim *
+                            #looseIsoMuonsForPbPbZMuSkim * 
                             tightMuonsForPbPbZMuSkim *
                             dimuonsForPbPbZMuSkim *
                             dimuonsFilterForPbPbZMuSkim 

Configuration/Skimming/python/autoSkim.py

@@ -45,8 +45,8 @@
     autoSkim[f'HIForward{i_split}'] = 'UPCMonopole+LogError+LogErrorMonitor'
 
 # For 2023 PbPb skims
-for i_split in range(32):
-    autoSkim[f'HIPhysicsRawPrime{i_split}'] = 'PbPbEMu+PbPbZEE+PbPbZMM+PbPbHighPtJets+LogError+LogErrorMonitor'
+for i_split in range(60):
+    autoSkim[f'HIPhysicsRawPrime{i_split}'] = 'PbPbEMu+PbPbZEE+PbPbZMu+PbPbHighPtJets+LogError+LogErrorMonitor'
 
 autoSkimRunII = {
  'BTagCSV' : 'LogError+LogErrorMonitor',