# selectedComponents = selectedComponents[-1:] tbd = [ 'DY2JetsToLL_M50_LO', 'DY3JetsToLL_M50_LO', 'DY4JetsToLL_M50_LO', 'HiggsGGH125', 'WZTo3LNu_amcatnlo', ] selectedComponents = [s for s in selectedComponents if s.name in tbd] if doSUSY: from CMGTools.RootTools.samples.autoAAAconfig import autoAAA autoAAA(selectedComponents) preprocessor = None if cmssw: sequence.append(fileCleaner) cfg_name = "$CMSSW_BASE/src/CMGTools/H2TauTau/prod/h2TauTauMiniAOD_ditau_data_cfg.py" if data else "$CMSSW_BASE/src/CMGTools/H2TauTau/prod/h2TauTauMiniAOD_ditau_cfg.py" preprocessor = CmsswPreprocessor(cfg_name, addOrigAsSecondary=False) # the following is declared in case this cfg is used in input to the # heppy.py script from PhysicsTools.HeppyCore.framework.eventsfwlite import Events config = cfg.Config(components=selectedComponents, sequence=sequence, services=outputService, preprocessor=preprocessor,
for pd, triggers in DatasetsAndTriggers: for comp in byCompName(allData, [pd]): comp.triggers = triggers[:] comp.vetoTriggers = vetoTriggers[:] dataSamples.append(comp) vetoTriggers += triggers[:] selectedComponents = mcSamples + dataSamples if getHeppyOption('selectComponents'): if getHeppyOption('selectComponents')=='MC': selectedComponents = mcSamples elif getHeppyOption('selectComponents')=='DATA': selectedComponents = dataSamples else: selectedComponents = byCompName(selectedComponents, getHeppyOption('selectComponents').split(",")) autoAAA(selectedComponents, quiet=not(getHeppyOption("verboseAAA",False)), redirectorAAA="xrootd-cms.infn.it") # must be done before mergeExtensions configureSplittingFromTime(mcSamples,250 if preprocessor else 10,10) configureSplittingFromTime(dataSamples,80 if preprocessor else 10,10) selectedComponents, _ = mergeExtensions(selectedComponents, verbose=True) # create and set preprocessor if requested if getHeppyOption("nanoPreProcessor"): from CMGTools.Production.nanoAODPreprocessor import nanoAODPreprocessor preproc_cfg = {2016: ("mc94X2016","data94X2016"), 2017: ("mc94Xv2","data94Xv2"), 2018: ("mc102X","data102X_ABC","data102X_D")} preproc_cmsswArea = "/afs/cern.ch/work/v/vtavolar/SusySOSSW_2_clean/nanoAOD/clean/CMSSW_10_2_16_UL"#/afs/cern.ch/user/v/vtavolar/work/SusySOSSW_2_clean/nanoAOD/CMSSW_10_2_15" #MODIFY ACCORDINGLY preproc_mc = nanoAODPreprocessor(cfg='%s/src/PhysicsTools/NanoAOD/test/%s_NANO.py'%(preproc_cmsswArea,preproc_cfg[year][0]),cmsswArea=preproc_cmsswArea,keepOutput=True) if year==2018: preproc_data_ABC = nanoAODPreprocessor(cfg='%s/src/PhysicsTools/NanoAOD/test/%s_NANO.py'%(preproc_cmsswArea,preproc_cfg[year][1]),cmsswArea=preproc_cmsswArea,keepOutput=True, injectTriggerFilter=True, injectJSON=True) preproc_data_D = nanoAODPreprocessor(cfg='%s/src/PhysicsTools/NanoAOD/test/%s_NANO.py'%(preproc_cmsswArea,preproc_cfg[year][2]),cmsswArea=preproc_cmsswArea,keepOutput=True, injectTriggerFilter=True, injectJSON=True)
import PhysicsTools.HeppyCore.framework.config as cfg from PhysicsTools.HeppyCore.framework.config import printComps from PhysicsTools.HeppyCore.framework.heppy_loop import getHeppyOption from PhysicsTools.HeppyCore.framework.eventsfwlite import Events from CMGTools.H2TauTau.proto.analyzers.TauIsoAnalyzer import TauIsoAnalyzer from CMGTools.H2TauTau.proto.analyzers.TauIsoTreeProducer import TauIsoTreeProducer from CMGTools.H2TauTau.proto.samples.fall17.htt_common import DYJetsToLL_M50_LO, QCD_Pt_15to7000_TuneCP5_Flat2017, QCD_Pt_15to7000_TuneCP5_Flat from CMGTools.RootTools.samples.autoAAAconfig import autoAAA production = getHeppyOption('production', False) selected_components = [DYJetsToLL_M50_LO] selected_components = [QCD_Pt_15to7000_TuneCP5_Flat2017] selected_components = [QCD_Pt_15to7000_TuneCP5_Flat] autoAAA(selected_components) tau_iso_ana = cfg.Analyzer( TauIsoAnalyzer, name='TauIsoAnalyzer', ) tau_iso_tree = cfg.Analyzer( TauIsoTreeProducer, name='TauIsoTreeProducer', defaultFloatType='F', #save storage space reject_gen_leptons=True) sequence = cfg.Sequence([tau_iso_ana, tau_iso_tree]) if not production:
from CMGTools.RootTools.samples.samples_13TeV_RunIIFall17MiniAOD import samples as mcSamples_ from CMGTools.RootTools.samples.samples_13TeV_DATA2017 import dataSamples_31Mar2018 as allData elif year == 2016: from CMGTools.RootTools.samples.samples_13TeV_RunIISummer16MiniAODv3 import samples as mcSamples_ from CMGTools.RootTools.samples.samples_13TeV_DATA2016 import dataSamples_17Jul2018 as allData else: if year == 2018: from CMGTools.RootTools.samples.samples_13TeV_RunIIAutumn18NanoAODv4 import samples as mcSamples_ from CMGTools.RootTools.samples.samples_13TeV_DATA2018_NanoAOD import dataSamples_1June2019 as allData elif year == 2017: from CMGTools.RootTools.samples.samples_13TeV_RunIIFall17NanoAODv4 import samples as mcSamples_ from CMGTools.RootTools.samples.samples_13TeV_DATA2017_NanoAOD import dataSamples_1June2019 as allData elif year == 2016: from CMGTools.RootTools.samples.samples_13TeV_RunIISummer16NanoAODv4 import samples as mcSamples_ from CMGTools.RootTools.samples.samples_13TeV_DATA2016_NanoAOD import dataSamples_1June2019 as allData autoAAA(mcSamples_+allData, quiet=not(getHeppyOption("verboseAAA",False))) # must be done before mergeExtensions mcSamples_, _ = mergeExtensions(mcSamples_) # Triggers if year == 2018: from CMGTools.RootTools.samples.triggers_13TeV_DATA2018 import all_triggers as triggers elif year == 2017: from CMGTools.RootTools.samples.triggers_13TeV_DATA2017 import all_triggers as triggers triggers["FR_1mu_iso"] = [] # they probably existed but we didn't use them in 2017 elif year == 2016: from CMGTools.RootTools.samples.triggers_13TeV_DATA2016 import all_triggers as triggers triggers["FR_1mu_noiso_smpd"] = [] DatasetsAndTriggers = [] if analysis == "main":
# make data dataSamples = [] vetoTriggers = [] for pd, triggers in DatasetsAndTriggers: for comp in byCompName(allData, [pd]): comp.triggers = triggers[:] comp.vetoTriggers = vetoTriggers[:] dataSamples.append(comp) vetoTriggers += triggers[:] selectedComponents = mcSamples + dataSamples if getHeppyOption('selectComponents'): selectedComponents = byCompName( selectedComponents, getHeppyOption('selectComponents').split(",")) autoAAA(selectedComponents, quiet=not (getHeppyOption("verboseAAA", False))) configureSplittingFromTime(mcSamples, 250 if preprocessor else 10, 10) configureSplittingFromTime(dataSamples, 80 if preprocessor else 10, 10) selectedComponents, _ = mergeExtensions(selectedComponents) # create and set preprocessor if requested if getHeppyOption("nanoPreProcessor"): from CMGTools.Production.nanoAODPreprocessor import nanoAODPreprocessor preproc_cfg = { 2016: ("mc94X2016", "data94X2016"), 2017: ("mc94Xv2", "data94Xv2"), 2018: ("mc102X", "data102X_ABC", "data102X_D") } preproc_cmsswArea = "/afs/cern.ch/user/v/vtavolar/work/SusySOSSW_2_clean/nanoAOD/CMSSW_10_2_15" #MODIFY ACCORDINGLY preproc_mc = nanoAODPreprocessor( cfg='%s/src/PhysicsTools/NanoAOD/test/%s_NANO.py' %