Skip to content

Commit

Permalink
Changed style of plot before Final Reading
Browse files Browse the repository at this point in the history
  • Loading branch information
@mgrippo
mgrippo committed Aug 6, 2015
1 parent 5ce9959 commit b5687c2
Show file tree
Hide file tree
Showing 6 changed files with 85 additions and 51 deletions.
24 changes: 12 additions & 12 deletions Analysis/config/sources.cfg
View file
Original file line number Diff line number Diff line change
Expand Up @@ -118,7 +118,7 @@ draw_sf: 50
[ggHhh260]
type: SIGNAL
type: LIMITS
title: hh#rightarrow#tau#taubb(m_{H}=260 tan#beta=2)
title: H#rightarrowhh(m_{H}=260,tan#beta=2)
datacard: ggHTohhTo2Tau2B260
color: blue
file: ggH_hh_bbtautau_260.root 0.0358936
Expand All @@ -128,7 +128,7 @@ draw_sf: 10
[ggHhh270]
type: SIGNAL
type: LIMITS
title: hh#rightarrow#tau#taubb(m_{H}=270 tan#beta=2)
title: H#rightarrowhh(m_{H}=270,tan#beta=2)
datacard: ggHTohhTo2Tau2B270
color: blue
file: ggH_hh_bbtautau_270.root 0.03057
Expand All @@ -138,7 +138,7 @@ draw_sf: 10
[ggHhh280]
type: SIGNAL
type: LIMITS
title: hh#rightarrow#tau#taubb(m_{H}=280 tan#beta=2)
title: H#rightarrowhh(m_{H}=280,tan#beta=2)
datacard: ggHTohhTo2Tau2B280
color: blue
file: ggH_hh_bbtautau_280.root 0.0284073
Expand All @@ -148,7 +148,7 @@ draw_sf: 10
[ggHhh290]
type: SIGNAL
type: LIMITS
title: hh#rightarrow#tau#taubb(m_{H}=290 tan#beta=2)
title: H#rightarrowhh(m_{H}=290,tan#beta=2)
datacard: ggHTohhTo2Tau2B290
color: blue
file: ggH_hh_bbtautau_290.root 0.0201117
Expand All @@ -158,7 +158,7 @@ draw_sf: 10
[ggHhh300]
type: SIGNAL
type: LIMITS
title: hh#rightarrow#tau#taubb(m_{H}=300 tan#beta=2)
title: H#rightarrowhh(m_{H}=300,tan#beta=2)
datacard: ggHTohhTo2Tau2B300
color: blue
file: ggH_hh_bbtautau_300.root 0.0053075
Expand All @@ -168,7 +168,7 @@ draw_sf: 10
[ggHhh310]
type: SIGNAL
type: LIMITS
title: hh#rightarrow#tau#taubb(m_{H}=310 tan#beta=2)
title: H#rightarrowhh(m_{H}=310,tan#beta=2)
datacard: ggHTohhTo2Tau2B310
color: blue
file: ggH_hh_bbtautau_310.root 0.0221424
Expand All @@ -178,7 +178,7 @@ draw_sf: 10
[ggHhh320]
type: SIGNAL
type: LIMITS
title: hh#rightarrow#tau#taubb(m_{H}=320 tan#beta=2)
title: H#rightarrowhh(m_{H}=320,tan#beta=2)
datacard: ggHTohhTo2Tau2B320
color: blue
file: ggH_hh_bbtautau_320.root 0.0215317
Expand All @@ -188,7 +188,7 @@ draw_sf: 10
[ggHhh330]
type: SIGNAL
type: LIMITS
title: hh#rightarrow#tau#taubb(m_{H}=330 tan#beta=2)
title: H#rightarrowhh(m_{H}=330,tan#beta=2)
datacard: ggHTohhTo2Tau2B330
color: blue
file: ggH_hh_bbtautau_330.root 0.0199818
Expand All @@ -198,7 +198,7 @@ draw_sf: 10
[ggHhh340]
type: SIGNAL
type: LIMITS
title: hh#rightarrow#tau#taubb(m_{H}=340 tan#beta=2)
title: H#rightarrowhh(m_{H}=340,tan#beta=2)
datacard: ggHTohhTo2Tau2B340
color: blue
file: ggH_hh_bbtautau_340.root 0.0207792
Expand All @@ -208,7 +208,7 @@ draw_sf: 10
[ggHhh350]
type: SIGNAL
type: LIMITS
title: hh#rightarrow#tau#taubb(m_{H}=350 tan#beta=2)
title: H#rightarrowhh(m_{H}=350,tan#beta=2)
datacard: ggHTohhTo2Tau2B350
color: blue
file: ggH_hh_bbtautau_350.root 0.0179311
Expand Down Expand Up @@ -600,7 +600,7 @@ uncertainty: zttNorm
type: BACKGROUND
type: ZTT
type: LIMITS
title: Z #rightarrow tautau
title: Z#rightarrow#tau#tau
datacard: ZTT
draw: true
color: yellow_custom
Expand All @@ -611,7 +611,7 @@ type: BACKGROUND
type: COMPOSIT
type: LIMITS
datacard: ZLL
title: Z #rightarrow ll
title: Z#rightarrowll
draw: true
color: azure_custom
subcategory: ZL
Expand Down
2 changes: 1 addition & 1 deletion Analysis/include/BaseFlatTreeAnalyzer.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -849,7 +849,7 @@ class BaseFlatTreeAnalyzer {

static const std::map<std::string, size_t> histogramsToBlind = {
{ FlatAnalyzerData_semileptonic::m_sv_Name(), 1 }, { FlatAnalyzerData::m_vis_Name(), 1 },
{ FlatAnalyzerData::m_bb_Name(), 1 }, { FlatAnalyzerData::m_ttbb_Name(), 2 },
{ FlatAnalyzerData_semileptonic::m_bb_Name(), 1 }, { FlatAnalyzerData::m_ttbb_Name(), 2 },
{ FlatAnalyzerData::m_ttbb_kinfit_Name(), 2 }, { FlatAnalyzerData::m_bb_slice_Name(), 3 }
};

Expand Down
43 changes: 27 additions & 16 deletions Analysis/include/FlatAnalyzerData.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -37,20 +37,20 @@ class FlatAnalyzerData : public root_ext::AnalyzerData {
public:
TH1D_ENTRY_CUSTOM_EX(m_vis, M_tt_Bins(), "M_{vis}(GeV)", "Events", false, 1.1, false, SaveAll)
TH1D_ENTRY_CUSTOM_EX(m_bb_slice, M_tt_bbSlice_Bins(), "2DM_{sv}(GeV)", "Events", false, 1.1, false, SaveAll)
TH1D_ENTRY_CUSTOM_EX(m_ttbb, M_ttbb_Bins(), "M_{#tau#tau+jj} (GeV)", "dN/dM_{#tau#tau+jj} (1/GeV)", false, 1.3, true, SaveAll)
TH1D_ENTRY_CUSTOM_EX(m_ttbb_kinfit, M_ttbb_Bins(), "M_{H}^{kinfit} (GeV)", "dN/dM_{H}^{kinfit} (1/GeV)", false, 1.3, true, true)
TH1D_ENTRY_CUSTOM_EX(m_ttbb, M_ttbb_Bins(), "M_{#tau#tau+jj} (GeV)", "dN/dm_{#tau#tau+jj} (1/GeV)", false, 1.5, true, SaveAll)
TH1D_ENTRY_CUSTOM_EX(m_ttbb_kinfit, M_ttbb_Bins(), "M_{H}^{kinfit} (GeV)", "dN/dm_{H}^{kinfit} (1/GeV)", false, 1.4, true, true)


TH1D_ENTRY_EX(pt_b1, 20, 0, 200, "Leading selected jet p_{T} (GeV)", "Events", false, 1.3, false, SaveAll)
TH1D_ENTRY_EX(pt_b1, 20, 0, 200, "Leading selected jet p_{T} (GeV)", "Events", false, 1.4, false, SaveAll)
TH1D_ENTRY_EX(eta_b1, 25, -2.5, 2.5, "Leading selected jet #eta", "Events", false, 1.8, false, SaveAll)
TH1D_ENTRY_EX(csv_b1, 25, 0, 1, "Leading selected jet CSV", "Events", false, 1.3, false, SaveAll)
TH1D_ENTRY_EX(pt_b2, 20, 0, 200, "Subleading selected jet p_{T} (GeV)", "Events", false, 1.3, false, SaveAll)
TH1D_ENTRY_EX(csv_b1, 25, 0, 1, "Leading selected jet CSV", "Events", false, 1.4, false, SaveAll)
TH1D_ENTRY_EX(pt_b2, 20, 0, 200, "Subleading selected jet p_{T} (GeV)", "Events", false, 1.4, false, SaveAll)
TH1D_ENTRY_EX(eta_b2, 25, -2.5, 2.5, "Subleading selected jet #eta", "Events", false, 1.8, false, SaveAll)
TH1D_ENTRY_EX(csv_b2, 25, 0, 1, "Subleading selected jet CSV", "Events", false, 1.3, false, SaveAll)
TH1D_ENTRY_EX(csv_b2, 25, 0, 1, "Subleading selected jet CSV", "Events", false, 1.4, false, SaveAll)
TH1D_ENTRY_EX(pt_H_tt, 20, 0, 300, "P_{T}(GeV)", "Events", false, 1.1, false, SaveAll)
TH1D_ENTRY_EX(pt_H_bb, 20, 0, 300, "P_{T}(GeV)", "Events", false, 1.1, false, SaveAll)
TH1D_ENTRY_EX(pt_H_hh, 20, 0, 300, "P_{T}(GeV)", "Events", false, 1.1, false, SaveAll)
TH1D_ENTRY_EX(m_bb, 30, 0, 600, "M_{jj} (GeV)", "dN/dM_{jj} (1/GeV)", false, 1.25, true, SaveAll)
//TH1D_ENTRY_EX(m_bb, 20, 0, 400, "M_{jj} (GeV)", "dN/dm_{jj} (1/GeV)", false, 1.4, true, SaveAll)
TH1D_ENTRY_EX(DeltaPhi_tt, 22, 0., 3.3, "#Delta#Phi_{#tau#tau}[rad]", "Events", false, 1.3, false, SaveAll)
TH1D_ENTRY_EX(DeltaPhi_bb, 22, 0., 3.3, "#Delta#Phi_{bb}[rad]", "Events", false, 1.8, false, SaveAll)
TH1D_ENTRY_EX(DeltaPhi_bb_MET, 22, 0., 3.3, "#Delta#Phi_{bb,MET}[rad]", "Events", false, 1.5, false, SaveAll)
Expand All @@ -67,7 +67,7 @@ class FlatAnalyzerData : public root_ext::AnalyzerData {
TH1D_ENTRY_EX(pull_balance, 20, -10, 10, "pull_balance", "Events", false, 2, false, SaveAll)
TH1D_ENTRY_EX(pull_balance_1, 100, -10, 10, "pull_balance_1", "Events", false, 1.1, false, SaveAll)
TH1D_ENTRY_EX(pull_balance_2, 100, -10, 10, "pull_balance_1", "Events", false, 1.1, false, SaveAll)
TH1D_ENTRY_EX(MET, 20, 0, 100, "E_{T}^{miss} (GeV)", "Events", false, 1.3, false, SaveAll)
TH1D_ENTRY_EX(MET, 20, 0, 100, "E_{T}^{miss} (GeV)", "Events", false, 1.4, false, SaveAll)
TH1D_ENTRY_EX(nJets_Pt30, 20, -0.5, 19.5, "nJets", "Events", false, 1.1, false, SaveAll)
TH2D_ENTRY_EX(csv_b1_vs_ptb1, 20, 0, 200, 25, 0, 1, "P_{T}(GeV)(leading_jet)", "CSV(leading_jet)", false, 1, SaveAll)
TH2D_ENTRY_EX(chi2_vs_ptb1, 20, 0, 200, 20, 0, 100, "P_{T}(GeV)(leading_jet)", "#chi^{2}", false, 1, SaveAll)
Expand Down Expand Up @@ -139,7 +139,7 @@ class FlatAnalyzerData : public root_ext::AnalyzerData {
DeltaR_bb().Fill(eventInfo.bjet_momentums.at(eventInfo.selected_bjets.first).DeltaR(
eventInfo.bjet_momentums.at(eventInfo.selected_bjets.second)), weight);
pt_H_bb().Fill(eventInfo.Hbb.Pt(),weight);
m_bb().Fill(eventInfo.Hbb.M(), weight);
// m_bb().Fill(eventInfo.Hbb.M(), weight);
DeltaPhi_bb_MET().Fill(std::abs(eventInfo.Hbb.DeltaPhi(eventInfo.MET)), weight);
DeltaPhi_hh().Fill(std::abs(eventInfo.Htt.DeltaPhi(eventInfo.Hbb)), weight);
DeltaR_hh().Fill(eventInfo.Htt.DeltaR(eventInfo.Hbb), weight);
Expand All @@ -165,7 +165,7 @@ class FlatAnalyzerData : public root_ext::AnalyzerData {
virtual void CreateAll()
{
m_vis(); m_bb_slice(); m_ttbb(); m_ttbb_kinfit(); pt_b1(); eta_b1(); csv_b1(); pt_b2(); eta_b2();
csv_b2(); pt_H_tt(); pt_H_bb(); pt_H_hh(); m_bb(); DeltaPhi_tt(); DeltaPhi_bb(); DeltaPhi_bb_MET();
csv_b2(); pt_H_tt(); pt_H_bb(); pt_H_hh(); /*m_bb();*/ DeltaPhi_tt(); DeltaPhi_bb(); DeltaPhi_bb_MET();
DeltaPhi_tt_MET(); DeltaPhi_hh(); DeltaR_tt(); DeltaR_bb(); DeltaR_hh(); mt_2(); pt_H_tt_MET(); convergence();
chi2(); fit_probability(); pull_balance(); pull_balance_1(); pull_balance_2(); MET(); nJets_Pt30();
csv_b1_vs_ptb1(); chi2_vs_ptb1(); mH_vs_chi2();
Expand Down Expand Up @@ -210,13 +210,14 @@ class FlatAnalyzerData : public root_ext::AnalyzerData {

class FlatAnalyzerData_semileptonic : public FlatAnalyzerData {
public:
TH1D_ENTRY_CUSTOM_EX(m_sv, M_tt_Bins(), "M_{#tau#tau} (GeV)", "dN/dM_{#tau#tau} (1/GeV)", false, 1.5, true, true)
TH1D_ENTRY_CUSTOM_EX(m_sv, M_tt_Bins(), "M_{#tau#tau} (GeV)", "dN/dm_{#tau#tau} (1/GeV)", false, 1.5, true, true)
TH1D_ENTRY_EX(pt_1, 20, 0, 200, "P_{T}(leading#tau_{h})(GeV)", "Events", false, 1.1, false, SaveAll)
TH1D_ENTRY_EX(eta_1, 25, -2.5, 2.5, "#eta(leading#tau_{h})", "Events", false, 2, false, SaveAll)
TH1D_ENTRY_EX(pt_2, 20, 0, 200, "P_{T}(subleading#tau_{h})(GeV)", "Events", false, 1, false, SaveAll)
TH1D_ENTRY_EX(eta_2, 25, -2.5, 2.5, "#eta(subleading#tau_{h})", "Events", false, 2, false, SaveAll)

TH1D_ENTRY_EX(mt_1, 50, 0, 50, "M_{T}(GeV)", "Events", false, 1.1, false, SaveAll)
TH1D_ENTRY_EX(m_bb, 20, 0, 400, "M_{jj} (GeV)", "Events/ 20 GeV", false, 1.4, true, SaveAll)

explicit FlatAnalyzerData_semileptonic(bool _fill_all) : FlatAnalyzerData(_fill_all) {}

Expand All @@ -234,6 +235,8 @@ class FlatAnalyzerData_semileptonic : public FlatAnalyzerData {
pt_2().Fill(event.pt_2, weight);
eta_2().Fill(event.eta_2, weight);
mt_1().Fill(event.mt_1, weight);
if(!eventInfo.has_bjet_pair) return;
m_bb().Fill(eventInfo.Hbb.M(), weight);
}

virtual root_ext::SmartHistogram<TH1D>& m_sv_base() override { return m_sv(); }
Expand All @@ -248,7 +251,7 @@ class FlatAnalyzerData_semileptonic : public FlatAnalyzerData {
virtual void CreateAll() override
{
FlatAnalyzerData::CreateAll();
m_sv(); pt_1(); eta_1(); pt_2(); eta_2(); mt_1();
m_sv(); pt_1(); eta_1(); pt_2(); eta_2(); mt_1();m_bb();
}
};

Expand All @@ -269,9 +272,9 @@ class FlatAnalyzerData_semileptonic_2tag : public FlatAnalyzerData_semileptonic

class FlatAnalyzerData_tautau : public FlatAnalyzerData {
public:
TH1D_ENTRY_EX(pt_1, 20, 0, 200, "Leading tau p_{T} (GeV)", "Events", false, 1.3, false, SaveAll)
TH1D_ENTRY_EX(pt_1, 20, 0, 200, "Leading tau p_{T} (GeV)", "Events", false, 1.4, false, SaveAll)
TH1D_ENTRY_EX(eta_1, 25, -2.5, 2.5, "Leading tau #eta", "Events", false, 1.8, false, SaveAll)
TH1D_ENTRY_EX(pt_2, 20, 0, 200, "Subleading tau p_{T} (GeV)", "Events", false, 1.3, false, SaveAll)
TH1D_ENTRY_EX(pt_2, 20, 0, 200, "Subleading tau p_{T} (GeV)", "Events", false, 1.4, false, SaveAll)
TH1D_ENTRY_EX(eta_2, 25, -2.5, 2.5, "Subleading tau #eta", "Events", false, 1.8, false, SaveAll)

TH1D_ENTRY_EX(mt_1, 20, 0, 200, "M_{T}(GeV)", "Events", false, 1.1, false, SaveAll)
Expand Down Expand Up @@ -312,7 +315,8 @@ class FlatAnalyzerData_tautau : public FlatAnalyzerData {

class FlatAnalyzerData_tautau_other_tag : public FlatAnalyzerData_tautau {
public:
TH1D_ENTRY_CUSTOM_EX(m_sv, M_tt_Bins(), "M_{#tau#tau} (GeV)", "dN/dM_{#tau#tau} (1/GeV)", false, 1.3, true, true)
TH1D_ENTRY_CUSTOM_EX(m_sv, M_tt_Bins(), "M_{#tau#tau} (GeV)", "dN/dm_{#tau#tau} (1/GeV)", false, 1.5, true, true)
TH1D_ENTRY_EX(m_bb, 20, 0, 400, "M_{jj} (GeV)", "Events/ 20 GeV", false, 1.4, true, SaveAll)

explicit FlatAnalyzerData_tautau_other_tag(bool _fill_all) : FlatAnalyzerData_tautau(_fill_all) {}

Expand All @@ -323,6 +327,8 @@ class FlatAnalyzerData_tautau_other_tag : public FlatAnalyzerData_tautau {
{
FlatAnalyzerData_tautau::Fill(eventInfo, weight);
m_sv().Fill(eventInfo.event->m_sv_MC, weight);
if(!eventInfo.has_bjet_pair) return;
m_bb().Fill(eventInfo.Hbb.M(), weight);
}

virtual root_ext::SmartHistogram<TH1D>& m_sv_base() override { return m_sv(); }
Expand All @@ -331,13 +337,15 @@ class FlatAnalyzerData_tautau_other_tag : public FlatAnalyzerData_tautau {
{
FlatAnalyzerData_tautau::CreateAll();
m_sv();
m_bb();
}
};


class FlatAnalyzerData_tautau_2tag : public FlatAnalyzerData_tautau {
public:
TH1D_ENTRY_CUSTOM_EX(m_sv, M_tt_Bins(), "M_{#tau#tau} (GeV)", "dN/dM_{#tau#tau} (1/GeV)", false, 1.7, true, true)
TH1D_ENTRY_CUSTOM_EX(m_sv, M_tt_Bins(), "M_{#tau#tau} (GeV)", "dN/dm_{#tau#tau} (1/GeV)", false, 1.7, true, true)
TH1D_ENTRY_EX(m_bb, 10, 0, 400, "M_{jj} (GeV)", "Events/ 40 GeV", false, 1.4, true, SaveAll)

explicit FlatAnalyzerData_tautau_2tag(bool _fill_all) : FlatAnalyzerData_tautau(_fill_all) {}

Expand All @@ -348,6 +356,8 @@ class FlatAnalyzerData_tautau_2tag : public FlatAnalyzerData_tautau {
{
FlatAnalyzerData_tautau::Fill(eventInfo, weight);
m_sv().Fill(eventInfo.event->m_sv_MC, weight);
if(!eventInfo.has_bjet_pair) return;
m_bb().Fill(eventInfo.Hbb.M(), weight);
}

virtual root_ext::SmartHistogram<TH1D>& m_sv_base() override { return m_sv(); }
Expand All @@ -362,6 +372,7 @@ class FlatAnalyzerData_tautau_2tag : public FlatAnalyzerData_tautau {
{
FlatAnalyzerData_tautau::CreateAll();
m_sv();
m_bb();
}
};

Expand Down
Loading

0 comments on commit b5687c2

Please sign in to comment.