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]
