Python AnalysisCategoryDef.AnalysisCategoryの例

コード例 #1

ファイル: InclusiveTauPt30to50_NJETS_Category.py プロジェクト: aloeliger/DatacardCreator

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInLowTauPtCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary["TauPt"] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and theEventDictionary.eventDictionary["TauPt"] <= 50.0
            and not (theEventDictionary.eventDictionary['Njets'] == 0
                     and theEventDictionary.eventDictionary['DeltaR'] < 2.0)):
        return True
    else:
        return False
    return False


LowTauPtCategory = AnalysisCategoryDef.AnalysisCategory()
LowTauPtCategory.name = 'htt_NJ_mt_LowTauPt'
LowTauPtCategory.IsInCategory = IsInLowTauPtCategory
LowTauPtCategory.rollingVariable = 'Njets'
LowTauPtCategory.rollingBins = [0, 1, 2, 3, 4, 20]
LowTauPtCategory.reconstructionVariable = 'M_sv'
LowTauPtCategory.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 290.0
]
#LowTauPtCategory.reconstructionBins = [75.0,90.0,105.0,120.0,135.0,150.0,170.0,190.0,210,250.0]

コード例 #2

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInInclusiveCategory(theAnalysisCategory, theEventDictionary):
    isGoodEvent = False
    if (theEventDictionary.eventDictionary["TauPt"] > 30.0
            and theEventDictionary.eventDictionary["MT"] < 50.0):
        isGoodEvent = True
    return isGoodEvent


MTInclusive = AnalysisCategoryDef.AnalysisCategory()
MTInclusive.name = "mtt"
MTInclusive.IsInCategory = IsInInclusiveCategory
MTInclusive.rollingVariable = "TauPt"
MTInclusive.rollingBins = [0.0, 9000.0]
MTInclusive.reconstructionVariable = "M_sv"
MTInclusive.reconstructionBins = [
    50.0,
    60.0,
    70.0,
    80.0,
    90.0,
    100.0,
    110.0,
    120.0,
    130.0,
    140.0,
    150.0,
    160.0,

コード例 #3

ファイル: Boosted_1J_Config.py プロジェクト: samhiggie/DatacardCreator

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInBoostedOneJetCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary['TauPt'] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and (theEventDictionary.eventDictionary['Njets'] == 1.0 or
                 (theEventDictionary.eventDictionary['Njets'] >= 2.0
                  and theEventDictionary.eventDictionary['mjj'] <= 350.0))
            and theEventDictionary.eventDictionary['Njets'] == 1.0):
        return True
    else:
        return False
    return False


Boosted1J = AnalysisCategoryDef.AnalysisCategory()
Boosted1J.name = 'mt_boosted_1J'
Boosted1J.IsInCategory = IsInBoostedOneJetCategory
Boosted1J.rollingVariable = 'HiggsPt'
Boosted1J.rollingBins = [0.0, 60.0, 120.0, 200.0, 250.0, 300.0, 10000.0]
Boosted1J.reconstructionVariable = 'M_sv'
Boosted1J.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 9000.0
]

コード例 #4

ファイル: InclusiveTauPt50to70_HiggsPt_Category.py プロジェクト: aloeliger/DatacardCreator

def IsInIntermediateTauPtCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary["TauPt"] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and theEventDictionary.eventDictionary["TauPt"] >= 50.0
            and theEventDictionary.eventDictionary["TauPt"] <= 70.0
            and not (theEventDictionary.eventDictionary['Njets'] == 0
                     and theEventDictionary.eventDictionary['DeltaR'] < 2.0)
            and
            not (abs(theEventDictionary.eventDictionary['TauEta']) > 0.2
                 and abs(theEventDictionary.eventDictionary['TauEta']) < 0.3)
            #and theEventDictionary.eventDictionary['MET'] < 100.0
        ):
        return True
    else:
        return False
    return False


IntermediateTauPtCategory = AnalysisCategoryDef.AnalysisCategory()
IntermediateTauPtCategory.name = 'htt_PTH_mt_IntermediateTauPt'
IntermediateTauPtCategory.IsInCategory = IsInIntermediateTauPtCategory
IntermediateTauPtCategory.rollingVariable = 'differentialHiggsPt'
IntermediateTauPtCategory.rollingBins = [0, 45, 80, 120, 200, 350, 100000000]
IntermediateTauPtCategory.reconstructionVariable = 'M_sv'
IntermediateTauPtCategory.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 290.0
]
#IntermediateTauPtCategory.reconstructionBins = [75.0,90.0,105.0,120.0,135.0,150.0,170.0,190.0,210,250.0]

コード例 #5

ファイル: VBF_PTH_0_200_Config.py プロジェクト: aloeliger/DatacardCreator

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInVBFLowCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary['Njets'] >= 2.0
            and theEventDictionary.eventDictionary['mjj'] > 350.0
            and theEventDictionary.eventDictionary['TauPt'] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and theEventDictionary.eventDictionary['HiggsPt'] < 200.0
            and theEventDictionary.eventDictionary['HiggsPt'] > 0.0
            and (theEventDictionary.eventDictionary['M_sv'] <= 90.0
                 or theEventDictionary.eventDictionary['M_sv'] >= 150.0)):
        return True
    else:
        return False
    return False


VBFLow = AnalysisCategoryDef.AnalysisCategory()
VBFLow.name = 'mt_vbf_PTH_0_200'
VBFLow.IsInCategory = IsInVBFLowCategory
VBFLow.rollingVariable = 'mjj'
VBFLow.rollingBins = [350.0, 700.0, 1000.0, 1500.0, 1800.0, 10000.0]
VBFLow.reconstructionVariable = 'M_sv'
VBFLow.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 290.0
]

コード例 #6

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInBoostedGTETwoJetCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary['TauPt'] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and (theEventDictionary.eventDictionary['Njets'] == 1.0 or
                 (theEventDictionary.eventDictionary['Njets'] >= 2.0
                  and theEventDictionary.eventDictionary['mjj'] <= 350.0))
            and theEventDictionary.eventDictionary['Njets'] >= 2.0):
        return True
    else:
        return False
    return False


BoostedGTE2J = AnalysisCategoryDef.AnalysisCategory()
BoostedGTE2J.name = 'mt_boosted_GE2J'
BoostedGTE2J.IsInCategory = IsInBoostedGTETwoJetCategory
BoostedGTE2J.rollingVariable = 'HiggsPt'
BoostedGTE2J.rollingBins = [0.0, 60.0, 120.0, 200.0, 250.0, 300.0, 10000.0]
BoostedGTE2J.reconstructionVariable = 'M_sv'
BoostedGTE2J.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 9000.0
]

コード例 #7

ファイル: ZeroJetCategory.py プロジェクト: aloeliger/DatacardCreator

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInZeroJetCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary['TauPt'] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and theEventDictionary.eventDictionary['Njets'] == 0
            and theEventDictionary.eventDictionary['DeltaR'] > 2.0):
        return True
    else:
        return False
    return False


ZeroJet = AnalysisCategoryDef.AnalysisCategory()
ZeroJet.name = "mt_0jet"
ZeroJet.IsInCategory = IsInZeroJetCategory
#ZeroJet.rollingVariable = 'HiggsPt'
#ZeroJet.rollingBins = [0,45,80,120,200,350,10000]
ZeroJet.rollingVariable = 'TauPt'
ZeroJet.rollingBins = [30.0, 40.0, 50.0, 10000.0]
ZeroJet.reconstructionVariable = 'M_sv'
ZeroJet.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 9000.0
]

コード例 #8

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInTwoJetLowCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary['TauPt'] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and theEventDictionary.eventDictionary['Njets'] == 2
            and theEventDictionary.eventDictionary['mjj'] < 450.0):
        return True
    else:
        return False
    return False


TwoJetLow = AnalysisCategoryDef.AnalysisCategory()
TwoJetLow.name = "mt_2jet_low"
TwoJetLow.IsInCategory = IsInTwoJetLowCategory
TwoJetLow.rollingVariable = 'HiggsPt'
TwoJetLow.rollingBins = [0, 45, 80, 120, 200, 350, 10000]
TwoJetLow.reconstructionVariable = 'M_sv'
TwoJetLow.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 9000.0
]

