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_I810030 : public Analysis {
10  public:
11
12    RIVET_DEFAULT_ANALYSIS_CTOR(STAR_2008_I810030);
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_I810030, STAR_2008_S7993412);
74
75}