def main():
PreselectionRunner17 = ScriptRunner(year='2017', analysisname=analysisname, t3_base_path=t3_base_path, result_base_path=result_base_path, selectionstage='Studies', selectionname='LQDMHLTCheck', backgrounds=['LQLQToBTau_MLQ1000_MPS117_MC1100_L1p0', 'LQLQToBTauPsiChi_MLQ1000_MPS117_MC1100_L1p0', 'LQLQToPsiChi_MLQ1000_MPS117_MC1100_L1p0', 'PsiPsiToLQChi_MLQ1000_MPS2221_MC12089_L1p0'], signals=[])
# PreselectionRunner17.CalculateSelectionEfficiencies(num_and_denom=[('singlemu', 'input'), ('doublemu', 'input'), ('mutau', 'input'), ('singleele', 'input'), ('eletau', 'input'), ('elejet', 'input'), ('ditau', 'input'), ('met', 'input'), ('triggercategories', 'input')])
PreselectionRunner17.CalculateSelectionEfficiencies(num_and_denom=[('muon', 'input'), ('electron', 'input'), ('tau', 'input'), ('MET', 'input'), ('combined', 'input')])
def main():
# print 'hello1'
# StudiesRunner17 = ScriptRunner(year='2017', analysisname=analysisname, t3_base_path=t3_base_path, result_base_path=result_base_path, selectionstage='Studies', selectionname='PsiPsi_LeptonIDEfficiencies', backgrounds=[], signals=['PsiPsiToLQChi_MLQ1000_MPS117_MC1100_L1p0', 'PsiPsiToLQChi_MLQ1000_MPS244_MC1214_L1p0', 'PsiPsiToLQChi_MLQ1000_MPS567_MC1457_L1p0', 'PsiPsiToLQChi_MLQ1000_MPS1051_MC1977_L1p0', 'PsiPsiToLQChi_MLQ1000_MPS2221_MC12089_L1p0'])
# StudiesRunner17.CalculateLeptonReconstructionEfficiencies()
# StudiesRunner17.CalculateLeptonIDEfficiencies()
# PreselectionRunner17 = ScriptRunner(year='2017', analysisname=analysisname, t3_base_path=t3_base_path, result_base_path=result_base_path, selectionstage='Preselection', selectionname='PsiPsi_Initial', backgrounds=['TT', 'DYJets', 'ST', 'WJets', 'QCD_Had'], signals=['PsiPsiToLQChi_MLQ1000_MPS117_MC1100_L1p0', 'PsiPsiToLQChi_MLQ1000_MPS244_MC1214_L1p0', 'PsiPsiToLQChi_MLQ1000_MPS567_MC1457_L1p0', 'PsiPsiToLQChi_MLQ1000_MPS1051_MC1977_L1p0', 'PsiPsiToLQChi_MLQ1000_MPS2221_MC12089_L1p0'])
# PreselectionRunner17.CalculateSelectionEfficiencies(num_and_denom=[('trigger', 'input'), ('base', 'input'), ('plateau', 'input'), ('tach_leptonveto', 'input'), ('tach_plateau', 'input'), ('tach_tau', 'input'), ('tach_met', 'input'), ('final', 'input')])
# LQ pair-production
# systematics17 = ['nominal', 'lumi', 'rate_tt', 'rate_dyjets', 'rate_st']
# # backgrounds17 = ['TT', 'DYJets', 'ST', 'WJets', 'QCD_Had']
# backgrounds17 = ['TT', 'DYJets', 'ST']
# signals17 = ['LQLQToBTau_MLQ1000_MPS117_MC1100_L1p0', 'LQLQToBTau_MLQ1360_MPS117_MC1100_L1p0', 'LQLQToBTau_MLQ1810_MPS117_MC1100_L1p0', 'LQLQToBTau_MLQ2170_MPS117_MC1100_L1p0', 'LQLQToBTau_MLQ2620_MPS117_MC1100_L1p0']
# channels17 = ['elch', 'much', 'elmuch', 'tach']
# histnames_in_out_per_category17 = {'cat1': ('stmet3', 'STMET'), 'cat2': ('stmet3', 'STMET')}
# limithisttag17 = 'Limits'
#
# processes_per_systematic17 = {'lumi': 'all', 'rate_tt': 'TT', 'rate_dyjets': 'DYJets', 'rate_st': 'ST', 'rate_wjets': 'WJets', 'rate_qcd': 'QCD_Had', 'pu': 'all', 'pdf': 'all', 'muid': 'all', 'muiso': 'all', 'mutrigger': 'all', 'eleid': 'all', 'elereco': 'all', 'btag_bc': 'all', 'btag_udsg': 'all', 'scale_tt': 'TT', 'scale_dyjets': 'DYJets', 'scale_st': 'ST'}
# pdf_per_systematic17 = {'lumi': 'lnN', 'pu': 'shape', 'pdf': 'shape', 'muid': 'shape', 'muiso': 'shape', 'mutrigger': 'shape', 'eleid': 'shape', 'elereco': 'shape', 'btag_bc': 'shape', 'btag_udsg': 'shape', 'rate_tt': 'lnN', 'rate_dyjets': 'lnN', 'rate_st': 'lnN', 'rate_wjets': 'lnN', 'rate_qcd': 'lnN', 'scale_tt': 'shape', 'scale_dyjets': 'shape', 'scale_st': 'shape'}
# value_per_systematic17 = {'lumi': 1.025, 'pu': 1, 'pdf': 1, 'muid': 1, 'muiso': 1, 'mutrigger': 1, 'eleid': 1, 'elereco': 1, 'btag_bc': 1, 'btag_udsg': 1, 'rate_tt': 1.056, 'rate_dyjets': 1.1, 'rate_st': 1.1, 'rate_wjets': 1.25, 'rate_qcd': 1.5, 'scale_tt': 1, 'scale_dyjets': 1, 'scale_st': 1}
# Analyzer17 = PostAnalyzer('2017', analysisname=analysisname, t3_base_path=t3_base_path, result_base_path=result_base_path, selectionname='TauTau_ElTau_MuTau_ElMu_2Jets_12LooseBTag_met50', systematics=systematics17, backgrounds=backgrounds17, signals=signals17, channels=channels17, histnames_in_out_per_category=histnames_in_out_per_category17, limithisttag=limithisttag17, processes_per_systematic=processes_per_systematic17, pdf_per_systematic=pdf_per_systematic17, value_per_systematic=value_per_systematic17)
# Analyzer17.ProduceCombineHistograms(signal_scaled_by=0.01)
# Analyzer17.CreateDatacards()
# Analyzer17.CombineChannels()
# Analyzer17.ExecuteCombineCombination()
# Analyzer17.PlotLimits(signal_scaled_by=0.01, draw_observed=False)
# LQ t-channel
PreselectionRunner17_tach = ScriptRunner(year='2017', analysisname=analysisname, t3_base_path=t3_base_path, result_base_path=result_base_path, selectionstage='Preselection', selectionname='TChannelTauNu_tau', backgrounds=['TT', 'VV', 'DYJets', 'ST', 'WJets', 'QCD_Had'], signals=['LQTChannelTauNu_DynamicScale_MLQ1000_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ1360_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ1810_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ2620_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ2980_MPS117_MC1100_L1p0'])
# PreselectionRunner17_tach.CalculateSelectionEfficiencies(num_and_denom=[('trigger', 'input'), ('base', 'input'), ('plateau', 'input'), ('tach_leptonveto', 'input'), ('tach_plateau', 'input'), ('tach_tau', 'input'), ('tach_met', 'input'), ('final', 'input')])
PreselectionRunner17_much = ScriptRunner(year='2017', analysisname=analysisname, t3_base_path=t3_base_path, result_base_path=result_base_path, selectionstage='Preselection', selectionname='TChannelTauNu_muon', backgrounds=['TT', 'VV', 'DYJets', 'ST', 'WJets', 'QCD_Had'], signals=['LQTChannelTauNu_DynamicScale_MLQ1000_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ1360_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ1810_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ2620_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ2980_MPS117_MC1100_L1p0'])
# PreselectionRunner17_much.CalculateSelectionEfficiencies(num_and_denom=[('trigger', 'input'), ('plateau', 'input'), ('much_leptonveto', 'input'), ('much_lepton', 'input'), ('much_plateau', 'input'), ('much_dphi', 'input'), ('much_ptratio', 'input'), ('final', 'input')])
FullselectionRunner17_tach = ScriptRunner(year='2017', analysisname=analysisname, t3_base_path=t3_base_path, result_base_path=result_base_path, selectionstage='Fullselection', selectionname='TChannelTauNu_tau', backgrounds=['TT', 'VV', 'DYJets', 'ST', 'WJets', 'QCD_Had'], signals=['LQTChannelTauNu_DynamicScale_MLQ1000_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ1360_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ1810_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ2620_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ2980_MPS117_MC1100_L1p0'])
# FullselectionRunner17_tach.CalculateSelectionEfficiencies(num_and_denom=[('tach_trigger', 'input'), ('tach_taucleaner', 'input'), ('tach_tauselection', 'input'), ('tach_dphi', 'input'), ('tach_ptratio', 'input'), ('tach_cat1_nominal', 'input')])
# FullselectionRunner17_tach.CalculateSelectionEfficiencies(num_and_denom=[('tach_trigger', 'input'), ('tach_taucleaner', 'input'), ('tach_tauselection', 'input'), ('tach_cat1jets_initial', 'input'), ('tach_cat1jets_dphi', 'input'), ('tach_cat1jets_ptratio', 'input'), ('tach_cat1jets_1btagloose', 'input'), ('tach_cat1jets_1btagloose_nominal', 'input'), ('tach_cat1jets_0btagloose', 'input'), ('tach_cat1jets_0btagloose_nominal', 'input'), ('tach_cat1jets_nominal', 'input'), ('tach_cat0jets_initial', 'input'), ('tach_cat0jets_dphi', 'input'), ('tach_cat0jets_ptratio', 'input'), ('tach_cat0jets_nominal', 'input'), ('tach_cat1_nominal', 'input')])
FullselectionRunner17_much = ScriptRunner(year='2017', analysisname=analysisname, t3_base_path=t3_base_path, result_base_path=result_base_path, selectionstage='Fullselection', selectionname='TChannelTauNu_muon_wprime', backgrounds=['TT', 'VV', 'DYJets', 'ST', 'WJets', 'QCD_Had'], signals=['LQTChannelTauNu_DynamicScale_MLQ1000_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ1360_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ1810_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ2620_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ2980_MPS117_MC1100_L1p0'])
# FullselectionRunner17_much.CalculateSelectionEfficiencies(num_and_denom=[('much_jetveto5', 'input'), ('much_btagveto1', 'input'), ('much_dphi25', 'input'), ('much_ptratio0415', 'input'), ('much_ptmu50', 'input'), ('much_cat1_nominal', 'input')])
systematics17 = ['nominal', 'lumi', 'rate_tt', 'rate_dyjets', 'rate_st', 'rate_vv', 'rate_wjets', 'rate_qcd']
backgrounds17 = ['TT', 'DYJets', 'ST', 'WJets', 'QCD_Had']
# backgrounds17 = ['TT', 'DYJets', 'ST', 'WJets']
signals17 = ['LQTChannelTauNu_DynamicScale_MLQ1000_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ1360_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ1810_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ2620_MPS117_MC1100_L1p0', 'LQTChannelTauNu_DynamicScale_MLQ2980_MPS117_MC1100_L1p0']
channels17 = ['tach']
histnames_in_out_per_category17 = {
'cat1': ('mt_taumet_rebinlimit2', 'MT'),
# 'cat1jets_0btagloose': ('mt_taumet_rebinlimit2', 'MT'),
# 'cat1jets_1btagloose': ('mt_taumet_rebinlimit2', 'MT'),
# 'cat0jets': ('mt_taumet_rebinlimit2', 'MT'),
}
# histnames_in_out_per_category17 = {'cat1': ('mt_mumet_rebinlimit2', 'MT')}
limithisttag17 = 'General'
processes_per_systematic17 = {'lumi': 'all', 'rate_tt': 'TT', 'rate_dyjets': 'DYJets', 'rate_st': 'ST', 'rate_wjets': 'WJets', 'rate_qcd': 'QCD_Had', 'rate_vv': 'VV', 'pu': 'all', 'pdf': 'all', 'muid': 'all', 'muiso': 'all', 'mutrigger': 'all', 'eleid': 'all', 'elereco': 'all', 'btag_bc': 'all', 'btag_udsg': 'all', 'scale_tt': 'TT', 'scale_dyjets': 'DYJets', 'scale_st': 'ST'}
pdf_per_systematic17 = {'lumi': 'lnN', 'pu': 'shape', 'pdf': 'shape', 'muid': 'shape', 'muiso': 'shape', 'mutrigger': 'shape', 'eleid': 'shape', 'elereco': 'shape', 'btag_bc': 'shape', 'btag_udsg': 'shape', 'rate_tt': 'lnN', 'rate_dyjets': 'lnN', 'rate_st': 'lnN', 'rate_vv': 'lnN', 'rate_wjets': 'lnN', 'rate_qcd': 'lnN', 'scale_tt': 'shape', 'scale_dyjets': 'shape', 'scale_st': 'shape'}
value_per_systematic17 = {'lumi': 1.025, 'pu': 1, 'pdf': 1, 'muid': 1, 'muiso': 1, 'mutrigger': 1, 'eleid': 1, 'elereco': 1, 'btag_bc': 1, 'btag_udsg': 1, 'rate_tt': 1.056, 'rate_dyjets': 1.1, 'rate_st': 1.1, 'rate_vv': 1.2, 'rate_wjets': 1.25, 'rate_qcd': 1.5, 'scale_tt': 1, 'scale_dyjets': 1, 'scale_st': 1}
Analyzer17 = PostAnalyzer('2017', analysisname=analysisname, t3_base_path=t3_base_path, result_base_path=result_base_path, selectionstage='Fullselection', selectionname='TChannelTauNu_wprime_tau', systematics=systematics17, backgrounds=backgrounds17, signals=signals17, channels=channels17, histnames_in_out_per_category=histnames_in_out_per_category17, limithisttag=limithisttag17, processes_per_systematic=processes_per_systematic17, pdf_per_systematic=pdf_per_systematic17, value_per_systematic=value_per_systematic17)
