-
Notifications
You must be signed in to change notification settings - Fork 10
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge pull request #85 from kreczko/config_changes
Config changes
- Loading branch information
Showing
41 changed files
with
1,090 additions
and
2,610 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,262 @@ | ||
| from __future__ import division | ||
|
|
||
| __all__ = [ | ||
| 'XSectionConfig', | ||
| ] | ||
|
|
||
| class XSectionConfig(): | ||
| current_analysis_path = '/storage/TopQuarkGroup/results/histogramfiles/AN-14-071_2nd_draft/' | ||
| known_centre_of_mass_energies = [7, 8] | ||
| # has to be separate as many variables depend on it | ||
| luminosities = {7:5050, 8:19584} | ||
| parameters = ['SingleTop_category_templates', 'SingleTop_file', | ||
| 'VJets_category_templates', 'analysis_types', | ||
| 'categories_and_prefixes', 'central_general_template', | ||
| 'centre_of_mass_energy', 'current_analysis_path', | ||
| 'data_file_electron', 'data_file_muon', | ||
| 'data_muon_category_templates', 'electron_QCD_MC_file', | ||
| 'electron_control_region', | ||
| 'electron_control_region_systematic', | ||
| 'general_category_templates', | ||
| 'generator_systematic_ttbar_templates', | ||
| 'generator_systematic_vjets_templates', | ||
| 'generator_systematics', | ||
| 'higgs_category_templates', 'higgs_file', | ||
| 'histogram_path_templates', 'include_higgs', | ||
| 'k_values_combined', 'k_values_electron', 'k_values_muon', | ||
| 'known_centre_of_mass_energies', 'luminosities', | ||
| 'luminosity', 'luminosity_scale', 'met_systematics_suffixes', | ||
| 'muon_QCD_MC_category_templates', 'muon_QCD_MC_file', | ||
| 'muon_QCD_file', 'muon_control_region', | ||
| 'muon_control_region_systematic', 'new_luminosity', | ||
| 'parameters', 'path_to_files', 'path_to_unfolding_histograms', | ||
| 'pdf_uncertainty_template', 'rate_changing_systematics', | ||
| 'rebin', 'special_muon_histogram', 'translate_options', | ||
| 'ttbar_category_templates', 'ttbar_generator_systematics', | ||
| 'ttbar_theory_systematic_prefix', 'ttbar_xsection', | ||
| 'unfolding_madgraph', 'unfolding_madgraph_raw', | ||
| 'unfolding_matching_down', 'unfolding_matching_down_raw', | ||
| 'unfolding_matching_up', 'unfolding_matching_up_raw', | ||
| 'unfolding_mcatnlo', 'unfolding_mcatnlo_raw', | ||
| 'unfolding_powheg', 'unfolding_powheg_raw', | ||
| 'unfolding_scale_down', 'unfolding_scale_down_raw', | ||
| 'unfolding_scale_up', 'unfolding_scale_up_raw', | ||
| 'vjets_theory_systematic_prefix' | ||
| ] | ||
|
|
||
| def __init__( self, centre_of_mass_energy ): | ||
| if not centre_of_mass_energy in self.known_centre_of_mass_energies: | ||
| raise AttributeError( 'Uknown centre of mass energy' ) | ||
| self.centre_of_mass_energy = centre_of_mass_energy | ||
| self.__fill_defaults__() | ||
|
|
||
| def __fill_defaults__( self ): | ||
| self.path_to_files = self.current_analysis_path + str( self.centre_of_mass_energy ) + 'TeV/' | ||
| self.path_to_unfolding_histograms = self.path_to_files + 'unfolding/' | ||
| path_to_files = self.path_to_files | ||
| path_to_unfolding_histograms = self.path_to_unfolding_histograms | ||
|
|
||
| self.luminosity = self.luminosities[self.centre_of_mass_energy] | ||
|
|
||
| # general | ||
| self.met_systematics_suffixes = [ | ||
| "ElectronEnUp", | ||
| "ElectronEnDown", | ||
| "MuonEnUp", | ||
| "MuonEnDown", | ||
| "TauEnUp", | ||
| "TauEnDown", | ||
| "JetResUp", | ||
| "JetResDown", | ||
| "JetEnUp", | ||
| "JetEnDown", | ||
| "UnclusteredEnUp", | ||
| "UnclusteredEnDown" | ||
| ] | ||
|
|
||
| self.analysis_types = { | ||
| 'electron':'EPlusJets', | ||
| 'muon':'MuPlusJets' | ||
| } | ||
|
|
||
| # measurement script options | ||
| self.translate_options = { | ||
| '0':'0btag', | ||
| '1':'1btag', | ||
| '2':'2btags', | ||
| '3':'3btags', | ||
| '0m':'0orMoreBtag', | ||
| '1m':'1orMoreBtag', | ||
| '2m':'2orMoreBtags', | ||
| '3m':'3orMoreBtags', | ||
| '4m':'4orMoreBtags', | ||
| # mettype: | ||
| 'pf':'PFMET', | ||
| 'type1':'patType1CorrectedPFMet', | ||
| } | ||
|
|
||
| self.ttbar_theory_systematic_prefix = 'TTJets_' | ||
| self.vjets_theory_systematic_prefix = 'VJets_' | ||
| # files | ||
| self.middle = '_' + str( self.luminosity ) + 'pb_PFElectron_PFMuon_PF2PATJets_PFMET' | ||
| middle = self.middle | ||
|
|
||
| self.data_file_muon = path_to_files + 'central/SingleMu' + middle + '.root' | ||
| self.muon_QCD_file = path_to_files + 'QCD_data_mu.root' | ||
| self.SingleTop_file = path_to_files + 'central/SingleTop' + middle + '.root' | ||
| self.electron_QCD_MC_file = path_to_files + 'central/QCD_Electron' + middle + '.root' | ||
| self.muon_QCD_MC_file = path_to_files + 'central/QCD_Muon' + middle + '.root' | ||
| self.higgs_file = path_to_files + 'central/TTH_Inclusive_M-125' + middle + '.root' | ||
|
|
||
| self.categories_and_prefixes = { | ||
| 'central':'', | ||
| 'Electron_down':'_minusElectron', | ||
| 'Electron_up':'_plusElectron', | ||
| 'Muon_down':'_minusMuon', | ||
| 'Muon_up':'_plusMuon', | ||
| 'BJet_down':'_minusBJet', | ||
| 'BJet_up':'_plusBjet', | ||
| 'JES_down':'_minusJES', | ||
| 'JES_up':'_plusJES', | ||
| 'JER_down':'_minusJER', | ||
| 'JER_up':'_plusJER', | ||
| 'LightJet_down':'_minusLightJet', | ||
| 'LightJet_up':'_plusLightJet', | ||
| } | ||
|
|
||
| # now fill in the centre of mass dependent values | ||
| # this position is crucial | ||
| if self.centre_of_mass_energy == 7: | ||
| self.__fill_defaults_7TeV__() | ||
| if self.centre_of_mass_energy == 8: | ||
| self.__fill_defaults_8TeV__() | ||
|
|
||
| self.generator_systematics = [ 'matchingup', 'matchingdown', 'scaleup', 'scaledown'] | ||
| self.ttbar_generator_systematics = [ 'matchingup', 'matchingdown', 'scaleup', 'scaledown', 'mcatnlo', 'ptreweight'] | ||
| self.central_general_template = path_to_files + 'central/%s' + middle + '.root' | ||
| self.generator_systematic_ttbar_templates = { systematic: path_to_files + 'central/TTJets-%s_%dpb_PFElectron_PFMuon_PF2PATJets_PFMET.root' % ( systematic, self.luminosity ) for systematic in self.ttbar_generator_systematics} | ||
| self.generator_systematic_vjets_templates = { systematic: path_to_files + 'central/VJets-%s_%dpb_PFElectron_PFMuon_PF2PATJets_PFMET.root' % ( systematic, self.luminosity ) for systematic in self.generator_systematics} | ||
| self.pdf_uncertainty_template = path_to_files + 'PDFWeights/TTJet' + middle + '_PDFWeights_%d.root' | ||
|
|
||
| categories_and_prefixes = self.categories_and_prefixes | ||
|
|
||
| self.general_category_templates = {category: path_to_files + category + '/%s' + middle + prefix + '.root' for category, prefix in categories_and_prefixes.iteritems()} | ||
| self.ttbar_category_templates = {category: path_to_files + category + '/TTJet' + middle + prefix + '.root' for category, prefix in categories_and_prefixes.iteritems()} | ||
| self.SingleTop_category_templates = {category: path_to_files + category + '/SingleTop' + middle + prefix + '.root' for ( category, prefix ) in categories_and_prefixes.iteritems()} | ||
| self.VJets_category_templates = {category: path_to_files + category + '/VJets' + middle + prefix + '.root' for ( category, prefix ) in categories_and_prefixes.iteritems()} | ||
| self.higgs_category_templates = {category: path_to_files + category + '/TTH_Inclusive_M-125' + middle + prefix + '.root' for ( category, prefix ) in categories_and_prefixes.iteritems()} | ||
| self.muon_QCD_MC_category_templates = {category: path_to_files + category + '/QCD_Muon' + middle + prefix + '.root' for ( category, prefix ) in categories_and_prefixes.iteritems()} | ||
|
|
||
| self.data_muon_category_templates = { | ||
| 'central': self.data_file_muon, | ||
| 'JES_up': path_to_files + 'JES_up/SingleMu' + middle + self.categories_and_prefixes['JES_up'] + '.root', | ||
| 'JES_down': path_to_files + 'JES_down/SingleMu' + middle + self.categories_and_prefixes['JES_down'] + '.root' | ||
| } | ||
|
|
||
| self.unfolding_madgraph_raw = path_to_unfolding_histograms + 'unfolding_merged.root' | ||
| self.unfolding_powheg_raw = path_to_unfolding_histograms + 'unfolding_TTJets_%dTeV_powheg.root' % self.centre_of_mass_energy | ||
| self.unfolding_mcatnlo_raw = path_to_unfolding_histograms + 'unfolding_TTJets_%dTeV_mcatnlo.root' % self.centre_of_mass_energy | ||
|
|
||
| self.unfolding_scale_down_raw = path_to_unfolding_histograms + 'unfolding_TTJets_%dTeV_scaledown.root' % self.centre_of_mass_energy | ||
| self.unfolding_scale_up_raw = path_to_unfolding_histograms + 'unfolding_TTJets_%dTeV_scaleup.root' % self.centre_of_mass_energy | ||
| self.unfolding_matching_down_raw = path_to_unfolding_histograms + 'unfolding_TTJets_%dTeV_matchingdown.root' % self.centre_of_mass_energy | ||
| self.unfolding_matching_up_raw = path_to_unfolding_histograms + 'unfolding_TTJets_%dTeV_matchingup.root' % self.centre_of_mass_energy | ||
|
|
||
| self.unfolding_madgraph = self.unfolding_madgraph_raw.replace( '.root', '_asymmetric.root' ) | ||
| self.unfolding_powheg = self.unfolding_powheg_raw.replace( '.root', '_asymmetric.root' ) | ||
| self.unfolding_mcatnlo = self.unfolding_mcatnlo_raw.replace( '.root', '_asymmetric.root' ) | ||
|
|
||
| self.unfolding_scale_down = self.unfolding_scale_down_raw.replace( '.root', '_asymmetric.root' ) | ||
| self.unfolding_scale_up = self.unfolding_scale_up_raw.replace( '.root', '_asymmetric.root' ) | ||
| self.unfolding_matching_down = self.unfolding_matching_down_raw.replace( '.root', '_asymmetric.root' ) | ||
| self.unfolding_matching_up = self.unfolding_matching_up_raw.replace( '.root', '_asymmetric.root' ) | ||
|
|
||
| self.histogram_path_templates = { | ||
| 'MET' : 'TTbar_plus_X_analysis/%s/Ref selection/Binned_MET_Analysis/%s_bin_%s/%s_absolute_eta', | ||
| 'HT' : 'TTbar_plus_X_analysis/%s/Ref selection/Binned_HT_Analysis/HT_bin_%s/%s_absolute_eta', | ||
| 'ST': 'TTbar_plus_X_analysis/%s/Ref selection/Binned_ST_Analysis/ST_with_%s_bin_%s/%s_absolute_eta', | ||
| 'MT': 'TTbar_plus_X_analysis/%s/Ref selection/Binned_MT_Analysis/MT_with_%s_bin_%s/%s_absolute_eta', | ||
| 'WPT': 'TTbar_plus_X_analysis/%s/Ref selection/Binned_WPT_Analysis/WPT_with_%s_bin_%s/%s_absolute_eta' | ||
| } | ||
|
|
||
| # optimal regularisation parameters | ||
| self.k_values_electron = { | ||
| 'MET' : 3, | ||
| 'HT' : 7, | ||
| 'ST' : 6, | ||
| 'MT' : 2, | ||
| 'WPT' : 3 | ||
| } | ||
|
|
||
| self.k_values_muon = { | ||
| 'MET' : 3, | ||
| 'HT' : 6, | ||
| 'ST' : 6, | ||
| 'MT' : 2, | ||
| 'WPT' : 6 | ||
| } | ||
|
|
||
| self.k_values_combined = { | ||
| 'MET' : 3, | ||
| 'HT' : 3, | ||
| 'ST' : 3, | ||
| 'MT' : 2, | ||
| 'WPT' : 4 | ||
| } | ||
|
|
||
| self.electron_control_region = 'QCDConversions' | ||
| self.electron_control_region_systematic = 'QCD non iso e+jets' | ||
|
|
||
| self.muon_control_region = 'QCD non iso mu+jets ge3j' | ||
| self.muon_control_region_systematic = 'QCD non iso mu+jets ge3j' # no systematic yet | ||
| self.rebin = 2 | ||
|
|
||
| self.include_higgs = False | ||
|
|
||
| self.luminosity_scale = self.new_luminosity / self.luminosity | ||
|
|
||
| def __fill_defaults_7TeV__( self ): | ||
| middle = self.middle | ||
| path_to_files = self.path_to_files | ||
|
|
||
| self.new_luminosity = self.luminosity # pb^-1 | ||
| self.ttbar_xsection = 164 # pb | ||
|
|
||
| self.data_file_electron = path_to_files + 'central/ElectronHad' + middle + '.root' | ||
| self.rate_changing_systematics = { | ||
| 'luminosity': 0.022, # https://twiki.cern.ch/twiki/bin/viewauth/CMS/PileupSystematicErrors | ||
| 'SingleTop_cross_section': 0.3, | ||
| 'TTJet_cross_section': 0.15 | ||
| } | ||
| self.categories_and_prefixes['PU_down'] = '_PU_64600mb' | ||
| self.categories_and_prefixes['PU_up'] = '_PU_71400mb' | ||
|
|
||
| self.special_muon_histogram = 'etaAbs_ge2j_data' | ||
|
|
||
| self.data_electron_category_templates = {'central': self.data_file_electron, | ||
| 'JES_up': path_to_files + 'JES_up/ElectronHad' + middle + self.categories_and_prefixes['JES_up'] + '.root', | ||
| 'JES_down': path_to_files + 'JES_down/ElectronHad' + middle + self.categories_and_prefixes['JES_down'] + '.root' | ||
| } | ||
|
|
||
| def __fill_defaults_8TeV__( self ): | ||
| middle = self.middle | ||
| path_to_files = self.path_to_files | ||
|
|
||
| self.new_luminosity = self.luminosity # pb^-1 | ||
| self.ttbar_xsection = 245.8 # pb | ||
|
|
||
| self.data_file_electron = path_to_files + 'central/SingleElectron' + middle + '.root' | ||
| self.rate_changing_systematics = { | ||
| 'luminosity': 0.026, # https://hypernews.cern.ch/HyperNews/CMS/get/physics-announcements/2526.html | ||
| 'SingleTop_cross_section': 0.034, # https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSectionsat8TeV | ||
| 'TTJet_cross_section': 0.043 | ||
| } | ||
| self.categories_and_prefixes['PU_down'] = '_PU_65835mb' | ||
| self.categories_and_prefixes['PU_up'] = '_PU_72765mb' | ||
|
|
||
| self.special_muon_histogram = 'muon_AbsEta_0btag' | ||
|
|
||
| self.data_electron_category_templates = {'central': self.data_file_electron, | ||
| 'JES_up': path_to_files + 'JES_up/SingleElectron' + middle + self.categories_and_prefixes['JES_up'] + '.root', | ||
| 'JES_down': path_to_files + 'JES_down/SingleElectron' + middle + self.categories_and_prefixes['JES_down'] + '.root' | ||
| } |
Oops, something went wrong.