-
Notifications
You must be signed in to change notification settings - Fork 7
/
haddoutput2017.py
140 lines (116 loc) · 6.64 KB
/
haddoutput2017.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
import os,sys,datetime,time, subprocess, math
from ROOT import *
execfile("/uscms_data/d3/jmanagan/EOSSafeUtils.py")
start_time = time.time()
#IO directories must be full paths
input = sys.argv[1]
output = sys.argv[2]
inDir='/store/user/escharni/FWLJMET102X_1lep2017Dnn_Mar2020_step1/'
outDir='/store/user/escharni/FWLJMET102X_1lep2017Dnn_Mar2020_step1hadds/'
scratchDir='/uscmst1b_scratch/lpc1/3DayLifetime/escharni/'
if not os.path.exists(scratchDir): os.system('mkdir '+scratchDir)
os.system('eos root://cmseos.fnal.gov/ mkdir -p '+outDir)
dirList = [
'BprimeBprime_M-1000_TuneCP5_13TeV-madgraph-pythia8',
'BprimeBprime_M-1100_TuneCP5_13TeV-madgraph-pythia8',
'BprimeBprime_M-1200_TuneCP5_13TeV-madgraph-pythia8',
'BprimeBprime_M-1300_TuneCP5_13TeV-madgraph-pythia8',
'BprimeBprime_M-1400_TuneCP5_13TeV-madgraph-pythia8',
'BprimeBprime_M-1500_TuneCP5_13TeV-madgraph-pythia8',
'BprimeBprime_M-1600_TuneCP5_13TeV-madgraph-pythia8',
'BprimeBprime_M-1700_TuneCP5_13TeV-madgraph-pythia8',
'BprimeBprime_M-1800_TuneCP5_13TeV-madgraph-pythia8',
'BprimeBprime_M-900_TuneCP5_13TeV-madgraph-pythia8',
'DYJetsToLL_M-50_HT-1200to2500_TuneCP5_13TeV-madgraphMLM-pythia8',
'DYJetsToLL_M-50_HT-2500toInf_TuneCP5_13TeV-madgraphMLM-pythia8',
'DYJetsToLL_M-50_HT-400to600_TuneCP5_13TeV-madgraphMLM-pythia8',
'DYJetsToLL_M-50_HT-600to800_TuneCP5_13TeV-madgraphMLM-pythia8',
'DYJetsToLL_M-50_HT-800to1200_TuneCP5_13TeV-madgraphMLM-pythia8',
'QCD_HT1000to1500_TuneCP5_13TeV-madgraph-pythia8',
'QCD_HT1500to2000_TuneCP5_13TeV-madgraph-pythia8',
'QCD_HT2000toInf_TuneCP5_13TeV-madgraph-pythia8',
'QCD_HT300to500_TuneCP5_13TeV-madgraph-pythia8',
'QCD_HT500to700_TuneCP5_13TeV-madgraph-pythia8',
'QCD_HT700to1000_TuneCP5_13TeV-madgraph-pythia8',
'ST_s-channel_antitop_leptonDecays_13TeV-PSweights_powheg-pythia',
'ST_s-channel_top_leptonDecays_13TeV-PSweights_powheg-pythia',
'ST_t-channel_antitop_4f_InclusiveDecays_TuneCP5_PSweights_13TeV-powheg-pythia8',
'ST_t-channel_top_4f_InclusiveDecays_TuneCP5_PSweights_13TeV-powheg-pythia8',
'ST_tW_antitop_5f_inclusiveDecays_TuneCP5_PSweights_13TeV-powheg-pythia8',
'ST_tW_top_5f_inclusiveDecays_TuneCP5_PSweights_13TeV-powheg-pythia8',
'SingleElectron',
'SingleMuon',
'TTTo2L2Nu_TuneCP5_PSweights_13TeV-powheg-pythia8',
'TTToHadronic_TuneCP5_PSweights_13TeV-powheg-pythia8',
'TTToSemiLeptonic_TuneCP5_PSweights_13TeV-powheg-pythia8',
'TTWJetsToLNu_TuneCP5_PSweights_13TeV-amcatnloFXFX-madspin-pythia8',
'TTZToLLNuNu_M-10_TuneCP5_PSweights_13TeV-amcatnlo-pythia8',
'TTZToLL_M-1to10_TuneCP5_13TeV-amcatnlo-pythia8',
'TT_Mtt-1000toInf_TuneCP5_PSweights_13TeV-powheg-pythia8',
'TT_Mtt-700to1000_TuneCP5_PSweights_13TeV-powheg-pythia8',
'TprimeTprime_M-1000_TuneCP5_13TeV-madgraph-pythia8',
'TprimeTprime_M-1100_TuneCP5_13TeV-madgraph-pythia8',
'TprimeTprime_M-1200_TuneCP5_13TeV-madgraph-pythia8',
'TprimeTprime_M-1300_TuneCP5_13TeV-madgraph-pythia8',
'TprimeTprime_M-1400_TuneCP5_13TeV-madgraph-pythia8',
'TprimeTprime_M-1500_TuneCP5_13TeV-madgraph-pythia8',
'TprimeTprime_M-1600_TuneCP5_13TeV-madgraph-pythia8',
'TprimeTprime_M-1700_TuneCP5_13TeV-madgraph-pythia8',
'TprimeTprime_M-1800_TuneCP5_13TeV-madgraph-pythia8',
'WJetsToLNu_HT-1200To2500_TuneCP5_13TeV-madgraphMLM-pythia8',
'WJetsToLNu_HT-2500ToInf_TuneCP5_13TeV-madgraphMLM-pythia8',
'WJetsToLNu_HT-400To600_TuneCP5_13TeV-madgraphMLM-pythia8',
'WJetsToLNu_HT-600To800_TuneCP5_13TeV-madgraphMLM-pythia8',
'WJetsToLNu_HT-800To1200_TuneCP5_13TeV-madgraphMLM-pythia8',
'WW_TuneCP5_13TeV-pythia8',
'WZ_TuneCP5_13TeV-pythia8',
'ZZ_TuneCP5_13TeV-pythia8',
'ttHToNonbb_M125_TuneCP5_13TeV-powheg-pythia8',
'ttHTobb_M125_TuneCP5_13TeV-powheg-pythia8',
]
for sample in dirList:
outList = ['none']
if 'Tprime' in sample: outList = ['BWBW','TZBW','THBW','TZTH','TZTZ','THTH']
elif 'Bprime' in sample: outList = ['TWTW','BZTW','BHTW','BZBH','BZBZ','BHBH']
elif 'TTTo' in sample: outList = ['Mtt0to700','Mtt700to1000','Mtt1000toInf']
for outlabel in outList:
outsample = sample+'_'+outlabel
if outlabel == 'none': outsample = sample
rootfiles = EOSlist_root_files(inDir+'/'+outsample)
print "------------ hadding Sample:",outsample,"---------------"
print 'N root files in',outsample,'=',len(rootfiles)
nFilesPerHadd = 900
onefile = ' root://cmseos.fnal.gov/'+inDir+'/'+outsample+'/'+rootfiles[-1]
manyfiles = nFilesPerHadd*onefile
lengthcheck = len('hadd -f root://cmseos.fnal.gov/'+outDir+'/'+outsample+'_hadd.root '+manyfiles)
if lengthcheck > 126000.:
toolong = lengthcheck - 126000.
num2remove = math.ceil(toolong/len(onefile))
nFilesPerHadd = int(nFilesPerHadd - num2remove)
print 'Length estimate reduced from',lengthcheck,'by',toolong,'via removing',num2remove,'files for nFilesPerHadd of',nFilesPerHadd
if len(rootfiles) < nFilesPerHadd:
haddcommand = 'hadd -f '+scratchDir+'/'+outsample+'_hadd.root '
for file in rootfiles:
haddcommand+=' root://cmseos.fnal.gov/'+inDir+'/'+outsample+'/'+file
print 'Length of hadd command =',len(haddcommand)
subprocess.call(haddcommand,shell=True)
xrdcpcommand = 'xrdcp '+scratchDir+'/'+outsample+'_hadd.root root://cmseos.fnal.gov/'+outDir+'/'+outsample+'_hadd.root'
subprocess.call(xrdcpcommand,shell=True)
if bool(EOSisfile(outDir+'/'+outsample+'_hadd.root')) != True:
print haddcommand
else:
for i in range(int(math.ceil(len(rootfiles)/float(nFilesPerHadd)))):
haddcommand = 'hadd -f '+scratchDir+'/'+outsample+'_'+str(i+1)+'_hadd.root '
begin=i*nFilesPerHadd
end=begin+nFilesPerHadd
if end > len(rootfiles): end=len(rootfiles)
print 'begin:',begin,'end:',end
for j in range(begin,end):
haddcommand+=' root://cmseos.fnal.gov/'+inDir+'/'+outsample+'/'+rootfiles[j]
print 'Length of hadd command =',len(haddcommand)
subprocess.call(haddcommand,shell=True)
xrdcpcommand = 'xrdcp '+scratchDir+'/'+outsample+'_'+str(i+1)+'_hadd.root root://cmseos.fnal.gov/'+outDir+'/'+outsample+'_'+str(i+1)+'_hadd.root'
subprocess.call(xrdcpcommand,shell=True)
if bool(EOSisfile(outDir+'/'+outsample+'_'+str(i+1)+'_hadd.root')) != True:
print haddcommand
print("--- %s minutes ---" % (round(time.time() - start_time, 2)/60))