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