| 1 |
#! /usr/bin/env python
|
| 2 |
|
| 3 |
|
| 4 |
import sys
|
| 5 |
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]
|
| 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]
|
| 12 |
|
| 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];
|
| 16 |
|
| 17 |
xSec = xSec8ZH
|
| 18 |
|
| 19 |
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"]
|
| 20 |
|
| 21 |
|
| 22 |
files=[
|
| 23 |
#"vhbb_DC_BDTRMed_M110_Zee_8TeV.txt",
|
| 24 |
#"vhbb_DC_BDTRMed_M110_Zmm_8TeV.txt",
|
| 25 |
#"vhbb_DC_BDTRMed_M115_Zee_8TeV.txt",
|
| 26 |
#"vhbb_DC_BDTRMed_M115_Zmm_8TeV.txt",
|
| 27 |
#"vhbb_DC_BDTRMed_M120_Zee_8TeV.txt",
|
| 28 |
#"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"
|
| 53 |
]
|
| 54 |
|
| 55 |
|
| 56 |
|
| 57 |
def Process(file, toMass, fromMass):
|
| 58 |
newcard = file.replace( '.txt', "_.txt" )
|
| 59 |
newcard = newcard.replace( '.txt', mass[toMass] )
|
| 60 |
newcard += ".txt"
|
| 61 |
os.system('more %s | grep rate > a' % file)
|
| 62 |
os.system("more a | awk '{print $2}' > aVH1")
|
| 63 |
os.system("more a | awk '{print $11}' > aVH2")
|
| 64 |
os.system("more a | awk '{print $20}' > aVH3")
|
| 65 |
os.system("more a | awk '{print $29}' > aVH4")
|
| 66 |
|
| 67 |
fVH1 = open("aVH1", "r")
|
| 68 |
fVH2 = open("aVH2", "r")
|
| 69 |
fVH3 = open("aVH3", "r")
|
| 70 |
fVH4 = open("aVH4", "r")
|
| 71 |
|
| 72 |
|
| 73 |
sVH1 = fVH1.read().rstrip('\n')
|
| 74 |
sVH2 = fVH2.read().rstrip('\n')
|
| 75 |
sVH3 = fVH3.read().rstrip('\n')
|
| 76 |
sVH4 = fVH4.read().rstrip('\n')
|
| 77 |
|
| 78 |
|
| 79 |
VH1 = float(sVH1)
|
| 80 |
VH2 = float(sVH2)
|
| 81 |
VH3 = float(sVH3)
|
| 82 |
VH4 = float(sVH4)
|
| 83 |
|
| 84 |
|
| 85 |
bVH1 = (xSec[toMass]/xSec[fromMass])*VH1
|
| 86 |
bVH2 = (xSec[toMass]/xSec[fromMass])*VH2
|
| 87 |
bVH3 = (xSec[toMass]/xSec[fromMass])*VH3
|
| 88 |
bVH4 = (xSec[toMass]/xSec[fromMass])*VH4
|
| 89 |
|
| 90 |
os.system("sed 's/%s/%f/g' a > b" % (sVH1,bVH1))
|
| 91 |
os.system("sed -i 's/%s/%f/g' b" % (sVH2,bVH2))
|
| 92 |
os.system("sed -i 's/%s/%f/g' b" % (sVH3,bVH3))
|
| 93 |
os.system("sed -i 's/%s/%f/g' b" % (sVH4,bVH4))
|
| 94 |
|
| 95 |
fRO = open("a", "r")
|
| 96 |
sRO = fRO.read().rstrip('\n')
|
| 97 |
fRF = open("b", "r")
|
| 98 |
sRF = fRF.read().rstrip('\n')
|
| 99 |
os.system("sed 's/%s/%s/g' %s > %s" % (sRO,sRF,file,newcard))
|
| 100 |
|
| 101 |
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))
|
| 102 |
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))
|
| 103 |
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))
|
| 104 |
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))
|
| 105 |
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))
|
| 125 |
os.system("echo %s written" % (newcard))
|
| 126 |
|
| 127 |
for item in files:
|
| 128 |
|
| 129 |
if "110" in item:
|
| 130 |
Process(item,0,0);
|
| 131 |
Process(item,1,0);
|
| 132 |
Process(item,2,0);
|
| 133 |
Process(item,3,0);
|
| 134 |
Process(item,4,0);
|
| 135 |
|
| 136 |
if "115" in item:
|
| 137 |
for x in range(5,15):
|
| 138 |
Process(item, x , 10);
|
| 139 |
|
| 140 |
if "120" in item:
|
| 141 |
for x in range(15,25):
|
| 142 |
Process(item, x , 20);
|
| 143 |
|
| 144 |
|
| 145 |
if "125" in item:
|
| 146 |
for x in range(25,35):
|
| 147 |
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);
|
