Rivet analyses
Measurement of proton-proton elastic scattering and total cross section at $\sqrt{s} = 7~\TeV$
Experiment: TOTEM (LHC)
Inspire ID: 1220862
Status: VALIDATED
Authors: - Sercan Sen - Peter Skands
References: - CERN-PH-EP-2012-239 - http://cds.cern.ch/record/1472948
Beams: p+ p+
Beam energies: (3500.0, 3500.0)GeV
Run details: - Elastic events only.
Measurement of the elastic differential cross-section in proton-proton interactions at a centre-of-mass energy $\sqrt{s} = 7~\TeV$ at the LHC. The data, which cover the |t| range $0.005-0.2~\GeV^2$, were collected using Roman Pot detectors very close to the outgoing beam in October 2011, allowing the precise extrapolation down to the optical point, t = 0, and hence the derivation of the elastic as well as the total cross-section via the optical theorem.
Source
code:TOTEM_2012_I1220862.cc
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
namespace Rivet {
/// TOTEM elastic and total cross-section measurement
class TOTEM_2012_I1220862 : public Analysis {
public:
TOTEM_2012_I1220862()
: Analysis("TOTEM_2012_I1220862")
{ }
void init() {
declare(ChargedFinalState(), "CFS");
book(_hist_tlow ,1, 1, 1);
book(_hist_thigh ,2, 1, 1);
book(_hist_sigma ,3, 1, 1);
}
void analyze(const Event& event) {
const ChargedFinalState& cfs = apply<ChargedFinalState>(event, "CFS");
if (cfs.size() > 2) MSG_DEBUG("Final state includes more than two charged particles!");
_hist_sigma->fill(sqrtS()/GeV);
for (const Particle& p : cfs.particles(Cuts::eta > 0)) { // && Cuts::pid == PID::PROTON)) {
if (p.pid() != PID::PROTON) continue;
const double t = sqr(p.pT());
_hist_tlow->fill(t);
_hist_thigh->fill(t);
}
}
void finalize() {
normalize(_hist_tlow, crossSection()/millibarn);
normalize(_hist_thigh, crossSection()/millibarn);
normalize(_hist_sigma, crossSection()/millibarn);
}
private:
Histo1DPtr _hist_tlow, _hist_thigh, _hist_sigma;
};
RIVET_DECLARE_ALIASED_PLUGIN(TOTEM_2012_I1220862, TOTEM_2012_002);
}Aliases: - TOTEM_2012_002