File indexing completed on 2024-04-06 12:09:15
0001
0002
0003
0004
0005
0006 import FWCore.ParameterSet.Config as cms
0007
0008 from Configuration.Eras.Era_Run3_cff import Run3
0009
0010 process = cms.Process('HARVESTING',Run3)
0011
0012 import FWCore.ParameterSet.VarParsing as VarParsing
0013
0014 options = VarParsing.VarParsing()
0015 options.register('maxEvents',
0016 -1,
0017 VarParsing.VarParsing.multiplicity.singleton,
0018 VarParsing.VarParsing.varType.int,
0019 "maximum events")
0020 options.register('globalTag',
0021 'auto:phase1_2022_realistic',
0022 VarParsing.VarParsing.multiplicity.singleton,
0023 VarParsing.VarParsing.varType.string,
0024 "conditions")
0025 options.register('resonance',
0026 'Z',
0027 VarParsing.VarParsing.multiplicity.singleton,
0028 VarParsing.VarParsing.varType.string,
0029 "resonance type")
0030 options.parseArguments()
0031
0032
0033 process.load('Configuration.StandardSequences.Services_cff')
0034 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0035 process.load('FWCore.MessageService.MessageLogger_cfi')
0036 process.load('Configuration.EventContent.EventContent_cff')
0037 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0038 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0039 process.load('Configuration.StandardSequences.MagneticField_cff')
0040 process.load('Configuration.StandardSequences.DQMSaverAtRunEnd_cff')
0041 process.load('Configuration.StandardSequences.Harvesting_cff')
0042 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0043
0044 process.maxEvents = cms.untracked.PSet(
0045 input = cms.untracked.int32(options.maxEvents),
0046 output = cms.optional.untracked.allowed(cms.int32,cms.PSet)
0047 )
0048
0049
0050 process.source = cms.Source("DQMRootSource",
0051 fileNames = cms.untracked.vstring('file:step3_inDQM_'+options.resonance+'.root')
0052 )
0053
0054 process.options = cms.untracked.PSet(
0055 IgnoreCompletely = cms.untracked.vstring(),
0056 Rethrow = cms.untracked.vstring('ProductNotFound'),
0057 accelerators = cms.untracked.vstring('*'),
0058 allowUnscheduled = cms.obsolete.untracked.bool,
0059 canDeleteEarly = cms.untracked.vstring(),
0060 deleteNonConsumedUnscheduledModules = cms.untracked.bool(True),
0061 dumpOptions = cms.untracked.bool(False),
0062 emptyRunLumiMode = cms.obsolete.untracked.string,
0063 eventSetup = cms.untracked.PSet(
0064 forceNumberOfConcurrentIOVs = cms.untracked.PSet(
0065 allowAnyLabel_=cms.required.untracked.uint32
0066 ),
0067 numberOfConcurrentIOVs = cms.untracked.uint32(0)
0068 ),
0069 fileMode = cms.untracked.string('FULLMERGE'),
0070 forceEventSetupCacheClearOnNewRun = cms.untracked.bool(False),
0071 makeTriggerResults = cms.obsolete.untracked.bool,
0072 numberOfConcurrentLuminosityBlocks = cms.untracked.uint32(0),
0073 numberOfConcurrentRuns = cms.untracked.uint32(1),
0074 numberOfStreams = cms.untracked.uint32(0),
0075 numberOfThreads = cms.untracked.uint32(1),
0076 printDependencies = cms.untracked.bool(False),
0077 sizeOfStackForThreadsInKB = cms.optional.untracked.uint32,
0078 throwIfIllegalParameter = cms.untracked.bool(True),
0079 wantSummary = cms.untracked.bool(False)
0080 )
0081
0082
0083 process.configurationMetadata = cms.untracked.PSet(
0084 annotation = cms.untracked.string('DiMuonVertex nevts:-1'),
0085 name = cms.untracked.string('Applications'),
0086 version = cms.untracked.string('$Revision: 1.19 $')
0087 )
0088
0089
0090
0091
0092
0093
0094 from Configuration.AlCa.GlobalTag import GlobalTag
0095 process.GlobalTag = GlobalTag(process.GlobalTag, options.globalTag, '')
0096
0097 process.dqmsave_step = cms.Path(process.DQMSaver)
0098
0099 if (options.resonance == 'Z'):
0100 print('',30*"#",'\n # will harvest Z file \n',30*"#")
0101 _folderName = cms.string('AlCaReco/TkAlDiMuonAndVertex')
0102 _doBkgFit = cms.bool(True)
0103 _useRooCBShape = cms.bool(False)
0104 _useRooCMSShape = cms.bool(True)
0105 _fitPar = cms.PSet(mean_par = cms.vdouble(90.,80.,100.),
0106 width_par = cms.vdouble(2.49,2.48,2.50),
0107 sigma_par = cms.vdouble(1.5,0.,10.))
0108 elif (options.resonance == 'Jpsi'):
0109 print('',30*"#",'\n # will harvest J/psi file \n',30*"#")
0110 _folderName = cms.string('AlCaReco/TkAlJpsiMuMu')
0111 _doBkgFit = cms.bool(True)
0112 _useRooCBShape = cms.bool(True)
0113 _useRooCMSShape = cms.bool(True)
0114 _fitPar = cms.PSet(mean_par = cms.vdouble(3.09, 2.7, 3.4),
0115 width_par = cms.vdouble(1.0, 0.0, 5.0),
0116 sigma_par = cms.vdouble(0.01, 0.0, 5.0))
0117 elif (options.resonance == 'Upsilon'):
0118 print('',30*"#",'\n # will harvest Upsilon file \n',30*"#")
0119 _folderName = cms.string('AlCaReco/TkAlUpsilonMuMu')
0120 _doBkgFit = cms.bool(True)
0121 _useRooCBShape = cms.bool(True)
0122 _useRooCMSShape = cms.bool(False)
0123 _fitPar = cms.PSet(mean_par = cms.vdouble(9.46, 8.9, 9.9),
0124 width_par = cms.vdouble(1.0, 0.0, 5.0),
0125 sigma_par = cms.vdouble(1.0, 0.0, 5.0))
0126
0127
0128 process.DiMuonMassBiasClient = cms.EDProducer("DiMuonMassBiasClient",
0129 FolderName = _folderName,
0130 fitBackground = _doBkgFit,
0131 debugMode = cms.bool(True),
0132 fit_par = _fitPar,
0133 useRooCBShape = _useRooCBShape,
0134 useRooCMSShape = _useRooCMSShape,
0135 MEtoHarvest = cms.vstring(
0136 'DiMuMassVsMuMuPhi',
0137 'DiMuMassVsMuMuEta',
0138 'DiMuMassVsMuPlusPhi',
0139 'DiMuMassVsMuPlusEta',
0140 'DiMuMassVsMuMinusPhi',
0141 'DiMuMassVsMuMinusEta',
0142 'DiMuMassVsMuMuDeltaEta',
0143 'DiMuMassVsCosThetaCS'
0144 )
0145 )
0146
0147 process.diMuonBiasClient = cms.Sequence(process.DiMuonMassBiasClient)
0148
0149 process.alcaHarvesting.insert(1,process.diMuonBiasClient)
0150
0151
0152 process.schedule = cms.Schedule(process.alcaHarvesting,process.dqmsave_step)
0153 from PhysicsTools.PatAlgos.tools.helpers import associatePatAlgosToolsTask
0154 associatePatAlgosToolsTask(process)
0155
0156
0157
0158
0159 from Configuration.StandardSequences.earlyDeleteSettings_cff import customiseEarlyDelete
0160 process = customiseEarlyDelete(process)
0161