Rivet analyses
title: CDF_2001_I564673
Field \& Stuart Run I underlying event analysis.
Experiment: CDF (Tevatron Run 1)
Inspire ID: 564673
Status: VALIDATED
Authors: - Andy Buckley
References: - Phys.Rev.D65:092002,2002 - FNAL-PUB 01/211-E
Beams: p- p+
Beam energies: (900.0, 900.0)GeV
Run details: - $p\bar{p}$ QCD interactions at 1800 GeV. The leading jet is binned from 0--49 GeV, and histos can usually can be filled with a single generator run without kinematic sub-samples.
The original CDF underlying event analysis, based on decomposing each event into a transverse structure with toward'',away'' and `transverse'' regions defined relative to the azimuthal direction of the leading jet in the event. Since the toward region is by definition dominated by the hard process, as is the away region by momentum balance in the matrix element, the transverse region is most sensitive to multi-parton interactions. The transverse regions occupy $|\phi| \in [60\degree, 120\degree]$ for $|\eta| < 1$. The $\pT$ ranges for the leading jet are divided experimentally into themin-bias' sample from 0--20 GeV, and the `JET20' sample from 18--49 GeV.
Source code:CDF_2001_I564673.cc
```c++ // -- C++ --
include "Rivet/Analysis.hh"
include "Rivet/Projections/ChargedFinalState.hh"
include "Rivet/Projections/ConstLossyFinalState.hh"
include "Rivet/Projections/FastJets.hh"
include "Rivet/Projections/TriggerCDFRun0Run1.hh"
//#include "Rivet/Projections/SmearedParticles.hh"
namespace Rivet {
/// @brief Field-Stuart CDF Run I track-jet underlying event analysis /// /// @author Andy Buckley /// /// The "original" underlying event analysis, using a non-standard track-jet algorithm. /// /// @par Run conditions /// /// @arg \f$ \sqrt{s} = \f$ 1800 GeV /// @arg Run with generic QCD events. /// @arg Several \f$ p_\perp^\text{min} \f$ cutoffs are probably required to fill the profile histograms: /// @arg \f$ p_\perp^\text{min} = \f$ 0 (min bias), 10, 20 GeV class CDF_2001_I564673 : public Analysis { public:
RIVET_DEFAULT_ANALYSIS_CTOR(CDF_2001_I564673);
/// @name Analysis methods
/// @{
// Book histograms
void init() {
declare(TriggerCDFRun0Run1(), "Trigger");
// Randomly discard 8% of charged particles as a kind of hacky detector correction.
const ChargedFinalState cfs(Cuts::abseta < 1.0 && Cuts::pT > 0.5 * GeV);
/// @todo Replace ConstLossyFinalState with SmearedParticles
const ConstLossyFinalState lossyfs(cfs, 0.08);
//const SmearedParticles lossyfs(cfs, [](const Particle&){ return 0.92; });
declare(lossyfs, "FS");
declare(FastJets(lossyfs, JetAlg::TRACKJET, 0.7), "TrackJet");
book(_numvsDeltaPhi2, 1, 1, 1);
book(_numvsDeltaPhi5, 1, 1, 2);
book(_numvsDeltaPhi30, 1, 1, 3);
book(_pTvsDeltaPhi2, 2, 1, 1);
book(_pTvsDeltaPhi5, 2, 1, 2);
book(_pTvsDeltaPhi30, 2, 1, 3);
book(_numTowardMB, 3, 1, 1);
book(_numTransMB, 3, 1, 2);
book(_numAwayMB, 3, 1, 3);
book(_numTowardJ20, 4, 1, 1);
book(_numTransJ20, 4, 1, 2);
book(_numAwayJ20, 4, 1, 3);
book(_ptsumTowardMB, 5, 1, 1);
book(_ptsumTransMB, 5, 1, 2);
book(_ptsumAwayMB, 5, 1, 3);
book(_ptsumTowardJ20, 6, 1, 1);
book(_ptsumTransJ20, 6, 1, 2);
book(_ptsumAwayJ20, 6, 1, 3);
book(_ptTrans2, 7, 1, 1);
book(_ptTrans5, 7, 1, 2);
book(_ptTrans30, 7, 1, 3);
book(_totalNumTrans2, "totalNumTrans2");
book(_totalNumTrans5, "totalNumTrans5");
book(_totalNumTrans30, "totalNumTrans30");
book(_sumWeightsPtLead2, "sumWeightsPtLead2");
book(_sumWeightsPtLead5, "sumWeightsPtLead5");
book(_sumWeightsPtLead30, "sumWeightsPtLead30");
_htmp_num_dphi_2 = std::make_shared<YODA::Profile1D>(refData(1, 1, 1).binning());
_htmp_num_dphi_5 = std::make_shared<YODA::Profile1D>(refData(1, 1, 2).binning());
_htmp_num_dphi_30 = std::make_shared<YODA::Profile1D>(refData(1, 1, 3).binning());
_htmp_pt_dphi_2 = std::make_shared<YODA::Profile1D>(refData(2, 1, 1).binning());
_htmp_pt_dphi_5 = std::make_shared<YODA::Profile1D>(refData(2, 1, 2).binning());
_htmp_pt_dphi_30 = std::make_shared<YODA::Profile1D>(refData(2, 1, 3).binning());
}
/// Do the analysis
void analyze(const Event& event) {
// Trigger
const bool trigger = apply<TriggerCDFRun0Run1>(event, "Trigger").minBiasDecision();
if (!trigger) vetoEvent;
// Get jets, sorted by pT
const Jets jets = apply<JetFinder>(event, "TrackJet").jetsByPt();
if (jets.empty()) vetoEvent;
const Jet jet1 = jets.front();
const double ptLead = jet1.pT();
// Cut on highest pT jet: combined 0.5 GeV < pT(lead) < 50 GeV
if (ptLead / GeV < 0.5) vetoEvent;
if (ptLead / GeV > 50.0) vetoEvent;
// Count sum of all event weights in three pT_lead regions
if (ptLead / GeV > 2.0) _sumWeightsPtLead2->fill();
if (ptLead / GeV > 5.0) _sumWeightsPtLead5->fill();
if (ptLead / GeV > 30.0) _sumWeightsPtLead30->fill();
// Run over tracks
double ptSumToward(0.0), ptSumAway(0.0), ptSumTrans(0.0);
size_t numToward(0), numTrans(0), numAway(0);
_htmp_num_dphi_2->reset();
_htmp_num_dphi_5->reset();
_htmp_num_dphi_30->reset();
_htmp_pt_dphi_2->reset();
_htmp_pt_dphi_5->reset();
_htmp_pt_dphi_30->reset();
// Final state charged particles
/// @todo Non-trackjet track efficiencies are corrected?
const Particles& tracks = apply<FinalState>(event, "FS").particles();
for (const Particle& p : tracks) {
const double dPhi = deltaPhi(p, jet1);
const double pT = p.pT();
if (dPhi < PI / 3.0) {
ptSumToward += pT;
++numToward;
}
else if (dPhi < 2 * PI / 3.0) {
ptSumTrans += pT;
++numTrans;
// Fill transverse pT distributions
if (ptLead / GeV > 2.0) {
_ptTrans2->fill(pT / GeV);
_totalNumTrans2->fill();
}
if (ptLead / GeV > 5.0) {
_ptTrans5->fill(pT / GeV);
_totalNumTrans5->fill();
}
if (ptLead / GeV > 30.0) {
_ptTrans30->fill(pT / GeV);
_totalNumTrans30->fill();
}
}
else {
ptSumAway += pT;
++numAway;
}
// Fill tmp histos to bin event's track Nch & pT in dphi
const double dPhideg = 180 * dPhi / M_PI;
if (ptLead / GeV > 2.0) {
_htmp_num_dphi_2->fill(dPhideg, 1);
_htmp_pt_dphi_2->fill(dPhideg, pT / GeV);
}
if (ptLead / GeV > 5.0) {
_htmp_num_dphi_5->fill(dPhideg, 1);
_htmp_pt_dphi_5->fill(dPhideg, pT / GeV);
}
if (ptLead / GeV > 30.0) {
_htmp_num_dphi_30->fill(dPhideg, 1);
_htmp_pt_dphi_30->fill(dPhideg, pT / GeV);
}
}
// Update the "proper" dphi profile histograms
for (int i = 1; i < 51;
i++) { ///< @todo Should really explicitly iterate over nbins for each temp histo
if (ptLead / GeV > 2.0) {
const double x2 = _htmp_pt_dphi_2->bin(i).xMid();
const double num2 = (_htmp_num_dphi_2->bin(i).numEntries()) ? _htmp_num_dphi_2->bin(i).yMean()
: 0.0;
const double pt2 = (_htmp_num_dphi_2->bin(i).numEntries()) ? _htmp_pt_dphi_2->bin(i).yMean() : 0.0;
_numvsDeltaPhi2->fill(x2, num2);
_pTvsDeltaPhi2->fill(x2, pt2);
}
if (ptLead / GeV > 5.0) {
const double x5 = _htmp_pt_dphi_5->bin(i).xMid();
const double num5 = (_htmp_num_dphi_5->bin(i).numEntries()) ? _htmp_num_dphi_5->bin(i).yMean()
: 0.0;
const double pt5 = (_htmp_num_dphi_5->bin(i).numEntries()) ? _htmp_pt_dphi_5->bin(i).yMean() : 0.0;
_numvsDeltaPhi5->fill(x5, num5);
_pTvsDeltaPhi5->fill(x5, pt5);
}
if (ptLead / GeV > 30.0) {
const double x30 = _htmp_pt_dphi_30->bin(i).xMid();
const double num30 = (_htmp_num_dphi_30->bin(i).numEntries()) ? _htmp_num_dphi_30->bin(i).yMean()
: 0.0;
const double pt30 = (_htmp_num_dphi_30->bin(i).numEntries()) ? _htmp_pt_dphi_30->bin(i).yMean()
: 0.0;
_numvsDeltaPhi30->fill(x30, num30);
_pTvsDeltaPhi30->fill(x30, pt30);
}
}
// Log some event details about pT
MSG_DEBUG("pT [lead; twd, away, trans] = [" << ptLead << "; " << ptSumToward << ", " << ptSumAway
<< ", " << ptSumTrans << "]");
// Update the pT profile histograms
_ptsumTowardMB->fill(ptLead / GeV, ptSumToward / GeV);
_ptsumTowardJ20->fill(ptLead / GeV, ptSumToward / GeV);
_ptsumTransMB->fill(ptLead / GeV, ptSumTrans / GeV);
_ptsumTransJ20->fill(ptLead / GeV, ptSumTrans / GeV);
_ptsumAwayMB->fill(ptLead / GeV, ptSumAway / GeV);
_ptsumAwayJ20->fill(ptLead / GeV, ptSumAway / GeV);
// Log some event details about Nch
MSG_DEBUG("N [twd, away, trans] = [" << ptLead << "; " << numToward << ", " << numTrans << ", "
<< numAway << "]");
// Update the N_track profile histograms
_numTowardMB->fill(ptLead / GeV, numToward);
_numTowardJ20->fill(ptLead / GeV, numToward);
_numTransMB->fill(ptLead / GeV, numTrans);
_numTransJ20->fill(ptLead / GeV, numTrans);
_numAwayMB->fill(ptLead / GeV, numAway);
_numAwayJ20->fill(ptLead / GeV, numAway);
}
/// Normalize histos
void finalize() {
normalize(_ptTrans2, *_totalNumTrans2 / *_sumWeightsPtLead2);
normalize(_ptTrans5, *_totalNumTrans5 / *_sumWeightsPtLead5);
normalize(_ptTrans30, *_totalNumTrans30 / *_sumWeightsPtLead30);
}
/// @}
private:
/// Sum total number of charged particles in the trans region, in 3 \f$ p_\perp^\text{lead} \f$ bins.
CounterPtr _totalNumTrans2, _totalNumTrans5, _totalNumTrans30;
/// Sum the total number of events in 3 \f$ p_\perp^\text{lead} \f$ bins.
CounterPtr _sumWeightsPtLead2, _sumWeightsPtLead5, _sumWeightsPtLead30;
/// @name Histogram collections
/// @{
// The sumpt vs. dphi and Nch vs. dphi histos
Profile1DPtr _numvsDeltaPhi2, _numvsDeltaPhi5, _numvsDeltaPhi30;
Profile1DPtr _pTvsDeltaPhi2, _pTvsDeltaPhi5, _pTvsDeltaPhi30;
/// Profile histograms, binned in the \f$ p_T \f$ of the leading jet, for
/// the \f$ p_T \f$ sum in the toward, transverse and away regions.
Profile1DPtr _ptsumTowardMB, _ptsumTransMB, _ptsumAwayMB;
Profile1DPtr _ptsumTowardJ20, _ptsumTransJ20, _ptsumAwayJ20;
/// Profile histograms, binned in the \f$ p_T \f$ of the leading jet, for
/// the number of charged particles per jet in the toward, transverse and
/// away regions.
Profile1DPtr _numTowardMB, _numTransMB, _numAwayMB;
Profile1DPtr _numTowardJ20, _numTransJ20, _numAwayJ20;
/// Histogram of \f$ p_T \f$ distribution for 3 different \f$ p_{T1} \f$ IR cutoffs.
Histo1DPtr _ptTrans2, _ptTrans5, _ptTrans30;
// Temporary histos that bin N and pT in dphi
YODA::Profile1DPtr _htmp_num_dphi_2, _htmp_num_dphi_5, _htmp_num_dphi_30;
YODA::Profile1DPtr _htmp_pt_dphi_2, _htmp_pt_dphi_5, _htmp_pt_dphi_30;
/// @}
};
RIVET_DECLARE_ALIASED_PLUGIN(CDF_2001_I564673, CDF_2001_S4751469);
} ```
Aliases: - CDF_2001_S4751469