1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/ImpactParameterProjection.hh"
 4#include "Rivet/Analyses/AliceCommon.hh"
 5
 6namespace Rivet {
 7
 8
 9  class ALICE_2015_CENT_PP : public Analysis {
10  public:
11
12    RIVET_DEFAULT_ANALYSIS_CTOR(ALICE_2015_CENT_PP);
13
14    // Book histograms, initialize projections.
15    void init() {
16      declare(ALICE::V0AndTrigger(),"V0-AND");
17      declare(ALICE::V0MMultiplicity(),"V0M");
18      declare(ImpactParameterProjection(), "IMP");
19      book(_v0m, "V0M",100,0,200);
20      book(_imp, "V0M_IMP",100,0,20);
21    }
22
23
24    // Per-event analysis
25    void analyze(const Event& event) {
26      _imp->fill(apply<SingleValueProjection>(event,"IMP")());
27
28      // Check if we have any hit in either V0-A or -C.  If not, the
29      // event is not selected and we get out.
30      if (!apply<ALICE::V0AndTrigger>(event,"V0-AND")()) return;
31
32      // Fill in the V0 multiplicity for this event
33      _v0m->fill(apply<ALICE::V0MMultiplicity>(event,"V0M")());
34    }
35
36
37    void finalize() {
38      _v0m->normalize();
39      _imp->normalize();
40    }
41
42
43    Histo1DPtr _v0m;
44    Histo1DPtr _imp;
45  };
46
47
48  RIVET_DECLARE_PLUGIN(ALICE_2015_CENT_PP);
49
50}