File indexing completed on 2023-03-17 10:57:31
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 FailPath = cms.untracked.vstring(),
0056 IgnoreCompletely = cms.untracked.vstring(),
0057 Rethrow = cms.untracked.vstring('ProductNotFound'),
0058 SkipEvent = cms.untracked.vstring(),
0059 accelerators = cms.untracked.vstring('*'),
0060 allowUnscheduled = cms.obsolete.untracked.bool,
0061 canDeleteEarly = cms.untracked.vstring(),
0062 deleteNonConsumedUnscheduledModules = cms.untracked.bool(True),
0063 dumpOptions = cms.untracked.bool(False),
0064 emptyRunLumiMode = cms.obsolete.untracked.string,
0065 eventSetup = cms.untracked.PSet(
0066 forceNumberOfConcurrentIOVs = cms.untracked.PSet(
0067 allowAnyLabel_=cms.required.untracked.uint32
0068 ),
0069 numberOfConcurrentIOVs = cms.untracked.uint32(0)
0070 ),
0071 fileMode = cms.untracked.string('FULLMERGE'),
0072 forceEventSetupCacheClearOnNewRun = cms.untracked.bool(False),
0073 makeTriggerResults = cms.obsolete.untracked.bool,
0074 numberOfConcurrentLuminosityBlocks = cms.untracked.uint32(0),
0075 numberOfConcurrentRuns = cms.untracked.uint32(1),
0076 numberOfStreams = cms.untracked.uint32(0),
0077 numberOfThreads = cms.untracked.uint32(1),
0078 printDependencies = cms.untracked.bool(False),
0079 sizeOfStackForThreadsInKB = cms.optional.untracked.uint32,
0080 throwIfIllegalParameter = cms.untracked.bool(True),
0081 wantSummary = cms.untracked.bool(False)
0082 )
0083
0084
0085 process.configurationMetadata = cms.untracked.PSet(
0086 annotation = cms.untracked.string('DiMuonVertex nevts:-1'),
0087 name = cms.untracked.string('Applications'),
0088 version = cms.untracked.string('$Revision: 1.19 $')
0089 )
0090
0091
0092
0093
0094
0095
0096 from Configuration.AlCa.GlobalTag import GlobalTag
0097 process.GlobalTag = GlobalTag(process.GlobalTag, options.globalTag, '')
0098
0099 process.dqmsave_step = cms.Path(process.DQMSaver)
0100
0101 if (options.resonance == 'Z'):
0102 print('',30*"#",'\n # will harvest Z file \n',30*"#")
0103 _folderName = cms.string('AlCaReco/TkAlDiMuonAndVertex')
0104 _doBkgFit = cms.bool(True)
0105 _useRooCBShape = cms.bool(False)
0106 _useRooCMSShape = cms.bool(True)
0107 _fitPar = cms.PSet(mean_par = cms.vdouble(90.,80.,100.),
0108 width_par = cms.vdouble(2.49,2.48,2.50),
0109 sigma_par = cms.vdouble(1.5,0.,10.))
0110 elif (options.resonance == 'Jpsi'):
0111 print('',30*"#",'\n # will harvest J/psi file \n',30*"#")
0112 _folderName = cms.string('AlCaReco/TkAlJpsiMuMu')
0113 _doBkgFit = cms.bool(True)
0114 _useRooCBShape = cms.bool(True)
0115 _useRooCMSShape = cms.bool(True)
0116 _fitPar = cms.PSet(mean_par = cms.vdouble(3.09, 2.7, 3.4),
0117 width_par = cms.vdouble(1.0, 0.0, 5.0),
0118 sigma_par = cms.vdouble(0.01, 0.0, 5.0))
0119 elif (options.resonance == 'Upsilon'):
0120 print('',30*"#",'\n # will harvest Upsilon file \n',30*"#")
0121 _folderName = cms.string('AlCaReco/TkAlUpsilonMuMu')
0122 _doBkgFit = cms.bool(True)
0123 _useRooCBShape = cms.bool(True)
0124 _useRooCMSShape = cms.bool(False)
0125 _fitPar = cms.PSet(mean_par = cms.vdouble(9.46, 8.9, 9.9),
0126 width_par = cms.vdouble(1.0, 0.0, 5.0),
0127 sigma_par = cms.vdouble(1.0, 0.0, 5.0))
0128
0129
0130 process.DiMuonMassBiasClient = cms.EDProducer("DiMuonMassBiasClient",
0131 FolderName = _folderName,
0132 fitBackground = _doBkgFit,
0133 debugMode = cms.bool(True),
0134 fit_par = _fitPar,
0135 useRooCBShape = _useRooCBShape,
0136 useRooCMSShape = _useRooCMSShape,
0137 MEtoHarvest = cms.vstring(
0138 'DiMuMassVsMuMuPhi',
0139 'DiMuMassVsMuMuEta',
0140 'DiMuMassVsMuPlusPhi',
0141 'DiMuMassVsMuPlusEta',
0142 'DiMuMassVsMuMinusPhi',
0143 'DiMuMassVsMuMinusEta',
0144 'DiMuMassVsMuMuDeltaEta',
0145 'DiMuMassVsCosThetaCS'
0146 )
0147 )
0148
0149 process.diMuonBiasClient = cms.Sequence(process.DiMuonMassBiasClient)
0150
0151 process.alcaHarvesting.insert(1,process.diMuonBiasClient)
0152
0153
0154 process.schedule = cms.Schedule(process.alcaHarvesting,process.dqmsave_step)
0155 from PhysicsTools.PatAlgos.tools.helpers import associatePatAlgosToolsTask
0156 associatePatAlgosToolsTask(process)
0157
0158
0159
0160
0161 from Configuration.StandardSequences.earlyDeleteSettings_cff import customiseEarlyDelete
0162 process = customiseEarlyDelete(process)
0163