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.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 17 | Line 17 | xSec8ZH = [0.5869,0.5788,0.5708,0.5629,0
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 – xSec = xSec8ZH
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 < "WH_2bin_110.txt",
26 < "WH_2bin_115.txt",
27 < "WH_2bin_120.txt",
28 < "vhbb_Znn_7TeV.txt",
29 < "WH_2bin_130.txt",
30 < "WH_2bin_135.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  
# Line 36 | Line 35 | files=[
35   def Process(file, toMass, fromMass):
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")
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 <
56 <       sVH1 = fVH1.read().rstrip('\n')
57 <       sVH2 = fVH2.read().rstrip('\n')
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  
49 –       VH1 = float(sVH1)
50 –       VH2 = float(sVH2)
67  
52 –       bVH1 = (xSec[toMass]/xSec[fromMass])*VH1
53 –       bVH2 = (xSec[toMass]/xSec[fromMass])*VH2
68  
69 <       os.system("sed 's/%s/%f/g' a > b" % (sVH1,bVH1))
70 <       os.system("sed -i 's/%s/%f/g' b" % (sVH2,bVH2))
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 <       os.system("sed -i 's/WmunuH_110_NWS.root/WmunuH_110_%s.root/g' %s " % (mass[toMass],newcard))
97 <       os.system("sed -i 's/WmunuH_115_NWS.root/WmunuH_115_%s.root/g' %s " % (mass[toMass],newcard))
98 <       os.system("sed -i 's/WmunuH_120_NWS.root/WmunuH_120_%s.root/g' %s " % (mass[toMass],newcard))
99 <       os.system("sed -i 's/WmunuH_125_NWS.root/WmunuH_125_%s.root/g' %s " % (mass[toMass],newcard))
100 <       os.system("sed -i 's/WmunuH_130_NWS.root/WmunuH_130_%s.root/g' %s " % (mass[toMass],newcard))
101 <       os.system("sed -i 's/WmunuH_135_NWS.root/WmunuH_135_%s.root/g' %s " % (mass[toMass],newcard))
69 <       os.system("sed -i 's/WenuH_110_NWS.root/WenuH_110_%s.root/g' %s " % (mass[toMass],newcard))
70 <       os.system("sed -i 's/WenuH_115_NWS.root/WenuH_115_%s.root/g' %s " % (mass[toMass],newcard))
71 <       os.system("sed -i 's/WenuH_120_NWS.root/WenuH_120_%s.root/g' %s " % (mass[toMass],newcard))
72 <       os.system("sed -i 's/WenuH_125_NWS.root/WenuH_125_%s.root/g' %s " % (mass[toMass],newcard))
73 <       os.system("sed -i 's/WenuH_130_NWS.root/WenuH_130_%s.root/g' %s " % (mass[toMass],newcard))
74 <       os.system("sed -i 's/WenuH_135_NWS.root/WenuH_135_%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   Process(files[0],0,0);

Diff Legend

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