def xsec_kfact_effic(self, isample):
# use yellowhiggs for cross sections
xs, _ = yellowhiggs.xs(self.energy, self.masses[isample],
self.modes[isample])
log.debug("{0} {1} {2} {3} {4}".format(self.samples[isample],
self.masses[isample],
self.modes[isample],
self.energy, xs))
kfact = 1.
effic = 1.
return xs, kfact, effic
def xsec_kfact_effic(self, isample):
# use yellowhiggs for cross sections
xs, _ = yellowhiggs.xs(
self.energy, self.masses[isample], self.modes[isample])
log.debug("{0} {1} {2} {3} {4}".format(
self.samples[isample],
self.masses[isample],
self.modes[isample],
self.energy,
xs))
kfact = 1.
effic = 1.
return xs, kfact, effic
for mass in range(110, 165, 5):
datadefs['VHtautau_lepdecay_%i' % mass] = {
'analyses': ['VH'],
'datasetpath': "/WH_ZH_TTH_HToTauTau_M-%i_lepdecay_8TeV-pythia6-tauola/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM" % mass,
'pu' : 'S10',
'calibrationTarget':'Summer12_DR53X_HCP2012',
'x_sec' : xsbr(8,mass,'wh','tautau')[0]*br_w_leptons + xsbr(8,mass,'zh','tautau')[0]*br_z_leptons + xsbr(8,mass,'tth','tautau')[0]*br_w_leptons**2,
}
for mass in range(110, 150, 10):
datadefs['VHWW_lepdecay_%i' % mass] = {
'analyses': ['VH'],
'datasetpath': "/WH_ZH_TTH_HToWW_M-%i_lepdecay_8TeV-pythia6/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM" % mass,
'pu' : 'S10',
'calibrationTarget':'Summer12_DR53X_HCP2012',
'x_sec' : (xs(8,mass,'wh')[0]*br_w_leptons+xs(8,mass,'zh')[0]*br_z_leptons+xs(8,mass,'tth')[0])*br(mass,'WW')*br_w_leptons**2,
}
# We use the same name for the 53X lepdecay only samples (sigh)
for mass in range(90, 165, 5):
datadefs['VH_H2Tau_M-%s' % mass] = {}
datadefs['VH_H2Tau_M-%s' % mass]['x_sec'] = xsbr(8,mass,'wh','tautau')[0] + xsbr(8,mass,'zh','tautau')[0] + xsbr(8,mass,'tth','tautau')[0]
## datadefs['VH_H2Tau_M-120']['x_sec'] = (0.7966*br_w_leptons + 0.4483*br_z_leptons + 0.1470*square(br_w_leptons))*7.04E-02
## datadefs['VH_H2Tau_M-130']['x_sec'] = (0.6095*br_w_leptons + 0.3473*br_z_leptons + 0.1157*square(br_w_leptons))*5.48E-02
## datadefs['VH_H2Tau_M-140']['x_sec'] = (0.4713*br_w_leptons + 0.2728*br_z_leptons + 0.09207*square(br_w_leptons))*3.54E-02
# Add data files
for mass in range(110, 165, 5):
datadefs['VHtautau_lepdecay_%i' % mass] = {
'analyses': ['VH'],
'datasetpath': "/WH_ZH_TTH_HToTauTau_M-%i_lepdecay_8TeV-pythia6-tauola/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM" % mass,
'pu' : 'S10',
'calibrationTarget':'Summer12_DR53X_HCP2012',
'x_sec' : xsbr(8,mass,'wh','tautau')[0]*br_w_leptons + xsbr(8,mass,'zh','tautau')[0]*br_z_leptons + xsbr(8,mass,'tth','tautau')[0]*br_w_leptons**2,
}
for mass in range(110, 150, 10):
datadefs['VHWW_lepdecay_%i' % mass] = {
'analyses': ['VH'],
'datasetpath': "/WH_ZH_TTH_HToWW_M-%i_lepdecay_8TeV-pythia6/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM" % mass,
'pu' : 'S10',
'calibrationTarget':'Summer12_DR53X_HCP2012',
'x_sec' : (xs(8,mass,'wh')[0]*br_w_leptons+xs(8,mass,'zh')[0]*br_z_leptons+xs(8,mass,'tth')[0])*br(mass,'WW')*br_w_leptons**2,
}
# We use the same name for the 53X lepdecay only samples (sigh)
for mass in range(90, 165, 5):
datadefs['VH_H2Tau_M-%s' % mass] = {}
datadefs['VH_H2Tau_M-%s' % mass]['x_sec'] = xsbr(8,mass,'wh','tautau')[0] + xsbr(8,mass,'zh','tautau')[0] + xsbr(8,mass,'tth','tautau')[0]
## datadefs['VH_H2Tau_M-120']['x_sec'] = (0.7966*br_w_leptons + 0.4483*br_z_leptons + 0.1470*square(br_w_leptons))*7.04E-02
## datadefs['VH_H2Tau_M-130']['x_sec'] = (0.6095*br_w_leptons + 0.3473*br_z_leptons + 0.1157*square(br_w_leptons))*5.48E-02
## datadefs['VH_H2Tau_M-140']['x_sec'] = (0.4713*br_w_leptons + 0.2728*br_z_leptons + 0.09207*square(br_w_leptons))*3.54E-02
# Add data files
xsbr(8, mass, 'zh', 'tautau')[0] * br_z_leptons +
xsbr(8, mass, 'tth', 'tautau')[0] * br_w_leptons**2,
}
for mass in range(110, 150, 10):
datadefs['VHWW_lepdecay_%i' % mass] = {
'analyses': ['VH'],
'datasetpath':
"/WH_ZH_TTH_HToWW_M-%i_lepdecay_8TeV-pythia6/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM"
% mass,
'pu':
'S10',
'calibrationTarget':
'Summer12_DR53X_HCP2012',
'x_sec':
(xs(8, mass, 'wh')[0] * br_w_leptons +
xs(8, mass, 'zh')[0] * br_z_leptons + xs(8, mass, 'tth')[0]) *
br(mass, 'WW') * br_w_leptons**2,
}
# We use the same name for the 53X lepdecay only samples (sigh)
for mass in range(90, 165, 5):
datadefs['VH_H2Tau_M-%s' % mass] = {}
datadefs['VH_H2Tau_M-%s' % mass]['x_sec'] = xsbr(
8, mass, 'wh', 'tautau')[0] + xsbr(8, mass, 'zh', 'tautau')[0] + xsbr(
8, mass, 'tth', 'tautau')[0]
## datadefs['VH_H2Tau_M-120']['x_sec'] = (0.7966*br_w_leptons + 0.4483*br_z_leptons + 0.1470*square(br_w_leptons))*7.04E-02
## datadefs['VH_H2Tau_M-130']['x_sec'] = (0.6095*br_w_leptons + 0.3473*br_z_leptons + 0.1157*square(br_w_leptons))*5.48E-02
## datadefs['VH_H2Tau_M-140']['x_sec'] = (0.4713*br_w_leptons + 0.2728*br_z_leptons + 0.09207*square(br_w_leptons))*3.54E-02
#datadefs['VH_H2Tau_M-110']['x_sec'] = (1.060 + 0.5869 + 0.1887)*7.95E-02
#datadefs['VH_H2Tau_M-120']['x_sec'] = (0.7966 + 0.4483 + 0.1470)*7.04E-02
#datadefs['VH_H2Tau_M-130']['x_sec'] = (0.6095 + 0.3473 + 0.1157)*5.48E-02
#datadefs['VH_H2Tau_M-140']['x_sec'] = (0.4713 + 0.2728 + 0.09207)*3.54E-02
# Add the cross sections for WH->HWW samples. We use the 7TeV ones here,
# and then just change the xsec.
datadefs['WH_110_HWW3l'] = {'x_sec': 1.060 * cube(br_w_leptons) * 4.82E-02}
datadefs['WH_120_HWW3l'] = {'x_sec': 0.7966 * cube(br_w_leptons) * 1.43E-01}
datadefs['WH_130_HWW3l'] = {'x_sec': 0.6095 * cube(br_w_leptons) * 3.05E-01}
datadefs['WH_140_HWW3l'] = {'x_sec': 0.4713 * cube(br_w_leptons) * 5.03E-01}
#VH->HWW xsec: WH + ZH; ZH --> totalxsec * BR(ZtoLL) * BR(HtoWW) * BR( WtoLL )^2
for mass in range(90, 150, 10):
datadefs['VH_%s_HWW' % mass] = {
'x_sec': (xs(8, mass, 'wh')[0] * br_w_leptons + xs(8, mass, 'zh')[0] *
br_z_leptons + xs(8, mass, 'tth')[0] * br_w_leptons**2) *
br(mass, 'WW') * br_w_leptons**2
}
# Add data files
def build_data_set(pd, analyses):
subsample_dict = {
'data_%s_Run2012A_PromptReco_v1' % pd: {
'datasetpath': "/%s/Run2012A-PromptReco-v1/AOD" % pd,
'lumi_mask':
"FinalStateAnalysis/RecoTools/data/masks/Cert_190456-194479_8TeV_PromptReco_Collisions12_JSON.txt",
'firstRun': 190450,
'lastRun': 193686,
'analyses': analyses,
datadefs["VHtautau_lepdecay_%i" % mass] = {
"analyses": ["VH"],
"datasetpath": "/WH_ZH_TTH_HToTauTau_M-%i_lepdecay_8TeV-pythia6-tauola/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM"
% mass,
"pu": "S10",
"calibrationTarget": "Summer12_DR53X_HCP2012",
"x_sec": -999,
}
for mass in range(110, 150, 10):
datadefs["VHWW_lepdecay_%i" % mass] = {
"analyses": ["VH"],
"datasetpath": "/WH_ZH_TTH_HToWW_M-%i_lepdecay_8TeV-pythia6/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM" % mass,
"pu": "S10",
"calibrationTarget": "Summer12_DR53X_HCP2012",
"x_sec": (xs(8, mass, "wh")[0] * br_w_leptons + xs(8, mass, "zh")[0] * br_z_leptons + xs(8, mass, "tth")[0])
* br(mass, "WW")
* br_w_leptons ** 2,
}
# We use the same name for the 53X lepdecay only samples (sigh)
for mass in range(110, 160, 5):
datadefs["VH_H2Tau_M-%s" % mass] = {}
datadefs["VH_H2Tau_M-%s" % mass]["x_sec"] = (
xsbr(8, mass, "wh", "tautau")[0] * br_w_leptons
+ xsbr(8, mass, "zh", "tautau")[0] * br_z_leptons
+ xsbr(8, mass, "tth", "tautau")[0] * br_w_leptons ** 2
)
## datadefs['VH_H2Tau_M-120']['x_sec'] = (0.7966*br_w_leptons + 0.4483*br_z_leptons + 0.1470*square(br_w_leptons))*7.04E-02
# Add the only one we are currently interested int
#datadefs['VH_H2Tau_M-110']['x_sec'] = (1.060 + 0.5869 + 0.1887)*7.95E-02
#datadefs['VH_H2Tau_M-120']['x_sec'] = (0.7966 + 0.4483 + 0.1470)*7.04E-02
#datadefs['VH_H2Tau_M-130']['x_sec'] = (0.6095 + 0.3473 + 0.1157)*5.48E-02
#datadefs['VH_H2Tau_M-140']['x_sec'] = (0.4713 + 0.2728 + 0.09207)*3.54E-02
# Add the cross sections for WH->HWW samples. We use the 7TeV ones here,
# and then just change the xsec.
datadefs['WH_110_HWW3l'] = { 'x_sec' : 1.060*cube(br_w_leptons)* 4.82E-02 }
datadefs['WH_120_HWW3l'] = { 'x_sec' : 0.7966*cube(br_w_leptons)*1.43E-01 }
datadefs['WH_130_HWW3l'] = { 'x_sec' : 0.6095*cube(br_w_leptons)*3.05E-01 }
datadefs['WH_140_HWW3l'] = { 'x_sec' : 0.4713*cube(br_w_leptons)*5.03E-01 }
#VH->HWW xsec: WH + ZH; ZH --> totalxsec * BR(ZtoLL) * BR(HtoWW) * BR( WtoLL )^2
for mass in range(110, 150, 10):
datadefs['VH_%s_HWW' % mass] = {'x_sec' : (xs(8,mass,'wh')[0]*br_w_leptons+xs(8,mass,'zh')[0]*br_z_leptons+xs(8,mass,'tth')[0])*br(mass,'WW')*br_w_leptons**2}
# Add data files
def build_data_set(pd, analyses):
subsample_dict = {
'data_%s_Run2012A_PromptReco_v1' % pd : {
'datasetpath' : "/%s/Run2012A-PromptReco-v1/AOD" % pd,
'lumi_mask' : "FinalStateAnalysis/RecoTools/data/masks/Cert_190456-194479_8TeV_PromptReco_Collisions12_JSON.txt",
'firstRun' : 190450,
'lastRun' : 193686,
'analyses' : analyses,
'calibrationTarget':'ICHEP2012'
},
'data_%s_Run2012B_PromptReco_v1_a' % pd : {
'datasetpath' : "/%s/Run2012B-PromptReco-v1/AOD" % pd,