1#include <Rivet/Analysis.hh>
 2#include <Rivet/Projections/AliceCommon.hh>
 3#include <Rivet/Projections/HepMCHeavyIon.hh>
 4
 5namespace Rivet {
 6
 7  /// Dummy analysis for centrality calibration in Pb-Pb at 5.02TeV
 8  ///
 9  /// @author Christian Holm Christensen <cholm@nbi.dk>
10  class ALICE_2015_PBPBCentrality : public Analysis {
11  public:
12
13    /// Constructor 
14    ALICE_2015_PBPBCentrality()
15      : Analysis("ALICE_2015_PBPBCentrality")
16    {    }
17
18    /// Initialize this analysis. 
19    void init() {
20      ALICE::V0AndTrigger v0and;
21      declare<ALICE::V0AndTrigger>(v0and,"V0-AND");
22
23      ALICE::V0MMultiplicity v0m;
24      declare<ALICE::V0MMultiplicity>(v0m,"V0M");
25
26       // Access the HepMC heavy ion info
27      declare(HepMCHeavyIon(), "HepMC");
28
29      book(_v0m, "V0M", 500, -5.0, 39995.0);
30      book(_imp, "V0M_IMP",100,0,20);
31    }
32
33
34    /// Analyse a single event.
35    void analyze(const Event& event) {
36      // Get and fill in the impact parameter value if the information is valid.
37      _imp->fill(apply<HepMCHeavyIon>(event, "HepMC").impact_parameter());
38	  
39      // Check if we have any hit in either V0-A or -C.  If not, the
40      // event is not selected and we get out.
41      if (!apply<ALICE::V0AndTrigger>(event,"V0-AND")()) return;
42
43      // Fill in the V0 multiplicity for this event
44      _v0m->fill(apply<ALICE::V0MMultiplicity>(event,"V0M")());
45    }
46
47
48    /// Finalize this analysis
49    void finalize() {
50      _v0m->normalize();
51      _imp->normalize();
52    }
53
54    /// The distribution of V0M multiplicity
55    Histo1DPtr _v0m;
56    /// The distribution of impact parameters
57    Histo1DPtr _imp;
58
59  };
60
61
62  // The hook for the plugin system
63  RIVET_DECLARE_PLUGIN(ALICE_2015_PBPBCentrality);
64
65}