ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/VHbbAnalysis/additionalFiles/DataCards/scalerWH.py

Comparing UserCode/VHbbAnalysis/additionalFiles/DataCards/scalerWH.py (file contents):
Revision 1.2 by madfish, Sat Oct 27 20:34:50 2012 UTC vs.
Revision 1.3 by madfish, Sat Nov 3 13:02:28 2012 UTC

#	Line 23 | Line 23 | mass=["110","110.5","111","111.5","112",
23
24
25		files=[
26	<	"WH_2bin_110.txt",
27	<	"WH_2bin_115.txt",
28	<	"WH_2bin_120.txt",
29	<	"vhbb_Wln_7TeV.txt",
30	<	"WH_2bin_130.txt",
31	<	"WH_2bin_135.txt"
26	>	"combined_vhbb_DC_BDT_M110_8TeV.txt",
27	>	"combined_vhbb_DC_BDT_M115_8TeV.txt",
28	>	"combined_vhbb_DC_BDT_M120_8TeV.txt",
29	>	"combined_vhbb_DC_BDT_M125_8TeV.txt",
30	>	"combined_vhbb_DC_BDT_M130_8TeV.txt",
31	>	"combined_vhbb_DC_BDT_M135_8TeV.txt"
32		]
33
34
35	–
35		def Process(file, toMass, fromMass):
36	<	newcard = file.replace( '.txt', "_.txt"  )
37	<	newcard = newcard.replace( '.txt', mass[toMass]  )
38	<	newcard += ".txt"
36	>	file2 = file.replace("combined_","")
37	>	file2 = file2.replace("DC_BDT_M110","Wln")
38	>	file2 = file2.replace("DC_BDT_M115","Wln")
39	>	file2 = file2.replace("DC_BDT_M120","Wln")
40	>	file2 = file2.replace("DC_BDT_M125","Wln")
41	>	file2 = file2.replace("DC_BDT_M130","Wln")
42	>	file2 = file2.replace("DC_BDT_M135","Wln")
43	>	newcard = mass[toMass] + "/" + file2
44		os.system('more %s | grep rate > a' % file)
45	<	os.system("more a |  awk '{print $2}'  > aVH1")
46	<	os.system("more a |  awk '{print $11}'  > aVH2")
47	<	os.system("more a |  awk '{print $20}'  > aVH3")
48	<	os.system("more a |  awk '{print $29}'  > aVH4")
49	<
50	<	fVH1 = open("aVH1", "r")
51	<	fVH2 = open("aVH2", "r")
52	<	fVH3 = open("aVH3", "r")
53	<	fVH4 = open("aVH4", "r")
54	<
55	<	sVH1 = fVH1.read().rstrip('\n')
56	<	sVH2 = fVH2.read().rstrip('\n')
57	<	sVH3 = fVH3.read().rstrip('\n')
58	<	sVH4 = fVH4.read().rstrip('\n')
59	<
60	<	VH1 = float(sVH1)
61	<	VH2 = float(sVH2)
62	<	VH3 = float(sVH3)
63	<	VH4 = float(sVH4)
64	<
65	<	bVH1 = (xSec[toMass]/xSec[fromMass])*VH1
66	<	bVH2 = (xSec[toMass]/xSec[fromMass])*VH2
67	<	bVH3 = (xSec[toMass]/xSec[fromMass])*VH3
68	<	bVH4 = (xSec[toMass]/xSec[fromMass])*VH4
69	<
70	<	os.system("sed 's/%s/%f/g' a > b" % (sVH1,bVH1))
71	<	os.system("sed -i 's/%s/%f/g' b" % (sVH2,bVH2))
72	<	os.system("sed -i 's/%s/%f/g' b" % (sVH3,bVH3))
73	<	os.system("sed -i 's/%s/%f/g' b" % (sVH4,bVH4))
45	>	os.system("more a |  awk '{print $2} '  > aZH1")
46	>	os.system("more a |  awk '{print $12}'  > aZH2")
47	>	os.system("more a |  awk '{print $22}'  > aZH3")
48	>	os.system("more a |  awk '{print $32}'  > aZH4")
49	>	os.system("more a |  awk '{print $42}'  > aZH5")
50	>	os.system("more a |  awk '{print $52}'  > aZH6")
51	>	os.system("more a |  awk '{print $3}'   > aWH1")
52	>	os.system("more a |  awk '{print $13}'  > aWH2")
53	>	os.system("more a |  awk '{print $23}'  > aWH3")
54	>	os.system("more a |  awk '{print $33}'  > aWH4")
55	>	os.system("more a |  awk '{print $43}'  > aWH5")
56	>	os.system("more a |  awk '{print $53}'  > aWH6")
57	>
58	>	fZH1 = open("aZH1", "r")
59	>	fZH2 = open("aZH2", "r")
60	>	fZH3 = open("aZH3", "r")
61	>	fZH4 = open("aZH4", "r")
62	>	fZH5 = open("aZH5", "r")
63	>	fZH6 = open("aZH6", "r")
64	>	fWH1 = open("aWH1", "r")
65	>	fWH2 = open("aWH2", "r")
66	>	fWH3 = open("aWH3", "r")
67	>	fWH4 = open("aWH4", "r")
68	>	fWH5 = open("aWH5", "r")
69	>	fWH6 = open("aWH6", "r")
70	>
71	>	sZH1 = fZH1.read().rstrip('\n')
72	>	sZH2 = fZH2.read().rstrip('\n')
73	>	sZH3 = fZH3.read().rstrip('\n')
74	>	sZH4 = fZH4.read().rstrip('\n')
75	>	sZH5 = fZH5.read().rstrip('\n')
76	>	sZH6 = fZH6.read().rstrip('\n')
77	>	sWH1 = fWH1.read().rstrip('\n')
78	>	sWH2 = fWH2.read().rstrip('\n')
79	>	sWH3 = fWH3.read().rstrip('\n')
80	>	sWH4 = fWH4.read().rstrip('\n')
81	>	sWH5 = fWH5.read().rstrip('\n')
82	>	sWH6 = fWH6.read().rstrip('\n')
83	>
84	>
85	>	ZH1 = float(sZH1)
86	>	ZH2 = float(sZH2)
87	>	ZH3 = float(sZH3)
88	>	ZH4 = float(sZH4)
89	>	ZH5 = float(sZH5)
90	>	ZH6 = float(sZH6)
91	>
92	>	WH1 = float(sWH1)
93	>	WH2 = float(sWH2)
94	>	WH3 = float(sWH3)
95	>	WH4 = float(sWH4)
96	>	WH5 = float(sWH5)
97	>	WH6 = float(sWH6)
98	>
99	>
100	>	bZH1 = (xSec[toMass]/xSec[fromMass])*ZH1
101	>	bZH2 = (xSec[toMass]/xSec[fromMass])*ZH2
102	>	bZH3 = (xSec[toMass]/xSec[fromMass])*ZH3
103	>	bZH4 = (xSec[toMass]/xSec[fromMass])*ZH4
104	>	bZH5 = (xSec[toMass]/xSec[fromMass])*ZH5
105	>	bZH6 = (xSec[toMass]/xSec[fromMass])*ZH6
106	>
107	>	bWH1 = (xSec[toMass]/xSec[fromMass])*WH1
108	>	bWH2 = (xSec[toMass]/xSec[fromMass])*WH2
109	>	bWH3 = (xSec[toMass]/xSec[fromMass])*WH3
110	>	bWH4 = (xSec[toMass]/xSec[fromMass])*WH4
111	>	bWH5 = (xSec[toMass]/xSec[fromMass])*WH5
112	>	bWH6 = (xSec[toMass]/xSec[fromMass])*WH6
113	>
114	>
115	>
116	>	os.system("sed 's/%s/%f/g' a > b" %(sZH1,bZH1))
117	>	os.system("sed -i 's/%s/%f/g' b" % (sZH2,bZH2))
118	>	os.system("sed -i 's/%s/%f/g' b" % (sZH3,bZH3))
119	>	os.system("sed -i 's/%s/%f/g' b" % (sZH4,bZH4))
120	>	os.system("sed -i 's/%s/%f/g' b" % (sZH5,bZH5))
121	>	os.system("sed -i 's/%s/%f/g' b" % (sZH6,bZH6))
122	>
123	>	os.system("sed -i 's/%s/%f/g' b" % (sWH1,bWH1))
124	>	os.system("sed -i 's/%s/%f/g' b" % (sWH2,bWH2))
125	>	os.system("sed -i 's/%s/%f/g' b" % (sWH3,bWH3))
126	>	os.system("sed -i 's/%s/%f/g' b" % (sWH4,bWH4))
127	>	os.system("sed -i 's/%s/%f/g' b" % (sWH5,bWH5))
128	>	os.system("sed -i 's/%s/%f/g' b" % (sWH6,bWH6))
129	>
130
131		fRO = open("a", "r")
132		sRO = fRO.read().rstrip('\n')
133		fRF = open("b", "r")
134		sRF = fRF.read().rstrip('\n')
135		os.system("sed 's/%s/%s/g' %s > %s" % (sRO,sRF,file,newcard))
136	<	os.system("sed -i 's/WmunuH_110_NWS.root/WmunuH_110_%s.root/g' %s " % (mass[toMass],newcard))
137	<	os.system("sed -i 's/WmunuH_115_NWS.root/WmunuH_115_%s.root/g' %s " % (mass[toMass],newcard))
138	<	os.system("sed -i 's/WmunuH_120_NWS.root/WmunuH_120_%s.root/g' %s " % (mass[toMass],newcard))
139	<	os.system("sed -i 's/WmunuH_125_NWS.root/WmunuH_125_%s.root/g' %s " % (mass[toMass],newcard))
140	<	os.system("sed -i 's/WmunuH_130_NWS.root/WmunuH_130_%s.root/g' %s " % (mass[toMass],newcard))
141	<	os.system("sed -i 's/WmunuH_135_NWS.root/WmunuH_135_%s.root/g' %s " % (mass[toMass],newcard))
142	<	os.system("sed -i 's/WenuH_110_NWS.root/WenuH_110_%s.root/g' %s " % (mass[toMass],newcard))
143	<	os.system("sed -i 's/WenuH_115_NWS.root/WenuH_115_%s.root/g' %s " % (mass[toMass],newcard))
144	<	os.system("sed -i 's/WenuH_120_NWS.root/WenuH_120_%s.root/g' %s " % (mass[toMass],newcard))
145	<	os.system("sed -i 's/WenuH_125_NWS.root/WenuH_125_%s.root/g' %s " % (mass[toMass],newcard))
146	<	os.system("sed -i 's/WenuH_130_NWS.root/WenuH_130_%s.root/g' %s " % (mass[toMass],newcard))
147	<	os.system("sed -i 's/WenuH_135_NWS.root/WenuH_135_%s.root/g' %s " % (mass[toMass],newcard))
136	>
137	>	os.system("sed -i 's/vhbb_WS_BDT_M110_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
138	>	os.system("sed -i 's/vhbb_WS_BDT_M110_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
139	>	os.system("sed -i 's/vhbb_WS_BDT_M110_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
140	>	os.system("sed -i 's/vhbb_WS_BDT_M110_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
141	>	os.system("sed -i 's/vhbb_WS_BDT_M110_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
142	>	os.system("sed -i 's/vhbb_WS_BDT_M110_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
143	>
144	>	os.system("sed -i 's/vhbb_WS_BDT_M115_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
145	>	os.system("sed -i 's/vhbb_WS_BDT_M115_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
146	>	os.system("sed -i 's/vhbb_WS_BDT_M115_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
147	>	os.system("sed -i 's/vhbb_WS_BDT_M115_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
148	>	os.system("sed -i 's/vhbb_WS_BDT_M115_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
149	>	os.system("sed -i 's/vhbb_WS_BDT_M115_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
150	>
151	>	os.system("sed -i 's/vhbb_WS_BDT_M120_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
152	>	os.system("sed -i 's/vhbb_WS_BDT_M120_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
153	>	os.system("sed -i 's/vhbb_WS_BDT_M120_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
154	>	os.system("sed -i 's/vhbb_WS_BDT_M120_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
155	>	os.system("sed -i 's/vhbb_WS_BDT_M120_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
156	>	os.system("sed -i 's/vhbb_WS_BDT_M120_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
157	>
158	>	os.system("sed -i 's/vhbb_WS_BDT_M125_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
159	>	os.system("sed -i 's/vhbb_WS_BDT_M125_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
160	>	os.system("sed -i 's/vhbb_WS_BDT_M125_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
161	>	os.system("sed -i 's/vhbb_WS_BDT_M125_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
162	>	os.system("sed -i 's/vhbb_WS_BDT_M125_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
163	>	os.system("sed -i 's/vhbb_WS_BDT_M125_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
164	>
165	>	os.system("sed -i 's/vhbb_WS_BDT_M130_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
166	>	os.system("sed -i 's/vhbb_WS_BDT_M130_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
167	>	os.system("sed -i 's/vhbb_WS_BDT_M130_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
168	>	os.system("sed -i 's/vhbb_WS_BDT_M130_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
169	>	os.system("sed -i 's/vhbb_WS_BDT_M130_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
170	>	os.system("sed -i 's/vhbb_WS_BDT_M130_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
171	>
172	>	os.system("sed -i 's/vhbb_WS_BDT_M135_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
173	>	os.system("sed -i 's/vhbb_WS_BDT_M135_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
174	>	os.system("sed -i 's/vhbb_WS_BDT_M135_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
175	>	os.system("sed -i 's/vhbb_WS_BDT_M135_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
176	>	os.system("sed -i 's/vhbb_WS_BDT_M135_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
177	>	os.system("sed -i 's/vhbb_WS_BDT_M135_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
178	>
179	>
180	>
181		os.system("echo %s written" % (newcard))
182
183		Process(files[0],0,0);

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines
>	Changed lines