Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-03-17 10:57:31

0001 # Auto generated configuration file
0002 # using: 
0003 # Revision: 1.19 
0004 # Source: /local/reps/CMSSW/CMSSW/Configuration/Applications/python/ConfigBuilder.py,v 
0005 # with command line options: DiMuonVertex -s HARVESTING:alcaHarvesting --conditions auto:phase1_2022_realistic --mc --filetype DQM --scenario pp --era Run3 -n -1 --filein file:step3_inDQM.root
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 # import of standard configurations
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 # Input source
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 # Production Info
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 # Output definition
0092 
0093 # Additional output definition
0094 
0095 # Other statements
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 # the module to harvest
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 # trick to run the harvester module
0151 process.alcaHarvesting.insert(1,process.diMuonBiasClient)
0152 
0153 # Schedule definition
0154 process.schedule = cms.Schedule(process.alcaHarvesting,process.dqmsave_step)
0155 from PhysicsTools.PatAlgos.tools.helpers import associatePatAlgosToolsTask
0156 associatePatAlgosToolsTask(process)
0157 
0158 # Customisation from command line
0159 
0160 # Add early deletion of temporary data products to reduce peak memory need
0161 from Configuration.StandardSequences.earlyDeleteSettings_cff import customiseEarlyDelete
0162 process = customiseEarlyDelete(process)
0163 # End adding early deletion