|
|
[Rivet-svn] r3357 - in trunk: data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Wed Sep 14 16:01:27 BST 2011
Author: buckley
Date: Wed Sep 14 16:01:27 2011
New Revision: 3357
Log:
more MC_TTBAR tweaks
Modified:
trunk/data/plotinfo/MC_TTBAR.plot
trunk/src/Analyses/MC_TTBAR.cc
Modified: trunk/data/plotinfo/MC_TTBAR.plot
==============================================================================
--- trunk/data/plotinfo/MC_TTBAR.plot Wed Sep 14 15:09:30 2011 (r3356)
+++ trunk/data/plotinfo/MC_TTBAR.plot Wed Sep 14 16:01:27 2011 (r3357)
@@ -1,78 +1,72 @@
# BEGIN PLOT /MC_TTBAR/jet_1_pT
-Title=Transverse momentum distribution for Jet 1
-XLabel=$p_\perp$
-YLabel=$\frac{dN}{d p_\perp}$ [GeV]
+Title=Transverse momentum distribution for jet 1
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}p_\perp$ [GeV$^-1$]
# END PLOT
# BEGIN PLOT /MC_TTBAR/jet_2_pT
-Title=Transverse momentum distribution for Jet 2
-XLabel=$p_\perp$
-YLabel=$\frac{dN}{d p_\perp}$ [GeV]
+Title=Transverse momentum distribution for jet 2
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}p_\perp$ [GeV$^-1$]
# END PLOT
# BEGIN PLOT /MC_TTBAR/jet_3_pT
-Title=Transverse momentum distribution for Jet 3
-XLabel=$p_\perp$
-YLabel=$\frac{dN}{d p_\perp}$ [GeV]
+Title=Transverse momentum distribution for jet 3
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}p_\perp$ [GeV$^-1$]
# END PLOT
# BEGIN PLOT /MC_TTBAR/jet_4_pT
-Title=Transverse momentum distribution for Jet 4
-XLabel=$p_\perp$
-YLabel=$\frac{dN}{d p_\perp}$ [GeV]
+Title=Transverse momentum distribution for jet 4
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}p_\perp$ [GeV$^-1$]
+# END PLOT
+
+# BEGIN PLOT /MC_TTBAR/jet_HT
+Title=$H_T$ distribution for all jets
+XLabel=$H_T$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}H_T$ [GeV$^-1$]
+LogX=1
# END PLOT
# BEGIN PLOT /MC_TTBAR/jetb_1_pT
-Title=Transverse momentum distribution for b-jet 1
-XLabel=$p_\perp$
-YLabel=$\frac{dN}{d p_\perp}$ [GeV]
+Title=Transverse momentum distribution for $b$-jet 1
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}p_\perp$ [GeV$^-1$]
# END PLOT
# BEGIN PLOT /MC_TTBAR/jetb_2_pT
-Title=Transverse momentum distribution for b-jet 2
-XLabel=$p_\perp$
-YLabel=$\frac{dN}{d p_\perp}$ [GeV]
+Title=Transverse momentum distribution for $b$-jet 2
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}p_\perp$ [GeV$^-1$]
# END PLOT
#
# BEGIN PLOT /MC_TTBAR/jetl_1_pT
Title=Transverse momentum distribution for light jet 1
-XLabel=$p_\perp$
-YLabel=$\frac{dN}{d p_\perp}$ [GeV]
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}p_\perp$ [GeV$^-1$]
# END PLOT
# BEGIN PLOT /MC_TTBAR/jetl_2_pT
Title=Transverse momentum distribution for light jet 2
-XLabel=$p_\perp$
-YLabel=$\frac{dN}{d p_\perp}$ [GeV]
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}p_\perp$ [GeV$^-1$]
# END PLOT
# BEGIN PLOT /MC_TTBAR/W_mass
-Title=Mass distribution for W bosons
-XLabel=$m_W$
-YLabel=$\frac{dN}{d m_W}$ [GeV]
+Title=Mass distribution for $W$ bosons
+XLabel=$m_{jj}$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}m_{jj}$ [GeV$^-1$]
# END PLOT
# BEGIN PLOT /MC_TTBAR/t_mass
-Title=Mass distribution for reconstructed top quark
-XLabel=$m_t$
-YLabel=$\frac{dN}{d m_t}$ [GeV]
+Title=Mass distribution for reconstructed top
+XLabel=$m_{q\bar{q}b$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}m_{q\bar{q}b}$ [GeV$^-1$]
# END PLOT
# BEGIN PLOT /MC_TTBAR/t_mass_W_cut
-Title=Mass distribution for reconstructed top quark after $m_W$ cut
-XLabel=$m_t$
-YLabel=$\frac{dN}{d m_t}$ [GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_TTBAR/W_comb_mass
-Title=Combinatoric W mass
-XLabel=$m_W$
-YLabel=$\frac{dN}{d m_W}$ [GeV]
+Title=Mass distribution for reconstructed top after $m_W$ cut
+XLabel=$m_{q\bar{q}b$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}m_{q\bar{q}b}$ [GeV$^-1$]
# END PLOT
-
-# BEGIN PLOT /MC_TTBAR/t_comb_mass
-Title=Combinatoric t mass
-XLabel=$m_t$
-YLabel=$\frac{dN}{d m_t}$ [GeV]
-# END PLOT
-
Modified: trunk/src/Analyses/MC_TTBAR.cc
==============================================================================
--- trunk/src/Analyses/MC_TTBAR.cc Wed Sep 14 15:09:30 2011 (r3356)
+++ trunk/src/Analyses/MC_TTBAR.cc Wed Sep 14 16:01:27 2011 (r3357)
@@ -43,19 +43,17 @@
_h_jet_2_pT = bookHistogram1D("jet_2_pT", 50, 0, 400);
_h_jet_3_pT = bookHistogram1D("jet_3_pT", 50, 0, 300);
_h_jet_4_pT = bookHistogram1D("jet_4_pT", 50, 0, 200);
- _h_jet_HT = bookHistogram1D("jet_HT", 50, 0, 1000);
+ _h_jet_HT = bookHistogram1D("jet_HT", logspace(0, 2000, 50));
//
- _h_bjet_1_pT = bookHistogram1D("jetb_1_pT", 50, 0, 250);
- _h_bjet_2_pT = bookHistogram1D("jetb_2_pT", 50, 0, 250);
+ _h_bjet_1_pT = bookHistogram1D("jetb_1_pT", 50, 0, 400);
+ _h_bjet_2_pT = bookHistogram1D("jetb_2_pT", 50, 0, 300);
//
- _h_ljet_1_pT = bookHistogram1D("jetl_1_pT", 50, 0, 250);
- _h_ljet_2_pT = bookHistogram1D("jetl_2_pT", 50, 0, 250);
+ _h_ljet_1_pT = bookHistogram1D("jetl_1_pT", 50, 0, 400);
+ _h_ljet_2_pT = bookHistogram1D("jetl_2_pT", 50, 0, 300);
//
_h_W_mass = bookHistogram1D("W_mass", 75, 30, 180);
_h_t_mass = bookHistogram1D("t_mass", 150, 130, 430);
_h_t_mass_W_cut = bookHistogram1D("t_mass_W_cut", 150, 130, 430);
- // _h_W_comb_mass = bookHistogram1D("W_comb_mass", 75, 30, 180);
- // _h_t_comb_mass = bookHistogram1D("t_comb_mass", 150, 130, 430);
}
@@ -184,21 +182,7 @@
_h_t_mass->fill(t1.mass(), weight);
_h_t_mass->fill(t2.mass(), weight);
- // // Fill histograms for all of the combinatoric 2-jet masses
- // _h_W_comb_mass->fill(mass(jets[0].momentum() + jets[1].momentum()), weight);
- // _h_W_comb_mass->fill(mass(jets[0].momentum() + jets[2].momentum()), weight);
- // _h_W_comb_mass->fill(mass(jets[0].momentum() + jets[3].momentum()), weight);
- // _h_W_comb_mass->fill(mass(jets[1].momentum() + jets[2].momentum()), weight);
- // _h_W_comb_mass->fill(mass(jets[1].momentum() + jets[3].momentum()), weight);
- // _h_W_comb_mass->fill(mass(jets[2].momentum() + jets[3].momentum()), weight);
- // // Fill histograms for all of the combinatoric 3-jet masses
- // _h_t_comb_mass->fill(mass(jets[0].momentum() + jets[1].momentum() + jets[2].momentum()), weight);
- // _h_t_comb_mass->fill(mass(jets[0].momentum() + jets[1].momentum() + jets[3].momentum()), weight);
- // _h_t_comb_mass->fill(mass(jets[0].momentum() + jets[2].momentum() + jets[3].momentum()), weight);
- // _h_t_comb_mass->fill(mass(jets[1].momentum() + jets[2].momentum() + jets[3].momentum()), weight);
-
- // Placing a cut on the W mass (which is well known) helps to reduce
- // backgrounds.
+ // Placing a cut on the well-known W mass helps to reduce backgrounds
if (inRange(W.mass()/GeV, 75, 85)) {
MSG_DEBUG("W found with mass " << W.mass()/GeV << " GeV");
_sumwPassedWMass += weight;
@@ -219,8 +203,6 @@
scale(_h_bjet_2_pT, 1/_sumwPassedJetID);
scale(_h_ljet_1_pT, 1/_sumwPassedJetID);
scale(_h_ljet_2_pT, 1/_sumwPassedJetID);
- // scale(_h_W_comb_mass, 1/_sumwPassedJetID);
- // scale(_h_t_comb_mass, 1/_sumwPassedJetID);
scale(_h_W_mass, 1/_sumwPassedJetID);
scale(_h_t_mass, 1/_sumwPassedJetID);
scale(_h_t_mass_W_cut, 1/_sumwPassedWMass);
@@ -242,10 +224,7 @@
AIDA::IHistogram1D *_h_bjet_1_pT, *_h_bjet_2_pT;
AIDA::IHistogram1D *_h_ljet_1_pT, *_h_ljet_2_pT;
AIDA::IHistogram1D *_h_W_mass;
- AIDA::IHistogram1D *_h_t_mass;
- // AIDA::IHistogram1D *_h_W_comb_mass;
- // AIDA::IHistogram1D *_h_t_comb_mass;
- AIDA::IHistogram1D *_h_t_mass_W_cut;
+ AIDA::IHistogram1D *_h_t_mass, *_h_t_mass_W_cut;
//@}
More information about the Rivet-svn
mailing list
|