コード例 #9

ファイル: TwoJetHighMjjCategory.py プロジェクト: aloeliger/DatacardCreator

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInTwoJetHighCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary['TauPt'] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and theEventDictionary.eventDictionary['Njets'] == 2
            and theEventDictionary.eventDictionary['mjj'] > 450.0):
        return True
    else:
        return False
    return False


TwoJetHigh = AnalysisCategoryDef.AnalysisCategory()
TwoJetHigh.name = "mt_2jet_high"
TwoJetHigh.IsInCategory = IsInTwoJetHighCategory
TwoJetHigh.rollingVariable = 'HiggsPt'
TwoJetHigh.rollingBins = [0, 45, 80, 120, 200, 350, 10000]
TwoJetHigh.reconstructionVariable = 'M_sv'
TwoJetHigh.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 9000.0
]

コード例 #10

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInVBFHighCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary["Njets"] >= 2.0
            and theEventDictionary.eventDictionary['mjj'] > 350.0
            and theEventDictionary.eventDictionary['TauPt'] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and theEventDictionary.eventDictionary['HiggsPt'] >= 200.0
            and (theEventDictionary.eventDictionary['M_sv'] <= 90.0
                 or theEventDictionary.eventDictionary['M_sv'] >= 150.0)):
        return True
    else:
        return False
    return False


VBFHigh = AnalysisCategoryDef.AnalysisCategory()
VBFHigh.name = 'mt_vbf_PTH_GE_200'
VBFHigh.IsInCategory = IsInVBFHighCategory
VBFHigh.rollingVariable = 'mjj'
VBFHigh.rollingBins = [350.0, 700.0, 1200.0, 10000.0]
VBFHigh.reconstructionVariable = 'M_sv'
VBFHigh.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 290.0
]

コード例 #11

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInOneJetCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary['TauPt'] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and theEventDictionary.eventDictionary['Njets'] == 1):
        return True
    else:
        return False
    return False


OneJet = AnalysisCategoryDef.AnalysisCategory()
OneJet.name = "mt_1jet"
OneJet.IsInCategory = IsInOneJetCategory
OneJet.rollingVariable = 'HiggsPt'
OneJet.rollingBins = [0, 45, 80, 120, 200, 350, 10000]
OneJet.reconstructionVariable = 'M_sv'
OneJet.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 9000.0
]

コード例 #12

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInFourJetCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary['TauPt'] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and theEventDictionary.eventDictionary['Njets'] >= 4):
        return True
    else:
        return False
    return False


FourJet = AnalysisCategoryDef.AnalysisCategory()
FourJet.name = "mt_4jet"
FourJet.IsInCategory = IsInFourJetCategory
FourJet.rollingVariable = 'HiggsPt'
FourJet.rollingBins = [0, 45, 80, 120, 200, 350, 10000]
FourJet.reconstructionVariable = 'M_sv'
FourJet.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 9000.0
]

コード例 #13

import ROOT
import AnalysisCategories.AnalysisCategoryDef as AnalysisCategoryDef


def IsInThreeJetHighCategory(theAnalysisCategory, theEventDictionary):
    if (theEventDictionary.eventDictionary['TauPt'] >= 30.0
            and theEventDictionary.eventDictionary['MT'] < 50.0
            and theEventDictionary.eventDictionary['Njets'] >= 3):
        return True
    else:
        return False
    return False


ThreeJetHigh = AnalysisCategoryDef.AnalysisCategory()
ThreeJetHigh.name = "mt_3jet"
ThreeJetHigh.IsInCategory = IsInThreeJetHighCategory
ThreeJetHigh.rollingVariable = 'HiggsPt'
ThreeJetHigh.rollingBins = [0, 45, 80, 120, 200, 350, 10000]
ThreeJetHigh.reconstructionVariable = 'M_sv'
ThreeJetHigh.reconstructionBins = [
    50.0, 70.0, 90.0, 110.0, 130.0, 150.0, 170.0, 210.0, 250.0, 9000.0
]