def config():
cfg = configtools.getConfig('mc', 2018, 'mm', bunchcrossing='25ns')
cfg["InputFiles"].set_input(
#ekppath='/storage/c/tberger/testfiles/skimming_output/mc/ZJet_DYJetsToLL_amcatnloFXFX-pythia8_RunIISummer16_testfile.root',
#bmspath='/storage/c/tberger/testfiles/skimming_output/mc/ZJet_DYJetsToLL_amcatnloFXFX-pythia8_RunIISummer16_localtestfile.root',
#ekppath='/storage/c/tberger/testfiles/skimming_output/mc/ZJet_DYJetsToLL_amcatnloFXFX-pythia8_RunIISummer16_inclPFcand_testfile.root',
#nafpath="/pnfs/desy.de/cms/tier2/store/user/tberger/Skimming/testfiles/Zll_DYJetsToLL_M-50_amcatnloFXFX-pythia8_RunIISummer16_testfile_noJTB.root",
#bmspath='/ceph/mschnepf/skimming/job_0_output.root'
bmspath='root://cmsxrootd-redirectors.gridka.de//pnfs/gridka.de/cms/disk-only/store/user/mschnepf/Skimming/ZJet_DYJetsToLL_M-50_TuneCP5_13TeV-amcatnloFXFX-pythia8/*.root'
)
cfg = configtools.expand(cfg,
['nocuts','leptoncuts','genleptoncuts', 'allleptoncuts','zjetcuts','genzjetcuts','allzjetcuts'],
['None','L1','L1L2L3'])#,'zcuts','genzcuts','allzcuts'
configtools.remove_quantities(cfg, ['jet1flavor','jet1rc'])
configtools.add_quantities(cfg, [ 'mu1IDSFWeight', 'mu2IDSFWeight', 'leptonIDSFWeight',
'mu1IsoSFWeight','mu2IsoSFWeight','leptonIsoSFWeight',
'mu1TrackingSFWeight','mu2TrackingSFWeight','leptonTrackingSFWeight',
'mu1TriggerSFWeight','mu2TriggerSFWeight','leptonTriggerSFWeight',
'jet1puidraw',
])
##### Add Producers: #####
cfg['Processors'] = [ 'producer:MuonTriggerMatchingProducer',
] + cfg['Processors'] + [
'producer:LeptonIDSFProducer',
'producer:LeptonIsoSFProducer',
# 'producer:LeptonTrackingSFProducer',
'producer:LeptonTriggerSFProducer',
'producer:LeptonSFProducer',
]
cfg['Processors'].insert(cfg['Processors'].index('producer:ValidMuonsProducer'), 'producer:MuonCorrectionsProducer',)
##### Specify input sources for Jets & Muons: #####
cfg['VertexSummary'] = 'offlinePrimaryVerticesSummary'
cfg['ValidMuonsInput'] = "corrected"
cfg['GenJets'] = 'ak4GenJetsNoNu' # JTB switched off
cfg['TaggedJets'] = 'ak4PFJetsCHS'
# cfg['Jec'] = os.path.join(configtools.getPath(), '../JECDatabase/textFiles/Summer16_07Aug2017_V11_MC/Summer16_07Aug2017_V11_MC') ## need to be updated when available
##### Change selection: (see also http://cms.cern.ch/iCMS/analysisadmin/cadilines?line=SMP-17-002&tp=an&id=1891&ancode=SMP-17-002) #####
cfg['MuonIso'] = 'loose_2018'
cfg['MuonID'] = 'tight'
cfg['CutMuonPtMin'] = 25.0
cfg['CutMuonEtaMax'] = 2.4
cfg['ZMassRange'] = 20.0
cfg['CutLeadingJetPtMin'] = 15.0
cfg['MinPUJetID'] = -9999
cfg['HltPaths'] = ['HLT_IsoMu24', 'HLT_IsoTkMu24']
cfg["MuonTriggerFilterNames"] = ['HLT_IsoMu24_v2:hltL3crIsoL1sMu22L1f0L2f10QL3f24QL3trkIsoFiltered0p09','HLT_IsoTkMu24_v3:hltL3fL1sMu22L1f0Tkf24QL3trkIsoFiltered0p09']
##### LeptonSF files: #####
## https://twiki.cern.ch/twiki/bin/view/CMS/MuonPOG ##
cfg['LeptonIDSFRootfile'] = os.path.join(configtools.getPath(),"data/scalefactors/2017/ID_EfficienciesAndSF_BCDEF_Nov17Nov2017.root")
cfg['LeptonIsoSFRootfile'] = os.path.join(configtools.getPath(),"data/scalefactors/2017/Iso_EfficienciesAndSF_BCDEF_Nov17Nov2017.root")
cfg['LeptonTriggerSFRootfile'] = os.path.join(configtools.getPath(),"data/scalefactors/2017/Trigger_EfficienciesAndSF_BCDEF_Nov17Nov2017.root")
# cfg['LeptonTrackingSFRootfile'] = os.path.join(configtools.getPath(),"data/scalefactors/2016latest/Tracking_EfficienciesAndSF_BCDEFGH.root")
##### MC specific properties: #####
cfg['NumberGeneratedEvents'] = 122055388 # from geteventsscript
cfg['GeneratorWeight'] = 0.670123731536
cfg['CrossSection'] = 1921.8*3
cfg['PileupWeightFile'] = os.path.join(configtools.getPath() , 'data/pileup/PUWeights_BCDEFGH_13TeV_23Sep2016ReReco_DYJetsToLL_M-50_amcatnloFXFX-pythia8_RunIISummer16.root')
return cfg
def modify_config(cfg):
"""
Add quantities from leading skim jets
"""
configtools.add_quantities(
cfg, ("skimjet1pt", "skimjet1phi", "skimjet1eta", "skimjet1validity"))
return cfg
def config():
cfg = configtools.getConfig('data', 2016, 'mm', JEC=JEC, IOV=RUN)
cfg["InputFiles"].set_input(
#path = '/portal/ekpbms1/home/tberger/Excalibur/CMSSW_9_4_12/src/Kappa/Skimming/zjet/output_numEvent100.root',
pathF ='root://cmsxrootd-redirectors.gridka.de//store/user/tberger/Skimming_94X/SingleMuon_Run2016F-17Jul2018-v1/*.root',
pathG ='root://cmsxrootd-redirectors.gridka.de//store/user/tberger/Skimming_94X/SingleMuon_Run2016G-17Jul2018-v1/*.root',
pathH ='root://cmsxrootd-redirectors.gridka.de//store/user/tberger/Skimming_94X/SingleMuon_Run2016H-17Jul2018-v1/*.root',
)
cfg = configtools.expand(cfg,
['nocuts','zjetcuts'],
['None','L1L2L3Res'])
configtools.remove_quantities(cfg, [
'jet1area','jet1l1', 'jet1rc', 'jet1l2','jet1ptraw', 'jet1ptl1','jet1ptl1l2l3', 'jet1res',
#'jet2pt', 'jet2eta', 'jet2y', 'jet2phi',
#'jet3pt', 'jet3eta', 'jet3y', 'jet3phi',
'mpf', 'rawmpf', 'met', 'metphi', 'rawmet', 'rawmetphi', 'sumet','mettype1vecpt', 'mettype1pt',
])
configtools.add_quantities(cfg, ['lepton'+x+'SFWeight'+y for x in ['ID','Iso','Trigger'] for y in ['Up','Down']])
##### Add Producers: #####
cfg['Processors'].remove("producer:JetEtaPhiCleaner")
cfg['Processors'] = ['producer:MuonTriggerMatchingProducer'] + cfg['Processors']
cfg['Processors'].insert(cfg['Processors'].index('producer:ValidMuonsProducer'), 'producer:MuonCorrectionsProducer',)
cfg['Processors'].insert(cfg['Processors'].index('producer:ValidMuonsProducer'), 'producer:PFCandidatesProducer',)
cfg['Processors'].insert(cfg['Processors'].index('producer:ValidMuonsProducer')+1, 'producer:ZJetDressedMuonsProducer',)
##### Specify input sources for Jets & Muons: #####
cfg['PackedPFCandidates'] = 'pfCandidates'
cfg['MaxZJetDressedMuonDeltaR'] = 0.1
cfg['VertexSummary'] = 'offlinePrimaryVerticesSummary'
cfg['ValidMuonsInput'] = "corrected"
cfg['UseObjectJetYCut'] = True
cfg['TaggedJets'] = 'ak4PFJetsCHS'
cfg['ProvideJECUncertainties'] = True
cfg['JetEnergyCorrectionUncertaintyShift'] = -1
cfg['JetID'] = 'loose'
cfg['PUJetID'] = 'medium'
##### Change selection: (see also http://cms.cern.ch/iCMS/analysisadmin/cadilines?line=SMP-17-002&tp=an&id=1891&ancode=SMP-17-002) #####
cfg['MuonIso'] = 'loose_2016'
cfg['MuonID'] = 'tight'
cfg['CutMuonPtMin'] = 25.0
cfg['CutMuonEtaMax'] = 2.4
cfg['ZMassRange'] = 20.0
cfg['CutLeadingJetPtMin'] = 10.0
cfg['MinPUJetID'] = -0.2
cfg['HltPaths'] = ['HLT_IsoMu24', 'HLT_IsoTkMu24']
cfg["MuonTriggerFilterNames"] = ['HLT_IsoMu24_v2:hltL3crIsoL1sMu22L1f0L2f10QL3f24QL3trkIsoFiltered0p09','HLT_IsoTkMu24_v3:hltL3fL1sMu22L1f0Tkf24QL3trkIsoFiltered0p09']
##### LeptonSF files: #####
cfg['LeptonSFVariation'] = True
cfg['LeptonIDSFRootfile'] = os.path.join(configtools.getPath(),"data/scalefactors/2016/RunGH_SF_ID.root")
cfg['LeptonIsoSFRootfile'] = os.path.join(configtools.getPath(),"data/scalefactors/2016/RunGH_SF_ISO.root")
cfg['LeptonTriggerSFRootfile'] = os.path.join(configtools.getPath(),"data/scalefactors/2016/RunGH_SF_Trigger.root")
##### Json & JEC #####
cfg['JsonFiles'] = [os.path.join(configtools.getPath(),'data/json/Cert_GH_13TeV_23Sep2016ReReco_Collisions16_JSON.txt')]
return cfg
def config():
cfg = configtools.getConfig('mc', 2016, 'mm', JEC=JEC)
cfg["InputFiles"].set_input(
path=
"root://cmsxrootd-redirectors.gridka.de//store/user/tberger/Skimming_94X/DYJetsToLL_madgraphMLM-herwigpp_RunIISummer16/*.root",
)
cfg = configtools.expand(cfg, ['nocuts', 'zjetcuts', 'genzjetcuts'],
['None', 'L1L2L3'])
configtools.remove_quantities(
cfg,
[
'jet1area',
'jet1l1',
'jet1rc',
'jet1l2',
'jet1ptraw',
'jet1ptl1',
#'jet2pt', 'jet2eta', 'jet2y', 'jet2phi',
#'jet3pt', 'jet3eta', 'jet3y', 'jet3phi',
'mpf',
'rawmpf',
'met',
'metphi',
'rawmet',
'rawmetphi',
'sumet',
'mettype1vecpt',
'mettype1pt',
#'genjet2pt','genjet2eta','genjet2y','genjet2phi',
#'genjet3pt','genjet3eta','genjet3y','genjet3phi',
])
configtools.add_quantities(cfg, [
'parton1flavour',
'parton1pt',
'parton1y',
'parton1phi',
'parton1mass',
'parton2flavour',
'parton2pt',
'parton2y',
'parton2phi',
'parton2mass',
])
##### Add Producers: #####
cfg['Processors'] = [
'producer:MuonTriggerMatchingProducer',
] + cfg['Processors']
cfg['Processors'].insert(
cfg['Processors'].index('producer:ValidMuonsProducer'),
'producer:MuonCorrectionsProducer',
)
cfg['Processors'].insert(
cfg['Processors'].index('producer:ValidMuonsProducer'),
'producer:PFCandidatesProducer',
)
cfg['Processors'].insert(
cfg['Processors'].index('producer:ValidMuonsProducer') + 1,
'producer:ZJetDressedMuonsProducer',
)
cfg['Processors'].insert(
cfg['Processors'].index('producer:GenZmmProducer'),
'producer:ZJetGenPhotonsProducer',
)
cfg['Processors'].insert(
cfg['Processors'].index('producer:GenZmmProducer'),
'producer:ZJetDressedGenMuonsProducer',
)
#cfg['Processors'].insert(cfg['Processors'].index('producer:GenZmmProducer'), 'producer:ZJetTrueGenMuonsProducer',)
##### Specify input sources for Jets & Muons: #####
cfg['PackedPFCandidates'] = 'pfCandidates'
cfg['MaxZJetDressedMuonDeltaR'] = 0.1
cfg['VertexSummary'] = 'offlinePrimaryVerticesSummary'
cfg['ValidMuonsInput'] = "corrected"
cfg['GenJets'] = 'ak4GenJets'
cfg['TaggedJets'] = 'ak4PFJetsCHS'
cfg['UseObjectJetYCut'] = True
cfg['Jec'] = os.path.join(
configtools.getPath(),
'../JECDatabase/textFiles/Summer16_07Aug2017_V11_MC/Summer16_07Aug2017_V11_MC'
)
cfg['JetID'] = 'loose'
cfg['PUJetID'] = 'medium'
cfg['DeltaRMatchingRecoJetGenJet'] = 0.2
##### Change selection: (see also http://cms.cern.ch/iCMS/analysisadmin/cadilines?line=SMP-17-002&tp=an&id=1891&ancode=SMP-17-002) #####
cfg['MuonIso'] = 'loose_2016'
cfg['MuonID'] = 'tight'
cfg['CutMuonPtMin'] = 25.0
cfg['CutMuonEtaMax'] = 2.4
cfg['ZMassRange'] = 20.0
cfg['CutLeadingJetPtMin'] = 10.0
cfg['MinPUJetID'] = -0.2
cfg['HltPaths'] = ['HLT_IsoMu24', 'HLT_IsoTkMu24']
cfg["MuonTriggerFilterNames"] = [
'HLT_IsoMu24_v2:hltL3crIsoL1sMu22L1f0L2f10QL3f24QL3trkIsoFiltered0p09',
'HLT_IsoTkMu24_v3:hltL3fL1sMu22L1f0Tkf24QL3trkIsoFiltered0p09'
]
##### MC specific properties: #####
cfg['NumberGeneratedEvents'] = 19554896
cfg['GeneratorWeight'] = 1.0
cfg['CrossSection'] = 358.6
cfg['PileupWeightFile'] = os.path.join(
configtools.getPath(),
'data/pileup/PUWeights_BCDEFGH_13TeV_23Sep2016ReReco_DYJetsToLL_M-50_madgraphMLM-herwigpp_RunIISummer16.root'
)
return cfg
