def xsec_kfact_effic(self, isample):
# use yellowhiggs for cross sections
xs, _ = yellowhiggs.xsbr(self.energy, self.masses[isample],
Higgs.MODES_DICT[self.modes[isample]][0],
'tautau')
log.debug("{0} {1} {2} {3} {4} {5}".format(
self.samples[isample], self.masses[isample], self.modes[isample],
Higgs.MODES_DICT[self.modes[isample]][0], self.energy, xs))
if not self.inclusive_decays:
xs *= TAUTAUHADHADBR
kfact = 1.
effic = 1.
return xs, kfact, effic
def xsec_kfact_effic(self, isample):
# use yellowhiggs for cross sections
xs, _ = yellowhiggs.xsbr(
self.energy, self.masses[isample],
Higgs.MODES_DICT[self.modes[isample]][0], 'tautau')
log.debug("{0} {1} {2} {3} {4} {5}".format(
self.samples[isample],
self.masses[isample],
self.modes[isample],
Higgs.MODES_DICT[self.modes[isample]][0],
self.energy,
xs))
xs *= TAUTAUHADHADBR
kfact = 1.
effic = 1.
return xs, kfact, effic
ver),
'pu' : 'S10',
'calibrationTarget':'Summer12_DR53X_HCP2012',
'x_sec' : cs_hzg[mass][ch]*hbr_hzg[mass]*zbr_hzg[ch]*picobarns,
'analyses' : ['HZG']
}
# Add VH 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):
ver),
'pu' : 'S10',
'calibrationTarget':'Summer12_DR53X_HCP2012',
'x_sec' : cs_hzg[mass][ch]*hbr_hzg[mass]*zbr_hzg[ch]*picobarns,
'analyses' : ['HZG']
}
# Add VH 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):
-999,
}
# Add WH TauTau signal samples
for mass in range(90, 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,
}
datadefs['VH_H2Tau_M-%i' % mass] = {
'analyses': ['HTT', 'VH'],
'datasetpath':
'/WH_ZH_TTH_HToTauTau_M-%i_8TeV-pythia6-tauola/Summer12-PU_S7_START52_V9-v2/AODSIM'
% mass,
'pu':
'S7',
'calibrationTarget':
'Summer12',
'x_sec':
xsbr(8, mass, 'wh', 'tautau')[0] + xsbr(8, mass, 'zh', 'tautau')[0] +
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
## 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
datadefs["data_DoubleMu_Run2012B_13Jul2012_v4"] = {
"datasetpath": "/DoubleMu/Run2012B-13Jul2012-v4/AOD",
"lumi_mask": "FinalStateAnalysis/RecoTools/data/masks/Cert_190456-196531_8TeV_13Jul2012ReReco_Collisions12_JSON_v2.txt",
"firstRun": 193834,
"lastRun": 196531,
-999,
}
# Add WH TauTau signal samples
for mass in range(110, 165, 5):
datadefs['VH_H2Tau_M-%i' % mass] = {
'analyses': ['HTT', 'VH'],
'datasetpath':
'/WH_ZH_TTH_HToTauTau_M-%i_8TeV-pythia6-tauola/Summer12-PU_S7_START52_V9-v2/AODSIM'
% mass,
'pu':
'S7',
'calibrationTarget':
'Summer12',
'x_sec':
xsbr(8, mass, 'wh', 'tautau')[0] + xsbr(8, mass, 'zh', 'tautau')[0] +
xsbr(8, mass, 'tth', 'tautau')[0],
}
if mass == 110:
# Special case use v3 instead of v2, which doesn't exist
datadefs['VH_H2Tau_M-110']['datasetpath'] = datadefs['VH_H2Tau_M-110'][
'datasetpath'].replace('V9-v2', 'V9-v3')
## datadefs['VH_H2Tau_M-110']['x_sec'] = (1.060*br_w_leptons + 0.5869*br_z_leptons + 0.1887*square(br_w_leptons))*7.95E-02
## 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
# fix me
for mass in range(110, 150, 10):
datadefs['VHWW_lepdecay_%i' % mass] = {
datadefs['VBF_HZZ4L_M-%i' % mass] = {
'analyses': ['4L'],
'datasetpath': '/VBF_HToZZTo4L_M-%i_8TeV-powheg-pythia6/Summer12-PU_S7_START52_V9-v1/AODSIM' % mass,
'pu': 'S7',
'calibrationTarget': 'Summer12',
'x_sec': -999,
}
# Add WH TauTau signal samples
for mass in range(110, 165, 5):
datadefs['VH_H2Tau_M-%i' % mass] = {
'analyses': ['HTT', 'VH'],
'datasetpath': '/WH_ZH_TTH_HToTauTau_M-%i_8TeV-pythia6-tauola/Summer12-PU_S7_START52_V9-v2/AODSIM' % mass,
'pu': 'S7',
'calibrationTarget': 'Summer12',
'x_sec': xsbr(8,mass,'wh','tautau')[0] + xsbr(8,mass,'zh','tautau')[0] + xsbr(8,mass,'tth','tautau')[0],
}
if mass == 110:
# Special case use v3 instead of v2, which doesn't exist
datadefs['VH_H2Tau_M-110']['datasetpath'] = datadefs['VH_H2Tau_M-110']['datasetpath'].replace(
'V9-v2', 'V9-v3')
## datadefs['VH_H2Tau_M-110']['x_sec'] = (1.060*br_w_leptons + 0.5869*br_z_leptons + 0.1887*square(br_w_leptons))*7.95E-02
## 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
# fix me
for mass in range(110, 150, 10):
datadefs['VHWW_lepdecay_%i' % mass] = {
'analyses': ['VH'],
datadefs['VBF_HZZ4L_M-%i' % mass] = {
'analyses': ['4L'],
'datasetpath': '/VBF_HToZZTo4L_M-%i_8TeV-powheg-pythia6/Summer12-PU_S7_START52_V9-v1/AODSIM' % mass,
'pu': 'S7',
'calibrationTarget': 'Summer12',
'x_sec': -999,
}
# Add WH TauTau signal samples
for mass in range(90, 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,
}
datadefs['VH_H2Tau_M-%i' % mass] = {
'analyses': ['HTT', 'VH'],
'datasetpath': '/WH_ZH_TTH_HToTauTau_M-%i_8TeV-pythia6-tauola/Summer12-PU_S7_START52_V9-v2/AODSIM' % mass,
'pu': 'S7',
'calibrationTarget': 'Summer12',
'x_sec': xsbr(8,mass,'wh','tautau')[0] + xsbr(8,mass,'zh','tautau')[0] + xsbr(8,mass,'tth','tautau')[0],
}
if mass == 110:
# Special case use v3 instead of v2, which doesn't exist
datadefs['VH_H2Tau_M-110']['datasetpath'] = datadefs['VH_H2Tau_M-110']['datasetpath'].replace(
'V9-v2', 'V9-v3')
}
############################################################################
#### Signal datasets ####################################
############################################################################
datadefs['ggHiggsToETau'] = {'analyses': ['ETau'],'datasetpath': '/hdfs/store/user/mcepeda/LFV_GG_H2ETau-MuTauMC', 'pu': 'S10', 'calibrationTarget': 'Summer12_DR53X_HCP2012', 'x_sec': 0.1927 } #if BR=1%, 0.171503 if BR 0.89%
datadefs['vbfHiggsToETau'] = { 'analyses': ['ETau'],'datasetpath': '/hdfs/store/user/mcepeda/LFV_VBF_H2ETau-MuTauMC','pu': 'S10','calibrationTarget': 'Summer12_DR53X_HCP2012','x_sec': 0.01579} #if BR=1%, 0.0140531 if BR0.89%
datadefs['ggHiggsToMuTau'] = {'analyses': ['MuTau'],'datasetpath': '/hdfs/store/user/mcepeda/LFV_GG_H2MuTau_LONG-MuTauMC', 'pu': 'S10', 'calibrationTarget': 'Summer12_DR53X_HCP2012', 'x_sec': 0.1927 }
datadefs['vbfHiggsToMuTau'] = { 'analyses': ['MuTau'],'datasetpath': '/hdfs/store/user/mcepeda/LFV_VBF_H2MuTau_LONG-MuTauMC','pu': 'S10','calibrationTarget': 'Summer12_DR53X_HCP2012','x_sec':0.01579}
datadefs['GluGluToHToTauTau_M-125_8TeV-powheg-pythia6'] = {
'analyses': ['HTT'],
'datasetpath': "/GluGluToHToTauTau_M-125_8TeV-powheg-pythia6/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM",
'pu' : 'S10',
'x_sec' : xsbr(8,125,'ggf','tautau')[0],# 19.27 ,
'calibrationTarget':'Summer12_DR53X_HCP2012'
}
datadefs['VBF_HToTauTau_M-125_8TeV-powheg-pythia6'] = {
'analyses': ['HTT'],
'datasetpath': "/VBF_HToTauTau_M-%i_8TeV-powheg-pythia6/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM",
'pu' : 'S10',
'x_sec' : xsbr(8,125,'vbf','tautau')[0],#1.578, # to be checked
'calibrationTarget':'Summer12_DR53X_HCP2012'
}
for mass in range(80,150, 10) + range(160, 220, 20) + range(250, 550, 50) + range(600, 1100, 100) :
datadefs['SUSYGluGluToHToTauTau_M-%i_8TeV-pythia6-tauola' % mass] = {
'analyses': ['HTT'],
'datasetpath': "/SUSYGluGluToHToTauTau_M-%i_8TeV-pythia6-tauola/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM" % mass,