def makeMyPrivateMCComponent(self, name, dataset, user, pattern, dbsInstance, xSec=1, useAAA=False): if self.__class__.useAAA is not None: useAAA = self.__class__.useAAA component = cfg.MCComponent( dataset=dataset, name=name, files=self.getMyFiles(dataset, user, pattern, dbsInstance, useAAA=useAAA), xSection=xSec, nGenEvents=1, triggers=[], effCorrFactor=1, ) component.splitFactor = 100 return component
def makePrivateMCComponent(self, name, dataset, files, xSec=1, prefix="auto"): if len(files) == 0: raise RuntimeError, "Trying to make a component %s with no files" % name dprefix = dataset + "/" if files[0][0] != "/" else "" if prefix == "auto": if (dprefix + files[0]).startswith("/store"): prefix = "root://eoscms.cern.ch//eos/cms" else: prefix = "" # prefix filenames with dataset unless they start with "/" component = cfg.MCComponent( dataset=dataset, name=name, files=[''.join([prefix, dprefix, f]) for f in files], xSection=xSec, nGenEvents=1, triggers=[], effCorrFactor=1, ) component.splitFactor = 100 return component
def makeMCComponent(self, name, dataset, user, pattern, xSec=1, useAAA=False, unsafe=False): component = cfg.MCComponent( dataset=dataset, name=name, files=self.getFiles(dataset, user, pattern, useAAA=useAAA, unsafe=unsafe), xSection=xSec, nGenEvents=1, triggers=[], effCorrFactor=1, ) component.dataset_entries = self.getPrimaryDatasetEntries( dataset, user, pattern, useAAA=useAAA) return component
def makeMCComponent(self, name, dataset, user, pattern, xSec=1, useAAA=False, unsafe=False, fracNegWeights=None, dbsInstance=None): if self.__class__.useAAA is not None: useAAA = self.__class__.useAAA component = cfg.MCComponent( dataset=dataset, name=name, files=self.getFiles(dataset, user, pattern, useAAA=useAAA, unsafe=unsafe, dbsInstance=dbsInstance), xSection=xSec, nGenEvents=1, triggers=[], effCorrFactor=1, ) component.splitFactor = 100 component.fracNegWeights = fracNegWeights component.dataset_entries = self.getPrimaryDatasetEntries( dataset, user, pattern, useAAA=useAAA) return component
def makePrivateMCComponentFromDir(self, name, dataset, filedir, pattern, xSec=1): import glob globpath = filedir + dataset + pattern files = glob.glob(globpath) if len(files) == 0: raise RuntimeError, "Trying to make a component %s with no files" % name # prefix filenames with dataset unless they start with "/" #dprefix = dataset +"/" if files[0][0] != "/" else "" component = cfg.MCComponent( dataset=dataset, name=name, files=files, #[filedir+'%s%s' % (dprefix,f) for f in files], xSection=xSec, nGenEvents=1, triggers=[], effCorrFactor=1, ) return component
def makeMCComponentFromEOS(self,name,dataset,path,pattern=".*root",xSec=1): component = cfg.MCComponent( dataset=dataset, name = name, files = self.getFilesFromEOS(name,dataset,path,pattern), xSection = xSec, nGenEvents = 1, triggers = [], effCorrFactor = 1, ) return component
def makeMyPrivateMCComponent(self,name,dataset,user,pattern,dbsInstance, xSec=1,useAAA=False): component = cfg.MCComponent( dataset=dataset, name = name, files = self.getMyFiles(dataset, user, pattern, dbsInstance, useAAA=useAAA), xSection = xSec, nGenEvents = 1, triggers = [], effCorrFactor = 1, ) return component
def makeMyPrivateMCComponentFromDPM(name, dataset, xSec=1): component = cfg.MCComponent( dataset=dataset, name = name, #files = [ 'root://hephyse.oeaw.ac.at/%s'%f for f in getFilesFromDAS(dataset, instance='phys03') ], files = getFilesOnDPM( dataset, pattern=".*root"), xSection = xSec, nGenEvents = 1, triggers = [], effCorrFactor = 1, ) return component
def makePrivateMCComponent(self,name,dataset,files,xSec=1): if len(files) == 0: raise RuntimeError, "Trying to make a component %s with no files" % name # prefix filenames with dataset unless they start with "/" dprefix = dataset +"/" if files[0][0] != "/" else "" component = cfg.MCComponent( dataset=dataset, name = name, files = ['root://eoscms.cern.ch//eos/cms%s%s' % (dprefix,f) for f in files], xSection = xSec, nGenEvents = 1, triggers = [], effCorrFactor = 1, ) return component
output_service2 = cfg.Service(TFileService, 'genonlyfile', name="genonlyfile", fname='treegen.root', option='recreate') # the following two lines are just for automatic testing # they are not needed for running on your own samples from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles testfiles = miniAodFiles() print "Running on test file %s" % testfiles sample = cfg.MCComponent( #specify the file you want to run on # files = ["/scratch/arizzi/Hbb/CMSSW_7_2_2_patch2/src/VHbbAnalysis/Heppy/test/ZLL-8A345C56-6665-E411-9C25-1CC1DE04DF20.root"], files=testfiles, name="SingleSample", isMC=True, isEmbed=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 selectedComponents = [sample] config = cfg.Config(components=selectedComponents, sequence=sequence, services=[output_service, output_service2], events_class=Events) # and the following runs the process directly if running as with python filename.py if __name__ == '__main__': from PhysicsTools.HeppyCore.framework.looper import Looper
comp.files = comp.files[:1] comp.splitFactor = 1 comp.fineSplitFactor = 4 elif test == '5': for comp in selectedComponents: comp.files = comp.files[:5] comp.splitFactor = 1 comp.fineSplitFactor = 5 elif test == '21': for comp in selectedComponents: comp.files = comp.files[:7] comp.splitFactor = 1 comp.fineSplitFactor = 3 elif test == "tau-sync": comp = cfg.MCComponent(files=[ "root://eoscms.cern.ch//store/mc/RunIISpring16MiniAODv2/TTWJetsToLNu_TuneCUETP8M1_13TeV-amcatnloFXFX-madspin-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_miniAODv2_v0-v1/50000/8E84F4BB-B620-E611-BBD8-B083FECFF2BF.root" ], name="TTW_Tau") comp.triggers = [] comp.splitFactor = 1 comp.fineSplitFactor = 6 selectedComponents = [comp] sequence.remove(jsonAna) ttHLepSkim.minLeptons = 0 elif test == '80X-MC': what = getHeppyOption("sample", "TTLep") if what == "TTLep": TTLep_pow = kreator.makeMCComponent( "TTLep_pow", "/TTTo2L2Nu_13TeV-powheg/RunIISpring16MiniAODv1-PUSpring16_80X_mcRun2_asymptotic_2016_v3_ext1-v1/MINIAODSIM", "CMS", ".*root", 831.76 * ((3 * 0.108)**2)) selectedComponents = [TTLep_pow]
VertexAna, LepAna, PhoAna, TauAna, JetAna, METAna, METPuppiAna, METNoHFAna, PdfAna, HighMassVHbb, treeProducer ] #,sh] from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles sample = cfg.MCComponent( files=[ #"root://xrootd.unl.edu//store/mc/RunIISpring15MiniAODv2/WprimeToWhToWlephbb_narrow_M-4000_13TeV-madgraph/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/60000/E4EB03F2-4170-E511-B5C4-00259073E466.root" "root://xrootd.unl.edu//store/mc/RunIISpring15MiniAODv2/WprimeToMuNu_M-3000_TuneCUETP8M1_13TeV-pythia8/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/60000/10EFFE5D-4374-E511-905E-003048F0113E.root" #"root://xrootd.unl.edu//store/mc/RunIISpring15MiniAODv2/WprimeToWhToWlephbb_narrow_M-1200_13TeV-madgraph/MINIAODSIM/Asympt50ns_74X_mcRun2_asymptotic50ns_v0-v1/30000/6248D38C-8D76-E511-B243-20CF305B058C.root" #'file:68790C43-4971-E511-AFD7-782BCB20E307.root' # "root://cms-xrd-global.cern.ch//store/mc/RunIISpring15MiniAODv2/DYJetsToLL_M-50_HT-100to200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/50000/009AE141-CA6D-E511-A060-002590A3716C.root" #'root://xrootd.unl.edu//store/mc/RunIISpring15MiniAODv2/WprimeToWhToWlephbb_narrow_M-1200_13TeV-madgraph/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/10000/DE7A5FB2-DC71-E511-9539-001517F8083C.root' #'file:DE7A5FB2-DC71-E511-9539-001517F8083C.root' #"root://xrootd.unl.edu//store/mc/RunIISpring15MiniAODv2/ZprimeToZhToZinvhbb_narrow_M-1200_13TeV-madgraph/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/30000/8C77EFF8-4771-E511-97BD-B083FED40671.root" #"/store/mc/RunIISpring15MiniAODv2/ZprimeToZhToZinvhbb_narrow_M-1200_13TeV-madgraph/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/30000/8C77EFF8-4771-E511-97BD-B083FED40671.root" # "/scratch/arizzi/00B6C8DE-E76E-E511-AEDE-008CFA000BB8.root" ##ttbar ], #files = ["226BB247-A565-E411-91CF-00266CFF0AF4.root"], name="ZHLL125", isEmbed=False, puFileMC="puMC.root", puFileData="puData.root", splitFactor=5) sample.isMC = True # the following is declared in case this cfg is used in input to the heppy.py script selectedComponents = [sample] from PhysicsTools.HeppyCore.framework.eventsfwlite import Events
comp.files = comp.files[:1] comp.splitFactor = 1 comp.fineSplitFactor = 1 elif test == '3': for comp in selectedComponents: comp.files = comp.files[:1] comp.splitFactor = 1 comp.fineSplitFactor = 4 elif test == '5': for comp in selectedComponents: comp.files = comp.files[:5] comp.splitFactor = 1 comp.fineSplitFactor = 5 elif test == "ewkinosync": comp = cfg.MCComponent(files=[ "root://eoscms.cern.ch//store/mc/RunIIFall15MiniAODv2/TTWJetsToLNu_TuneCUETP8M1_13TeV-amcatnloFXFX-madspin-pythia8/MINIAODSIM/PU25nsData2015v1_76X_mcRun2_asymptotic_v12-v1/60000/14C51DB0-D6B8-E511-8D9B-8CDCD4A9A484.root" ], name="TTW_EWK_sync") comp.triggers = [] comp.splitFactor = 1 comp.fineSplitFactor = 1 selectedComponents = [comp] sequence.remove(jsonAna) elif test == "ra5-sync-mc": comp = cfg.MCComponent(files=[ "root://eoscms.cern.ch//store/mc/RunIISpring16MiniAODv1/TTWJetsToLNu_TuneCUETP8M1_13TeV-amcatnloFXFX-madspin-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_v3-v1/00000/6E02CA07-BA02-E611-A59E-14187741208F.root" ], name="TTW_RA5_sync") comp.triggers = [] comp.splitFactor = 1 comp.fineSplitFactor = 1 selectedComponents = [comp]
prefetch=True) elif test == "prompt2015D": comp = cfg.DataComponent(files=[ "root://eoscms.cern.ch//store/data/Run2015D/DoubleEG/MINIAOD/PromptReco-v3/000/256/676/00000/E6FC61D6-C55F-E511-88D2-02163E011BE0.root" ], name="Prompt_2015D", intLumi=1) comp.triggers = [] comp.json = None selectedComponents = [comp] ttHLepSkim.minLeptons = 0 sequence.remove(jsonAna) if is50ns: raise RuntimeError, 'Incorrect is50ns configuration' elif test == 'miniAODv2': comp = cfg.MCComponent(files=[ "root://eoscms.cern.ch//store/mc/RunIISpring15MiniAODv2/TTJets_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/00000/0014DC94-DC5C-E511-82FB-7845C4FC39F5.root" ], name="TTJets_miniAODv2") comp.triggers = [] selectedComponents = [comp] ttHLepSkim.minLeptons = 0 sequence.remove(jsonAna) if is50ns: raise RuntimeError, 'Incorrect is50ns configuration' elif test == "ttH-sync": # TTHnobb_mWCutfix_ch0.files = ["/data/p/peruzzi/088378DB-3D24-E511-8B0E-20CF3027A589.root"] TTHnobb_mWCutfix_ch0.files = [ "root://eoscms.cern.ch//store/mc/RunIISpring15DR74/ttHJetToNonbb_M125_13TeV_amcatnloFXFX_madspin_pythia8_mWCutfix/MINIAODSIM/Asympt25ns_MCRUN2_74_V9-v1/30000/088378DB-3D24-E511-8B0E-20CF3027A589.root" ] selectedComponents = [TTHnobb_mWCutfix_ch0] for comp in selectedComponents: comp.splitFactor = 1 comp.fineSplitFactor = 30
import PhysicsTools.HeppyCore.framework.config as cfg from DMPD.Heppy.samples.Spring15.fileLists import mcsamples sampleDYJetsToLL_M50_HT100to200_madgraphMLM_pythia8_v2 = cfg.MCComponent( files=mcsamples["DYJetsToLL_M50_HT100to200_madgraphMLM_pythia8_v2"] ["files"], name="DYJetsToLL_M50_HT100to200_madgraphMLM_pythia8_v2", isMC=True, isEmbed=False, splitFactor=20, ) sampleTTbarDMJets_scalar_Mchi_1_Mphi_50_madgraphMLM_pythia8_v1 = cfg.MCComponent( files=mcsamples["TTbarDMJets_scalar_Mchi_1_Mphi_50_madgraphMLM_pythia8_v1"] ["files"], name="TTbarDMJets_scalar_Mchi_1_Mphi_50_madgraphMLM_pythia8_v1", isMC=True, isEmbed=False, splitFactor=5, ) sampleST_tW_top_5f_inclusiveDecays_powheg_pythia8_v1 = cfg.MCComponent( files=mcsamples["ST_tW_top_5f_inclusiveDecays_powheg_pythia8_v1"]["files"], name="ST_tW_top_5f_inclusiveDecays_powheg_pythia8_v1", isMC=True, isEmbed=False, splitFactor=5, ) sampleTTbarDMJets_scalar_Mchi_150_Mphi_500_madgraphMLM_pythia8_v3 = cfg.MCComponent( files=mcsamples[
import copy import re import PhysicsTools.HeppyCore.framework.config as cfg from CMGTools.RootTools.yellowreport.YRParser import yrparser13TeV # from CMGTools.H2TauTau.proto.samples.sampleShift import sampleShift HiggsVBF125 = cfg.MCComponent(name='HiggsVBF125', files=[], xSection=None, nGenEvents=0, triggers=[], effCorrFactor=1) ############# Gluon fusion ############### HiggsGGH125 = cfg.MCComponent(name='HiggsGGH125', files=[], xSection=None, nGenEvents=0, triggers=[], effCorrFactor=1) ############# VH ############### HiggsVH125 = cfg.MCComponent(name='HiggsVH125', files=[], xSection=None, nGenEvents=0, triggers=[], effCorrFactor=1)
from __future__ import print_function import os import PhysicsTools.HeppyCore.framework.config as cfg from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles # set to True if you want several parallel processes multi_thread = False # input component # several input components can be declared, # and added to the list of selected components inputSample = cfg.MCComponent( 'test_component', files=miniAodFiles(), # returns a list of file for this release # a list of local or xrootd files can be specified by hand. ) if 'RelValZMM' not in inputSample.files[0]: print('''WARNING: this tutorial is supposed to run on Z->mumu events. Do not expect meaningful results for this sample:''') print(inputSample) if multi_thread: inputSample.splitFactor = len(inputSample.files) selectedComponents = [inputSample] # a very simple muon analyzer # read miniaod muons and wrap them in python muons from PhysicsTools.Heppy.analyzers.examples.SimpleMuonAnalyzer import SimpleMuonAnalyzer muons = cfg.Analyzer(
files=datasamples[i]['files'], name=i, json= '%s/src/DMPD/Heppy/python/tools/JSON/Cert_246908-260627_13TeV_PromptReco_Collisions15_25ns_JSON.txt' % os.environ['CMSSW_BASE'], splitFactor=int(datasamples[i]['nevents'] / (maxlsftime * 3600 * eventspersec)), ) for i in mcsamples: if int(mcsamples[i]['nevents']) < 1: continue #print i, " - ", int(mcsamples[i]['nevents']/(maxlsftime*3600*eventspersec)) sample[i] = cfg.MCComponent( files=mcsamples[i]['files'], name=i, isMC=True, isEmbed=False, splitFactor=int(mcsamples[i]['nevents'] / (maxlsftime * 3600 * eventspersec)), ) testDataCompontent = cfg.Component( files=[ 'dcap://t2-srm-02.lnl.infn.it/pnfs/lnl.infn.it/data/cms//store/mc/RunIISpring15MiniAODv2/ST_t-channel_top_4f_leptonDecays_13TeV-powheg-pythia8_TuneCUETP8M1/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/50000/548FB0B9-D072-E511-84CA-0025908653C4.root', ], name="test", json= '%s/src/DMPD/Heppy/python/tools/JSON/Cert_246908-260627_13TeV_PromptReco_Collisions15_25ns_JSON.txt' % os.environ['CMSSW_BASE'], splitFactor=1, )
comp.files = comp.files[:1] dmCoreSequence.remove(jsonAna) elif test == 'simone': ##dmCoreSequence.remove(monoXSubJetPuppiAna) ##dmCoreSequence.remove(monoXSubJetSoftDropAna) from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles sample = cfg.MCComponent( files = [ #"root://cms-xrd-global.cern.ch//store/data/Run2016B/SingleElectron/RAW/v2/000/273/150/00000/A8A803BC-CC17-E611-BA8C-02163E014550.root", # "root://eoscms//eos/cms/store/data/Run2016B/DoubleMuon/MINIAOD/PromptReco-v2/000/273/448/00000/7ED72389-581C-E611-BF09-02163E011BF0.root" #"root://cms-xrd-global.cern.ch//store/data/Run2016B/SingleElectron/MINIAOD/PromptReco-v2/000/273/158/00000/06277EC1-181A-E611-870F-02163E0145E5.root" #"root://xrootd.unl.edu//store/mc/RunIISpring16DR80/TT_TuneCUETP8M1_13TeV-powheg-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_v3_ext3-v1/60000/08078977-2F1F-E611-AF79-001E675053A5.root" #"root://cms-xrd-global.cern.ch//store/mc/RunIISpring16MiniAODv2/TT_TuneCUETP8M1_13TeV-powheg-pythia8/MINIAODSIM/PUSpring16RAWAODSIM_80X_mcRun2_asymptotic_2016_miniAODv2_v0_ext3-v1/00000/00A10CC4-4227-E611-BBF1-C4346BBCD528.root" #"root://cms-xrd-global.cern.ch//store/mc/RunIISpring16MiniAODv1/TT_TuneCUETP8M1_13TeV-powheg-pythia8-evtgen/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_v3-v1/00000/4603CC0B-D012-E611-972B-90B11C06E1A0.root" #"root://cms-xrd-global.cern.ch//store/mc/RunIISpring16MiniAODv2/TTJets_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_miniAODv2_v0-v1/00000/04FB4BAA-3A33-E611-BC64-008CFA197A90.root", #"root://cms-xrd-global.cern.ch//store/mc/RunIISpring16MiniAODv2/TTGJets_TuneCUETP8M1_13TeV-amcatnloFXFX-madspin-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_miniAODv2_v0-v1/00000/9C6425C8-081C-E611-9F71-001E67E71A56.root", #"file:A8A803BC-CC17-E611-BA8C-02163E014550.root", "root://cms-xrd-global.cern.ch//store/mc/RunIISpring16MiniAODv2/GluGluSpin0ToZGamma_ZToQQ_W_0-p-014_M_1000_TuneCUEP8M1_13TeV_pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_miniAODv2_v0-v1/10000/3E743410-C82B-E611-8948-0025907FD40C.root", #"root://cms-xrd-global.cern.ch//store/data/Run2016B/SinglePhoton/MINIAOD/PromptReco-v2/000/273/158/00000/00DD3222-261A-E611-9FD2-02163E011E34.root", ], name="ZHLL125", isEmbed=False, puFileMC="puMC.root", puFileData="puData.root", splitFactor = 5 ) sample.isMC=True#False selectedComponents = [sample] elif test== 'simoneComponent': comp = kreator.makeMCComponent("ZGamma_Signal_1000TeV","/GluGluSpin0ToZGamma_ZToQQ_W_0-p-014_M_1000_TuneCUEP8M1_13TeV_pythia8/RunIISpring16MiniAODv2-PUSpring16_80X_mcRun2_asymptotic_2016_miniAODv2_v0-v1/MINIAODSIM", "CMS", ".*root", 1.) selectedComponents = [ comp ]
nEvents = 100 elif os.environ.get("SAMPLE", "tth") == "ttjets": fns = map(getSitePrefix, [ "/store/user/jpata/tth/Sep29_v1/TTToSemilepton_TuneCUETP8M2_ttHtranche3_13TeV-powheg-pythia8/Sep29_v1/161005_125253/0000/tree_1.root" ]) dataset = "TTToSemilepton_TuneCUETP8M2_ttHtranche3_13TeV-powheg-pythia8" firstEvent = 0 nEvents = 1000 print "files", fns config = cfg.Config( #Run across these inputs components=[ cfg.MCComponent( dataset, files=fns, tree_name="vhbb/tree", ) ], sequence=sequence, services=[output_service], events_class=Events, ) from PhysicsTools.HeppyCore.framework.looper import Looper looper = Looper( 'Loop', config, nPrint=0, firstEvent=firstEvent, nEvents=nEvents, )
import os #####COMPONENT CREATOR from CMGTools.RootTools.samples.ComponentCreator import ComponentCreator kreator = ComponentCreator() # WToTauTo3Mu = kreator.makeMCComponent("DoubleMuonLowMass_Run2016B_23Sep2016" , "/DoubleMuonLowMass/Run2016B-23Sep2016-v3/MINIAOD" , "CMS", ".*root", json) WToTauTo3Mu = cfg.MCComponent( dataset='WToTauTo3Mu', name='WToTauTo3Mu', files=[ '/afs/cern.ch/work/m/mverzett/public/perRic/t3mMINIAODSIM/t3mu_MINIAODSIM_0.root', '/afs/cern.ch/work/m/mverzett/public/perRic/t3mMINIAODSIM/t3mu_MINIAODSIM_1.root', '/afs/cern.ch/work/m/mverzett/public/perRic/t3mMINIAODSIM/t3mu_MINIAODSIM_2.root', '/afs/cern.ch/work/m/mverzett/public/perRic/t3mMINIAODSIM/t3mu_MINIAODSIM_3.root', '/afs/cern.ch/work/m/mverzett/public/perRic/t3mMINIAODSIM/t3mu_MINIAODSIM_4.root', '/afs/cern.ch/work/m/mverzett/public/perRic/t3mMINIAODSIM/t3mu_MINIAODSIM_5.root', '/afs/cern.ch/work/m/mverzett/public/perRic/t3mMINIAODSIM/t3mu_MINIAODSIM_6.root', '/afs/cern.ch/work/m/mverzett/public/perRic/t3mMINIAODSIM/t3mu_MINIAODSIM_7.root', '/afs/cern.ch/work/m/mverzett/public/perRic/t3mMINIAODSIM/t3mu_MINIAODSIM_8.root', '/afs/cern.ch/work/m/mverzett/public/perRic/t3mMINIAODSIM/t3mu_MINIAODSIM_9.root', ], xSection=20508.9 * 1.e-7, # W to lep nu / 3.[pb] x BR # xSection = 20508.9 * 0.005, # W to lep nu / 3.[pb] x BR nGenEvents=20000., triggers=['HLT_DoubleMu3_Trk_Tau3mu_v%d' % i for i in range(1, 5)], effCorrFactor=1, )
for comp in selectedComponents: comp.json = json comp.splitFactor = 7 comp.fineSplitFactor = 1 if getHeppyOption("single") else 8 if not getHeppyOption("all"): comp.files = comp.files[:1] dmCoreSequence.remove(jsonAna) elif test == 'simone': from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles sample = cfg.MCComponent( files=[ # "root://eoscms//eos/cms/store/data/Run2016B/DoubleMuon/MINIAOD/PromptReco-v2/000/273/448/00000/7ED72389-581C-E611-BF09-02163E011BF0.root" #"root://cms-xrd-global.cern.ch//store/data/Run2016B/SingleElectron/MINIAOD/PromptReco-v2/000/273/158/00000/06277EC1-181A-E611-870F-02163E0145E5.root" #"root://xrootd.unl.edu//store/mc/RunIISpring16DR80/TT_TuneCUETP8M1_13TeV-powheg-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_v3_ext3-v1/60000/08078977-2F1F-E611-AF79-001E675053A5.root" #"root://cms-xrd-global.cern.ch//store/mc/RunIISpring16MiniAODv2/TT_TuneCUETP8M1_13TeV-powheg-pythia8/MINIAODSIM/PUSpring16RAWAODSIM_80X_mcRun2_asymptotic_2016_miniAODv2_v0_ext3-v1/00000/00A10CC4-4227-E611-BBF1-C4346BBCD528.root" "root://cms-xrd-global.cern.ch//store/mc/RunIISpring16MiniAODv1/TT_TuneCUETP8M1_13TeV-powheg-pythia8-evtgen/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_v3-v1/00000/4603CC0B-D012-E611-972B-90B11C06E1A0.root" ], name="ZHLL125", isEmbed=False, puFileMC="puMC.root", puFileData="puData.root", splitFactor=5) sample.isMC = True selectedComponents = [sample] ## output histogram outputService = [] from PhysicsTools.HeppyCore.framework.services.tfile import TFileService output_service = cfg.Service(TFileService,
jsonAna, LHEAna, LHEWeightAna, FlagsAna, hbheAna, GenAna, VHGenAna, PUAna, TrigAna, VertexAna, LepAna, PhoAna, TauAna, JetAna, ttHLeptonMVA, METAna, METPuppiAna, PdfAna, VHbb, TTHtoTauTau, TTHtoTauTauGen, TriggerObjectsAna, treeProducer ] #,sh] from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles sample = cfg.MCComponent( files=[ #"root://xrootd.ba.infn.it//store/mc/RunIIFall15MiniAODv1/TT_TuneCUETP8M1_13TeV-powheg-scaledown-pythia8/MINIAODSIM/PU25nsData2015v1_76X_mcRun2_asymptotic_v12-v1/30000/045996FE-A19D-E511-B76D-D4AE526A0B47.root" ##ttbar #"045996FE-A19D-E511-B76D-D4AE526A0B47.root" ##ttbar #"/scratch/bianchi/06EC3B62-39B8-E511-9957-0002C94CD0C0.root" # "027D7153-29BF-E511-A2BC-0025B3E025B6.root" "root://xrootd.ba.infn.it//store/mc/RunIIFall15MiniAODv2/DYJetsToLL_M-100to200_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/MINIAODSIM/PU25nsData2015v1_76X_mcRun2_asymptotic_v12_ext1-v1/60000/027D7153-29BF-E511-A2BC-0025B3E025B6.root" ], #files = ["226BB247-A565-E411-91CF-00266CFF0AF4.root"], name="ZHLL125", isEmbed=False, puFileMC="puMC.root", puFileData="puData.root", puFileDataPlus="puDataPlus.root", puFileDataMinus="puDataMinus.root", splitFactor=5) sample.isMC = True # the following is declared in case this cfg is used in input to the heppy.py script selectedComponents = [sample] from PhysicsTools.HeppyCore.framework.eventsfwlite import Events config = cfg.Config(components=selectedComponents, sequence=sequence,
comp.files = comp.files[:1] comp.splitFactor = 1 comp.fineSplitFactor = 1 elif test == '3': for comp in selectedComponents: comp.files = comp.files[:1] comp.splitFactor = 1 comp.fineSplitFactor = 4 elif test == '5': for comp in selectedComponents: comp.files = comp.files[:5] comp.splitFactor = 1 comp.fineSplitFactor = 5 elif test == "ra5-sync-mc": comp = cfg.MCComponent(files=[ "root://eoscms.cern.ch//store/mc/RunIISpring16MiniAODv1/TTWJetsToLNu_TuneCUETP8M1_13TeV-amcatnloFXFX-madspin-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_v3-v1/00000/6E02CA07-BA02-E611-A59E-14187741208F.root" ], name="TTW_RA5_sync") comp.triggers = [] comp.splitFactor = 1 comp.fineSplitFactor = 1 selectedComponents = [comp] sequence.remove(jsonAna) elif test == '80X-MC': what = getHeppyOption("sample", "TTLep") if what == "TTLep": TTLep_pow = kreator.makeMCComponent( "TTLep_pow", "/TTTo2L2Nu_13TeV-powheg/RunIISpring16MiniAODv1-PUSpring16_80X_mcRun2_asymptotic_2016_v3_ext1-v1/MINIAODSIM", "CMS", ".*root", 831.76 * ((3 * 0.108)**2)) selectedComponents = [TTLep_pow] comp = selectedComponents[0]