Add muon skim for 2025 light ion run

Author: flodamas

Configuration/Skimming/python/Ion_MuonSkim_cff.py

@@ -0,0 +1,59 @@
+import FWCore.ParameterSet.Config as cms
+
+# single high-pT muon skim sequence
+
+HighPtMuonSelection = "(isTrackerMuon || isGlobalMuon) && abs(eta) <= 2.4 && pt > 10."
+
+highPtMuonSelectorForMuonIon = cms.EDFilter("PATMuonRefSelector",
+    src = cms.InputTag("slimmedMuons"),
+    cut = cms.string(HighPtMuonSelection),
+    filter = cms.bool(True)
+    )
+
+highPtMuonCountFilterForMuonIon = cms.EDFilter("MuonRefPatCount",
+    src = cms.InputTag("slimmedMuons"),
+    cut = cms.string(HighPtMuonSelection),
+    minNumber = cms.uint32(1)
+    )
+
+HighPtMuonIonSkimSequence = cms.Sequence(
+    highPtMuonSelectorForMuonIon *
+    highPtMuonCountFilterForMuonIon
+    )
+
+
+# loose dimuon skim sequence
+
+LooseMuonSelection = "(isTrackerMuon || isGlobalMuon) && ((abs(eta) <= 1.0 && pt > 3.3) || (1.0 < abs(eta) <= 2.4 && pt > 1.0))"
+
+looseMuonSelectorForMuonIon = cms.EDFilter("PATMuonRefSelector",
+    src = cms.InputTag("slimmedMuons"),
+    cut = cms.string(LooseMuonSelection),
+    filter = cms.bool(True)
+    )
+
+looseMuonCountFilterForMuonIon = cms.EDFilter("MuonRefPatCount",
+    src = cms.InputTag("slimmedMuons"),
+    cut = cms.string(LooseMuonSelection),
+    minNumber = cms.uint32(2)
+    )
+
+
+dimuonMassCutForMuonIon = cms.EDProducer("CandViewShallowCloneCombiner",
+    checkCharge = cms.bool(True),
+    cut = cms.string(''),
+    decay = cms.string("looseMuonSelectorForMuonIon@+ looseMuonSelectorForMuonIon@-")
+    )
+
+dimuonCountFilterForMuonIon = cms.EDFilter("CandViewCountFilter",
+    src = cms.InputTag("dimuonMassCutForMuonIon"),
+    minNumber = cms.uint32(1)
+    )
+
+# dimuon skim sequence
+DimuonIonSkimSequence = cms.Sequence(
+    looseMuonSelectorForMuonIon *
+    looseMuonCountFilterForMuonIon *
+    dimuonMassCutForMuonIon *
+    dimuonCountFilterForMuonIon
+    )

Configuration/Skimming/python/Skims_Ion_cff.py

@@ -0,0 +1,31 @@
+import FWCore.ParameterSet.Config as cms
+
+from Configuration.EventContent.EventContent_cff import MINIAODEventContent
+
+skimMINIAODEventContent = MINIAODEventContent.clone()
+skimMINIAODEventContent.outputCommands.append("drop *_MEtoEDMConverter_*_*")
+skimMINIAODEventContent.outputCommands.append("drop *_*_*_SKIM")
+
+from Configuration.Skimming.Ion_MuonSkim_cff import *
+
+HighPtMuonIonSkimPath = cms.Path( HighPtMuonIonSkimSequence )
+SKIMStreamIonHighPtMuon = cms.FilteredStream(
+    responsible = 'HI PAG',
+    name = 'IonHighPtMuon',
+    paths = (HighPtMuonIonSkimPath),
+    content = skimMINIAODEventContent.outputCommands,
+    selectEvents = cms.untracked.PSet(),
+    dataTier = cms.untracked.string('USER')
+    )
+
+DimuonIonSkimPath = cms.Path( DimuonIonSkimSequence )
+SKIMStreamIonDimuon = cms.FilteredStream(
+    responsible = 'HI PAG',
+    name = 'IonDimuon',
+    paths = (DimuonIonSkimPath),
+    content = skimMINIAODEventContent.outputCommands,
+    selectEvents = cms.untracked.PSet(),
+    dataTier = cms.untracked.string('USER')
+    )
+
+#####################

Configuration/Skimming/python/autoSkim.py

@@ -40,6 +40,10 @@
  #'SingleMuon': 'LogError+LogErrorMonitor',
 }
 
+# For 2025 Ion skims
+for i_split in range(60):
+    autoSkim[f'IonPhysics{i_split}'] = 'IonHighPtMuon+IonDimuon+LogError+LogErrorMonitor'
+    
 # For 2024 PbPb skims
 for i_split in range(20):
     autoSkim[f'HIForward{i_split}'] = 'UPCMonopole+LogError+LogErrorMonitor'

Configuration/StandardSequences/python/Skims_cff.py

@@ -49,3 +49,4 @@ def getSkimDataTier(skimname):
 from Configuration.Skimming.Skims_PDWG_cff import *
 from Configuration.Skimming.Skims_PA_cff import *
 from Configuration.Skimming.Skims_PbPb_cff import *
+from Configuration.Skimming.Skims_Ion_cff import *