TOTEM_2012_I1220862.cc

Go to the documentation of this file.

// -*- 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");
      _hist_tlow  = bookHisto1D(1, 1, 1);
      _hist_thigh = bookHisto1D(2, 1, 1);
      _hist_sigma = bookHisto1D(3, 1, 1);
    }


    void analyze(const Event& event) {
      const double weight = event.weight();

      const ChargedFinalState& cfs = apply<ChargedFinalState>(event, "CFS");
      if (cfs.size() > 2) MSG_WARNING("Final state includes more than two charged particles!");
      _hist_sigma->fill(sqrtS()/GeV, weight);

      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, weight);
        _hist_thigh->fill(t, weight);
      }
    }


    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;

  };


  DECLARE_ALIASED_RIVET_PLUGIN(TOTEM_2012_I1220862, TOTEM_2012_002);

}