Page MenuHomeHEPForge
Files F10881346

No OneTemporary
Actions

View Options

diff --git a/analyses/pluginMC/MC_PARTONICTOPS.cc b/analyses/pluginMC/MC_PARTONICTOPS.cc
--- a/analyses/pluginMC/MC_PARTONICTOPS.cc
+++ b/analyses/pluginMC/MC_PARTONICTOPS.cc
@@ -1,108 +1,129 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/PartonicTops.hh"
namespace Rivet {
/// Find and plot partonic top properties (requires tops in event record)
class MC_PARTONICTOPS : public Analysis {
public:
/// Constructor
DEFAULT_RIVET_ANALYSIS_CTOR(MC_PARTONICTOPS);
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
// Initialise and register projections
declare(PartonicTops(PartonicTops::ALL), "AllTops");
declare(PartonicTops(PartonicTops::ALL, true, false, Cuts::OPEN, PartonicTops::FIRST), "AllTopsFirst"); ///< @todo API ick!
declare(PartonicTops(PartonicTops::E_MU), "LeptonicTops");
declare(PartonicTops(PartonicTops::HADRONIC), "HadronicTops");
// Book histograms
_h_tall_n = bookHisto1D("t_all_n", linspace(5, -0.5, 4.5));
_h_tall_pt = bookHisto1D("t_all_pT", logspace(50, 1, 500));
_h_tall_y = bookHisto1D("t_all_y", linspace(50, -5, 5));
_h_tall_n_first = bookHisto1D("t_all_n_firsttop", linspace(5, -0.5, 4.5));
_h_tall_pt_first = bookHisto1D("t_all_pT_firsttop", logspace(50, 1, 500));
_h_tall_y_first = bookHisto1D("t_all_y_firsttop", linspace(50, -5, 5));
+ _h_tall_pt_dfirstlast = bookHisto1D("t_all_pT_dfirstlast", linspace(100, -100, 100));
+ _p_tall_pt_dfirstlast = bookProfile1D("t_all_pT_dfirstlast_prof", logspace(50, 1, 500));
+ // _h_tall_y_dfirstlast = bookHisto1D("t_all_y_dfirstlast", linspace(50, -2, 2));
+ // _p_tall_y_dfirstlast = bookProfile1D("t_all_y_dfirstlast_prof", linspace(50, -5, 5));
+
_h_tlep_n = bookHisto1D("t_lep_n", linspace(5, -0.5, 4.5));
_h_tlep_pt = bookHisto1D("t_lep_pT", logspace(50, 1, 500));
_h_tlep_y = bookHisto1D("t_lep_y", linspace(50, -5, 5));
_h_thad_n = bookHisto1D("t_had_n", linspace(5, -0.5, 4.5));
_h_thad_pt = bookHisto1D("t_had_pT", logspace(50, 1, 500));
_h_thad_y = bookHisto1D("t_had_y", linspace(50, -5, 5));
}
/// Perform the per-event analysis
void analyze(const Event& event) {
+ const double w = event.weight();
+ // Last tops (standard)
const Particles& alltops = apply<PartonicTops>(event, "AllTops").particlesByPt();
- _h_tall_n->fill(alltops.size(), event.weight());
+ _h_tall_n->fill(alltops.size(), w);
for (const Particle& t : alltops) {
- _h_tall_pt->fill(t.pT()/GeV, event.weight());
- _h_tall_y->fill(t.rap(), event.weight());
+ _h_tall_pt->fill(t.pT()/GeV, w);
+ _h_tall_y->fill(t.rap(), w);
}
+ // First tops
const Particles& alltops_first = apply<PartonicTops>(event, "AllTopsFirst").particlesByPt();
- _h_tall_n_first->fill(alltops_first.size(), event.weight());
+ _h_tall_n_first->fill(alltops_first.size(), w);
for (const Particle& t : alltops_first) {
- _h_tall_pt_first->fill(t.pT()/GeV, event.weight());
- _h_tall_y_first->fill(t.rap(), event.weight());
+ _h_tall_pt_first->fill(t.pT()/GeV, w);
+ _h_tall_y_first->fill(t.rap(), w);
}
- const Particles& leptops = apply<PartonicTops>(event, "LeptonicTops").particlesByPt();
- _h_tlep_n->fill(leptops.size(), event.weight());
- for (const Particle& t : leptops) {
- _h_tlep_pt->fill(t.pT()/GeV, event.weight());
- _h_tlep_y->fill(t.rap(), event.weight());
+ // Match first and last tops
+ for (const Particle& tf : alltops_first) {
+ for (const Particle& tl : alltops) {
+ //if (deltaR(tf, tl) > 1) continue;
+ if (tf.pid() != tl.pid()) continue;
+ const double dpt = tl.pT() - tf.pT(); //< defined as change due to PS
+ _h_tall_pt_dfirstlast->fill(dpt/GeV, w);
+ _p_tall_pt_dfirstlast->fill(tf.pT()/GeV, fabs(dpt)/GeV, w);
+ }
}
- const Particles& hadtops = apply<PartonicTops>(event, "HadronicTops").particlesByPt();
- _h_thad_n->fill(hadtops.size(), event.weight());
- for (const Particle& t : hadtops) {
- _h_thad_pt->fill(t.pT()/GeV, event.weight());
- _h_thad_y->fill(t.rap(), event.weight());
+ // Leptonic (last) tops
+ const Particles& leptops = apply<PartonicTops>(event, "LeptonicTops").particlesByPt();
+ _h_tlep_n->fill(leptops.size(), w);
+ for (const Particle& t : leptops) {
+ _h_tlep_pt->fill(t.pT()/GeV, w);
+ _h_tlep_y->fill(t.rap(), w);
}
+ // Hadronic (last) tops
+ const Particles& hadtops = apply<PartonicTops>(event, "HadronicTops").particlesByPt();
+ _h_thad_n->fill(hadtops.size(), w);
+ for (const Particle& t : hadtops) {
+ _h_thad_pt->fill(t.pT()/GeV, w);
+ _h_thad_y->fill(t.rap(), w);
+ }
}
/// Normalise histograms etc., after the run
void finalize() {
normalize({_h_tall_n, _h_tall_n_first, _h_tlep_n, _h_thad_n});
normalize({_h_tall_pt, _h_tall_pt_first, _h_tlep_pt, _h_thad_pt});
normalize({_h_tall_y, _h_tall_y_first, _h_tlep_y, _h_thad_y});
+ normalize(_h_tall_pt_dfirstlast);
}
//@}
/// @name Histograms
//@{
Histo1DPtr _h_tall_n, _h_tall_n_first, _h_tlep_n, _h_thad_n;
Histo1DPtr _h_tall_pt, _h_tall_pt_first, _h_tlep_pt, _h_thad_pt;
Histo1DPtr _h_tall_y, _h_tall_y_first, _h_tlep_y, _h_thad_y;
+ Histo1DPtr _h_tall_pt_dfirstlast;
+ Profile1DPtr _p_tall_pt_dfirstlast;
//@}
};
- // The hook for the plugin system
DECLARE_RIVET_PLUGIN(MC_PARTONICTOPS);
-
}
diff --git a/analyses/pluginMC/MC_PARTONICTOPS.plot b/analyses/pluginMC/MC_PARTONICTOPS.plot
--- a/analyses/pluginMC/MC_PARTONICTOPS.plot
+++ b/analyses/pluginMC/MC_PARTONICTOPS.plot
@@ -1,45 +1,59 @@
# BEGIN PLOT /MC_PARTONICTOPS/t_.*_n
XLabel=$n$
YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}n$
# END PLOT
# BEGIN PLOT /MC_PARTONICTOPS/t_.*_pT
XLabel=$p_\perp$~[GeV]
YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}p_\perp$~[GeV$^{-1}$]
# END PLOT
# BEGIN PLOT /MC_PARTONICTOPS/t_.*_y
XLabel=$y$
YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}y$
# END PLOT
# BEGIN PLOT /MC_PARTONICTOPS/t_all_n
Title=Top-quark per-event multiplicity
# END PLOT
# BEGIN PLOT /MC_PARTONICTOPS/t_all_pT
Title=Top-quark transverse momentum spectrum
# END PLOT
# BEGIN PLOT /MC_PARTONICTOPS/t_all_y
Title=Top-quark rapidity distribution
# END PLOT
+# BEGIN PLOT /MC_PARTONICTOPS/t_all_pT_dfirstlast$
+Title=First-to-last top-quark transverse momentum difference
+XLabel=$\Delta p_\perp = p_\perp^\mathrm{last} - p_\perp^\mathrm{first}$~[GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}\Delta p_\perp$~[GeV$^{-1}$]
+# END PLOT
+
+# BEGIN PLOT /MC_PARTONICTOPS/t_all_pT_dfirstlast_prof
+Title=First-to-last top-quark $p_\perp$ difference vs $p_\perp^\mathrm{first}$
+XLabel=$p_\perp^\mathrm{first}$~[GeV]
+YLabel=$\langle \Delta p_\perp \rangle$~[GeV]
+LogY=0
+LegendYPos=0.2
+# END PLOT
+
# BEGIN PLOT /MC_PARTONICTOPS/t_lep_n
Title=Leptonic top-quark per-event multiplicity
# END PLOT
# BEGIN PLOT /MC_PARTONICTOPS/t_lep_pT
Title=Leptonic top-quark transverse momentum spectrum
# END PLOT
# BEGIN PLOT /MC_PARTONICTOPS/t_lep_y
Title=Leptonic top-quark rapidity distribution
# END PLOT
# BEGIN PLOT /MC_PARTONICTOPS/t_had_n
Title=Hadronic top-quark per-event multiplicity
# END PLOT
# BEGIN PLOT /MC_PARTONICTOPS/t_had_pT
Title=Hadronic top-quark transverse momentum spectrum
# END PLOT
# BEGIN PLOT /MC_PARTONICTOPS/t_had_y
Title=Hadronic top-quark rapidity distribution
# END PLOT
diff --git a/analyses/pluginMC/MC_TTBAR.cc b/analyses/pluginMC/MC_TTBAR.cc
--- a/analyses/pluginMC/MC_TTBAR.cc
+++ b/analyses/pluginMC/MC_TTBAR.cc
@@ -1,300 +1,297 @@
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
#include "Rivet/Projections/ChargedLeptons.hh"
#include "Rivet/Projections/MissingMomentum.hh"
#include "Rivet/Projections/FastJets.hh"
#include "Rivet/AnalysisLoader.hh"
namespace Rivet {
-
-
class MC_TTBAR : public Analysis {
public:
/// Minimal constructor
- MC_TTBAR() : Analysis("MC_TTBAR")
- {
- }
+ DEFAULT_RIVET_ANALYSIS_CTOR(MC_TTBAR);
/// @name Analysis methods
//@{
/// Set up projections and book histograms
void init() {
// A FinalState is used to select particles within |eta| < 4.2 and with pT
// > 30 GeV, out of which the ChargedLeptons projection picks only the
// electrons and muons, to be accessed later as "LFS".
- ChargedLeptons lfs(FinalState(-4.2, 4.2, 30*GeV));
+ ChargedLeptons lfs(FinalState(Cuts::abseta < 4.2 && Cuts::pT > 30*GeV));
declare(lfs, "LFS");
+
// A second FinalState is used to select all particles in |eta| < 4.2,
// with no pT cut. This is used to construct jets and measure missing
// transverse energy.
- VetoedFinalState fs(FinalState(-4.2, 4.2, 0*GeV));
+ VetoedFinalState fs(FinalState(Cuts::abseta < 4.2));
fs.addVetoOnThisFinalState(lfs);
declare(FastJets(fs, FastJets::ANTIKT, 0.6), "Jets");
declare(MissingMomentum(fs), "MissingET");
// Booking of histograms
_h_njets = bookHisto1D("jet_mult", 11, -0.5, 10.5);
//
_h_jet_1_pT = bookHisto1D("jet_1_pT", logspace(50, 20.0, 500.0));
_h_jet_2_pT = bookHisto1D("jet_2_pT", logspace(50, 20.0, 400.0));
_h_jet_3_pT = bookHisto1D("jet_3_pT", logspace(50, 20.0, 300.0));
_h_jet_4_pT = bookHisto1D("jet_4_pT", logspace(50, 20.0, 200.0));
_h_jet_HT = bookHisto1D("jet_HT", logspace(50, 100.0, 2000.0));
//
_h_bjet_1_pT = bookHisto1D("jetb_1_pT", logspace(50, 20.0, 400.0));
_h_bjet_2_pT = bookHisto1D("jetb_2_pT", logspace(50, 20.0, 300.0));
//
_h_ljet_1_pT = bookHisto1D("jetl_1_pT", logspace(50, 20.0, 400.0));
_h_ljet_2_pT = bookHisto1D("jetl_2_pT", logspace(50, 20.0, 300.0));
//
_h_W_mass = bookHisto1D("W_mass", 75, 30, 180);
_h_t_mass = bookHisto1D("t_mass", 150, 130, 430);
_h_t_mass_W_cut = bookHisto1D("t_mass_W_cut", 150, 130, 430);
//
_h_jetb_1_jetb_2_dR = bookHisto1D("jetb_1_jetb_2_dR", 20, 0.0, 7.0);
_h_jetb_1_jetb_2_deta = bookHisto1D("jetb_1_jetb_2_deta", 20, 0.0, 7.0);
_h_jetb_1_jetb_2_dphi = bookHisto1D("jetb_1_jetb_2_dphi", 20, 0.0, M_PI);
_h_jetb_1_jetl_1_dR = bookHisto1D("jetb_1_jetl_1_dR", 20, 0.0, 7.0);
_h_jetb_1_jetl_1_deta = bookHisto1D("jetb_1_jetl_1_deta", 20, 0.0, 7.0);
_h_jetb_1_jetl_1_dphi = bookHisto1D("jetb_1_jetl_1_dphi", 20, 0.0, M_PI);
_h_jetl_1_jetl_2_dR = bookHisto1D("jetl_1_jetl_2_dR", 20, 0.0, 7.0);
_h_jetl_1_jetl_2_deta = bookHisto1D("jetl_1_jetl_2_deta", 20, 0.0, 7.0);
_h_jetl_1_jetl_2_dphi = bookHisto1D("jetl_1_jetl_2_dphi", 20, 0.0, M_PI);
_h_jetb_1_W_dR = bookHisto1D("jetb_1_W_dR", 20, 0.0, 7.0);
_h_jetb_1_W_deta = bookHisto1D("jetb_1_W_deta", 20, 0.0, 7.0);
_h_jetb_1_W_dphi = bookHisto1D("jetb_1_W_dphi", 20, 0.0, M_PI);
_h_jetb_1_l_dR = bookHisto1D("jetb_1_l_dR", 20, 0.0, 7.0);
_h_jetb_1_l_deta = bookHisto1D("jetb_1_l_deta", 20, 0.0, 7.0);
_h_jetb_1_l_dphi = bookHisto1D("jetb_1_l_dphi", 20, 0.0, M_PI);
_h_jetb_1_l_mass = bookHisto1D("jetb_1_l_mass", 40, 0.0, 500.0);
}
void analyze(const Event& event) {
const double weight = event.weight();
// Use the "LFS" projection to require at least one hard charged
// lepton. This is an experimental signature for the leptonically decaying
// W. This helps to reduce pure QCD backgrounds.
const ChargedLeptons& lfs = apply<ChargedLeptons>(event, "LFS");
MSG_DEBUG("Charged lepton multiplicity = " << lfs.chargedLeptons().size());
foreach (const Particle& lepton, lfs.chargedLeptons()) {
MSG_DEBUG("Lepton pT = " << lepton.pT());
}
if (lfs.chargedLeptons().empty()) {
MSG_DEBUG("Event failed lepton multiplicity cut");
vetoEvent;
}
// Use a missing ET cut to bias toward events with a hard neutrino from
// the leptonically decaying W. This helps to reduce pure QCD backgrounds.
const MissingMomentum& met = apply<MissingMomentum>(event, "MissingET");
MSG_DEBUG("Vector ET = " << met.vectorEt().mod() << " GeV");
if (met.vectorEt().mod() < 30*GeV) {
MSG_DEBUG("Event failed missing ET cut");
vetoEvent;
}
// Use the "Jets" projection to check that there are at least 4 jets of
// any pT. Getting the jets sorted by pT ensures that the first jet is the
// hardest, and so on. We apply no pT cut here only because we want to
// plot all jet pTs to help optimise our jet pT cut.
const FastJets& jetpro = apply<FastJets>(event, "Jets");
const Jets alljets = jetpro.jetsByPt();
if (alljets.size() < 4) {
MSG_DEBUG("Event failed jet multiplicity cut");
vetoEvent;
}
// Update passed-cuts counter and fill all-jets histograms
_h_jet_1_pT->fill(alljets[0].pT()/GeV, weight);
_h_jet_2_pT->fill(alljets[1].pT()/GeV, weight);
_h_jet_3_pT->fill(alljets[2].pT()/GeV, weight);
_h_jet_4_pT->fill(alljets[3].pT()/GeV, weight);
// Insist that the hardest 4 jets pass pT hardness cuts. If we don't find
// at least 4 such jets, we abandon this event.
const Jets jets = jetpro.jetsByPt(30*GeV);
_h_njets->fill(jets.size(), weight);
double ht = 0.0;
foreach (const Jet& j, jets) { ht += j.pT(); }
_h_jet_HT->fill(ht/GeV, weight);
if (jets.size() < 4 ||
jets[0].pT() < 60*GeV ||
jets[1].pT() < 50*GeV ||
jets[3].pT() < 30*GeV) {
MSG_DEBUG("Event failed jet cuts");
vetoEvent;
}
// Sort the jets into b-jets and light jets. We expect one hard b-jet from
// each top decay, so our 4 hardest jets should include two b-jets. The
// Jet::bTagged() method is equivalent to perfect experimental
// b-tagging, in a generator-independent way.
Jets bjets, ljets;
foreach (const Jet& jet, jets) {
// // Don't count jets that overlap with the hard leptons
bool isolated = true;
foreach (const Particle& lepton, lfs.chargedLeptons()) {
if (deltaR(jet.momentum(), lepton.momentum()) < 0.3) {
isolated = false;
break;
}
}
if (!isolated) {
MSG_DEBUG("Jet failed lepton isolation cut");
break;
}
if (jet.bTagged()) {
bjets.push_back(jet);
} else {
ljets.push_back(jet);
}
}
MSG_DEBUG("Number of b-jets = " << bjets.size());
MSG_DEBUG("Number of l-jets = " << ljets.size());
if (bjets.size() != 2) {
MSG_DEBUG("Event failed post-lepton-isolation b-tagging cut");
vetoEvent;
}
if (ljets.size() < 2) {
MSG_DEBUG("Event failed since not enough light jets remaining after lepton-isolation");
vetoEvent;
}
// Plot the pTs of the identified jets.
_h_bjet_1_pT->fill(bjets[0].pT(), weight);
_h_bjet_2_pT->fill(bjets[1].pT(), weight);
_h_ljet_1_pT->fill(ljets[0].pT(), weight);
_h_ljet_2_pT->fill(ljets[1].pT(), weight);
// Construct the hadronically decaying W momentum 4-vector from pairs of
// non-b-tagged jets. The pair which best matches the W mass is used. We start
// with an always terrible 4-vector estimate which should always be "beaten" by
// a real jet pair.
FourMomentum W(10*(sqrtS()>0.?sqrtS():14000.), 0, 0, 0);
for (size_t i = 0; i < ljets.size()-1; ++i) {
for (size_t j = i + 1; j < ljets.size(); ++j) {
const FourMomentum Wcand = ljets[i].momentum() + ljets[j].momentum();
MSG_TRACE(i << "," << j << ": candidate W mass = " << Wcand.mass()/GeV
<< " GeV, vs. incumbent candidate with " << W.mass()/GeV << " GeV");
if (fabs(Wcand.mass() - 80.4*GeV) < fabs(W.mass() - 80.4*GeV)) {
W = Wcand;
}
}
}
MSG_DEBUG("Candidate W mass = " << W.mass() << " GeV");
// There are two b-jets with which this can be combined to make the
// hadronically decaying top, one of which is correct and the other is
// not... but we have no way to identify which is which, so we construct
// both possible top momenta and fill the histograms with both.
const FourMomentum t1 = W + bjets[0].momentum();
const FourMomentum t2 = W + bjets[1].momentum();
_h_W_mass->fill(W.mass(), weight);
_h_t_mass->fill(t1.mass(), weight);
_h_t_mass->fill(t2.mass(), weight);
// Placing a cut on the well-known W mass helps to reduce backgrounds
if (inRange(W.mass()/GeV, 75.0, 85.0)) {
MSG_DEBUG("W found with mass " << W.mass()/GeV << " GeV");
_h_t_mass_W_cut->fill(t1.mass(), weight);
_h_t_mass_W_cut->fill(t2.mass(), weight);
_h_jetb_1_jetb_2_dR->fill(deltaR(bjets[0].momentum(), bjets[1].momentum()),weight);
_h_jetb_1_jetb_2_deta->fill(fabs(bjets[0].eta()-bjets[1].eta()),weight);
_h_jetb_1_jetb_2_dphi->fill(deltaPhi(bjets[0].momentum(),bjets[1].momentum()),weight);
_h_jetb_1_jetl_1_dR->fill(deltaR(bjets[0].momentum(), ljets[0].momentum()),weight);
_h_jetb_1_jetl_1_deta->fill(fabs(bjets[0].eta()-ljets[0].eta()),weight);
_h_jetb_1_jetl_1_dphi->fill(deltaPhi(bjets[0].momentum(),ljets[0].momentum()),weight);
_h_jetl_1_jetl_2_dR->fill(deltaR(ljets[0].momentum(), ljets[1].momentum()),weight);
_h_jetl_1_jetl_2_deta->fill(fabs(ljets[0].eta()-ljets[1].eta()),weight);
_h_jetl_1_jetl_2_dphi->fill(deltaPhi(ljets[0].momentum(),ljets[1].momentum()),weight);
_h_jetb_1_W_dR->fill(deltaR(bjets[0].momentum(), W),weight);
_h_jetb_1_W_deta->fill(fabs(bjets[0].eta()-W.eta()),weight);
_h_jetb_1_W_dphi->fill(deltaPhi(bjets[0].momentum(),W),weight);
FourMomentum l=lfs.chargedLeptons()[0].momentum();
_h_jetb_1_l_dR->fill(deltaR(bjets[0].momentum(), l),weight);
_h_jetb_1_l_deta->fill(fabs(bjets[0].eta()-l.eta()),weight);
_h_jetb_1_l_dphi->fill(deltaPhi(bjets[0].momentum(),l),weight);
_h_jetb_1_l_mass->fill(FourMomentum(bjets[0].momentum()+l).mass(), weight);
}
}
void finalize() {
normalize(_h_njets);
normalize(_h_jet_1_pT);
normalize(_h_jet_2_pT);
normalize(_h_jet_3_pT);
normalize(_h_jet_4_pT);
normalize(_h_jet_HT);
normalize(_h_bjet_1_pT);
normalize(_h_bjet_2_pT);
normalize(_h_ljet_1_pT);
normalize(_h_ljet_2_pT);
normalize(_h_W_mass);
normalize(_h_t_mass);
normalize(_h_t_mass_W_cut);
normalize(_h_jetb_1_jetb_2_dR);
normalize(_h_jetb_1_jetb_2_deta);
normalize(_h_jetb_1_jetb_2_dphi);
normalize(_h_jetb_1_jetl_1_dR);
normalize(_h_jetb_1_jetl_1_deta);
normalize(_h_jetb_1_jetl_1_dphi);
normalize(_h_jetl_1_jetl_2_dR);
normalize(_h_jetl_1_jetl_2_deta);
normalize(_h_jetl_1_jetl_2_dphi);
normalize(_h_jetb_1_W_dR);
normalize(_h_jetb_1_W_deta);
normalize(_h_jetb_1_W_dphi);
normalize(_h_jetb_1_l_dR);
normalize(_h_jetb_1_l_deta);
normalize(_h_jetb_1_l_dphi);
normalize(_h_jetb_1_l_mass);
}
//@}
private:
// @name Histogram data members
//@{
Histo1DPtr _h_njets;
Histo1DPtr _h_jet_1_pT, _h_jet_2_pT, _h_jet_3_pT, _h_jet_4_pT;
Histo1DPtr _h_jet_HT;
Histo1DPtr _h_bjet_1_pT, _h_bjet_2_pT;
Histo1DPtr _h_ljet_1_pT, _h_ljet_2_pT;
Histo1DPtr _h_W_mass;
Histo1DPtr _h_t_mass, _h_t_mass_W_cut;
Histo1DPtr _h_jetb_1_jetb_2_dR, _h_jetb_1_jetb_2_deta, _h_jetb_1_jetb_2_dphi;
Histo1DPtr _h_jetb_1_jetl_1_dR, _h_jetb_1_jetl_1_deta, _h_jetb_1_jetl_1_dphi;
Histo1DPtr _h_jetl_1_jetl_2_dR, _h_jetl_1_jetl_2_deta, _h_jetl_1_jetl_2_dphi;
Histo1DPtr _h_jetb_1_W_dR, _h_jetb_1_W_deta, _h_jetb_1_W_dphi;
Histo1DPtr _h_jetb_1_l_dR, _h_jetb_1_l_deta, _h_jetb_1_l_dphi,_h_jetb_1_l_mass;
//@}
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(MC_TTBAR);
}

File Metadata

Mime Type
text/x-diff
Expires
Sat, May 3, 6:13 AM (22 h, 27 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4968060
Default Alt Text
(20 KB)

Event Timeline