|
|
[Rivet-svn] r1893 - in trunk: bin data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Thu Oct 8 16:04:20 BST 2009
Author: buckley
Date: Thu Oct 8 16:04:19 2009
New Revision: 1893
Log:
Various plot info improvements
Modified:
trunk/bin/make-plots
trunk/data/plotinfo/MC_LHC_SUSY.plot
trunk/data/plotinfo/MC_LHC_WANALYSIS.plot
trunk/data/plotinfo/MC_LHC_ZANALYSIS.plot
trunk/data/plotinfo/UA1_1990_S2044935.plot
trunk/src/Analyses/MC_LHC_SUSY.cc
Modified: trunk/bin/make-plots
==============================================================================
--- trunk/bin/make-plots Thu Oct 8 15:40:30 2009 (r1892)
+++ trunk/bin/make-plots Thu Oct 8 16:04:19 2009 (r1893)
@@ -94,6 +94,7 @@
else:
self.description['RatioPlotSizeY'] = 3.
+ ## TODO: Also accept "true", "True", "TRUE", "yes", etc.
self.description['LogX'] = self.description.has_key('LogX') and self.description['LogX']=='1'
self.description['LogY'] = self.description.has_key('LogY') and self.description['LogY']=='1'
Modified: trunk/data/plotinfo/MC_LHC_SUSY.plot
==============================================================================
--- trunk/data/plotinfo/MC_LHC_SUSY.plot Thu Oct 8 15:40:30 2009 (r1892)
+++ trunk/data/plotinfo/MC_LHC_SUSY.plot Thu Oct 8 16:04:19 2009 (r1893)
@@ -65,6 +65,14 @@
# END PLOT
# BEGIN PLOT
+/MC_LHC_SUSY/mll-ossf-mumu
+Title=$m_{\mu\mu}$ distribution for clean $\mu^+ \mu^-$
+XLabel=$m_{\mu\mu}$
+YLabel=$\mathrm{d}N_\text{evt}/\mathrm{d}m_{\mu\mu}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT
/MC_LHC_SUSY/mll-all-ossf-mumu
Title=$m_{\mu\mu}$ distribution for combinatoric $\mu^+ \mu^-$
XLabel=$m_{\mu\mu}$
@@ -81,6 +89,14 @@
# END PLOT
# BEGIN PLOT
+/MC_LHC_SUSY/mll-ossf-ee
+Title=$m_{ee}$ distribution for clean $e^+ e^-$
+XLabel=$m_{ee}$
+YLabel=$\mathrm{d}N_\text{evt}/\mathrm{d}m_{ee}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT
/MC_LHC_SUSY/mll-all-ossf-ee
Title=$m_{ee}$ distribution for combinatoric $e^+ e^-$
XLabel=$m_{ee}$
@@ -97,6 +113,14 @@
# END PLOT
# BEGIN PLOT
+/MC_LHC_SUSY/mll-osof-emu
+Title=$m_{e\mu}$ distribution for clean $e^\pm \mu^\mp$
+XLabel=$m_{e\mu}$
+YLabel=$\mathrm{d}N_\text{evt}/\mathrm{d}m_{e\mu}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT
/MC_LHC_SUSY/mll-all-osof-emu
Title=$m_{e\mu}$ distribution for combinatoric $e^\pm \mu^\mp$
XLabel=$m_{e\mu}$
Modified: trunk/data/plotinfo/MC_LHC_WANALYSIS.plot
==============================================================================
--- trunk/data/plotinfo/MC_LHC_WANALYSIS.plot Thu Oct 8 15:40:30 2009 (r1892)
+++ trunk/data/plotinfo/MC_LHC_WANALYSIS.plot Thu Oct 8 16:04:19 2009 (r1893)
@@ -19,25 +19,25 @@
# BEGIN PLOT /MC_LHC_WANALYSIS/pt-ch
Title=Charged track $p_\perp$
XLabel=$p_\perp$ / GeV
-YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}p_\perp^\text{ch}$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}p_\perp^\text{ch}$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_WANALYSIS/ptavg-ch
Title=Average charged track $p_\perp$
XLabel=$\langle p_\perp^\text{ch} \rangle$ / GeV
-YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}\langle p_\perp^\text{ch} \langle$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}\langle p_\perp^\text{ch} \langle$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_WANALYSIS/ptrms-ch
Title=Width of charged track $p_\perp$
XLabel=$\hat{\sigma}_{p_\perp^\text{ch}}$ / GeV
-YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}\hat{\sigma}_{p_\perp^\text{ch}}$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}\hat{\sigma}_{p_\perp^\text{ch}}$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_WANALYSIS/pt-w
Title=W candidate $p_\perp$
XLabel=$p_\perp(\text{W})$ / GeV
-YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}p_\perp(\text{W})$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}p_\perp(\text{W})$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_WANALYSIS/logpt-w
@@ -49,7 +49,7 @@
# BEGIN PLOT /MC_LHC_WANALYSIS/pt-jet
Title=Jet $p_\perp$
XLabel=$p_\perp(\text{jet})$ / GeV
-YLabel=$\mathrm{d}N_\text{jet}/\mathrm{d}p_\perp(\text{jet})$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{jet}/\mathrm{d}p_\perp(\text{jet})$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_WANALYSIS/logpt-jet
@@ -73,7 +73,7 @@
# BEGIN PLOT /MC_LHC_WANALYSIS/m-w
Title=W candidate mass
XLabel=$m(\text{W})$ / GeV
-YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}m(\text{W})$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{W}/\mathrm{d}m(\text{W})$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_WANALYSIS/m-w
Modified: trunk/data/plotinfo/MC_LHC_ZANALYSIS.plot
==============================================================================
--- trunk/data/plotinfo/MC_LHC_ZANALYSIS.plot Thu Oct 8 15:40:30 2009 (r1892)
+++ trunk/data/plotinfo/MC_LHC_ZANALYSIS.plot Thu Oct 8 16:04:19 2009 (r1893)
@@ -19,25 +19,25 @@
# BEGIN PLOT /MC_LHC_ZANALYSIS/pt-ch
Title=Charged track $p_\perp$
XLabel=$p_\perp$ / GeV
-YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}p_\perp^\text{ch}$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}p_\perp^\text{ch}$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_ZANALYSIS/ptavg-ch
Title=Average charged track $p_\perp$
XLabel=$\langle p_\perp^\text{ch} \rangle$ / GeV
-YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}\langle p_\perp^\text{ch} \langle$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}\langle p_\perp^\text{ch} \langle$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_ZANALYSIS/ptrms-ch
Title=Width of charged track $p_\perp$
XLabel=$\hat{\sigma}_{p_\perp^\text{ch}}$ / GeV
-YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}\hat{\sigma}_{p_\perp^\text{ch}}$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{ch}/\mathrm{d}\hat{\sigma}_{p_\perp^\text{ch}}$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_ZANALYSIS/pt-z
Title=Z candidate $p_\perp$
XLabel=$p_\perp(\text{Z})$ / GeV
-YLabel=$\mathrm{d}N_\text{Z}/\mathrm{d}p_\perp(\text{Z})$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{Z}/\mathrm{d}p_\perp(\text{Z})$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_ZANALYSIS/logpt-z
@@ -49,7 +49,7 @@
# BEGIN PLOT /MC_LHC_ZANALYSIS/pt-jet
Title=Jet $p_\perp$
XLabel=$p_\perp(\text{jet})$ / GeV
-YLabel=$\mathrm{d}N_\text{jet}/\mathrm{d}p_\perp(\text{jet})$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{jet}/\mathrm{d}p_\perp(\text{jet})$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_ZANALYSIS/logpt-jet
@@ -73,7 +73,7 @@
# BEGIN PLOT /MC_LHC_ZANALYSIS/m-z
Title=Z candidate mass
XLabel=$m(\text{Z})$ / GeV
-YLabel=$\mathrm{d}N_\text{Z}/\mathrm{d}m(\text{Z})$ / GeV$^-1$
+YLabel=$\mathrm{d}N_\text{Z}/\mathrm{d}m(\text{Z})$ / GeV$^{-1}$
# END PLOT
# BEGIN PLOT /MC_LHC_ZANALYSIS/m-z
Modified: trunk/data/plotinfo/UA1_1990_S2044935.plot
==============================================================================
--- trunk/data/plotinfo/UA1_1990_S2044935.plot Thu Oct 8 15:40:30 2009 (r1892)
+++ trunk/data/plotinfo/UA1_1990_S2044935.plot Thu Oct 8 16:04:19 2009 (r1893)
@@ -1,107 +1,107 @@
# BEGIN PLOT /UA1_1990_S2044935/d01-x01-y01
Title=Cross-section distribution at $\sqrt{s} = \unit{200}{\GeV}$
-XLabel=$n_{ch}$
-YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{ n_{ch}}$
+XLabel=$n_\text{ch}$
+YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{ n_\text{ch}}$ / mb
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d01-x01-y02
Title=Cross-section distribution at $\sqrt{s} = \unit{500}{\GeV}$
-XLabel=$n_{ch}$
-YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{ n_{ch}}$
+XLabel=$n_\text{ch}$
+YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{ n_\text{ch}}$ / mb
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d01-x01-y03
Title=Cross-section distribution at $\sqrt{s} = \unit{900}{\GeV}$
-XLabel=$n_{ch}$
-YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{ n_{ch}}$
+XLabel=$n_\text{ch}$
+YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{ n_\text{ch}}$ / mb
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d02-x01-y01
Title=Invariant cross section $\mathrm{E} \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ at $\eta =0$ and $\sqrt{s} = \unit{200}{\GeV}$
-XLabel=$p_{T}$
-YLabel=$\mathrm{E} \frac{d^{3}\sigma}{dp^{3}}$
+XLabel=$p_{T}$ / GeV
+YLabel=$E \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ / mb\,GeV$^{-2}
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d02-x01-y02
Title=Invariant cross section $\mathrm{E} \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ at $\eta =0$ and $\sqrt{s} = \unit{500}{\GeV}$
-XLabel=$p_{T}$
-YLabel=$\mathrm{E} \frac{d^{3}\sigma}{dp^{3}}$
+XLabel=$p_{T}$ / GeV
+YLabel=$E \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ / mb\,GeV$^{-2}
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d02-x01-y03
Title=Invariant cross section $\mathrm{E} \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ at $\eta =0$ and $\sqrt{s} = \unit{900}{\GeV}$
-XLabel=$p_{T}$
-YLabel=$\mathrm{E} \frac{d^{3}\sigma}{dp^{3}}$
+XLabel=$p_{T}$ / GeV
+YLabel=$E \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ / mb\,GeV$^{-2}
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d03-x01-y01
-Title=Cross section $\mathrm{E} \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ at $\eta =0$ and $\sqrt{s} = \unit{900}{\GeV}$ ($\frac{dn_{ch}}{d\eta} = 0.8 \dots 4 $)
-XLabel=$p_{T}$
-YLabel=$\mathrm{E} \frac{d^{3}\sigma}{dp^{3}}$
+Title=Cross section $\mathrm{E} \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ at $\eta =0$ and $\sqrt{s} = \unit{900}{\GeV}$ ($\mathrm{d}dn_\text{ch}/\mathrm{d}\eta = 0.8 \dots 4$)
+XLabel=$p_{T}$ / GeV
+YLabel=$E \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ / $\mu$b\,GeV$^{-2}
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d04-x01-y01
-Title=Cross section $\mathrm{E} \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ at $\eta =0$ and $\sqrt{s} = \unit{900}{\GeV}$ ($\frac{dn_{ch}}{d\eta} = 4 \dots 8 $)
-XLabel=$p_{T}$
-YLabel=$\mathrm{E} \frac{d^{3}\sigma}{dp^{3}}$
+Title=Cross section $\mathrm{E} \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ at $\eta =0$ and $\sqrt{s} = \unit{900}{\GeV}$ ($\mathrm{d}dn_\text{ch}/\mathrm{d}\eta = 4 \dots 8$)
+XLabel=$p_{T}$ / GeV
+YLabel=$E \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ / $\mu$b\,GeV$^{-2}
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d05-x01-y01
-Title=Cross section $\mathrm{E} \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ at $\eta =0$ and $\sqrt{s} = \unit{900}{\GeV}$ ($\frac{dn_{ch}}{d\eta} > 8 $)
-XLabel=$p_{T}$
-YLabel=$\mathrm{E} \frac{d^{3}\sigma}{dp^{3}}$
+Title=Cross section $\mathrm{E} \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ at $\eta =0$ and $\sqrt{s} = \unit{900}{\GeV}$ ($\mathrm{d}dn_\text{ch}/\mathrm{d}\eta > 8$)
+XLabel=$p_{T}$ / GeV
+YLabel=$E \mathrm{d}^{3}\sigma/\mathrm{d}p^{3}$ / $\mu$b\,GeV$^{-2}
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d06-x01-y01
-Title=$\langle p_{T} \rangle$ vs. $n_{ch}$ at $\sqrt{s} = \unit{200}{\GeV}$
-XLabel=$n_{ch}$
-YLabel=$\langle p_{T} \rangle$
+Title=$\langle p_{T} \rangle$ vs. $n_\text{ch}$ at $\sqrt{s} = \unit{200}{\GeV}$
+XLabel=$n_\text{ch}$
+YLabel=$\langle p_{T} \rangle$ / GeV
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d07-x01-y01
-Title=$\langle p_{T} \rangle$ vs. $n_{ch}$ at $\sqrt{s} = \unit{900}{\GeV}$
-XLabel=$n_{ch}$
-YLabel=$\langle p_{T} \rangle$
+Title=$\langle p_{T} \rangle$ vs. $n_\text{ch}$ at $\sqrt{s} = \unit{900}{\GeV}$
+XLabel=$n_\text{ch}$
+YLabel=$\langle p_{T} \rangle$ / GeV
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d08-x01-y01
-Title=$\langle p_{T} \rangle$ vs. $n_{ch}$ at $\sqrt{s} = \unit{63}{\GeV}$
-XLabel=$n_{ch}$
-YLabel=$\langle p_{T} \rangle$
+Title=$\langle p_{T} \rangle$ vs. $n_\text{ch}$ at $\sqrt{s} = \unit{63}{\GeV}$
+XLabel=$n_\text{ch}$
+YLabel=$\langle p_{T} \rangle$ / GeV
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d09-x01-y01
Title=Transverse energy cross section at $\sqrt{s} = \unit{200}{\GeV}$ and $|\eta|<6$
-XLabel=$\Sigma E_T$
-YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{\Sigma E_{t}}$
+XLabel=$\Sigma E_T$ / GeV
+YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{\Sigma E_{t}}$ / mb\,GeV$-1$
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d10-x01-y01
Title=Transverse energy cross section at $\sqrt{s} = \unit{500}{\GeV}$ and $|\eta|<6$
-XLabel=$\Sigma E_T$
-YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{\Sigma E_{t}}$
+XLabel=$\Sigma E_T$ / GeV
+YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{\Sigma E_{t}}$ / mb\,GeV$-1$
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d11-x01-y01
Title=Transverse energy cross section at $\sqrt{s} = \unit{900}{\GeV}$ and $|\eta|<6$
-XLabel=$\Sigma E_T$
-YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{\Sigma E_{t}}$
+XLabel=$\Sigma E_T$ / GeV
+YLabel=$\mathrm{d}{\sigma}/\mathrm{d}{\Sigma E_{t}}$ / mb\,GeV$-1$
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d12-x01-y01
-Title=$\langle \Sigma E_T \rangle$ vs. $n_{ch}$ at $\sqrt{s} = \unit{200}{\GeV}$ and $|\eta|<2.5$
-XLabel=$n_{ch}$
-YLabel=$\mathrm{d}\sigma/\mathrm{d}P_{t}}$
+Title=$\langle \Sigma E_T \rangle$ vs. $n_\text{ch}$ at $\sqrt{s} = \unit{200}{\GeV}$ and $|\eta|<2.5$
+XLabel=$n_\text{ch}$
+YLabel=$\Sigma E_T$ / GeV
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d12-x01-y02
-Title=$\langle \Sigma E_T \rangle$ vs. $n_{ch}$ at $\sqrt{s} = \unit{500}{\GeV}$ and $|\eta|<2.5$
-XLabel=$n_{ch}$
-YLabel=$\mathrm{d}\sigma/\mathrm{d}\phi}$
+Title=$\langle \Sigma E_T \rangle$ vs. $n_\text{ch}$ at $\sqrt{s} = \unit{500}{\GeV}$ and $|\eta|<2.5$
+XLabel=$n_\text{ch}$
+YLabel=$\Sigma E_T$ / GeV
# END PLOT
# BEGIN PLOT /UA1_1990_S2044935/d12-x01-y03
-Title=$\langle \Sigma E_T \rangle$ vs. $n_{ch}$ at $\sqrt{s} = \unit{900}{\GeV}$ and $|\eta|<2.5$
-XLabel=$n_{ch}$
-YLabel=$\mathrm{d}\sigma/\mathrm{d}\eta}$
+Title=$\langle \Sigma E_T \rangle$ vs. $n_\text{ch}$ at $\sqrt{s} = \unit{900}{\GeV}$ and $|\eta|<2.5$
+XLabel=$n_\text{ch}$
+YLabel=$\Sigma E_T$ / GeV
# END PLOT
Modified: trunk/src/Analyses/MC_LHC_SUSY.cc
==============================================================================
--- trunk/src/Analyses/MC_LHC_SUSY.cc Thu Oct 8 15:40:30 2009 (r1892)
+++ trunk/src/Analyses/MC_LHC_SUSY.cc Thu Oct 8 16:04:19 2009 (r1893)
@@ -97,6 +97,10 @@
_hist_mll_all_ossf_mumu = bookHistogram1D("mll-all-ossf-mumu", 50, 0.0, 500);
_hist_mll_all_osof_emu = bookHistogram1D("mll-all-osof-emu", 50, 0.0, 500);
+ _hist_mll_2_ossf_ee = bookHistogram1D("mll-2-ossf-ee", 50, 0.0, 500);
+ _hist_mll_2_ossf_mumu = bookHistogram1D("mll-2-ossf-mumu", 50, 0.0, 500);
+ _hist_mll_2_osof_emu = bookHistogram1D("mll-2-osof-emu", 50, 0.0, 500);
+
/// @todo LSP eta, pT, phi, mass: no reliable cross-scenario LSP PID but
/// maybe plot for all of chi^0_1, gravitino, sneutrino, gluino, ... or
/// identify the LSP as any PID::isSUSY (?) particle with status = 1?
@@ -232,20 +236,29 @@
if (eminus_ok && eplus_ok) {
const double m_ee = FourMomentum(peplus + peminus).mass();
_hist_mll_ossf_ee->fill(m_ee/GeV, weight);
+ if (epluses.size() == 1 && eminuses.size() == 1)
+ _hist_mll_2_ossf_ee->fill(m_ee/GeV, weight);
}
// m_mumu
if (muminus_ok && muplus_ok) {
const double m_mumu = FourMomentum(pmuplus + pmuminus).mass();
_hist_mll_ossf_mumu->fill(m_mumu/GeV, weight);
+ if (mupluses.size() == 1 && muminuses.size() == 1)
+ _hist_mll_2_ossf_mumu->fill(m_mumu/GeV, weight);
}
// m_emu (both configurations)
if (eminus_ok && muplus_ok) {
const double m_emu = FourMomentum(pmuplus + peminus).mass();
_hist_mll_osof_emu->fill(m_emu/GeV, weight);
+ if (mupluses.size() == 1 && eminuses.size() == 1)
+ _hist_mll_2_ossf_emu->fill(m_emu/GeV, weight);
+
}
if (muminus_ok && eplus_ok) {
const double m_mue = FourMomentum(peplus + pmuminus).mass();
_hist_mll_osof_emu->fill(m_mue/GeV, weight);
+ if (epluses.size() == 1 && muminuses.size() == 1)
+ _hist_mll_2_ossf_emu->fill(m_mue/GeV, weight);
}
@@ -297,6 +310,7 @@
AIDA::IHistogram1D *_hist_n_gamma, *_hist_phi_gamma, *_hist_eta_gamma, *_hist_pt_gamma;
AIDA::IHistogram1D *_hist_n_gammaiso, *_hist_phi_gammaiso, *_hist_eta_gammaiso, *_hist_pt_gammaiso;
AIDA::IHistogram1D *_hist_met;
+ AIDA::IHistogram1D *_hist_mll_2_ossf_ee, *_hist_mll_2_ossf_mumu, *_hist_mll_2_osof_emu;
AIDA::IHistogram1D *_hist_mll_ossf_ee, *_hist_mll_ossf_mumu, *_hist_mll_osof_emu;
AIDA::IHistogram1D *_hist_mll_all_ossf_ee, *_hist_mll_all_ossf_mumu, *_hist_mll_all_osof_emu;
};
More information about the Rivet-svn
mailing list
|