additionalFiles/DataCards/scalerWH.py

#! /usr/bin/env python


import sys
import os
import commands
import string

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]


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]

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]
          
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]


xSec = xSec8WH

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"];


files=[
"combined_vhbb_DC_BDT_M110_8TeV.txt",
"combined_vhbb_DC_BDT_M115_8TeV.txt",
"combined_vhbb_DC_BDT_M120_8TeV.txt",
"combined_vhbb_DC_BDT_M125_8TeV.txt",
"combined_vhbb_DC_BDT_M130_8TeV.txt",
"combined_vhbb_DC_BDT_M135_8TeV.txt"
]


def Process(file, toMass, fromMass):
       file2 = file.replace("combined_","")
       file2 = file2.replace("DC_BDT_M110","Wln")
       file2 = file2.replace("DC_BDT_M115","Wln")
       file2 = file2.replace("DC_BDT_M120","Wln")
       file2 = file2.replace("DC_BDT_M125","Wln")
       file2 = file2.replace("DC_BDT_M130","Wln")
       file2 = file2.replace("DC_BDT_M135","Wln")
       newcard = mass[toMass] + "/" + file2
       os.system('more %s | grep rate > a' % file)
       os.system("more a |  awk '{print $2} '  > aZH1")
       os.system("more a |  awk '{print $12}'  > aZH2")
       os.system("more a |  awk '{print $22}'  > aZH3")
       os.system("more a |  awk '{print $32}'  > aZH4")
       os.system("more a |  awk '{print $42}'  > aZH5")
       os.system("more a |  awk '{print $52}'  > aZH6")
       os.system("more a |  awk '{print $3}'   > aWH1")
       os.system("more a |  awk '{print $13}'  > aWH2")
       os.system("more a |  awk '{print $23}'  > aWH3")
       os.system("more a |  awk '{print $33}'  > aWH4")
       os.system("more a |  awk '{print $43}'  > aWH5")
       os.system("more a |  awk '{print $53}'  > aWH6")
 
       fZH1 = open("aZH1", "r")
       fZH2 = open("aZH2", "r")
       fZH3 = open("aZH3", "r")
       fZH4 = open("aZH4", "r")
       fZH5 = open("aZH5", "r")
       fZH6 = open("aZH6", "r")
       fWH1 = open("aWH1", "r")
       fWH2 = open("aWH2", "r")
       fWH3 = open("aWH3", "r")
       fWH4 = open("aWH4", "r")
       fWH5 = open("aWH5", "r")
       fWH6 = open("aWH6", "r")

       sZH1 = fZH1.read().rstrip('\n')
       sZH2 = fZH2.read().rstrip('\n')
       sZH3 = fZH3.read().rstrip('\n')
       sZH4 = fZH4.read().rstrip('\n')
       sZH5 = fZH5.read().rstrip('\n')
       sZH6 = fZH6.read().rstrip('\n')
       sWH1 = fWH1.read().rstrip('\n')
       sWH2 = fWH2.read().rstrip('\n')
       sWH3 = fWH3.read().rstrip('\n')
       sWH4 = fWH4.read().rstrip('\n')
       sWH5 = fWH5.read().rstrip('\n')
       sWH6 = fWH6.read().rstrip('\n')


       ZH1 = float(sZH1)
       ZH2 = float(sZH2)
       ZH3 = float(sZH3)
       ZH4 = float(sZH4)
       ZH5 = float(sZH5)
       ZH6 = float(sZH6)

       WH1 = float(sWH1)
       WH2 = float(sWH2)
       WH3 = float(sWH3)
       WH4 = float(sWH4)
       WH5 = float(sWH5)
       WH6 = float(sWH6)


       bZH1 = (xSec[toMass]/xSec[fromMass])*ZH1
       bZH2 = (xSec[toMass]/xSec[fromMass])*ZH2
       bZH3 = (xSec[toMass]/xSec[fromMass])*ZH3
       bZH4 = (xSec[toMass]/xSec[fromMass])*ZH4
       bZH5 = (xSec[toMass]/xSec[fromMass])*ZH5
       bZH6 = (xSec[toMass]/xSec[fromMass])*ZH6

       bWH1 = (xSec[toMass]/xSec[fromMass])*WH1
       bWH2 = (xSec[toMass]/xSec[fromMass])*WH2
       bWH3 = (xSec[toMass]/xSec[fromMass])*WH3
       bWH4 = (xSec[toMass]/xSec[fromMass])*WH4
       bWH5 = (xSec[toMass]/xSec[fromMass])*WH5
       bWH6 = (xSec[toMass]/xSec[fromMass])*WH6


       os.system("sed 's/%s/%f/g' a > b" %(sZH1,bZH1))
       os.system("sed -i 's/%s/%f/g' b" % (sZH2,bZH2))
       os.system("sed -i 's/%s/%f/g' b" % (sZH3,bZH3))
       os.system("sed -i 's/%s/%f/g' b" % (sZH4,bZH4))
       os.system("sed -i 's/%s/%f/g' b" % (sZH5,bZH5))
       os.system("sed -i 's/%s/%f/g' b" % (sZH6,bZH6))

       os.system("sed -i 's/%s/%f/g' b" % (sWH1,bWH1))
       os.system("sed -i 's/%s/%f/g' b" % (sWH2,bWH2))
       os.system("sed -i 's/%s/%f/g' b" % (sWH3,bWH3))
       os.system("sed -i 's/%s/%f/g' b" % (sWH4,bWH4))
       os.system("sed -i 's/%s/%f/g' b" % (sWH5,bWH5))
       os.system("sed -i 's/%s/%f/g' b" % (sWH6,bWH6))


       fRO = open("a", "r")
       sRO = fRO.read().rstrip('\n')
       fRF = open("b", "r")
       sRF = fRF.read().rstrip('\n')
       os.system("sed 's/%s/%s/g' %s > %s" % (sRO,sRF,file,newcard))

       os.system("sed -i 's/vhbb_WS_BDT_M110_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M110_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M110_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M110_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M110_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M110_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
 
       os.system("sed -i 's/vhbb_WS_BDT_M115_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M115_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M115_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M115_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M115_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M115_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
 
       os.system("sed -i 's/vhbb_WS_BDT_M120_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M120_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M120_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M120_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M120_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M120_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
 
       os.system("sed -i 's/vhbb_WS_BDT_M125_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M125_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M125_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M125_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M125_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M125_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
 
       os.system("sed -i 's/vhbb_WS_BDT_M130_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M130_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M130_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M130_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M130_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M130_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
 
       os.system("sed -i 's/vhbb_WS_BDT_M135_WenLowPt_8TeV.root/vhbb_WenLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M135_WmnLowPt_8TeV.root/vhbb_WmnLowPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M135_WenHighPt_8TeV.root/vhbb_WenHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M135_WmnHighPt_8TeV.root/vhbb_WmnHighPt_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M135_WenHighPtLooseCSV_8TeV.root/vhbb_WenHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
       os.system("sed -i 's/vhbb_WS_BDT_M135_WmnHighPtLooseCSV_8TeV.root/vhbb_WmnHighPtLooseCSV_8TeV.root/g' %s " % (newcard))
 

       os.system("echo %s written" % (newcard))

Process(files[0],0,0);
Process(files[0],1,0);
Process(files[0],2,0);
Process(files[0],3,0);
Process(files[0],4,0);
  
for x in range(5,15):
  Process(files[1], x , 10);
for x in range(15,25):
  Process(files[2], x , 20);
for x in range(25,35):
  Process(files[3], x , 30);
for x in range(35,45):
  Process(files[4], x , 40);
for x in range(45,51):
  Process(files[5], x , 50);
for x in range(51,67):
  Process(files[3], x , 30);
Revision:	1.3
Committed:	Sat Nov 3 13:02:28 2012 UTC (12 years, 6 months ago) by madfish
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	HEAD
Changes since 1.2:	+144 -51 lines
Log Message:	* empty log message *
#	User	Rev	Content
1	madfish	1.1	#! /usr/bin/env python
2
3
4			import sys
5			import os
6			import commands
7			import string
8
9	madfish	1.2	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	madfish	1.1
11
12
13	madfish	1.2	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			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	madfish	1.1
19
20	madfish	1.2	xSec = xSec8WH
21	madfish	1.1
22	madfish	1.2	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"];
23	madfish	1.1
24
25			files=[
26	madfish	1.3	"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	madfish	1.1	]
33
34
35			def Process(file, toMass, fromMass):
36	madfish	1.3	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	madfish	1.1	os.system('more %s \| grep rate > a' % file)
45	madfish	1.3	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	madfish	1.1
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	madfish	1.3
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	madfish	1.1	os.system("echo %s written" % (newcard))
182
183			Process(files[0],0,0);
184			Process(files[0],1,0);
185			Process(files[0],2,0);
186			Process(files[0],3,0);
187			Process(files[0],4,0);
188
189			for x in range(5,15):
190			Process(files[1], x , 10);
191			for x in range(15,25):
192			Process(files[2], x , 20);
193			for x in range(25,35):
194			Process(files[3], x , 30);
195			for x in range(35,45):
196			Process(files[4], x , 40);
197			for x in range(45,51):
198			Process(files[5], x , 50);
199	madfish	1.2	for x in range(51,67):
200			Process(files[3], x , 30);