ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/VHbbAnalysis/additionalFiles/DataCards/scalerZnn.py
(Generate patch)

Comparing UserCode/VHbbAnalysis/additionalFiles/DataCards/scalerZnn.py (file contents):
Revision 1.1 by madfish, Thu Oct 11 09:47:34 2012 UTC vs.
Revision 1.3 by madfish, Sat Nov 3 13:02:28 2012 UTC

# Line 6 | Line 6 | import os
6   import commands
7   import string
8  
9 < xSec7ZH = [0.4721,0.4655,0.4589,0.4525,0.4462,0.4400,0.4340,0.4280,0.4221,0.4164, 0.4107, 0.4052, 0.3998, 0.3945, 0.3893, 0.3842, 0.3791, 0.3742, 0.3693, 0.3645, 0.3598, 0.3551, 0.3505, 0.3459, 0.3414, 0.3370, 0.3326, 0.3283, 0.3241, 0.3199, 0.3158, 0.3117, 0.3077, 0.3038, 0.2999, 0.2961, 0.2923, 0.2886, 0.2849, 0.2813, 0.2778, 0.2743, 0.2709, 0.2675, 0.2642, 0.2609, 0.2577, 0.2545, 0.2514, 0.2483, 0.2453]
9 > xSec7ZH = [0.4721,0.4655,0.4589,0.4525,0.4462,0.4400,0.4340,0.4280,0.4221,0.4164, 0.4107, 0.4052, 0.3998, 0.3945, 0.3893, 0.3842, 0.3791, 0.3742, 0.3693, 0.3645, 0.3598, 0.3551, 0.3505, 0.3459, 0.3414, 0.3370, 0.3326, 0.3283, 0.3241, 0.3199, 0.3158, 0.3117, 0.3077, 0.3038, 0.2999, 0.2961, 0.2923, 0.2886, 0.2849, 0.2813, 0.2778, 0.2743, 0.2709, 0.2675, 0.2642, 0.2609, 0.2577, 0.2545, 0.2514, 0.2483, 0.2453, 0.31908, 0.31826, 0.31744, 0.31662, 0.314980, 0.314160, 0.313340, 0.312520, 0.310900, 0.310100, 0.309300, 0.308500, 0.306920, 0.306140, 0.305360, 0.304580]
10  
11 xSec7WH  = [0.8754,0.8623,0.8495,0.8368,0.8244, 0.8122, 0.8003, 0.7885, 0.7770, 0.7657, 0.7546, 0.7439, 0.7333, 0.7230, 0.7129, 0.7030, 0.6933, 0.6837, 0.6744, 0.6651, 0.6561, 0.6472, 0.6384, 0.6297, 0.6212, 0.6129, 0.6046, 0.5965, 0.5885, 0.5806, 0.5729, 0.5652, 0.5576, 0.5501, 0.5428, 0.5355, 0.5284, 0.5213, 0.5144, 0.5075, 0.5008, 0.4942, 0.4877, 0.4813, 0.4749, 0.4687, 0.4626, 0.4566, 0.4506, 0.4448, 0.4390]
11  
13 xSec8ZH = [0.5869,0.5788,0.5708,0.5629,0.5552,0.5476,0.5402,0.5329,0.5258,0.5187,0.5117,0.5049,0.4981,0.4916,0.4850,0.4787,0.4724,0.4662,0.4602,0.4542,0.4483,0.4426,0.4368,0.4312,0.4257,0.4203,0.4150,0.4096,0.4044,0.3993,0.3943,0.3893,0.3843,0.3794,0.3746,0.3699,0.3652,0.3606,0.3561,0.3516,0.3473,0.3430,0.3388,0.3347,0.3306,0.3266,0.3226,0.3188,0.3149,0.3112,0.3074];
14          
15 xSec8WH = [1.060 ,1.045 ,1.030 ,1.015 , 0.9998, 0.9852, 0.9709, 0.9570, 0.9432, 0.9297, 0.9165, 0.9035, 0.8907, 0.8782, 0.8659, 0.8538, 0.8420, 0.8303, 0.8187, 0.8075, 0.7966, 0.7859, 0.7753, 0.7649, 0.7547, 0.7446, 0.7347, 0.7249, 0.7154, 0.7060, 0.6966, 0.6873, 0.6782, 0.6691, 0.6602, 0.6515, 0.6429, 0.6344, 0.6260, 0.6177, 0.6095, 0.6015, 0.5936, 0.5859, 0.5783, 0.5708, 0.5634, 0.5562, 0.5491, 0.5420, 0.5351];
12  
13 < xSec = xSec8ZH
13 > xSec7WH  = [0.8754,0.8623,0.8495,0.8368,0.8244, 0.8122, 0.8003, 0.7885, 0.7770, 0.7657, 0.7546, 0.7439, 0.7333, 0.7230, 0.7129, 0.7030, 0.6933, 0.6837, 0.6744, 0.6651, 0.6561, 0.6472, 0.6384, 0.6297, 0.6212, 0.6129, 0.6046, 0.5965, 0.5885, 0.5806, 0.5729, 0.5652, 0.5576, 0.5501, 0.5428, 0.5355, 0.5284, 0.5213, 0.5144, 0.5075, 0.5008, 0.4942, 0.4877, 0.4813, 0.4749, 0.4687, 0.4626, 0.4566, 0.4506, 0.4448, 0.4390,0.579060,0.577520,0.575980,0.574440,0.571360,0.569820,0.568280,0.566740,0.563680,0.562160,0.560640,0.559120,0.556100,0.554600,0.553100,0.551600]
14  
15 < mass=["110","110_5","111","111_5","112","112_5","113","113_5","114","114_5","115","115_5","116","116_5","117","117_5","118","118_5","119","119_5","120","120_5","121","121_5","122","122_5","123","123_5","124","124_5","125","125_5","126","126_5","127","127_5","128","128_5","129","129_5","130","130_5","131","131_5","132","132_5","133","133_5","134","134_5","135"]
15 > xSec8ZH = [0.5869,0.5788,0.5708,0.5629,0.5552,0.5476,0.5402,0.5329,0.5258,0.5187,0.5117,0.5049,0.4981,0.4916,0.4850,0.4787,0.4724,0.4662,0.4602,0.4542,0.4483,0.4426,0.4368,0.4312,0.4257,0.4203,0.4150,0.4096,0.4044,0.3993,0.3943,0.3893,0.3843,0.3794,0.3746,0.3699,0.3652,0.3606,0.3561,0.3516,0.3473,0.3430,0.3388,0.3347,0.3306,0.3266,0.3226,0.3188,0.3149,0.3112,0.3074,0.39830,0.39730,0.39630,0.39530,0.39330,0.39230,0.39130,0.39030,0.38830,0.38730,0.38630,0.38530,0.38332,0.38234,0.38136,0.38038]
16 >          
17 > xSec8WH = [1.060 ,1.045 ,1.030 ,1.015 , 0.9998, 0.9852, 0.9709, 0.9570, 0.9432, 0.9297, 0.9165, 0.9035, 0.8907, 0.8782, 0.8659, 0.8538, 0.8420, 0.8303, 0.8187, 0.8075, 0.7966, 0.7859, 0.7753, 0.7649, 0.7547, 0.7446, 0.7347, 0.7249, 0.7154, 0.7060, 0.6966, 0.6873, 0.6782, 0.6691, 0.6602, 0.6515, 0.6429, 0.6344, 0.6260, 0.6177, 0.6095, 0.6015, 0.5936, 0.5859, 0.5783, 0.5708, 0.5634, 0.5562, 0.5491, 0.5420, 0.5351, 0.70412,0.70224,0.70036,0.69848,0.69474,0.69288,0.69102,0.68916,0.68548,0.68366,0.68184,0.68002,0.67638,0.67456,0.67274,0.67092]
18 >
19 >
20 >
21 > mass=["110","110.5","111","111.5","112","112.5","113","113.5","114","114.5","115","115.5","116","116.5","117","117.5","118","118.5","119","119.5","120","120.5","121","121.5","122","122.5","123","123.5","124","124.5","125","125.5","126","126.5","127","127.5","128","128.5","129","129.5","130","130.5","131","131.5","132","132.5","133","133.5","134","134.5","135","124.6","124.7","124.8","124.9","125.1","125.2","125.3","125.4","125.6","125.7","125.8","125.9","126.1","126.2","126.3","126.4"];
22  
23  
24   files=[
25 < #"vhbb_DC_BDTRMed_M110_Zee_8TeV.txt",
26 < #"vhbb_DC_BDTRMed_M110_Zmm_8TeV.txt",
27 < #"vhbb_DC_BDTRMed_M115_Zee_8TeV.txt",
28 < #"vhbb_DC_BDTRMed_M115_Zmm_8TeV.txt",
29 < #"vhbb_DC_BDTRMed_M120_Zee_8TeV.txt",
30 < #"vhbb_DC_BDTRMed_M120_Zmm_8TeV.txt",
29 < #"vhbb_DC_BDTRMed_M125_Zee_8TeV.txt",
30 < #"vhbb_DC_BDTRMed_M125_Zmm_8TeV.txt",
31 < #"vhbb_DC_BDTRMed_M130_Zee_8TeV.txt",
32 < #"vhbb_DC_BDTRMed_M130_Zmm_8TeV.txt",
33 < #"vhbb_DC_BDTRMed_M135_Zee_8TeV.txt",
34 < #"vhbb_DC_BDTRMed_M135_Zmm_8TeV.txt",
35 < #"vhbb_DC_BDTRTight_M110_Zee_8TeV.txt",
36 < #"vhbb_DC_BDTRTight_M110_Zmm_8TeV.txt",
37 < #"vhbb_DC_BDTRTight_M115_Zee_8TeV.txt",
38 < #"vhbb_DC_BDTRTight_M115_Zmm_8TeV.txt",
39 < #"vhbb_DC_BDTRTight_M120_Zee_8TeV.txt",
40 < #"vhbb_DC_BDTRTight_M120_Zmm_8TeV.txt",
41 < #"vhbb_DC_BDTRTight_M125_Zee_8TeV.txt",
42 < #"vhbb_DC_BDTRTight_M125_Zmm_8TeV.txt",
43 < #"vhbb_DC_BDTRTight_M130_Zee_8TeV.txt",
44 < #"vhbb_DC_BDTRTight_M130_Zmm_8TeV.txt",
45 < #"vhbb_DC_BDTRTight_M135_Zee_8TeV.txt",
46 < #"vhbb_DC_BDTRTight_M135_Zmm_8TeV.txt",
47 < "vhbb_DC_BDT_M110_Zll_8TeV.txt",
48 < "vhbb_DC_BDT_M115_Zll_8TeV.txt",
49 < "vhbb_DC_BDT_M120_Zll_8TeV.txt",
50 < "vhbb_DC_BDT_M125_Zll_8TeV.txt",
51 < "vhbb_DC_BDT_M130_Zll_8TeV.txt",
52 < "vhbb_DC_BDT_M135_Zll_8TeV.txt"
25 > "vhbb_DC_BDT_M110_Znn_8TeV.txt",
26 > "vhbb_DC_BDT_M115_Znn_8TeV.txt",
27 > "vhbb_DC_BDT_M120_Znn_8TeV.txt",
28 > "vhbb_DC_BDT_M125_Znn_8TeV.txt",
29 > "vhbb_DC_BDT_M130_Znn_8TeV.txt",
30 > "vhbb_DC_BDT_M135_Znn_8TeV.txt"
31   ]
32  
33  
34  
35   def Process(file, toMass, fromMass):
36 <       newcard = file.replace( '.txt', "_.txt"  )
37 <       newcard = newcard.replace( '.txt', mass[toMass]  )
38 <       newcard += ".txt"
36 >       newcard = "./%s/aaa" % mass[toMass]
37 >       newcard = newcard.replace( 'aaa', file)
38 >       newcard = newcard.replace('_DC_BDT_M110','')
39 >       newcard = newcard.replace('_DC_BDT_M115','')
40 >       newcard = newcard.replace('_DC_BDT_M120','')
41 >       newcard = newcard.replace('_DC_BDT_M125','')
42 >       newcard = newcard.replace('_DC_BDT_M130','')
43 >       newcard = newcard.replace('_DC_BDT_M135','')
44 >
45         os.system('more %s | grep rate > a' % file)
46 <       os.system("more a |  awk '{print $2}'  > aVH1")
47 <       os.system("more a |  awk '{print $11}'  > aVH2")
48 <       os.system("more a |  awk '{print $20}'  > aVH3")
49 <       os.system("more a |  awk '{print $29}'  > aVH4")
46 >       os.system("more a |  awk '{print $2}'  > aZH1")
47 >       os.system("more a |  awk '{print $12}'  > aZH2")
48 >       os.system("more a |  awk '{print $22}'  > aZH3")
49 >       os.system("more a |  awk '{print $3}'  > aWH1")
50 >       os.system("more a |  awk '{print $13}'  > aWH2")
51 >       os.system("more a |  awk '{print $23}'  > aWH3")
52  
53 <       fVH1 = open("aVH1", "r")
54 <       fVH2 = open("aVH2", "r")
55 <       fVH3 = open("aVH3", "r")
56 <       fVH4 = open("aVH4", "r")
57 <
53 >       fZH1 = open("aZH1", "r")
54 >       fZH2 = open("aZH2", "r")
55 >       fZH3 = open("aZH3", "r")
56 >       fWH1 = open("aWH1", "r")
57 >       fWH2 = open("aWH2", "r")
58 >       fWH3 = open("aWH3", "r")
59 >
60 >       sZH1 = fZH1.read().rstrip('\n')
61 >       sZH2 = fZH2.read().rstrip('\n')
62 >       sZH3 = fZH3.read().rstrip('\n')
63 >       sWH1 = fWH1.read().rstrip('\n')
64 >       sWH2 = fWH2.read().rstrip('\n')
65 >       sWH3 = fWH3.read().rstrip('\n')
66  
73       sVH1 = fVH1.read().rstrip('\n')
74       sVH2 = fVH2.read().rstrip('\n')
75       sVH3 = fVH3.read().rstrip('\n')
76       sVH4 = fVH4.read().rstrip('\n')
67  
68  
69 <       VH1 = float(sVH1)
70 <       VH2 = float(sVH2)
71 <       VH3 = float(sVH3)
72 <       VH4 = float(sVH4)
73 <
74 <
75 <       bVH1 = (xSec[toMass]/xSec[fromMass])*VH1
76 <       bVH2 = (xSec[toMass]/xSec[fromMass])*VH2
77 <       bVH3 = (xSec[toMass]/xSec[fromMass])*VH3
78 <       bVH4 = (xSec[toMass]/xSec[fromMass])*VH4
79 <
80 <       os.system("sed 's/%s/%f/g' a > b" % (sVH1,bVH1))
81 <       os.system("sed -i 's/%s/%f/g' b" % (sVH2,bVH2))
82 <       os.system("sed -i 's/%s/%f/g' b" % (sVH3,bVH3))
83 <       os.system("sed -i 's/%s/%f/g' b" % (sVH4,bVH4))
69 >       ZH1 = float(sZH1)
70 >       ZH2 = float(sZH2)
71 >       ZH3 = float(sZH3)
72 >       WH1 = float(sWH1)
73 >       WH2 = float(sWH2)
74 >       WH3 = float(sWH3)
75 >
76 >       bZH1 = (xSec8ZH[toMass]/xSec8ZH[fromMass])*ZH1
77 >       bZH2 = (xSec8ZH[toMass]/xSec8ZH[fromMass])*ZH2
78 >       bZH3 = (xSec8ZH[toMass]/xSec8ZH[fromMass])*ZH3
79 >       bWH1 = (xSec8WH[toMass]/xSec8WH[fromMass])*WH1
80 >       bWH2 = (xSec8WH[toMass]/xSec8WH[fromMass])*WH2
81 >       bWH3 = (xSec8WH[toMass]/xSec8WH[fromMass])*WH3
82 >
83 >
84 >       os.system("sed 's/%s/%f/g' a > b" % (sWH1,bWH1))
85 >       os.system("sed -i 's/%s/%f/g' b" % (sWH2,bWH2))
86 >       os.system("sed -i 's/%s/%f/g' b" % (sWH3,bWH3))
87 >       os.system("sed -i 's/%s/%f/g' b" % (sZH1,bZH1))
88 >       os.system("sed -i 's/%s/%f/g' b" % (sZH2,bZH2))
89 >       os.system("sed -i 's/%s/%f/g' b" % (sZH3,bZH3))
90  
91         fRO = open("a", "r")
92         sRO = fRO.read().rstrip('\n')
93         fRF = open("b", "r")
94         sRF = fRF.read().rstrip('\n')
95         os.system("sed 's/%s/%s/g' %s > %s" % (sRO,sRF,file,newcard))
96 <
97 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M110_Zmm_8TeV.root/vhbb_WS_BDTRMed_M110_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
98 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M115_Zmm_8TeV.root/vhbb_WS_BDTRMed_M115_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
99 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M120_Zmm_8TeV.root/vhbb_WS_BDTRMed_M120_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
100 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M125_Zmm_8TeV.root/vhbb_WS_BDTRMed_M125_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
101 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M130_Zmm_8TeV.root/vhbb_WS_BDTRMed_M130_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
106 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M135_Zmm_8TeV.root/vhbb_WS_BDTRMed_M135_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
107 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M110_Zmm_8TeV.root/vhbb_WS_BDTRTight_M110_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
108 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M115_Zmm_8TeV.root/vhbb_WS_BDTRTight_M115_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
109 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M120_Zmm_8TeV.root/vhbb_WS_BDTRTight_M120_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
110 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M125_Zmm_8TeV.root/vhbb_WS_BDTRTight_M125_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
111 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M130_Zmm_8TeV.root/vhbb_WS_BDTRTight_M130_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
112 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M135_Zmm_8TeV.root/vhbb_WS_BDTRTight_M135_Zmm_8TeV%s.root/g' %s " % (mass[toMass],newcard))
113 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M110_Zee_8TeV.root/vhbb_WS_BDTRMed_M110_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
114 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M115_Zee_8TeV.root/vhbb_WS_BDTRMed_M115_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
115 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M120_Zee_8TeV.root/vhbb_WS_BDTRMed_M120_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
116 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M125_Zee_8TeV.root/vhbb_WS_BDTRMed_M125_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
117 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M130_Zee_8TeV.root/vhbb_WS_BDTRMed_M130_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
118 <       os.system("sed -i 's/vhbb_WS_BDTRMed_M135_Zee_8TeV.root/vhbb_WS_BDTRMed_M135_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
119 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M110_Zee_8TeV.root/vhbb_WS_BDTRTight_M110_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
120 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M115_Zee_8TeV.root/vhbb_WS_BDTRTight_M115_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
121 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M120_Zee_8TeV.root/vhbb_WS_BDTRTight_M120_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
122 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M125_Zee_8TeV.root/vhbb_WS_BDTRTight_M125_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
123 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M130_Zee_8TeV.root/vhbb_WS_BDTRTight_M130_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
124 <       os.system("sed -i 's/vhbb_WS_BDTRTight_M135_Zee_8TeV.root/vhbb_WS_BDTRTight_M135_Zee_8TeV%s.root/g' %s " % (mass[toMass],newcard))
96 >       os.system("sed -i 's/vhbb_WS_BDT_M110_Znn_8TeV.root/vhbb_Znn_8TeV.root/g' %s " % (newcard))
97 >       os.system("sed -i 's/vhbb_WS_BDT_M115_Znn_8TeV.root/vhbb_Znn_8TeV.root/g' %s " % (newcard))
98 >       os.system("sed -i 's/vhbb_WS_BDT_M120_Znn_8TeV.root/vhbb_Znn_8TeV.root/g' %s " % (newcard))
99 >       os.system("sed -i 's/vhbb_WS_BDT_M125_Znn_8TeV.root/vhbb_Znn_8TeV.root/g' %s " % (newcard))
100 >       os.system("sed -i 's/vhbb_WS_BDT_M130_Znn_8TeV.root/vhbb_Znn_8TeV.root/g' %s " % (newcard))
101 >       os.system("sed -i 's/vhbb_WS_BDT_M135_Znn_8TeV.root/vhbb_Znn_8TeV.root/g' %s " % (newcard))
102         os.system("echo %s written" % (newcard))
103  
104 < for item in files:
105 <
106 < if "110" in item:
107 <  Process(item,0,0);
108 <  Process(item,1,0);
132 <  Process(item,2,0);
133 <  Process(item,3,0);
134 <  Process(item,4,0);
104 > Process(files[0],0,0);
105 > Process(files[0],1,0);
106 > Process(files[0],2,0);
107 > Process(files[0],3,0);
108 > Process(files[0],4,0);
109    
110 < if "115" in item:
111 <  for x in range(5,15):
112 <    Process(item, x , 10);
113 <
114 < if "120" in item:
115 <  for x in range(15,25):
116 <    Process(item, x , 20);
117 <
118 <
119 < if "125" in item:
120 <  for x in range(25,35):
121 <    Process(item, x , 30);
148 <
149 <
150 < if "130" in item:
151 <  for x in range(35,45):
152 <    Process(item, x , 40);
153 <
154 <
155 < if "135" in item:
156 <  for x in range(45,51):
157 <    Process(item, x , 50);
110 > for x in range(5,15):
111 >  Process(files[1], x , 10);
112 > for x in range(15,25):
113 >  Process(files[2], x , 20);
114 > for x in range(25,35):
115 >  Process(files[3], x , 30);
116 > for x in range(35,45):
117 >  Process(files[4], x , 40);
118 > for x in range(45,51):
119 >  Process(files[5], x , 50);
120 > for x in range(51,67):
121 >  Process(files[3], x , 30);

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines