Page MenuHomeHEPForge
Files F8723790

CDF_1988_S1865951.cc
No OneTemporary
Actions

CDF_1988_S1865951.cc
View Options

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/RivetAIDA.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/TriggerCDFRun0Run1.hh"
namespace Rivet {
/// @brief CDF track \f$ p_\perp \f$ distributions at 630 and 1800 GeV
class CDF_1988_S1865951 : public Analysis {
public:
/// Constructor
CDF_1988_S1865951()
: Analysis("CDF_1988_S1865951")
{
setBeams(PROTON, ANTIPROTON);
setNeedsCrossSection(true);
_sumWTrig = 0;
}
/// @name Analysis methods
//@{
/// Book histograms and set up projections
void init() {
// Set up projections
addProjection(TriggerCDFRun0Run1(), "Trigger");
const ChargedFinalState cfs(-1.0, 1.0, 0.4*GeV);
addProjection(cfs, "CFS");
// Book histo
if (fuzzyEquals(sqrtS()/GeV, 1800, 1E-3)) {
_hist_pt = bookHistogram1D(1, 1, 1);
} else if (fuzzyEquals(sqrtS()/GeV, 630, 1E-3)) {
_hist_pt = bookHistogram1D(2, 1, 1);
}
}
/// Do the analysis
void analyze(const Event& event) {
// Trigger
const bool trigger = applyProjection<TriggerCDFRun0Run1>(event, "Trigger").minBiasDecision();
if (!trigger) vetoEvent;
const double weight = event.weight();
_sumWTrig += weight;
const FinalState& trackfs = applyProjection<ChargedFinalState>(event, "CFS");
foreach (Particle p, trackfs.particles()) {
const double pt = p.momentum().pT()/GeV;
// Effective weight for d3sig/dp3 = weight / ( Delta eta * 2pi * pt ), with Delta(eta) = 2
const double eff_weight = weight/(2*2*TWOPI*pt);
_hist_pt->fill(pt, eff_weight);
}
}
/// Scale histos
void finalize() {
scale(_hist_pt, crossSectionPerEvent()/millibarn);
}
//@}
private:
/// @name Counters
//@{
double _sumWTrig;
//@}
/// @name Histos
//@{
AIDA::IHistogram1D* _hist_pt;
//@}
};
// This global object acts as a hook for the plugin system
AnalysisBuilder<CDF_1988_S1865951> plugin_CDF_1988_S1865951;
}

File Metadata

Mime Type
text/x-c
Expires
Mon, Jan 20, 9:29 PM (1 d, 4 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4242475
Default Alt Text
CDF_1988_S1865951.cc (2 KB)

Event Timeline