Beispiel #1
0
 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
Beispiel #2
0
 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):
Beispiel #5
0
        -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,
Beispiel #7
0
        -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,