rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

STAR_2008_S7993412

Di-hadron correlations in $d$-Au at 200 GeV
Experiment: STAR (RHIC d-Au 200 GeV)
Inspire ID: 810030
Status: UNVALIDATED
Authors:
  • Christine Nattrass
  • Hendrik Hoeth
References: Beams: p+ p+
Beam energies: (100.0, 100.0) GeV
Run details:
  • $d$-Au at 200 GeV, but use $pp$ Monte Carlo! See description.

Correlation in $\eta$ and $\phi$ between the charged hadron with the highest pT (``trigger particle'') and the other charged hadrons in the event (``associated particles''). The data was collected in d-Au collisions at 200 GeV. Nevertheless, it is very proton--proton like and can therefore be compared to $pp$ Monte Carlo (not for tuning, but for qualitative studies.)

Source code: STAR_2008_S7993412.cc
 1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/ChargedFinalState.hh"
 4
 5namespace Rivet {
 6
 7
 8  /// STAR di-hadron correlations in d-Au at 200 GeV
 9  class STAR_2008_S7993412 : public Analysis {
10  public:
11
12    RIVET_DEFAULT_ANALYSIS_CTOR(STAR_2008_S7993412);
13
14
15    /// @name Analysis methods
16    /// @{
17
18    /// Book projections and histograms
19    void init() {
20      ChargedFinalState fs((Cuts::etaIn(-1.0, 1.0) && Cuts::pT >=  1.0*GeV));
21      declare(fs, "FS");
22
23      book(_h_Y_jet_trigger ,1, 1, 1);
24      book(_h_Y_jet_associated ,2, 1, 1);
25    }
26
27
28    /// Do the analysis
29    void analyze(const Event& event) {
30      // Skip if the event is empty
31      const FinalState& fs = apply<FinalState>(event, "FS");
32      if (fs.empty()) {
33        MSG_DEBUG("Skipping event " << numEvents() << " because no final state found ");
34        vetoEvent;
35      }
36
37      for (const Particle& tp : fs.particles()) {
38        const double triggerpT = tp.pT();
39        if (triggerpT >= 2.0 && triggerpT < 5.0) {
40          int n_associated = 0;
41          for (const Particle& ap : fs.particles()) {
42            if (!inRange(ap.pT()/GeV, 1.5, triggerpT)) continue;
43            if (deltaPhi(tp.phi(), ap.phi()) > 1) continue;
44            if (fabs(tp.eta() - ap.eta()) > 1.75) continue;
45            n_associated += 1;
46          }
47          //const double dPhidEta = 2 * 2*1.75;
48          //_h_Y_jet_trigger->fill(triggerpT, n_associated/dPhidEta);
49          _h_Y_jet_trigger->fill(triggerpT, n_associated);
50        }
51      }
52    }
53
54
55    /// Finalize
56    // void finalize() {    }
57
58    /// @}
59
60
61  private:
62
63    /// @name Histograms
64    /// @{
65    Profile1DPtr _h_Y_jet_trigger;
66    Profile1DPtr _h_Y_jet_associated;
67    /// @}
68
69  };
70
71
72
73  RIVET_DECLARE_ALIASED_PLUGIN(STAR_2008_S7993412, STAR_2008_I810030);
74
75}