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